JPH08247931A - 生化学画像解析装置 - Google Patents

生化学画像解析装置

Info

Publication number
JPH08247931A
JPH08247931A JP7222207A JP22220795A JPH08247931A JP H08247931 A JPH08247931 A JP H08247931A JP 7222207 A JP7222207 A JP 7222207A JP 22220795 A JP22220795 A JP 22220795A JP H08247931 A JPH08247931 A JP H08247931A
Authority
JP
Japan
Prior art keywords
data
coordinate
lane
image
image data
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
Application number
JP7222207A
Other languages
English (en)
Other versions
JP3571819B2 (ja
Inventor
Makiko Fukushima
眞喜子 福島
Yoshiko Shiimori
佳子 椎森
Makoto Sato
真 佐藤
Masato Some
真人 染
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.)
Fujifilm Holdings Corp
Original Assignee
Fuji Photo Film Co Ltd
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 Fuji Photo Film Co Ltd filed Critical Fuji Photo Film Co Ltd
Priority to JP22220795A priority Critical patent/JP3571819B2/ja
Priority to US08/581,261 priority patent/US6144758A/en
Publication of JPH08247931A publication Critical patent/JPH08247931A/ja
Application granted granted Critical
Publication of JP3571819B2 publication Critical patent/JP3571819B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/20Image preprocessing
    • G06V10/25Determination of region of interest [ROI] or a volume of interest [VOI]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/60Type of objects
    • G06V20/69Microscopic objects, e.g. biological cells or cellular parts

Abstract

(57)【要約】 【課題】 本発明は、操作者による煩雑な操作を必要と
することなく、迅速に、生化学測定画像中の定量化ある
いは定量解析すべき関心領域を画定することのできる生
化学画像解析装置を提供することを課題する。 【解決手段】 本発明の課題は、少なくとも一つの試料
を一次元的に展開して得られた試料中の標識物質の位置
情報を、電気信号に変換して得られたデータに基づき、
生成された生化学画像の画像データを記憶する画像デー
タ記憶手段50と、画像データ記憶手段に記憶された画
像データの中から画像データを選択して、所定の処理を
施すことにより、平面座標系に展開された表示画像デー
タを生成する表示画像データ生成手段102と、表示画
像データに基づき、画像を表示する表示手段60と、表
示画像データを、表示画像データが展開される平面座標
系の一方の座標軸方向に分割して、一方の座標軸方向の
定量すべき領域を画定するレーンを形成するための第1
の分割座標データを生成する第1の画像データ分割手段
106とを備えた生化学画像解析装置により達成され
る。

Description

【発明の詳細な説明】
【0001】
【産業上の利用分野】本発明は、生化学画像解析装置に
関するものであり、さらに詳細には、少なくとも一種の
試料を一次元的に展開して得られた試料中の標識物質の
位置情報を含む生化学測定画像上の定量解析すべき領域
を自動的に画定することのできる生化学画像解析装置に
関するものである。
【0002】
【従来の技術】放射性標識を付与した物質を、生物体に
投与した後、その生物体あるいはその生物体の組織の一
部を試料とし、この試料を、高感度X線フィルムなどの
放射線フィルムに一定時間重ね合わせることによって、
放射線フィルムを感光させ或いは露光し、放射線フィル
ムの感光された部位に基づき、試料中の放射性標識物質
の位置情報を得るようにしたオートラジオグラフィ検出
方法や、蛋白質、核酸配列などの固定された高分子を、
化学発光物質と接触して、化学発光を生じさせる標識物
質により、選択的に標識し、標識物質によって選択的に
標識された高分子と、化学発光物質とを接触させて、化
学発光物質と標識物質との接触によって生ずる可視光波
長域の化学発光を検出することによって、遺伝子情報な
どの高分子に関する情報を得るようにした化学発光検出
方法が知られている。これらの方法は、従来、検出材料
として、写真フイルムを用い、写真フイルム上に、放射
線画像や化学発光画像を記録し、目視によって、可視画
像を検出することによって、おこなわれていたが、検出
材料として、写真フイルムを用いる場合には、オートラ
ジオグラフィ検出方法にあっては、放射線フイルムの感
度が低く、画像記録に多大な時間を要するという問題が
あり、また、化学発光検出方法にあっては、微弱な化学
発光を確実に検出するために、γ値の高い高感度フイル
ムを用いる必要があるが、γ値の高い高感度フイルムを
用いるときは、確実に、特性曲線の直線部を用いて、露
光することが困難であって、露光ミスが多く、露光条件
を変えて、繰り返し、露光する必要があるという問題が
あり、さらには、いずれの方法にあっても、現像処理と
いう化学的処理が必要不可欠であって、操作が煩雑であ
るという問題を有している。
【0003】そこで、従来の写真フイルムに代えて、放
射線、可視光、電子線などが照射されると、そのエネル
ギーを吸収して、蓄積し、その後に、特定の波長域の電
磁波を用いて励起すると、照射された放射線、可視光、
電子線などのエネルギーの量に応じた光量の輝尽光を発
する特性を有する輝尽性蛍光体を、放射線、可視光、電
子線などの検出材料として用い、輝尽性蛍光体から発せ
られた輝尽光を、光電的に検出して、ディジタル信号に
変換し、得られた画像データに所定の画像処理を施した
後に、画像を、CRT画面などの表示手段あるいは写真
フイルム上に再生するようにしたオートラジオグラフィ
検出方法および化学発光検出方法が提案されている(た
とえば、特公平1−60784号公報、特公平1−60
782号公報、特公平4−3952号公報、米国特許第
5,028,793号、英国特許出願公開GB第2,2
46,197Aなど)。この輝尽性蛍光体を用いた検出
方法によれば、現像処理という化学的処理が不必要であ
るだけでなく、オートラジオグラフィ検出方法にあって
は、露光時間を大幅に短縮することができ、化学発光検
出方法にあっては、露光ミスが少なく、容易に、露光を
おこなうことができるという利点があり、さらには、デ
ィジタル信号に変換された後に、画像が再生されるの
で、画像データに、信号処理を施すことによって、所望
のように、画像を再生し、あるいは、コンピュータによ
る定量解析が可能になり、好ましい。
【0004】したがって、かかる方法を利用し、少なく
とも一種の試料を一次元的に展開して得られた試料中の
放射性標識物質の位置情報を含むオートラジオグラフィ
画像あるいは試料中の標識物質の位置情報を含む化学発
光画像を読み取って、画像データを生成し、解析する画
像解析装置が提案されており、薄層クロマトグラフィー
(TLC:Thin Layer Chromatography)法によって得ら
れた画像、サザン・ブロット・ハイブリタイゼーション
法を利用した遺伝子の画像、ポリアクリルアミドゲル電
気泳動法によって、蛋白質の分離、同定、あるいは、分
子量、特性の評価などをおこなうための画像の解析に利
用されている。これらのうち、薬物代謝研究の分野にお
いて、薬物代謝を調べるために、しばしば用いられてい
る薄層クロマトグラフィー法は、固定相として、薬物と
の相互作用をもつ物質を均一な薄層にして、その上で、
薬物あるいは代謝物を分離するものであり、オートラジ
オグラフィによる場合には、実験動物に放射性物質(放
射性同位体:RI)で標識した薬物を投与し、その薬物
が体内で変化した場合の成分を分析することができる。
たとえば、放射性物質(放射性同位体:RI)で標識化
した薬物を実験動物に投与してから一定時間経過後に、
実験動物の特定の部位より採取した尿、血液、組織など
に所定の処理を施し、それを、ガラス板の上に0.5mm
程度の薄いシリカゲルの粉を塗布したTLCプレート上
の所定の位置にスポットし、これを展開溶媒に浸して、
クロマトグラフィー展開すると、TLCプレート上に、
成分ごとに分離されたスポットが形成される。このTL
Cプレートを、蓄積性蛍光体シートに転写し、画像読み
取り装置を用いて、蓄積性蛍光体シートを電磁波によっ
て走査して、蓄積性蛍光体シートから放出される輝尽光
を、光電的に検出し、試料から蓄積性蛍光体シートに転
写された放射性標識物質の位置情報を電気信号に変換し
てデータを生成する。画像解析装置は、こうして得られ
たデータを受信し、これをハードディスクなどのデータ
記憶手段に、記憶することによって、蓄積性蛍光体シー
トに転写された放射性標識物質の位置情報を含む画像デ
ータを記憶し、さらに、データ記憶手段に記憶された画
像データから所望のデータを選択して、これに所望の信
号処理を施し、あるいは、所定の拡大率で拡大すること
により、放射性標識物質の位置情報を精度良く得るのに
適した画像をCRTの画面などの表示手段上に再生する
ことができる。
【0005】このような画像解析装置においては、画像
データの所望の領域を関心領域として画定し、蓄積性蛍
光体シートからの発光量を、画像データの関心領域に含
まれる画像を構成する画素の濃度を数値化して、その総
和を求めることにより、画像の所望の領域の濃度、すな
わち、蓄積性蛍光体シートからの発光量を定量化するこ
とができ、さらに、複数の関心領域をグループ化して、
グループに属する関心領域の各々の濃度の比、すなわ
ち、蓄積性蛍光体シートからの発光量の比などを算出し
て、画像の所望の領域間の濃度、すなわち、蓄積性蛍光
体シートからの発光量を比較し、定量解析ができるよう
に構成されている。たとえば、前述した薄層クロマトグ
ラフィー法を用いた薬物代謝研究の分野において、TL
Cプレート上に、同時に複数の物質をスポットして、ク
ロマトグラフィー展開する場合、たとえば、ある実験動
物に薬物を投与してから、30分後、1時間後、2時間
後、5時間後に、その実験動物から血液を採取し、それ
らをTLCプレート上に横方向にスポットし、クロマト
グラフィー展開する場合に、このように生成されたTL
Cプレートを転写した蓄積性蛍光体シートから得られた
画像データから、TLCプレート上のスポットが形成さ
れている位置に対応する領域を関心領域として画定し、
さらに、複数の所定の領域をグループ化し、各関心領域
の濃度を数値化して、グループに属する各領域の濃度の
比を調べることができる。
【0006】このような定量化処理および定量解析をす
るために、オートラジオグラフィ画像解析装置あるいは
化学発光画像解析装置は、表示装置の画面上に形成され
た画像中に、矩形などで囲まれた関心領域を画定し、関
心領域内に存在する画像の濃度データを得ることができ
るように構成されている。
【0007】
【発明の解決しようとする課題】しかしながら、従来の
オートラジオグラフィ画像解析装置あるいは化学発光画
像解析装置においては、操作者が定量しようとする領域
ごとに、関心領域を画定する必要があり、その操作が煩
雑であるという問題点があった。たとえば、前述したよ
うに、TLCプレート上に、同時に二種以上の試料をス
ポットして、クロマトグラフィー展開したものを、転写
した蓄積性蛍光体シートから画像データを得た場合に、
画像の所望の関心領域間の濃度、すなわち、蓄積性蛍光
体シートからの発光量を比較し、定量解析をするために
は、操作者は、CRTなどの表示手段に表示された画像
を観察して、TLCプレート上に形成されているスポッ
トごとに、複数の関心領域を設定する必要がある。しか
しながら、このような方法では、得られた画像データに
形成されたスポットを含む所望の関心領域の定量処理
を、迅速にすることができないという問題点があった。
さらに、たとえば、ある試料をクロマトグラフィー展開
することにより得られた特定のレーン内に含まれる複数
のスポットに対応する関心領域の濃度データの比などを
求めて、定量解析をしようとする場合、或いは、クロマ
トグラフィー展開された二以上の試料間で、対応する位
置に形成されたスポットに対応する関心領域の濃度デー
タの比などを求めて、定量解析をしようとする場合に
は、設定した関心領域に基づき、複数の関心領域をグル
ープ化する必要がある。しかし、画像を観察して、操作
者が、このような複数の関心領域をグループ化するので
は、複数の所望の領域の定量解析を、迅速にすることが
できないという問題点があった。
【0008】これらは、薄層クロマトグラフィー法によ
って得られた画像に限らず、少なくとも一種の試料を一
次元的に展開して得られた試料中の放射性標識物質の位
置情報を含むオートラジオグラフィ画像や試料中の標識
物質の位置情報を含む化学発光画像の解析をおこなう必
要があるサザン・ブロット・ハイブリタイゼーション法
を利用した遺伝子の画像、ポリアクリルアミドゲル電気
泳動法によって、蛋白質の分離、同定、あるいは、分子
量、特性の評価などをおこなうための画像の解析をする
場合にも同様に問題となっていた。また、このことは、
試料中の標識物質の位置情報を、従来のオートラジオグ
ラフィや化学発光検出方法と同様に、写真フィルムに記
録し、可視画像を、一旦、再生した後に、再生された可
視画像を光電的に読み取って、電気信号化し、得られた
画像データに、所望の信号処理を施すことにより、標識
物質の位置情報を、画像として、CRT画面などの表示
手段上に再生する場合にも、同様に問題になる。
【0009】
【発明の目的】本発明は、操作者による煩雑な操作を必
要とすることなく、迅速に、生化学測定画像中の定量化
あるいは定量解析すべき関心領域を画定することのでき
る生化学画像解析装置を提供することを目的とするもの
である。
【0010】
【発明の構成】本発明のかかる目的は、少なくとも一つ
の試料を一次元的に展開して得られた試料中の標識物質
の位置情報を、電気信号に変換して得られたデータに基
づき、生成された生化学画像の画像データを記憶する画
像データ記憶手段と、前記画像データ記憶手段に記憶さ
れた画像データの中から画像データを選択して、所定の
処理を施すことにより、平面座標系に展開された表示画
像データを生成する表示画像データ生成手段と、前記表
示画像データに基づき、画像を表示する表示手段と、前
記表示画像データを、前記表示画像データが展開される
前記平面座標系の一方の座標軸方向に分割して、前記一
方の座標軸方向の定量すべき領域を画定するレーンを形
成するための第1の分割座標データを生成する第1の画
像データ分割手段とを備えた生化学画像解析装置により
達成される。本発明の好ましい実施態様においては、生
化学画像解析装置は、さらに、前記一方の座標軸方向に
分割された前記表示画像データを、前記平面座標系の他
方の座標軸方向に分割し、前記レーンごとに、前記他方
の座標軸方向の定量すべき領域を画定するスポットを形
成するための第2の分割座標データを生成する第2の画
像データ分割手段を備えている。
【0011】本発明のさらに好ましい実施態様において
は、生化学画像解析装置は、さらに、前記表示画像デー
タの中から、前記表示手段に表示されるべき画像のノイ
ズに対応するデータを除去して、ノイズ除去画像データ
を生成する第1のノイズ除去手段を備え、前記第1の画
像データ分割手段が、前記第1のノイズ除去手段によっ
てノイズに対応するデータが除去されたノイズ除去デー
タを、前記表示画像データが展開される前記平面座標系
の一方の座標軸方向に分割して、前記一方の座標軸方向
の定量すべき領域を画定するレーンを形成するための第
1の分割座標データを生成するように構成されている。
本発明のさらに好ましい実施態様においては、生化学画
像解析装置は、さらに、前記表示画像データの中から、
前記表示手段に表示されるべき画像のノイズに対応する
データを除去して、ノイズ除去画像データを生成する第
2のノイズ除去手段を備え、前記第2の画像データ分割
手段が、前記第2のノイズ除去手段によりノイズに対応
するデータが除去されたノイズ除去データを、前記平面
座標系の他方の座標軸方向に分割し、前記レーンごと
に、前記他方の座標軸方向の定量すべき領域を画定する
スポットを形成するための第2の分割座標データを生成
するように構成されている。
【0012】本発明のさらに好ましい実施態様において
は、前記第1の画像データ分割手段が、前記一方の座標
の値を共通にする画素データの有する濃度データの値を
累算し、前記一方の座標方向の周辺分布に対応するデー
タを算出する第1の周辺分布データ算出手段と、前記第
1の周辺分布データ算出手段により算出された周辺分布
に対応するデータに基づいて、前記レーンの開始点の座
標を算出するレーン開始点座標算出手段と、前記第1の
周辺分布データ算出手段により算出された周辺分布に対
応するデータに基づいて、前記レーン終了点の座標を算
出するレーン終了点座標算出手段と、前記レーン開始点
座標算出手段および前記レーン終了点座標算出手段によ
り算出された前記開始点および前記終了点の座標に基づ
いて、前記レーンを画定する分割座標データを生成する
レーン分割座標データ生成手段とを備えている。本発明
のさらに好ましい実施態様においては、前記第2の画像
データ分割手段が、前記第1の画像データ分割手段によ
り生成された第1の分割座標データに基づいて画定され
たレーンごとに、前記他方の座標の値を共通にする画素
データの有する濃度データの値を累算し、前記他方の座
標方向の周辺分布データを算出する第2の周辺分布デー
タ算出手段と、前記第2の周辺分布算出手段により算出
された周辺分布に対応するデータに基づいて、前記レー
ンごとのスポットを画定する第2の分割座標データを生
成するスポット分割座標データ生成手段を備えている。
【0013】本発明のさらに好ましい実施態様において
は、前記第1のノイズ除去手段が、前記第1の周辺分布
データ算出手段により算出された周辺分布に対応するデ
ータに対して、ノイズ除去処理をおこない、前記第1の
ノイズ除去手段によりノイズに対応するデータが除去さ
れた周辺分布に対応するデータに基づいて、前記レーン
開始点座標算出手段が、前記レーン終了点の座標を算出
し、前記レーン終了点座標算出手段が、前記レーン終了
点の座標を算出するように構成されている。本発明のさ
らに好ましい実施態様においては、前記第2のノイズ除
去手段が、前記第2の周辺分布データ算出手段により算
出された周辺分布に対応するデータに対して、ノイズ除
去処理をおこない、前記第2のノイズ除去手段によりノ
イズに対応するデータが除去された周辺分布に対応する
データに基づいて、前記スポット分割座標データ生成手
段が、第2の分割座標データを生成するように構成され
ている。本発明のさらに好ましい実施態様においては、
前記第1の画像分割手段が、前記第1のノイズ除去手段
によりノイズに対応するデータが除去されたノイズ除去
データに一次元フーリエ変換を施して、濃度データのピ
ークの座標を推定するピーク推定手段と、該ピーク推定
手段により推定されたピークの座標を補正するピーク座
標補正手段と、該ピーク座標補正手段により得られたピ
ークの座標にしたがって、第1の分割座標データを生成
するレーン分割座標データ生成手段を備えている。
【0014】本発明のさらに好ましい実施態様において
は、生化学画像解析装置は、さらに、前記第1の分割座
標データと前記第2の分割座標データとに基づいて、前
記各レーンに形成されたスポットが属すべきグループを
示すグループデータを生成するグループデータ生成手段
を備えている。本発明のさらに好ましい実施態様におい
ては、前記グループデータ生成手段が、各レーンの前記
第2の分割座標データを対比して、前記各レーンに形成
されたスポットをグループ化するスポットグループ化手
段と、前記各レーンの前記第2の分割座標データを比較
して、各レーンをグループ化するレーングループ化手段
と、該レーングループ化手段による前記各レーンのグル
ープ化の結果に基づいて、前記スポットグループ化手段
によりグループ化されたスポットのグループを補正する
スポットグループ化補正手段を備えている。本発明のさ
らに好ましい実施態様においては、前記スポット分割座
標データ生成手段が、前記第2の周辺分布データ算出手
段により算出された周辺分布に対応するデータの値を所
定値と比較して、領域開始点および領域終了点を検出す
る領域開始点/終了点検出手段と、前記領域開始点と前
記領域終了点との間に存在する前記第2の周辺分布に対
応するデータの微分値の符号が変化するピーク点を検出
するピーク点検出手段と、前記領域開始点および領域終
了点ならびに前記ピーク点に基づいて、前記スポットを
形成するための第2の分割座標データを算出する分割座
標データ算出手段を備えている。
【0015】本発明のさらに好ましい実施態様において
は、生化学画像解析装置は、さらに、前記表示画像デー
タの中から、前記表示手段に表示されるべき画像の背景
ノイズに対応するデータを除去して、ノイズ除去画像デ
ータを生成する背景ノイズ除去手段を備え、前記第1の
画像データ分割手段が、前記ノイズ除去画像データを、
前記表示画像データが展開される前記平面座標系の一方
の座標軸方向に分割して、前記一方の座標軸方向の定量
すべき領域を画定するレーンを形成するための第1の分
割座標データを生成するように構成されている。本発明
のさらに好ましい実施態様においては、前記背景ノイズ
除去手段が、前記表示画像データを構成する複数の画素
を含む領域を画定し、前記領域の所定の位置にある画素
データの有する濃度データが、前記領域の他の位置にあ
る画素データの有する濃度データと異なる場合に、前記
所定の位置にある画素データの有する濃度データの値
を、所定の値に変換し、前記領域中の孤立した濃度デー
タを有する画素データを除去する孤立点除去手段を備え
ている。本発明のさらに好ましい実施態様においては、
前記背景ノイズ除去手段が、さらに、前記孤立点除去手
段から出力されたデータを二値化する二値化手段と、前
記二値化手段から出力されたデータから、複数の画素デ
ータを含む領域を画定し、前記領域の所定の位置にある
画素データの有する濃度データが、前記領域の他の位置
にある画素データの有する濃度データと異なる場合に、
前記所定の位置にある画素データの有する濃度データの
値を、所定の値に変換し、該領域中の孤立した濃度デー
タを有する画素データを除去する第2の孤立点除去手段
を備えている。
【0016】本発明のさらに好ましい実施態様において
は、前記濃度データの値が、大きくなるにしたがって、
前記濃度データを有する画素データに対応する画像が、
前記表示手段上に白く表示されるように構成され、前記
背景ノイズ除去手段が、さらに、前記表示画像データを
構成する画素データの有する濃度データを、所定の範囲
に拡張して、コントラストを強調するコントラスト強調
手段と、前記コントラストの強調された表示画像データ
の濃度ヒストグラムを生成し、前記濃度ヒストグラムの
面積の所定の割合を占める点より大きな値を有する濃度
データに、濃度データの最大値を付与する濃度変換手段
とを有し、前記濃度変換手段により濃度データの変換さ
れたデータが、前記孤立点除去手段に与えられるように
構成されている。本発明のさらに好ましい実施態様にお
いては、前記濃度データの値が、大きくなるにしたがっ
て、前記濃度データを有する画素データに対応する画像
が、前記表示手段上に白く表示されるように構成され、
前記背景ノイズ除去手段が、さらに、前記表示画像デー
タを構成する画素データの有する濃度データを、所定の
範囲に拡張して、コントラストを強調するコントラスト
強調手段と、前記表示画像データから選択された所定の
領域に含まれる画素データの有する濃度データの平均値
を算出し、前記平均値より大きな値を有する濃度データ
に、濃度データの最大値を付与する濃度変換手段を有
し、前記濃度変換手段により濃度データの変換されたデ
ータが、前記孤立点除去手段に与えられるように構成さ
れている。
【0017】本発明のさらに好ましい実施態様において
は、画像データが、蓄積性蛍光体シートを用いて生成さ
れている。本発明のさらに好ましい実施態様において
は、前記画像データが、オートラジオグラフィ画像デー
タおよび化学発光画像データよりなる群から選ばれる画
像データにより構成されている。本発明のさらに好まし
い実施態様においては、オートラジオグラフィ画像デー
タが、試料から発せられる放射線または電子線を、輝尽
性蛍光体に蓄積、吸収させ、しかる後に、前記輝尽性蛍
光体に、電磁波を照射して、該輝尽性蛍光体から発せら
れた光を光電変換することにより得られている。本発明
のさらに好ましい実施態様においては、化学発光画像デ
ータが、試料から発せられる可視光を、輝尽性蛍光体に
蓄積、吸収させ、しかる後に、前記輝尽性蛍光体に、電
磁波を照射して、該輝尽性蛍光体から発せられた光を光
電変換することにより得られている。本発明において、
位置情報とは、試料中における放射性標識物質もしくは
その集合体の位置を中心とした各種の情報、たとえば、
試料中に存在する放射性標識物質の集合体の存在位置と
形状、その位置における放射性標識物質の濃度、分布な
どからなる情報の一つもしくは任意の組み合わせとして
得られる各種の情報を意味するものである。
【0018】本発明において、オートラジオグラフィ画
像を生成するために使用することのできる輝尽性蛍光体
としては、放射線のエネルギーを蓄積可能で、電磁波に
よって励起され、蓄積している放射線のエネルギーを光
の形で放出可能なものであればよく、とくに限定される
ものではないが、可視光波長域の光によって励起可能で
あるものが好ましい。具体的には、たとえば、特開昭5
5−12145号公報に開示されたアルカリ土類金属弗
化ハロゲン化物系蛍光体(Ba1-x,2+ x )FX:yA
(ここに、M2+はMg、Ca、Sr、ZnおよびCdか
らなる群より選ばれる少なくとも一種のアルカリ土類金
属元素、XはCl、BrおよびIからなる群より選ばれ
る少なくとも一種のハロゲン、AはEu、Tb、Ce、
Tm、Dy、Pr、He、Nd、YbおよびErからな
る群より選ばれる少なくとも一種の3価金属元素、xは
0≦x≦0.6、yは0≦y≦0.2である。)、特開
平2−276997号公報に開示されたアルカリ土類金
属弗化ハロゲン化物系蛍光体SrFX:Z(ここに、X
はCl、BrおよびIからなる群より選ばれる少なくと
も一種のハロゲン、ZはEuまたはCeである。)、特
開昭59−56479号公報に開示されたユーロピウム
付活複合ハロゲン物系蛍光体BaFX・xNaX’:a
Eu2+(ここに、XおよびX’はいずれも、Cl、Br
およびIからなる群より選ばれる少なくとも一種のハロ
ゲンであり、xは0<x≦2、aは0<a≦0.2であ
る。)、特開昭58−69281号公報に開示されたセ
リウム付活三価金属オキシハロゲン物系蛍光体であるM
OX:xCe(ここに、MはPr、Nd、Pm、Sm、
Eu、Tb、Dy、Ho、Er、Tm、YbおよびBi
からなる群より選ばれる少なくとも一種の三価金属元
素、XはBrおよびIのうちの一方あるいは双方、x
は、0<x<0.1である。)、特開昭60−1011
79号公報および同60−90288号公報に開示され
たセリウム付活希土類オキシハロゲン物系蛍光体である
LnOX:xCe(ここに、LnはY、La、Gdおよ
びLuからなる群より選ばれる少なくとも一種の希土類
元素、XはCl、BrおよびIからなる群より選ばれる
少なくとも一種のハロゲン、xは、0<x≦0.1であ
る。)および特開昭59−75200号公報に開示され
たユーロピウム付活複合ハロゲン物系蛍光体MIIFX・
aMI X’・bM'II '' 2 ・cMIII ''' 3 ・x
A:yEu2+(ここに、MIIはBa、SrおよびCaか
らなる群より選ばれる少なくとも一種のアルカリ土類金
属元素、MI はLi、Na、K、RbおよびCsからな
る群より選ばれる少なくとも一種のアルカリ金属元素、
M' IIはBeおよびMgからなる群より選ばれる少なく
とも一種の二価金属元素、MIII はAl、Ga、Inお
よびTlからなる群より選ばれる少なくとも一種の三価
金属元素、Aは少なくとも一種の金属酸化物、XはC
l、BrおよびIからなる群より選ばれる少なくとも一
種のハロゲン、X’、X''およびX''' はF、Cl、B
rおよびIからなる群より選ばれる少なくとも一種のハ
ロゲンであり、aは、0≦a≦2、bは、0≦b≦10
-2、cは、0≦c≦10-2で、かつ、a+b+c≧10
-2であり、xは、0<x≦0.5で、yは、0<y≦
0.2である。)が、好ましく使用し得る。
【0019】本発明において、化学発光画像を生成する
ために、使用することのできる輝尽性蛍光体としては、
可視光波長域の光のエネルギーを蓄積可能で、電磁波に
よって励起され、蓄積している可視光波長域の光のエネ
ルギーを光の形で放出可能なものであればよく、とくに
限定されるものではないが、可視光波長域の光によって
励起可能であるものが好ましい。具体的には、たとえ
ば、特開平4−232864号公報に開示された金属ハ
ロリン酸塩系蛍光体、希土類元素付活蛍光体、アルミン
酸塩系蛍光体、珪酸塩系蛍光体、フッ化物系蛍光体が、
好ましく使用し得る。
【0020】
【発明の作用】本発明によれば、表示画像データから、
生化学画像上の一方の座標方向の定量すべき領域を画定
するレーンを形成するための第1の分割座標データが生
成されるため、操作者による煩雑な操作を必要とするこ
となく、生化学画像上の定量解析すべき領域を生成する
ことができる。また、本発明の好ましい実施態様によれ
ば、生化学画像上の他方の座標方向の定量すべき領域を
画定するレーンを形成するための第2の分割座標データ
が生成されるため、操作者による煩雑な操作を必要とす
ることなく、レーンごとに、定量解析すべき領域を生成
することができる。また、本発明のさらに好ましい実施
態様によれば、各々のレーンに形成されたスポットが属
すべきグループを示すグループデータが生成されるた
め、操作者が煩雑な操作をして、一つのグループに属す
べきスポットを画定する必要がなくなり、迅速に、定量
処理および定量解析を実行することが可能になる。
【0021】
【発明の実施の形態】以下、添付図面に基づいて、本発
明の実施態様につき詳細に説明を加える。図1は、本発
明の実施態様にかかるオートラジオグラフィ画像解析装
置のための画像データを生成する画像読み取り装置の一
例を示す略斜視図である。図1において、蓄積性蛍光体
シート1には、試料(図示せず)に含まれる放射性標識
物質の位置情報が、放射線エネルギーの形で、蓄積され
ている。本実施態様においては、成分の異なる複数の種
の薬物を、実験動物に投与し、所定時間経過後に、その
尿、胆汁およびフンを採取して、これを、TLCプレー
ト上にクロマトグラフィー展開し、TLCプレートから
転写された画像が、蓄積性蛍光体シート1に蓄積記録さ
れている。したがって、このTLCプレートから転写さ
れた蓄積性蛍光体シート1には、複数の実験動物ごと
に、これから採取した尿、胆汁およびフンをクロマトグ
ラフィー展開した一次元的なレーンが形成され、各レー
ン内には、それぞれ、各個体の状態などに応じた位置
に、スポットが形成されている。こうして試料中の放射
性標識物質の位置情報が蓄積記録された蓄積性蛍光体シ
ート1を、レーザ光2により、走査して、励起し、輝尽
光を発生させる。
【0022】レーザ光2は、レーザ光源3により発生さ
れ、フィルタ4を通過することにより、レーザ光2によ
る励起によって蓄積性蛍光体シート1から発生する輝尽
光の波長領域に対応する波長領域の部分がカットされ
る。次いで、レーザ光2は、ビーム・エクスパンダ5に
より、そのビーム径が正確に調整され、ガルバノミラー
等の光偏向器6に入射する。光偏向器6によって偏向さ
れたレーザ光2は、fθレンズ7を介して、平面反射鏡
8により反射され、蓄積性蛍光体シート1上に、一次元
的に入射する。fθレンズ7は、蓄積性蛍光体シート1
上を、レーザ光2により走査するときに、つねに、均一
のビーム速度で、走査がなされることを保証するもので
ある。このようなレーザ光2による走査と同期して、蓄
積性蛍光体シート1は、図1において、矢印Aの方向に
移動され、その全面が、レーザ光2によって走査される
ようになっている。蓄積性蛍光体シート1は、レーザ光
2が照射されると、蓄積記録していた放射線エネルギー
に比例する光量の輝尽光を発光し、発光した輝尽光は、
導光性シート9に入射する。
【0023】導光性シート9は、その受光端部が直線状
をなし、蓄積性蛍光体シート1上の走査線に対向するよ
うに近接して配置され、また、その射出端部は、円環状
をなし、フォトマルチプライアなどの光電変換型の光検
出器10の受光面に接続されている。この導光性シート
9は、アクリル系合成樹脂などの透明な熱可塑性樹脂シ
ートを加工して作られており、受光端部から入射した光
が、その内面で、全反射を繰り返しながら、射出端部を
経て、光検出器10の受光面に伝達されるように、その
形状が定められている。したがって、レーザ光2の照射
に応じて、蓄積性蛍光体シート1から発光した輝尽光
は、導光性シート9に入射し、その内部で、全反射を繰
り返しながら、射出端部を経て、光検出器10によって
受光される。光検出器10の受光面には、蓄積性蛍光体
シート1から発光される輝尽光の波長領域の光のみを透
過し、レーザ光2の波長領域の光をカットするフィルタ
が貼着されており、光検出器10は、蓄積性蛍光体シー
ト1から発光された輝尽光のみを光電的に検出するよう
に構成されている。光検出器10によって光電的に検出
された輝尽光は、電気信号に変換され、所定の増幅率を
有する増幅器11によって、所定のレベルの電気信号に
増幅された後、A/D変換器12に入力される。電気信
号は、A/D変換器12において、信号変動幅に適した
スケールファクタで、ディジタル信号に変換され、ライ
ンバッファ13に入力される。ラインバッファ13は、
走査線1列分の画像データを一時的に記憶するものであ
り、以上のようにして、走査線1列分の画像データが記
憶されると、そのデータを、ラインバッファ13の容量
よりも大きな容量を有する送信バッファ14に出力し、
送信バッファ14は、所定の容量の画像データが記憶さ
れると、画像データを、オートラジオグラフィ画像解析
装置に出力するように構成されている。
【0024】図2は、本発明の実施態様にかかるオート
ラジオグラフィ画像解析装置および画像読み取り装置の
ブロックダイアグラムである。図2において、オートラ
ジオグラフィ画像解析装置30は、蓄積性蛍光体シート
1に蓄積記録され、画像読み取り装置20により読み取
られて、ディジタル信号に変換された試料に含まれる放
射性標識物質の位置情報を含む画像データを受理し、画
像データに、たとえば、濃度、色調、コントラストなど
が適正で、観察解析特性に優れた可視画像を再生し得る
ように信号処理を施すとともに、画像データの一部分に
所望の領域を画定し、該領域内の濃度などを算出する信
号処理手段40と、画像読み取り装置20から信号処理
手段40に入力された画像データを一時的に記憶すると
ともに、信号処理後の画像データを記憶する画像データ
記憶手段50と、試料に含まれる放射性標識物質の位置
情報を含む画像データを画像として再生するCRT60
と、操作者により操作され、信号処理手段40に所定の
指示を与えるための入力データが入力される入力装置7
0とを備えている。オートラジオグラフィ画像解析装置
30は、受信バッファ41と制御部42とを備えてお
り、画像読み取り装置20の送信バッファ14に、一次
的に記憶された画像データは、オートラジオグラフィ画
像解析装置30の信号処理手段40の受信バッファ41
に入力されて、一時的に記憶される。この受信バッファ
41内に、所定量の画像データが記憶されると、制御部
42の指令に基づき、記憶された画像データが、画像デ
ータ記憶手段50の画像データ一時記憶部51に出力さ
れ、記憶される。このようにして、画像読み取り装置2
0の送信バッファ14から、信号処理手段40の受信バ
ッファ41に送られ、一時的に記憶された画像データ
は、さらに、受信バッファ41から、画像データ記憶手
段50の画像データ一時記憶部51に記憶される。こう
して、蓄積性蛍光体シート1の全面を、レーザ光2によ
って走査して得られた画像データが、画像データ記憶手
段50の画像データ一時記憶部51に記憶されると、信
号処理手段40の制御部42は、画像データ一時記憶部
51から所定の画像データを読み出し、画像データ記憶
部52に記憶する。
【0025】さらに、オートラジオグラフィ画像解析装
置30の信号処理手段40は、画像データ記憶部52に
記憶された画像データを読み出し、これに所定の処理を
施して、CRT60に出力するととも、画像データの定
量解析をおこなう画像生成/解析部43を備えている。
図3は、本実施態様にかかるオートラジオグラフィ画像
解析装置30の画像生成/解析部43およびその周辺回
路を詳細に示したブロック図である。図3に示されるよ
うに、画像生成/解析部43は、操作者が入力装置70
を操作することにより、画像データ記憶部52から読み
出された画像データに所定の処理を施して、CRT60
に表示すべき画像に対応する表示画像データを生成する
表示画像生成部102と、表示画像生成部102により
生成された表示画像データから背景部分の雑音に対応す
るデータを除去する背景雑音処理部104と、表示画像
データが展開される平面座標のx軸方向に、表示画像デ
ータを分割して、x軸方向の定量すべき領域を画定する
レーンを形成するレーン分割処理部106と、レーン分
割処理部106により形成されたレーンごとに、y軸方
向の定量すべき領域を画定するスポットを認識して、ス
ポット領域データを生成するスポット認識処理部108
と、入力装置70に与えられた入力データにしたがっ
て、スポット認識処理部108によって生成されたスポ
ット領域データを修正するデータ修正部110と、スポ
ット領域データに基づいて、スポットの各々に含まれる
画素の濃度データを算出する定量処理部112と、定量
処理部112により生成された濃度データなどを記憶す
る定量データ記憶部114と、定量データ記憶部114
に記憶された濃度データなどに基づいて、CRT60に
表示すべき表に対応する表データを生成する表データ生
成部116と、表示画像生成部102によって生成され
た表示画像データ、レーン分割処理部106により形成
されたレーンおよびスポット認識処理部108により生
成されスポット領域データ、ならびに、表データ生成部
116により生成された表データを合成する画像合成部
118と、画像合成部118により合成された画像デー
タを展開するウィンドウメモリ120と、スポット領域
データに基づいて、各レーンごとに、そのレーンに含ま
れる各スポットを、属すべきグループに分類するグルー
プデータを生成する代謝物分類処理部122とを備えて
いる。
【0026】図4は、本実施態様にかかる画像読み取り
装置20により読み取られ、画像データ記憶部52に記
憶された画像データの中から、入力装置70から入力さ
れたデータにしたがって、表示画像生成部102により
生成された表示画像データに基づき、CRT60の画面
に表示される画像を示す図である。図4に示されるよう
に、CRT60の画面400には、TLCプレート上に
滴下された10の試料をクロマトグラフィー展開するこ
とによって形成され、各々が複数のスポットを含んだレ
ーン401−1ないし401−10が形成されている。
ここに、レーン401−1、401−2および401−
3は、それぞれ、オスの実験動物に薬物を投与し、4時
間後、8時間後および24時間後に、その胆汁を採取し
て、これをTLCプレート上にクロマトグラフ展開した
レーンに対応し、レーン401−4は、オスの実験動物
に薬物を投与してから、24時間後に採取した糞をクロ
マトグラフ展開したレーンに対応し、また、レーン40
1−5は、メスの実験動物に薬物を投与してから、24
時間後に採取した糞をクロマトグラフ展開したレーンに
対応する。さらに、レーン401−6および401−7
は、オスの実験動物に薬物を投与してから、それぞれ、
4時間後および24時間後に採取した尿をクロマトグラ
フ展開したレーンに対応し、レーン401−8および4
01−9は、メスの実験動物に薬物を投与してから、そ
れぞれ、4時間後および8時間後に採取した尿をクロマ
トグラフ展開したレーンに対応する。レーン401−1
0は、実験動物に投与した薬物自体を、すなわち、薬物
を実験動物に投与することなく、クロマトグラフィー展
開することにより形成されたレーンに対応するレファレ
ンスである。
【0027】また、蓄積性蛍光体シート1を露光する際
に、宇宙線もしくは地中に含まれる放射線、または、た
とえば、クロマトグラフ展開するためのTLCプレート
が有する固有の放射線により、ほぼ均一に蓄積性蛍光体
シート1上に生じるノイズ成分に対応する背景雑音が、
画面400のスポットが形成されていない部分に現れて
いる。表示画像データは、複数の画素から構成されてお
り、各画素ごとに、その濃度を示す濃度データの値を有
している。本実施態様においては、濃度データの値の範
囲は、0ないし255であり、濃度データの値が0のと
きに、CRT60の画面上で、その画素は黒で表わさ
れ、濃度データの値が255のときに、その画素は白で
表わされるようになっている。図5は、背景雑音処理部
104の構成を示すブロックダイアグラムである。図5
に示されるように、背景雑音処理部104は、表示画像
データを構成する画素データの有する濃度データの偏り
を補正して、コントラスト強調画像データを得るための
コントラスト強調処理部130と、背景部分に対応する
画像データを構成する画素データの有する濃度データの
値の範囲を決定する濃度変換パラメータを算出する濃度
変換パラメータ算出部132と、各画素データの有する
濃度データの値を変換する濃度変換処理部134と、変
換された濃度データの値が所定範囲の画素データ中で、
孤立している画素データを捜し出して、その画素データ
の有する濃度データを所定の値に変更することにより、
濃度データの値が孤立している画素データを除去する濃
度孤立点除去処理部136と、濃度孤立点除去処理部1
36によって得られた画像データを二値化する二値化処
理部140と、二値化された画像データを構成する所定
範囲の画素データ中で、さらに、濃度データの値が孤立
している画素データを求め、その画素データの有する濃
度データの値を変更して、これを除去する孤立点除去処
理部142と、孤立点除去処理部142によって得られ
た二値化画像データの背景部分以外の領域の画素データ
に対応する画像データに、もとの濃度データを与える画
像間演算部144を備えている。
【0028】表示画像生成部102から出力された表示
画像データが、コントラスト強調処理部130に入力さ
れると、コントラスト強調処理部130により、この表
示画像データを構成する各画素データの有する濃度デー
タの値の範囲が求められる。次いで、コントラスト強調
処理部130は、こうして求められた濃度データの値の
範囲を、0ないし255の範囲に拡張するように、各画
素データの有する濃度データの値を変換する。すなわ
ち、表示画像データを構成する各画素データの有する濃
度データの最小値をa、最大値をbとした場合に、濃度
データの値がcであった画素データの有する濃度データ
は、次式にしたがって、濃度データc’に変換される。 c’=(c−a)×255/(b−a) このように、各画素データの有する濃度データの値を変
換することによって、TLCプレートを転写した蓄積性
蛍光体シートの露光条件や、蓄積性蛍光体シートの特性
などにより各画素データの有する濃度データの値の範囲
に偏りがある場合に、表示画像データのコントラストを
所望のように変換することができる。コントラスト強調
処理部130は、こうして、濃度データを変換して得た
コントラスト強調画像データを濃度変換パラメータ算出
部132および濃度変換処理部134に出力する。
【0029】濃度変換パラメータ算出部132は、コン
トラスト強調処理部130から入力されたコントラスト
強調画像データを構成する各画素データの有する濃度デ
ータの値のヒストグラムを生成し、このヒストグラムに
基づき、背景部分に対応する画像データを構成する画素
データの有する濃度データの値の範囲を決定する。図6
は、濃度変換パラメータ算出部132により生成された
濃度ヒストグラムを示す図である。濃度変換パラメータ
算出部132は、ある濃度データの値を有する画素デー
タの数を、濃度データの値ごとに算出して、図6に示さ
れるような濃度ヒストグラムを生成し、表示画像データ
を構成する画素データの総和のs%(sは、100より
小さな正の数である。)を占める領域601に関連する
濃度データの値dを、濃度パラメータとして、算出し、
濃度変換処理部134に出力される。この濃度パラメー
タを決定するための値sは、入力装置70を用いて任意
に設定することができる。濃度変換処理部134は、コ
ントラスト強調画像データを構成する画素データの有す
る濃度データの値eが、濃度パラメータdよりも大きい
場合、すなわち、この画素データが、CRT60の画面
上に、濃度パラメータdに対応する色よりも白く表示さ
れる場合に、濃度データの値を、255に変換する。す
なわち、濃度変換処理部134は、その入力された濃度
データの値eに基づき、以下のようにして、濃度変換画
像データを生成して、濃度孤立点除去処理部136に出
力する。
【0030】 濃度変換画像データ = 255 (d<e≦255) 濃度変換画像データ = e (0≦e≦d) このようにして生成された濃度画像変換画像データは、
たとえば、図4の領域402に示されるように、背景部
分、すなわち、スポットが形成されていない領域に対応
する画素であるにもかかわらず、濃度パラメータd以下
の濃度データの値を有する画素データ、すなわち、CR
T60上の画面上に部分的に黒く表示されるような濃度
データの値を有する画素データを、依然として含んでい
る。したがって、濃度孤立点除去処理部136は、濃度
変換画像データを構成する画素データから、(2m+1)×(2
m+1)個の矩形の画素データの領域を画定し、その領域に
含まれる所定の画素データの有する濃度データの値を、
矩形に含まれる画素データの濃度データの低い方からf
番目の値に変換する。本実施態様においては、濃度孤立
点除去処理部136は、mを1として、3×3個の矩形
の画素データの領域を画定し、画定された領域の中央に
存在する画素データの有する濃度データの値を、これら
の矩形領域に含まれる画素データの有する濃度データの
値のうち、低い方から3番目の値に変換するように構成
されている。したがって、ある(2m+1)×(2m+1)の画素デ
ータから構成される濃度変換画像データの各画素データ
の有する濃度データが、図7(a)に示すような値を有
している場合は、濃度孤立点除去処理部136は、図7
(b)に示すように、その中央に位置する画素データの
有する濃度データを、3番目に低い値、すなわち、25
5に変換する。なお、前述した矩形の領域を構成する画
素データの数に関するパラメータm、および、濃度デー
タの値を変換するために用いられるパラメータfは、入
力装置70によって任意に設定することが可能である。
このような処理をすることにより、背景部分であるにも
かかわらず、CRT60の画面上に黒く表示される画素
データを除去すること、すなわち、その画素データの有
する濃度データの値を変更して、CRT60の画面上に
表示するのに適した画像データを得ることが可能にな
る。
【0031】濃度孤立点除去処理部136からの出力デ
ータは、二値化処理部140に与えられる。二値化処理
部140は、濃度データの値が255である領域以外の
領域は、すべてスポットが形成されている領域と判断
し、その領域に属する画素データが有する濃度データの
値を0にする。すなわち、二値化処理部140は、入力
された濃度データの値gから、次式にしたがって、出力
濃度データ値を生成し、孤立点除去処理部142に出力
する。 出力濃度データ値 = 255 (g=255) 出力濃度データ値 = 0 (g≠255) 孤立点除去処理部142は、二値化処理部140から入
力された画像データを構成する画素データから、n×n
個の矩形の画素データの領域を画定し、画定された領域
内に、濃度データの値が0である画素データの数が、h
個以下であるときには、その領域内の画素データが有す
る濃度データの値をすべて255に変換する。本実施態
様においては、孤立点除去処理部142は、n=3とし
て、3×3個の矩形の画素データの領域を画定し、画定
された領域内に濃度データの値が0である画素データの
数が、2個以下である場合に、濃度データの値を255
に変換するように構成されている。したがって、あるn
×nの画素データから構成される濃度変換画像データの
各画素データが有する濃度データが、図8(a)に示す
ような値を有している場合に、孤立点除去処理部142
は、各画素データが有する濃度データの値をすべて、図
8(b)に示されるように、255に変換する。なお、
臨界画素データ数を示すパラメータhは、入力装置70
により任意に設定することが可能である。このような処
理をすることにより、背景部分であるにもかかわらず、
CRT60の画面上に黒く表示される孤立点の画素デー
タを、さらに、除去することが可能になる。
【0032】孤立点除去処理部142により、孤立点が
除去された二値化データは、画像間演算部144に与え
られる。画像間演算部144は、二値化データおよび表
示画像生成部102から入力された表示画像データに基
づいて、二値化データの値が255であった場合には、
その画素データの濃度データの値を255に保持し、他
方、二値化データの値が0であった場合には、その画素
データに、表示画像データの対応する濃度データを与
え、背景雑音の除去された雑音除去画像データを生成す
る。すなわち、二値化データの値が255であった画素
データに対応する雑音除去画像データは255となり、
二値化データの値が0であった画素データに対応する雑
音除去画像データは、表示画像データに一致する値とな
る。このように生成された雑音除去画像データは、レー
ン分割処理部106に送られる。図9は、x軸方向の定
量すべき領域を示すレーンを画定するレーン分割処理部
106の構成を示すブロックダイアグラムである。図9
に示されるように、レーン分割処理部106は、表示画
像データが展開される平面座標の各x座標に含まれる画
素データの有する濃度データの値の和、すなわち、x座
標の値を共通にする画素データのy軸方向の濃度データ
の和を、それぞれ算出し、領域のx軸方向のデータ値の
周辺分布を示す周辺分布データを生成するx軸方向周辺
分布算出部150と、x軸方向周辺分布算出部150に
より生成された周辺分布データに基づき、複数のスポッ
トが形成された定量すべきレーンを画定するレーンの開
始点を検出するレーン開始点検出部152と、周辺分布
データに基づき、レーンの終了点を検出するレーン終了
点検出部154と、レーン開始点検出部152およびレ
ーン終了点検出部154により検出されたデータに基づ
いて、スポットの形成されたレーンを画定するレーン分
割座標データを生成するレーン分割座標算出部156と
を備えている。
【0033】CRT60の画面400上に、図10
(a)に示されるような画像が表示されていた場合に、
x軸方向周辺分布算出部150は、この画像に対応する
雑音除去画像データ中のx座標の値を共通にする画素デ
ータの濃度データの和を、図10(b)に示すように算
出する。次いで、x軸方向周辺分布算出部150は、算
出された濃度データの和に含まれる高周波成分を除去す
るために、この濃度データの和に、ローパスフィルタ処
理を施す。さらに、x軸方向周辺分布算出部150は、
ローパスフィルタ処理が施され、高周波数成分が除去さ
れたデータの値と、予め定められたしきい値iとを比較
し、データ値が、しきい値iよりも大きいときは、これ
をしきい値iに置き換えて、x軸方向周辺分布データを
生成する。こうして生成されたx軸方向周辺分布データ
は、レーン開始点算出部152およびレーン終了点算出
部154に出力される。レーン開始点算出部152は、
x軸方向周辺分布データに対して、さらに、k点平滑微
分を施す。本実施態様においては、次式に示されるよう
に、7点平滑微分が施され、平滑化されたデータが生成
される。 b(n) =((−2)×a(n-2) + (−1)×a(n-1) +0
×a(n)+1×a(n+1) +2×a(n+2))/5 上式において、a(n) は、あるx座標nに対するx軸方
向周辺分布データの値であり、b(n) は、あるx座標n
に対する平滑化されたデータの値である。
【0034】次いで、レーン開始点算出部152は、平
滑化されたデータの値が減少し始める点のx座標を算出
する。このようにして、図11に示されるように、下降
開始点S(1)ないしS(11) が求められる。レーン終了点算
出部154は、レーン開始点算出部152と同様に、x
軸方向周辺分布データに対して、k点平滑微分を施し
て、平滑化されたデータの値の増大が終了する点、すな
わち、データの値が最高値に到達した点のx座標を算出
する。このようにして、図11に示すように、増大終了
点E(1)ないしE(11) が求められる。レーン開始点算出部
152およびレーン終了点算出部154により算出され
た減少開始点S(1)ないしS(11) および増大終了点E(1)な
いしE(11) に対応するx座標データは、レーン分割座標
算出部156に出力される。レーン分割座標算出部15
6は、これらデータに基づいて、それぞれ、レーン分割
座標データを生成する。ここに、レーンLnは、レーン
分割座標データS(n)およびE(n+1)によって画定されるよ
うになっている。すなわち、図11に示す場合には、た
とえば、レーンL1は、レーン分割座標データS(1)およ
びE(2)により、レーンL2は、レーン分割座標データS
(2)およびE(3)により、さらに、レーンL10は、レー
ン分割座標データS(10) およびE(11) により画定され
る。
【0035】レーン分割座標算出部156によりレーン
分割座標データが生成され、x軸方向の定量すべき領域
を示すレーンが画定されると、レーン分割座標算出部1
56は、レーン分割座標データを、画像合成部118に
出力する。画像合成部118は、表示画像データ生成部
102から入力された表示画像データとレーン分割座標
算出部156から入力されたレーン分割座標データとを
合成し、この合成されたデータを、ウィンドウメモリ1
20を介して、CRT60に出力する。したがって、C
RT60の画面400上には、図12に示されるような
画像が表示される。同時に、レーン分割座標算出部15
6は、レーン分割座標データを、スポット認識処理部1
08に出力する。図13は、スポット認識処理部108
の構成を示すブロックダイアグラムである。図13に示
されるように、スポット認識処理部108は、画定され
たレーンごとに、表示画像データが展開される平面座標
のy座標の値を共通にする画素データの有するx軸方向
の濃度データの和を、それぞれ算出し、y軸方向の濃度
データの値の周辺分布を示すy軸方向周辺分布データを
算出するy軸方向周辺分布算出部160と、y軸方向周
辺分布データから、y軸方向に定量すべき領域を画定す
るスポットが形成された領域の開始点と終了点とを検出
する領域開始点/終了点検出部162と、検出された領
域内に存在する極大点を検出する領域内ピーク値検出部
164と、開始点および終了点ならびに極大点に基づい
て、スポットが形成された領域の座標を算出するスポッ
ト認識座標算出部166とを備えている。
【0036】y軸方向周辺分布算出部160は、x軸方
向周辺分布算出部150が、雑音除去画像データに対し
ておこなったのと同様に、雑音除去画像データ中の各レ
ーン内で同一のy座標を有する画素データの有する濃度
データの値を累算し、この処理を、すべてのy座標につ
いて繰り返し、さらに、得られたデータに対して、ロー
パスフィルタ処理を施して、高周波数成分を除去するこ
とによって、y軸方向周辺分布データを生成する。図1
4は、図12に示すレーンL6のy軸方向の周辺分布デ
ータを示す図である。次いで、領域開始点/終了点検出
部162は、y軸方向周辺分布データの最大値のj%の
値を算出し、周辺分布データの値が、最大値×j%の値
よりも小さくなる点である領域開始点ASと、最大値×j
%の値よりも大きくなる点である領域終了点AEを算出す
る。たとえば、レーンL6の周辺分布データについて、
図14に示すように、AS(1) およびAE(1) が求められ
る。このパラメータjは、入力装置70により、任意に
設定することができる。本実施態様においては、jは9
5に設定されている。さらに、領域内ピーク値検出部1
64は、領域開始点ASと領域終了点AEとの間に含まれる
領域内に存在するy軸方向周辺分布データに対して、p
点平滑微分を施す。次いで、領域内ピーク値検出部16
4は、p点平滑微分が施されたデータに基づいて、この
領域内に存在する極大点を検出する。本実施態様におい
ては、7点平滑微分が施される。たとえば、レーンL6
のy軸方向周辺分布データについて、AS(1) およびAE
(1) が求められていた場合には、図14に示すように、
極大点P(1)ないしP(8)が検出される。
【0037】このように、領域開始点/終了点検出部1
62により得られた領域開始点ASおよび領域終了点AEに
対応するデータと、領域内ピーク値検出部164によっ
て得られた極大点AS(1) およびAE(1) に対応するデータ
は、スポット認識座標算出部166に与えられる。スポ
ット認識座標算出部166は、まず、ある領域開始点と
これに隣接する極大点との間、隣接する極大点の間、あ
るいは、極大点とこれに隣接する領域終了点との間のy
軸方向の長さを算出する。この長さが、入力装置70に
よって入力されたしきい値qよりも小さな場合には、ス
ポット認識座標算出部166は、この長さを有する領域
を、スポットとして認識しない。その一方、算出された
長さが、しきい値q以上である場合には、スポット認識
座標算出部166は、この長さを有する領域をスポット
として認識する。たとえば、図15に示されるようなy
軸方向周辺分布データから、領域開始点AS(1) および領
域終了点AE(1) と、これらの間の極大点P(1)およびP(2)
とが得られるとともに、領域開始点AS(2) および領域終
了点AE(2) と、これらの間の極大点P(1)およびP(2)とが
得られた場合、スポット認識座標算出部166は、これ
らの間のy軸方向の長さl1、l2、l3、l4、l5およびl6
算出する。本実施態様においては、しきい値qを20に
設定している。したがって、l1=11、l2=12、l3
16、l4=110、l5=100およびl6=80という値
が算出された場合には、領域開始点AS(1) と極大点P(1)
との間の長さ、極大点P(1)とP(2)との間の長さ、およ
び、極大点p(2)と領域終了点AE(1) との間の長さは、そ
れぞれ、しきい値qの値より小さいから、領域開始点AS
(1) と領域終了点AE(1) との間の領域は、領域内に極大
点P(1)およびP(2)が含まれているにもかかわらず、一つ
のスポットSP(1) としてのみ認識される。これに対し
て、領域開始点AS(2) と極大点P(1)との間の長さ、極大
点P(1)とP(2)との間の長さ、および、極大点p(2)と領域
終了点AE(2) との間の長さは、それぞれ、しきい値qの
値以上であるから、これらは、スポットSP(2) 、SP(3)
およびSP(4) として認識される。スポット認識座標算出
部166は、このように画定したスポットの両端のy座
標の値から構成されるスポット領域データを生成し、画
像合成部118に出力する。
【0038】画像合成部118は、表示画像データ生成
部102から入力された表示画像データと、レーン分割
座標算出部156から入力されたレーン分割座標データ
と、スポット認識座標算出部166から入力されたスポ
ット領域データとを合成し、この合成されたデータをウ
ィンドウメモリ120を介して、CRT60に出力す
る。したがって、CRT60の画面400上には、図1
6に示すような画像が表示される。また、データ修正部
110に、入力装置70から、操作者により入力された
データが与えられた場合には、データ修正部110は、
入力データにしたがって、スポット認識処理部108に
より生成されたスポット領域データを削除し、スポット
の両端のy座標の値の一方または双方を変更し、或い
は、スポットの両端のy座標を新たに画定してスポット
領域データを生成する。データ修正部110により、ス
ポット領域データが削除され、変更され、或いは新たに
生成された場合には、これらスポット領域データの修正
に応じた画像が、CRT60の画面上に表示される。さ
らに、スポット領域データは、スポット認識座標算出部
166から、定量処理部112および代謝物分類処理部
122に与えられる。定量処理部112は、レーン分割
座標データにより画定されたレーンごとに、スポット領
域データにより画定されたスポットに含まれる画素デー
タの有する濃度データ、すなわち、蓄積性蛍光体シート
1からの発光量(PSL値)を、それぞれ算出して、ス
ポットごとの濃度データを、定量データ記憶部114の
定量データベースの所定の領域に記憶する。
【0039】代謝物分類処理部122は、すべてのレー
ンL1〜L9について、スポット領域データによって画
定されたスポットの中心座標を示す中心座標データを生
成する。第iレーンの第j番目のスポットに対応するス
ポット領域データを構成する二つのy座標データをys
(i,j) 、ye(i,j) とすると、中心座標データyc(i,j) の
値は、yc(i,j) ={ys(i,j) −ye(i,j) }/2と表わす
ことができる。このように、すべてのレーンのすべての
スポットの中心座標データを算出した後に、代謝物分類
処理部122は、これら中で、最も小さい中心座標デー
タを抽出して、この値からr1 の範囲に含まれる他のレ
ーンに存在するスポットに対応する中心座標データを検
索する。次いで、代謝物分類処理部122は、この範囲
に含まれるスポットを同一の代謝物と判断し、これらの
スポットを、第1のグループとして分類し、グループデ
ータを生成する。さらに、代謝物分類処理部122は、
第1のグループに属さないスポットの中心座標データの
うち、その値の最も小さな中心座標データを抽出し、こ
の値からr2 の範囲に含まれる他のレーンに存在するス
ポットに対応する中心座標データを検索する。代謝物分
類処理部122は、この範囲に含まれるスポットを同一
の代謝物と判断し、これらのスポットを、第2のグルー
プとして分類し、グループデータを生成する。このよう
な処理を繰り返すことにより、第1のグループないし第
tのグループが形成され、これらに属するスポットを示
すグループデータが形成される。上述したように、所定
の中心座標データの値から所定の範囲に含まれる中心座
標データを検索することにより、TLCプレートが均一
でないことなどによって、展開溶媒の展開速度がレーン
により異なり、その結果、同一物質であるにもかかわら
ず、レーンによってスポットの位置がy軸方向にずれる
場合であっても、これらスポットを同一のグループに属
するように分類し、グループデータを生成することがで
きる。
【0040】図17は、CRT60の画面上に表示され
た画像を用いて、代謝物分類処理部122によりなされ
る処理を説明するための図である。図17に示すよう
に、4つのレーンL1ないしL4が形成され、各レーン
に、複数のスポットが形成されている場合に、代謝物分
類処理部122は、スポットの中心座標データの中で、
値の最も小さい中心座標データであるyc(3,1) から、r
1 の範囲に含まれる値を有する中心座標データを検索す
る。図17に示される場合には、中心座標データyc(1,
1) 、yc(2,1) 、yc(4,1) の三つの値が、yc(3,1) か
ら、r1 の範囲に含まれる値を有する中心座標データに
該当する。したがって、これら中心座標データに対応す
るスポットSP(1,1) 、SP(2,1) 、SP(3,1) およびSP(4,
1) のy座標の値から構成されるスポット領域データ
が、第1のグループに属するものとして、分類されるこ
とになる。次いで、代謝物処理部122は、yc(1,1) 、
yc(2,1) 、yc(3,1) 及びyc(4,1)以外の中心座標データ
のうちで、最も値の小さい中心座標データを検索し、さ
らに、この中心座標データの値からr2 の範囲に含まれ
る値を有する中心座標データを検索する。図17に示さ
れる場合には、yc(4,2) が、最も値の小さい中心座標デ
ータに該当し、中心座標データyc(1,2) が、yc(4,2) か
らr2 の範囲に含まれる値を有している。したがって、
これら中心座標データに対応するスポットSP(1,2) およ
びSP(4,2) のy座標の値から構成されるスポット領域デ
ータが、第2のグループに属するものとして、分類され
ることになる。
【0041】同様に、yc(2,2) からr3 の範囲に含まれ
る値を有する中心座標データを検索すると、yc(3,2)
が、これに含まれる値を有する中心座標データとなる。
したがって、これら中心座標データに対応するスポット
SP(2,2) およびSP(3,2) のy座標の値から構成されるス
ポット領域データが、第3のグループに属するものとし
て、分類されることになる。同様にして、SP(1,3) 、SP
(2,3) およびSP(4,3) のy座標の値から構成されるスポ
ット領域データが第4のグループに属することになり、
SP(1,4) 、SP(2,4)、SP(3,3) およびSP(4,4) のy座標
から構成されるスポット領域データが第5のグループに
属するものとして、分類されることになる。ここに、r
1 ないしr5 は経験的に定められ、入力装置70に入力
される。こうして、代謝物分類処理部122により生成
されたグループデータは、定量データ記憶部114の定
量データベースに与えられる。図18は、図17に示さ
れる画像に基づき形成された定量データ記憶部114の
定量データベースのメモリマップを示す図である。図1
8に示されるように、定量データベースには、第iレー
ンの第j番目のスポットを画定するための座標データys
(i,j) およびye(i,j) と、そのスポットに含まれる画素
の濃度データyp(i,j) と、代謝物分類処理部122によ
り生成されたスポットSP(i,j) が属すべきグループのグ
ループナンバーを示すグループデータとを有している。
【0042】前述したように、スポット認識処理部10
8により形成され、あるいは、データ修正部110によ
り修正されたスポットSP(i,j) の座標データおよび濃度
データは、定量処理部112を介して、定量データ記憶
部114の定量データベースに与えられ、また、このス
ポットスポットSP(i,j) が属すべきグループのグループ
ナンバーは、グループデータとして代謝物分類処理部1
22により、定量データベースに与えられる。表データ
生成部116は、定量データ記憶部114の定量データ
ベースに記憶されたデータに基づき、同一のグループナ
ンバーを有するスポットの濃度データの総和を100と
した場合に、各スポットの濃度データの値の割合を生成
し、あるいは、同一のレーンに形成されたスポットの濃
度データの総和を100とした場合に、各スポットの濃
度データの値の割合を生成し、これら割合を示す表を生
成し、これを画像合成部118に出力する。以上のよう
に構成された本実施態様にかかるオートラジオグラフィ
画像解析装置30の動作を以下に説明する。まず、操作
者が、入力装置70を操作して、所定のデータを入力す
ることにより、画像データ記憶手段50の画像データ記
憶部52から、所定領域の画像データが読み出される。
【0043】次いで、画像データ記憶部52から読み出
された画像は、操作者が入力したデータにしたがって、
画像形成/解析部43の表示画像生成部102により拡
大または縮小され、さらに、所定の処理が施され、その
結果、得られた表示画像データが、画像合成部118お
よび背景雑音処理部104に与えられる。表示画像デー
タは、コントラスト強調処理部130により、表示画像
データを構成する各画素データの有する濃度データの値
が変換されて、そのコントラストが強調され、入力装置
70から与えられた濃度パラメータを決定するためのデ
ータsに基づいて、濃度変換パラメータ算出部132に
より得られた濃度パラメータdにしたがって、濃度変換
画像データが形成される。濃度変換処理部134により
得られた濃度変換画像データは、濃度孤立点除去処理部
136に与えられ、入力装置70から与えられたデータ
に基づいて得られた矩形の領域を構成する画素データ数
を画定するパラメータmと、濃度データの値を変換する
ために用いられるパラメータfとにしたがって、背景部
分、すなわち、スポットが形成されていないと考えられ
る領域であるにもかかわらず、CRT60の画面上に黒
く表示される画素データ、すなわち、白く表示されるべ
きであるのに、白の濃度データの値255以外の値を有
する画素データの有する濃度データの値を255にする
ことができ、このような領域に存在する孤立点を除去す
ることができる。
【0044】さらに、濃度孤立点除去処理部136から
の出力は、二値化処理部140により二値化され、二値
化されたデータは、入力装置70から与えられたデータ
に基づいて得られた矩形の領域を構成する画素データ数
を画定するパラメータnと、臨界画素データ数を示すパ
ラメータhとにしたがって、さらに、背景部分であるに
もかかわらず、白く表示されないような濃度データ値を
有する画素データの濃度データの値を255にすること
ができ、このような領域に存在する孤立点をさらに除去
することができる。次いで、画像間演算部144は、孤
立点除去処理部142からの出力である孤立点が除去さ
れた二値化データと、表示画像データとに基づいて、二
値化データの値が255である画素データに対しては、
その値を225に保持し、他方、二値化データの値が0
である画素データに対しては、この画素データに表示画
像データの対応する濃度データを与える。このようにし
て、背景雑音処理部104により、背景部分に存在する
画素データを白く表示することができ、スポットの形成
されている領域を、より明確に、CRT60の画面40
0上に表示することが可能になる。
【0045】こうして、背景雑音処理部104により生
成された雑音除去画像データは、背景雑音処理部104
から、レーン分割処理部106に出力される。レーン分
割処理部106のx軸方向周辺分布算出部150は、C
RT60の画面上に表示された画像に対応する雑音除去
画像データ中のx座標を共通にする画素データの有する
濃度データの和を生成し、この生成された濃度データの
和にローパスフィルタ処理を施し、これに含まれる高周
波数成分を除去する。さらに、x軸方向周辺分布算出部
150は、高周波数成分が除去された濃度データの和の
値と、しきい値iとを比較し、データ値がしきい値iよ
りも大きいときは、この値をしきい値iに置き換えて、
x軸方向周辺分布データを生成する。レーン開始点算出
部152は、x軸方向周辺分布算出部150から出力さ
れたx軸方向周辺分布データに対して、さらに、k点平
滑微分を施し、微分されたx軸方向周辺分布データの値
の減少が開始される点のx座標を算出し、他方、レーン
終了点算出部154は、x軸方向周辺分布データに、k
点平滑微分を施し、微分されたデータの値の増大が終了
する点、すなわち、データの値が最高値に到達する点の
x座標を算出する。前述したように、図10(a)に示
されるような画像が、CRT60の画面上に表示されて
いた場合には、図11に示すように、伝承開始点S(1)な
いしS(11) および増大終了点E(1)およびE(11) が、それ
ぞれ得られる。
【0046】次いで、レーン開始点算出部152および
レーン終了点算出部154により、それぞれ得られた減
少開始点および増大終了点に対応するx座標に基づい
て、レーン分割座標算出部156は、レーン分割座標デ
ータを算出する。レーン分割座標データは、減少開始点
S(n)および増大終了点E(n+1)に対応するデータから構成
され、これらにデータにより、n番目のレーンLnが画
定される。図10(a)に示されるような画像が、CR
T60の画面上に表示されていた場合には、図12に示
されるように、10個のレーンL1ないしL10が形成
される。レーン分割処理部106により生成されたレー
ン分割座標データは、スポット認識処理部108に与え
られる。スポット認識処理部108のy軸方向周辺分布
算出部160は、雑音除去画像データに基づき、各レー
ン内のy座標を共通にする画素データの有する濃度デー
タの値を累算する。たとえば、図12に示されるよう
に、10個のレーンL1ないしL10が形成されている
場合には、まず、レーンL1内のy座標を共通にする画
素データの有する濃度データの値を、すべてのy座標に
ついて累算する。レーンL2ないしL10に対しても、
同様な処理を施すことにより、各レーンごとのy軸方向
周辺分布データが生成される。
【0047】各レーンごとのy軸方向周辺分布データ
は、領域開始点/終了点検出部162に与えられる。領
域開始点/終了点検出部162は、入力装置70から与
えられたデータに基づいて得られたパラメータjにした
がって、各レーンごとのy軸方向周辺分布データの値
が、最大値×j%よりも小さくなる点である領域開始点
ASと、最大値×j%よりも大きくなる点である領域終了
点AEとを算出する。図12に示されるように、10のレ
ーンL1ないしL10が形成されている場合には、各レ
ーンごとに、少なくとも一組の領域開始点ASおよび領域
終了点AEが算出される。各レーンの領域開始点ASおよび
領域終了点AEに対応するデータは、領域内ピーク値検出
部164に与えられる。領域内ピーク値検出部164
は、各レーンごとに、領域開始点ASと領域終了点AEとの
間に含まれる領域内に存在するy軸方向周辺分布データ
に対して、p点平滑微分を施し、さらに、この微分され
たデータから、この領域内に存在する極大点を検出す
る。領域開始点/終了点検出部162により検出された
領域開始点ASおよび領域終了点AEに対応するデータなら
びに領域内ピーク点検出部164により検出された極大
点に対応するデータは、スポット認識座標算出部166
に与えられる。
【0048】スポット認識座標算出部166は、これら
データに基づき、各レーンごとに、画定されたスポット
の両端のy座標の値から構成されるスポット領域データ
を算出する。このスポット領域データは、画像合成部1
18に与えられ、たとえば、図16に示すような画像
が、CRT60の画面400上に表示される。また、入
力装置70に入力されたデータにより、データ修正部1
10が作動される場合には、このデータに基づいて、ス
ポット領域データが修正される。さらに、スポット領域
データは、定量処理部112および代謝物処理部122
に与えられる。定量処理部112は、レーンごとに、各
スポットに含まれる画素データの有する濃度データを算
出し、これを定量データ記憶部114の定量データベー
スに与える。たとえば、第1のレーンの第1のスポット
sp(1,1) に含まれる画素データの有する濃度データは、
図18に示す領域181に記憶され、第1レーンの第2
のスポットsp(1,2) に含まれる画素データの有する濃度
データは、領域182に記憶される。また、代謝物分類
処理部122は、各レーンに形成された各スポットが、
どのグループに所属するかを示すグループデータ生成し
て、これを定量データ記憶部114の定量データベース
に与える。第1のレーンの第1のスポットsp(1,1) が属
すべきグループを示すグループデータは、図18に示さ
れる領域191に記憶され、第1レーンの第2のスポッ
トsp(1,2) が属すべきグループを示すグループデータ
は、領域192に記憶される。したがって、同一の代謝
物から得られたオートラジオグラフィ画像のスポット
を、同一のグループに属するようにグループ化すること
ができる。また、パラメータ l1 、l2・・・・を入力装
置70に入力されたデータに基づいて変更することによ
り、クロマトグラフ展開した際に、種々の原因によりス
ポットが形成される位置に生じる誤差を吸収することが
可能になる。
【0049】さらに、操作者が入力装置70に所定のデ
ータを入力することにより、表データ生成部116は、
定量データ記憶部114の定量データベースに記憶され
たデータを読み出し、前述したように、入力されたデー
タにしたがって、同一のグループナンバーを有するスポ
ットの濃度データの総和を100とした場合に、各スポ
ットの濃度データの値の割合を生成し、あるいは、同一
のレーンに形成されたスポットの濃度データの総和を1
00とした場合に、各スポットの濃度データの値の割合
を生成し、これら割合を示す表データを生成する。たと
えば、定量データベースに、図18に示されるようなデ
ータが記憶されていた場合に、表データ生成部116
は、グループナンバー1に属するスポットSP(1,1) 、SP
(2,1) 、SP(3,1) などの濃度データyp(1,1) 、yp(2,1)
、yp(3,1) などの値の総和をとり、この値を100と
した場合の、各濃度データyp(1,1) 、yp(2,1) 、yp(3,
1) などの割合を示す表データを生成する。また、入力
されたデータにしたがって、他のグループナンバー2な
いし5に属するスポットの各濃度データの割合を示す表
データを生成することもできる。表データ生成部116
により生成された表データは、画像合成部118に与え
られ、画像合成部118からウィンドウメモリ120を
介して、この表データはCRT60に与えられて、CR
T60の画面上に、表データに対応する表が表示され
る。
【0050】本実施態様によれば、表示画像生成部10
2により生成された表示画像データに基づき、背景雑音
処理部104により、背景部分の雑音に対応するデータ
が除去されて、雑音除去画像データが生成され、レーン
分割処理部106により、雑音除去画像データに基づい
て、表示画像データが展開される平面座標のx軸方向
に、表示画像データを分割して、x軸方向の定量すべき
領域を画定するレーンが決定され、さらに、スポットの
認識処理部108により、レーンごとに、表示画像デー
タが展開される平面座標のy軸方向に、表示画像データ
を分割して、y軸方向の定量すべき領域を画定するスポ
ットが決定される。したがって、操作者が煩雑な操作を
する必要なく、定量すべき領域を画定することが可能と
なる。図19は、本発明の他の実施態様にかかるオート
ラジオグラフィ画像解析装置の画像生成/解析部のブロ
ックダイアグラムである。本実施態様においては、演算
時間を短縮するため、背景雑音を除去する処理はなされ
ておらず、前記実施態様における背景雑音処理部104
は設けられていない。図20は、x軸方向の定量すべき
領域を示すレーンを画定するレーン分割処理部206の
構成を示すブロックダイアグラムである。
【0051】図20に示されるように、レーン分割処理
部206は、表示画像データが展開される平面座標のx
座標の値を共通にする画素データの有するy軸方向の濃
度データの和を、それぞれ算出し、領域のx軸方向のデ
ータ値の周辺分布を示すx軸方向周辺分布データを生成
するx軸方向周辺分布算出部250を備えている。本実
施態様においては、x軸方向周辺分布算出部250は、
前記実施態様とは正負を反転した形で、x軸方向の周辺
分布データを算出するように構成されている。レーン分
割処理部206は、さらに、x軸方向周辺分布算出部2
50により算出されたx軸方向周辺分布データに基づ
き、x軸方向周辺分布データの濃度データのピーク値を
与えるピーク位置を算出するピーク位置算出部252
と、ピーク位置算出部252により算出されたピーク位
置に基づいて、複数のスポットが形成され、x軸方向の
定量すべき領域を画定するレーンの開始点および終了点
を検出するレーン開始点/終了点検出部254と、レー
ン開始点/終了点検出部254により得られたデータに
基づき、スポットの形成されたレーンを画定するレーン
分割座標データを生成するレーン分割座標算出部256
とを備えている。図21は、ピーク位置算出部252の
ブロックダイアグラムである。
【0052】図21に示されるように、ピーク位置算出
部252は、x軸方向周辺分布算出部250によって算
出されたx軸方向の周辺分布データのノイズを除去する
ノイズ除去部300と、ノイズ除去部300によりノイ
ズが除去されたx軸方向の周辺分布データに一次元フー
リエ変換を施して、ピークの位置を推定するピーク位置
推定部302と、ピーク位置推定部302により推定さ
れたピークの位置に基づき、ピークの位置を決定するピ
ーク位置決定部304とを備えている。ノイズ除去部3
00においては、x軸方向周辺分布算出部250により
算出された周辺分布データにローパスフイルタ処理が施
されて、高周波のノイズ成分が除去される。一般に、画
像データ中のノイズは、小さな濃度データのゆらぎとし
て現れることが多いため、このような濃度データのゆら
ぎを除去するため、次いで、ヒステリシス平滑化処理が
施される。図22は、ヒステリシス平滑化処理の方法を
示すものである。図22に示されるように、ヒステリシ
ス平滑化処理においては、入力されたx軸方向の周辺分
布データの波形500に沿って、高さhのカーソル50
2が移動され、カーソル502の中心点の軌跡として、
周辺分布データの出力波形504が得られる。具体的に
は、カーソル502を1ステップだけ、横軸方向に移動
させ、入力波形500とカーソル502とが交わってい
るか否かを判定する。その結果、交わっていれば、カー
ソル502を、さらに、1ステップだけ、横軸方向に移
動させる。他方、入力波形500とカーソル502とが
交わっていない場合には、カーソル502の上端部より
も入力波形500が上側にあるときは、カーソル502
の上端部が入力波形500と一致するように、カーソル
502を縦軸方向に移動させ、カーソル502の下端部
よりも入力波形500が下側にあるときは、カーソル5
02の下端部が入力波形500と一致するように、カー
ソル502を縦軸方向に移動させて、さらに、1ステッ
プだけ、横軸方向に移動させる操作を繰り返す。こうし
て、カーソル502の中心点の軌跡を求めて、出力波形
とすることにより、hよりも小さな入力波形500のゆ
らぎが除去されて、滑らかな出力波形504を得ること
ができる。こうして得られた出力波形504は、図22
に示されるように、ヒステリシス平滑化処理によって、
入力波形500に対して、位相がずれているため、ノイ
ズ除去部300は、さらに、出力波形504に対して、
逆方向から、ローパスフイルタ処理を施し、位相のずれ
を補正する。
【0053】こうして、ノイズ除去部300によりノイ
ズが除去されたx軸方向の周辺分布データは、ピーク位
置推定部302に入力される。ピーク位置推定部302
においては、x軸方向の周辺分布データに一次元フーリ
エ変換が施され、複数の正弦波が得られる。こうして得
られた複数の正弦波のうち、最も振幅が大きい正弦波
が、x軸方向の周辺分布データの主たる波形成分と考え
られるから、ピーク位置推定部302は、得られた複数
の正弦波の中から、振幅が最も大きく、かつ、操作者に
より、x軸方向の周辺分布データが有しているであろう
と予想され、入力装置70に入力された入力ピーク数以
上のピーク数を有する正弦波を選択する。ピーク位置推
定部302は、さらに、こうして選択された正弦波の位
相および周波数に基づき、図23に示されるように、正
弦波のピークの位置p(1)〜p(i)を求め、ピーク位置決定
部304に出力する。こうして求められた正弦波のピー
クの位置p(1)〜p(i)は、あくまでも、x軸方向の周辺分
布データの主たる波形成分と考えられる正弦波のピーク
の位置にすぎないから、x軸方向の周辺分布データの実
際のピークの位置とは、必ずしも一致してはいない。し
たがって、ピーク位置決定部304は、図24に示され
るように、入力された正弦波のピークの位置p(1)〜p(i)
に対応するピークの位置を中心として、所定の幅の範囲
のx軸方向の周辺分布データの面積を求め、面積が最大
となる位置をx軸方向の周辺分布データのピークの位置
と暫定的に決定する。
【0054】x軸方向の周辺分布データの波形によって
は、こうして決定されたピークの位置がきわめて近接し
ている場合があり、そのような場合は、一方は、ピーク
ではないにもかかわらず、ピークと認識されたと認めら
れるから、ピーク位置決定部304は、所定周期、たと
えば、選択された正弦波の周期の1/4の範囲内に、2
以上のピークが認識されているときは、位相の最も小さ
いピークのみをピークとして認識する。前述のように、
ピーク位置決定部304は、所定の幅の範囲のx軸方向
の周辺分布データの面積に基づき、x軸方向の周辺分布
データのピークを求めているため、こうして得られた複
数のピークの間に、濃度データの小さいピークが存在し
ているにもかかわらず、ピークとして認識されていない
場合がおこり得る。そこで、ピーク位置決定部304
は、図25に示されるように、左右のピークの位置か
ら、所定数の画素データ、たとえば、5画素データにわ
たり、連続的に濃度データが増大する部分があるか否か
を判定し、そのような部分があるときは、濃度データが
連続的に増大を開始した点の中点に、ピークがあると判
定している。こうして、ピーク位置算出部252により
算出されたx軸方向の周辺分布データのピーク値は、レ
ーン開始点/終了点検出部254に入力される。
【0055】レーン開始点/終了点検出部254は、ピ
ーク位置算出部252から入力されたx軸方向の周辺分
布データのピーク値に基づいて、ピーク間の中点を求
め、図26に示されるように、隣接する中点の間で、濃
度データが連続的に増大を開始する位置を、それぞれ、
レーン開始点およびレーン終了点と決定して、レーン分
割座標算出部256に出力する。レーン分割座標算出部
256は、レーン開始点/終了点検出部254により得
られたデータに基づき、前記実施態様と同様にして、ス
ポットの形成されたレーンを画定するレーン分割座標デ
ータを生成する。レーン分割座標算出部256によりレ
ーン分割座標データが生成され、x軸方向の定量すべき
領域を示すレーンが画定されると、レーン分割座標デー
タは、画像合成部118に出力される。画像合成部11
8は、表示画像データとレーン分割座標データとを合成
し、この合成されたデータを、ウィンドウメモリ120
を介して、CRT60に与える。したがって、CRT6
0の画面400上には、図12に示すような画像が表示
される。また、レーン分割座標データは、スポット認識
処理部208に与えられる。
【0056】図27は、スポット認識処理部208の構
成を示すブロックダイアグラムである。図27に示され
るように、スポット認識処理部208は、画定されたレ
ーンごとに、表示画像データが展開される平面座標のy
座標の値を共通にする画素データの有するx軸方向の濃
度データの和を、それぞれ算出し、y軸方向の濃度デー
タの値の周辺分布を示すy軸方向周辺分布データを算出
するy軸方向周辺分布算出部260と、y軸方向周辺分
布算出部260により算出されたy軸方向の周辺分布デ
ータのノイズ成分を除去するノイズ除去部262と、y
軸方向に定量すべき領域を画定するスポットが形成され
た領域の開始点と終了点とを検出する領域開始点/終了
点検出部262と、検出された領域内に存在する境界値
を検出する領域内極大値検出部264と、開始点および
終了点ならびに境界値に基づき、スポットが形成された
領域内の隣接するスポットの境界の座標を算出するスポ
ット境界座標算出部266と、スポット境界座標算出部
266によって算出された隣接するスポットの境界の座
標に基づき、スポットが形成された領域内の隣接するス
ポットの境界の座標を決定するスポット境界座標決定部
268を備えている。
【0057】次いで、領域開始点/終了点検出部262
は、ノイズ成分が除去されたy軸方向周辺分布データに
基づいて、前記実施態様の領域開始点/終了点検出部1
62と全く同様にして、領域開始点ASおよび領域終了点
AEを算出する。さらに、領域内極大値検出部264は、
領域開始点ASと領域終了点AEとの間に含まれる領域内に
存在する周辺分布データに対して、p点平滑微分を施
し、p点平滑微分が施されたデータから、この領域内に
存在する極大値を検出する。前記実施態様と同様に、本
実施態様においては、7点平滑微分が施される。こうし
て、領域開始点/終了点検出部262によって得られた
領域開始点ASおよび領域終了点AEに対応するデータと、
領域内極大値検出部264によって得られた極大値に対
応するデータは、それぞれ、スポット境界座標算出部2
66に与えられる。スポット認識処理部208において
は、前記実施態様と同様に、濃度データの値が0のとき
に、CRT60の画面上で、その画素データは黒で表わ
され、濃度データの値が255のときに、その画素デー
タは白を表わされるようになっており、したがって、領
域開始点ASおよび領域終了点AEおよび極大値を与える座
標値は、スポットの境界線の座標値に対応するものと推
定される。そこで、スポット境界座標算出部266は、
領域開始点/終了点検出部262によって得られた領域
開始点ASおよび領域終了点AEを与える座標値ならびに領
域内強大値検出部264によって得られた極大値を与え
る座標値を、隣接するスポットの境界に対応する境界座
標として検出して、スポット境界座標決定部268に出
力する。
【0058】しかしながら、こうして求めた境界座標
は、領域開始点ASおよび領域終了点AEおよび極大値を与
える座標値として求めたものにすぎず、その間に、スポ
ットが実際に存在するとは限らず、場合によっては、ノ
イズによる濃度データの微小な変化が存在するにすぎな
い場合があるので、ただちに、スポットの境界座標と決
定することはできない。そこで、スポット境界座標決定
部268は、隣接する境界に対応する境界座標の間の領
域内のy実願昭方向周辺分布データに対して、p点平滑
微分を施して、その領域内に極小値が存在するか否かを
判定する。その結果、極小値が検出されなかったとき
は、スポット境界座標決定部268は、その領域を画定
する境界座標を、スポットの境界座標としては認識しな
い。本実施態様においては、7点平滑微分が施される。
さらに、ノイズによる濃度データの微小な変化が存在す
るにすぎない場合に、これを、スポットとして誤って認
識することを防止するために、スポット境界座標決定部
268は、2以上の境界座標が、所定数の画素データ、
たとえば、5画素データ内に含まれるときは、1つの境
界座標のみをスポットの境界座標として認識し、その他
の境界座標はスポットの境界座標としては認識しない。
【0059】スポット境界座標決定部268は、こうし
て得られたスポットの境界座標に基づいて、各レーンご
とに、各スポットを画定する境界のy座標の値から構成
されるスポット領域データを生成し、画像合成部218
に出力する。その結果、前記実施態様と同様にして、た
とえば、図16に示すように、レーン分割線によりレー
ン分割がされ、スポット分割線によりスポット分割され
た画像がCRT60の画面400上に表示される。操作
者が、入力装置70を操作して、画像データを修正する
ためのデータが、データ修正部210に入力されたとき
は、前記実施態様と同様に、データ修正部210によ
り、画像データが修正され、修正された画像データに応
じた画像が、CRT60の画面上に表示される。さら
に、スポット境界座標決定部268は、スポット領域デ
ータを、定量処理部212および代謝物分類処理部22
2に出力する。定量処理部212は、レーン分割座標デ
ータにより画定されたレーンごとに、スポット領域デー
タにより画定されたスポットに含まれる画素データの有
する濃度データ、すなわち、蓄積性蛍光体シート1から
の発光量(PSL値)を、それぞれ、算出し、スポット
ごとの濃度データを、定量データ記憶部214の定量デ
ータベースの所定の領域に記憶する。
【0060】図28は、代謝物分類処理部222のブロ
ックダイアグラムである。図28に示されるように、代
謝物分類処理部222は、入力されたスポット領域デー
タに基づいて、レーンを2本づつ、比較して、スポット
分割線の位置の一致、不一致にしたがって、スポットを
グループ化し、グループ化がされる毎に、各グループ
に、試料を展開した側から順に、グループナンバーを付
す処理を実行するスポットグループ化部223と、濃度
データの分布が類似したレーンを同じグループに分類
し、グループ化するレーングループ化部224と、レー
ングループ化部224によるレーンのグループ化に基づ
いて、各スポットの属すべきグループを補正するスポッ
トグループ補正部225を備えている。スポット領域デ
ータを受けると、代謝物分類処理部122のスポットグ
ループ化部223は、図29(a)に示されるように、
2本のレーン内の対応すると認められるスポットを画定
するスポット分割線の位置を比較し、対応するスポット
分割線700および701、702および703の位置
が一致した場合には、これらのスポットを同一グループ
に分類する。これに対して、図29(b)に示されるよ
うに、対応するスポット分割線700および701、7
02および703の位置が、いずれも一致しないとき
は、これらのスポットを別のグループに分類する。他
方、図29(c)に示されるように、対応するスポット
を画定するスポット分割線のうち、一方のスポット分割
線の位置は一致したが、他方のスポット分割線の位置が
一致しない場合には、スポットグループ化部223は、
それぞれのスポット内のピークの位置を比較して、これ
が一致したときは、これらのスポットを同一グループに
分類する。さらに、図29(d)に示されるように、2
本のレーン内に、位置がほぼ一致するスポット分割線7
00および701、702および703の存在は認めら
れるが、一方のレーン内のスポット分割線701および
703の間に、さらに他のスポット分割線705が存在
する場合には、スポットグループ化部223は、一対の
スポット分割線700および702ならびに701およ
び703により画定されるスポット内のピークの位置を
比較して、これが一致したときは、これらのスポットを
同一グループに分類し、一致しないときは、別のグルー
プに分類する。
【0061】こうして、図12に示されるレーンL1
を、レーンL2ないしL9と順次比較し、さらに、レー
ンL2を、レーンL3ないしL9と、レーンL3を、レ
ーンL4ないしL9と順次比較するというようにして、
すべてのレーンを他のレーンと比較し、レーンL1ない
しL9に含まれる各スポットをグループ化する。こうし
てなされたスポットのグループ化の結果、1つのレーン
内の各スポットは、必ず、異なるグループに属すること
になるはずであるが、グループ化は、あるレーンの各ス
ポットを他のすべてのレーンのすべてのスポットと比較
して、あるスポットがどのスポットと同一グループに分
類されているか否かを判定することによりなされ、同じ
レーン内のスポットを比較していないため、1つのレー
ン内に含まれる2以上のスポットが、同一グループに分
類されていることが起こり得る。そこで、本実施態様に
おいては、スポットグループ化部223は、1つのレー
ン内の2以上のスポットが、同一グループに分類されて
いるときは、各スポットの面積を比較して、面積が最大
のスポットのみを、そのグループに残し、他のスポット
は、そのグループから外す処理をおこなう。その結果、
1つのレーン内の各スポットは、すべて、異なるグルー
プに分類されることになるが、あるスポットが、2以上
のグループに分類されている場合があり得る。そこで、
代謝物分類処理部222のスポットグループ化部223
は、各スポットを1つのグループに分類するため、2以
上のグループに分類されているスポットと同じグループ
に属している他のレーン内のスポットが、そのレーンの
レーンナンバーよりもレーンナンバーが小さいレーン内
に存在するか否かを、レーンナンバーの大きい順に検索
する。その結果、そのスポットと同じグループに分類さ
れているスポットが、異なるレーン内で検索されたとき
は、スポットグループ化部223は、そのスポットを、
レーンナンバーが最も大きいレーン内で検索されたスポ
ットと同じグループに分類する。他方、そのスポットが
含まれているレーンよりもレーンナンバーが小さい同一
レーン内に、そのスポットが分類されているのと同じグ
ループに属しているスポットが、2以上、検索されたと
きは、スポットグループ化部223は、そのスポットの
面積と、検索された他のレーン内のスポットの面積との
比を算出して、比較し、そのスポットを、面積比の大き
いスポットと同じグループに分類する。これに対して、
そのレーンよりもレーンナンバーの小さいレーン内に
は、そのスポットと同じグループに属するスポットが検
索されなかったときは、スポットグループ化部223
は、そのスポットを、グループナンバーの最も小さいグ
ループに分類する。
【0062】次いで、代謝物分類処理部222のレーン
グループ化部224は、L1ないしL9のレーンを比較
して、同じグループに属するスポットの数が多いレー
ン、すなわち、一致するスポット分割線の割合の多いレ
ーンを選択し、レーンのグループ化をおこなう。この
際、3以上のレーン内に、2つのスポット分割線が互い
に一致してはいるが、その間に、スポット分割線が存在
するスポットがある場合には、一致したスポット分割線
の間にあるスポットの面積の大きい方を同じグループの
レーンとしてグループ化する。こうして、レーンのグル
ープ化が完了すると、代謝物分類処理部222のスポッ
トグループ補正部225は、同じグループ内に分類され
たレーン間の対応する位置に、異なるグループに属する
スポットがあるときは、これらのスポットを同じグルー
プに分類する処理をおこない、スポットグループ化デー
タを書き換えを実行する。以上のようにして、代謝物分
類処理部222により生成されたグループデータは、定
量データ記憶部214の定量データベースに与えられ、
定量データ記憶部214の定量データベースには、図1
8に示されるのと同様なデータが記憶される。
【0063】前記実施態様と同様にして、表データ生成
部216によって、所望の表が生成され、画像合成図2
18に出力され、さらに、ウィンドメモリ220を介し
て、CRT60の画面上に表示される。本実施態様によ
れば、背景雑音を除去する処理がなされていないため、
小さな容量のメモリを用い、短時間に、操作者が煩雑な
操作をすることを要せずして、定量すべき領域を画定す
ることが可能になる。また、平滑化微分ではなく、ノイ
ズ除去部300によりノイズが除去されたx軸方向の周
辺分布データに一次元フーリエ変換を施して、ピークの
位置を推定しているので、小さなピークも精度良く認識
することができる。さらに、領域開始点ASおよび領域終
了点AEおよび極大値を与える座標値を、ただちに、スポ
ットの境界座標として認識せず、これらの間に、極小値
が存在するか否かおよび2以上の境界座標が、所定数の
画素に含まれるか否かを判定して、スポットの境界座標
を決定しているので、ノイズ成分により濃度が微小に変
化する部分をスポットとして誤って認識することを防止
することができ、スポットの認識精度を向上させること
ができる。また、前記実施態様においては、スポットの
中心座標に基づいて、スポットのグループ化をしている
ため、スポットの面積が大きく、濃度の高い部分がスポ
ットの一方の境界の近くに存在するようなスポットがあ
る場合には、2つのスポットの濃度の高い部分の座標は
ほぼ一致していても、濃度の高い部分の座標とスポット
の中心座標が一致しないため、2つのスポットは異なる
グループに分類されてしまうおそれがあった。しかしな
がら、本実施態様によれば、一方のスポット分割線の位
置は一致したが、他方のスポット分割線の位置が一致し
ない場合には、スポットグループ化部223は、それぞ
れのスポット内のピークの位置を比較して、これが一致
したときは、これらのスポットを同一グループにグルー
プ化するようにしているので、スポットの面積が大き
く、濃度の高い部分がスポットの一方の境界の近くに存
在するようなスポットも、所望のように、グループ化す
ることが可能になる。
【0064】図30は、本発明の他の実施態様にかかる
スポット認識処理部の構成を示すブロックダイアグラム
である。図30に示されるように、本実施態様にかかる
スポット認識処理部208は、画定されたレーンごと
に、表示画像データが展開される平面座標の各y座標に
含まれる画素の濃度データの和、すなわち、y座標の値
を共通にする画素のx軸方向の濃度データの和を、それ
ぞれ算出し、y軸方向の濃度データの値の周辺分布を示
す周辺分布データを算出するy軸方向周辺分布算出部8
60と、y軸方向周辺分布算出部860により算出され
たy軸方向の周辺分布データのノイズ成分を除去するノ
イズ除去部861と、ノイズ成分が除去されたy軸方向
の周辺分布データの極小値を検出する極小値検出部86
2と、ノイズ成分が除去されたy軸方向の周辺分布デー
タの極大値を検出する極大値検出部863と、極小値検
出部862によって検出された極小値に基づいて、スポ
ットを決定するためのしきい値を決定するしきい値決定
部864と、しきい値決定部864により決定されたし
きい値に基づき、スポットを決定するスポット決定部8
65と、スポット決定部865により決定されたスポッ
トに基づいて、スポット境界座標を決定するスポット境
界座標決定部868を備えている。
【0065】y軸方向周辺分布算出部860は、前記実
施態様におけるy軸方向周辺分布算出部160と全く同
様に、各レーン内で同一のy座標を有する画素に対応す
る画像データの濃度データの値を累算し、この処理を、
すべてのy座標について繰り返し、y軸方向の周辺分布
データを生成して、ノイズ除去部861に出力する。ノ
イズ除去部861は、入力されたy軸方向の周辺分布デ
ータに対して、ローパスフィルタ処理およびメディアン
・フィルタ処理を施すことによって、y軸方向周辺分布
データのノイズ成分を除去し、ノイズ成分が除去された
y軸方向の周辺分布データを、極小値検出部862に出
力する。極小値算出部862は、ノイズ成分が除去され
たy軸方向の周辺分布データに対して、p点平滑微分を
施し、p点平滑微分が施されたデータから、極小値を検
出する。前記実施態様と同様に、本実施態様において
も、7点平滑微分が施される。こうして検出された極大
値は、しきい値決定部864に出力される。しきい値決
定部864は、入力された極小値を比較して、値の小さ
い順に、s個の極小値を選択して、選択された極小値の
間のy軸方向周辺分布データの最大値を求める。図31
に示されるように、本実施態様においては、s=3に設
定され、3つの極小値V1、V2、V3が選択され、そ
の間でのy軸方向周辺分布データの最大値dmax が求め
られている。図31においては、極小値V1とV2との
間の極大値が選択されている。次いで、しきい値決定部
864は、検出された最大値のu%をしきい値として、
決定し、スポット決定部865に出力する。本実施態様
においては、u=80%に設定されている。
【0066】スポット決定部865は、しきい値決定部
864から入力されたしきい値よりも小さい値を有する
極小値をスポットとして決定し、スポット座標決定部8
68に出力する。他方、極大値検出部863は、ノイズ
成分が除去されたy軸方向の周辺分布データに対して、
p点平滑微分を施し、p点平滑微分が施されたデータか
ら、極大値を検出して、スポット境界座標決定部868
に出力する。前記実施態様と同様に、本実施態様におい
ても、7点平滑微分が施される。スポット境界座標決定
部868は、極大値検出部863から入力された極大値
のうち、スポット決定部865により決定されたスポッ
トを挟んで、隣接する極大値の座標を、スポットの境界
を確定するスポット境界座標として決定する。スポット
境界座標決定部868は、こうして得られたスポットの
境界座標に基づいて、各レーンごとに、各スポットを画
定する境界のy座標の値から構成されるスポット領域デ
ータを生成して、前記実施態様と全く同様にして、画像
合成部218、定量処理部212および代謝物分類処理
部222に出力する。第1の実施態様においては、y軸
方向周辺分布データの最大値を基準として、しきい値を
定めて、領域開始点ASと領域終了点AEとを決定し、その
間で、スポット分割線を求めているので、そのレーン内
のy軸方向周辺分布データの最大値がきわめて大きいと
き、すなわち、そのレーン内に、濃度がきわめて低い部
分が存在し、レーン内のスポット分割線の部分のy軸方
向周辺分布データの値が、最大値に比して、小さいとき
には、領域開始点ASと領域終了点AEとの間に、スポット
分割線を検出することができない場合がある。このよう
な問題は、しきい値を定めるjの値を小さく設定すれ
ば、ある程度、防止することは可能になるが、いかなる
画像に対しても、適切に、jの値を決定することは困難
である。しかしながら、本実施態様によれば、実際に、
スポットが存在する部分におけるy軸方向周辺分布デー
タの最大値、すなわち、濃度が低く、スポット分割線に
相当する部分のy軸方向周辺分布データの値を基準とし
て、スポット分割線を求めるしきい値を定めているの
で、いかなる画像に対しても、適切に、しきい値を決定
し、スポット分割線を欠点することが可能になる。
【0067】本発明は、以上の実施態様に限定されるこ
となく、特許請求の範囲に記載された発明の範囲内で種
々の変更が可能であり、それらも本発明の範囲内に包含
されるものであることがいうまでもない。たとえば、前
記実施態様においては、複数の試料を、それぞれ、クロ
マトグラフィー展開する薄層クロマトグラフィ法につ
き、説明を加えたが、本発明は、単一の試料を一次元的
に展開する薄層クロマトグラフィ法にも適用することも
できる。さらに、前記実施態様においては、実験動物か
ら採取した尿、胆汁およびフンをクロマトグラフィー展
開し、蓄積性蛍光体シートを用いて得られたオートラジ
オグラフィ画像の領域を画定する場合につき、説明を加
えたが、少なくとも一つの試料を一次元的に展開して得
られた試料中の位置情報を含む画像であれば、かかるオ
ートラジオグラフィ画像に限定されることなく、サザン
・ブロット・ハイブリダイゼーション法を利用した遺伝
子のオートラジオグラフィ画像や、ポリアクリルアミド
ゲル電気泳動法による蛋白質のオートラジオグラフィ画
像などの他のオートラジオグラフィ画像の領域を画定す
る場合はもとより、化学発光画像の領域を画定する場合
にも適用することができる。
【0068】また、第1の実施態様においては、レーン
分割がなされた画像データを、さらに、スポット認識座
標算出部166により、スポット領域データを生成して
いるが、スポット領域を画定することは、操作者が、C
RT60の画面400上のオートラジオグラフィ画像を
観察しながら、おこなう方が、正確になし得ることがあ
るので、レーン分割座標算出部156により、レーン分
割座標データが生成され、CRT60の画面400上の
オートラジオグラフィ画像が、自動的に、レーン分割さ
れれば、大幅に、オートラジオグラフィ画像の解析を迅
速化することが可能になり、その意味で、スポット認識
座標算出部166により、スポット領域データを生成す
ることは必ずしも必要でない。さらに、第1の実施態様
においては、背景雑音処理部104は、コントラスト強
調処理部130、濃度変換パラメータ算出部132、濃
度変換処理部134、濃度孤立点除去処理部136、二
値化処理部140、孤立点除去処理部142および画像
間演算部144を有し、これらによって、背景部分の雑
音に対応するデータを除去しているが、背景雑音処理部
104は、これに限定されるものではない。たとえば、
図30に示すように、背景雑音処理部104には、濃度
変換パラメータ算出部132の代わりに、CRT60の
画面上に表示されたオートラジオグラフィ画像から、入
力装置70を用いて選択された所望の領域の濃度データ
を算出する濃度データ算出部800と、この領域の面積
を算出し、濃度データと面積とに基づいて、濃度データ
の平均値、すなわち、単位面積当たりの濃度データの値
を算出する平均値算出部802とを設け、平均値算出部
202により算出された値を、濃度変換パラメータdと
して、濃度変換処理部134が、その入力されたデータ
の値eから、以下のような値を有する出力データ値を生
成して出力するようにしてもよい。
【0069】 出力データ値 = 255 (d≦e≦255) 出力データ値 = e (0≦e≦d) また、第1の実施態様においては、背景雑音処理部10
4が、濃度孤立点除去処理部136、二値化処理部14
0および孤立点除去処理部142とを備え、背景部分の
雑音に対応するデータを除去しているが、二値化処理部
140および孤立点除去処理部142を省略しても、あ
るいは、孤立点除去処理部142を省略してもよい。さ
らに、第1の実施態様における背景雑音処理部104か
ら、コントラスト強調処理部130あるいは濃度変換処
理部134を省略してもよい。また、第1の実施態様に
おいては、レーン分割処理部106のx軸方向周辺分布
算出部150が、濃度データの和にローパスフィルタ処
理を施し、さらに、処理の施されたデータとしきい値i
とを比較しているが、これらの処理の一方或いは双方を
省略してもよい。さらに、第1の実施態様においては、
スポット認識座標算出部166は、ある領域開始点とこ
れに隣接する極大点との間、隣接する極大点の間、或い
は、極大点とこれに隣接する領域終了点との間のy軸方
向の長さが、しきい値q以上である場合にかぎって、そ
の領域をスポットとして認識しているが、領域開始点と
これに隣接する極大点との間、隣接する極大点の間、或
いは、極大点とこれに隣接する領域終了点との間のy軸
方向の長さいかんにかかわらず、スポットと認識するよ
うにしてもよい。
【0070】また、第1の実施態様においては、画像デ
ータを構成する画素の濃度データの値が大きい場合に、
対応する画像がより白く表示され、小さい場合に、対応
する画像がより黒く表示されるように構成されている
が、逆に、画像データを構成する画素の濃度データの値
が小さい場合に、対応する画像がより白く表示され、大
きい場合に、対応する画像がより黒く表示されるように
構成してもよい。さらに、前記実施態様においては、蓄
積性蛍光体シート1を用いて得た画像データを、CRT
60の画面上に、画像として表示しているが、蓄積性蛍
光体シート1に代えて、写真フィルムを用いて、一旦、
可視画像を形成し、この可視画像を光電的に読み取り、
電気信号に変換した画像データに対して、同様の処理を
おこなうことも可能である。また、第2の実施態様にお
いては、ノイズ除去部300により、ローパスフイルタ
処理が施されて、高周波のノイズ成分が除去された画像
データに対して、ヒステリシス平滑化処理が施され、画
像データ中の濃度レベルのゆらぎを除去しているが、ヒ
ステリシス平滑化処理に代えて、メディアン・フィルタ
処理や移動平均法による処理を用いることもできる。同
様に、ノイズ除去部261は、ローパスフィルタ処理に
加えて、メディアン・フィルタ処理を施すことによっ
て、y軸方向周辺分布データを生成しているが、メディ
アン・フィルタ処理に代えて、ヒステリシス平滑化処理
あるいは移動平均法による処理を用いることもできる。
【0071】さらに、第1の実施態様においては、背景
雑音を除去することにより、画像データのノイズを除去
し、第2の実施態様においては、x軸方向の周辺分布デ
ータに、ローパスフイルタ処理およびヒステリシス平滑
化処理を施し、ローパスフィルタ処理およびメディアン
・フィルタ処理を施して、y軸方向周辺分布データを生
成することにより、画像データのノイズを除去している
が、第1の実施態様において、第2の実施態様において
用いたノイズ除去方法を用い、第2の実施態様におい
て、第1の実施態様において用いたノイズ除去方法を用
いて、画像データのノイズを除去するようにしてもよ
い。また、第1の実施態様におけるスポット認識処理部
108でなされるスポット認識処理に代えて、第2の実
施態様におけるスポット認識処理部208でなされるス
ポット認識処理をおこなっても、さらには、第2の実施
態様におけるスポット認識処理部208でなされるスポ
ット認識処理に代えて、第1の実施態様におけるスポッ
ト認識処理部108でなされるスポット認識処理をおこ
なうこともできる。さらに、第2の実施態様におけるス
ポット境界座標決定部268における処理は任意であ
り、これを省略することもできるし、また、スポット境
界座標決定部268における処理を、第1の実施態様に
おいて、実行することもできる。
【0072】また、第1の実施態様における代謝物分類
処理部122でなされる処理に代えて、第2の実施態様
における代謝物分類処理部222でなされる処理をおこ
なっても、さらには、第2の実施態様における代謝物分
類処理部222でなされる処理に代えて、第1の実施態
様における代謝物分類処理部122でなされる処理をお
こなうこともできる。さらに、本明細書において、手段
とは、必ずしも物理的手段を意味するものではなく、各
手段の機能が、ソフトウエアによって実現される場合も
包含する。また、一つの手段の機能が二以上の物理的手
段により実現されても、二以上の手段の機能が一つの物
理的手段により実現されてもよい。
【0073】
【発明の効果】本発明によれば、操作者による煩雑な操
作を必要とすることなく、迅速に、生化学測定画像中の
定量化あるいは定量解析すべき関心領域を画定すること
のできる生化学画像解析装置を提供することが可能とな
る。
【図面の簡単な説明】
【図1】図1は、本発明の実施態様にかかるオートラジ
オグラフィ画像解析装置のための画像データを生成する
画像読み取り装置の一例を示す略斜視図である。
【図2】図2は、本発明の実施態様にかかるオートラジ
オグラフィ画像解析装置および画像読み取り装置のブロ
ックダイアグラムである。
【図3】図3は、本実施態様にかかるオートラジオグラ
フィ画像解析装置の画像形成/解析部およびその周辺回
路を詳細に示したブロック図である。
【図4】図4は、本実施態様にかかる画像読み取り装置
により読み取られ、画像データ記憶部に記憶された画像
データから生成された表示画像データに対応し、CRT
に表示された画像を示す図である。
【図5】図5は、背景雑音処理部の構成を示すブロック
ダイヤグラムである。
【図6】図6は、濃度変換パラメータ算出部により生成
された濃度ヒストグラムを示す図である。
【図7】図7は、濃度孤立点除去処理部による処理を説
明するための図である。
【図8】図8は、孤立点除去処理部による処理を説明す
るための図である。
【図9】図9は、レーン分割処理部の構成を示すブロッ
クダイヤグラムである。
【図10】図10は、x軸方向周辺分布算出部による処
理を説明するための図である。
【図11】図11は、レーン分割処理部によるレーンの
分割方法を説明するための図である。
【図12】図12は、レーン分割処理部により形成され
たレーンが表示されたCRTの画面を示す図である。
【図13】図13は、スポット認識処理部の構成を示す
ブロックダイヤグラムである。
【図14】図14は、図12に示されるレーンL6のy
軸方向周辺分布データを示す図である。
【図15】図15は、スポット座標算出部による処理を
説明するための図である。
【図16】図16は、スポット認識処理部より形成され
たスポット表示されたCRTの画面を示す図である。
【図17】図17は、CRTの画面上に表示された画像
を用いて、代謝物分類処理部による処理を説明するため
の図である。
【図18】図18は、定量データ記憶部の定量データベ
ースを説明するための図である。
【図19】図19は、本発明の他の実施態様にかかるオ
ートラジオグラフィ画像解析装置の画像生成/解析部の
ブロックダイアグラムである。
【図20】図20は、x軸方向の定量すべき領域を示す
レーンを画定するレーン分割処理部の構成を示すブロッ
クダイアグラムである。
【図21】図21は、ピーク位置算出部のブロックダイ
アグラムである。
【図22】図22は、ヒステリシス平滑化処理の方法を
説明するための図である。
【図23】図23は、ピーク位置推定部における処理を
説明するための図である。
【図24】図24は、ピーク位置決定部における処理を
説明するための図である。
【図25】図25は、ピーク位置決定部における処理を
説明するための図である。
【図26】図26は、レーン開始点/終了点検出部にお
ける処理を説明するための図である。
【図27】図27は、スポット認識処理部の構成を示す
ブロックダイアグラムである。
【図28】図28は、代謝物分類処理部のブロックダイ
アグラムである。
【図29】図29は、代謝物分類処理部における処理を
説明するための図である。
【図30】図30は、本発明の他の実施態様におけるス
ポット認識処理部の構成を示すブロックダイアグラムで
ある。
【図31】図31は、しきい値を決定するための方法を
示すグラフである。
【図32】図32は、本発明の他の実施態様にかかる背
景雑音処理部の構成を示すブロックダイアグラムであ
る。
【符号の説明】
1 蓄積性蛍光体シート 2 レーザ光 3 レーザ光源 4 フィルタ 5 ビーム・エクスパンダ 6 光偏向器 7 fθレンズ 8 平面反射鏡 9 導光性シート 10 光検出器 11 増幅器 12 A/D変換器 13 ラインバッファ 20 画像読み取り装置 30 オートラジオグラフィ画像解析装置 40 信号処理手段 41 ラインバッファ 42 制御部 43 画像生成/解析部 50 画像データ記憶手段 51 画像データ一時記憶部 52 画像データ記憶部 60 CRT 70 入力装置 102、202 表示画像生成部 104 背景雑音処理部 106、206 レーン分割処理部 108、208 スポット認識処理部 110、210 データ修正部 112、212 定量処理部 114、214 定量データ記憶部 116、216 表データ生成部 118、218 画像合成部 120 ウィンドウメモリ 122、222 代謝物分類処理部 130 コントラスト強調処理部 132 濃度変換パラメータ算出部 134 濃度変換処理部 136 濃度孤立点除去処理部 140 二値化処理部 142 孤立点除去処理部 144 画像間演算部 150、250 x軸方向周辺分布算出部 152 レーン開始点検出部 154 レーン終了点検出部 156 レーン分割座標算出部 160、260、860 y軸方向周辺分布算出部 162、262 領域開始点/終了点検出部 164 領域内ピーク値検出部 166 スポット認識座標算出部 223 スポットグループ化部 224 レーングループ化部 225 スポットグループ補正部 252 ピーク位置算出部 254 レーン開始点/終了点検出部 256 レーン分割座標算出部 261 ノイズ除去部 264 領域内極大値検出部 266 スポット境界座標算出部 268 スポット境界座標決定部 300 ノイズ除去部 302 ピーク位置推定部 304 ピーク位置決定部 400 CRTの画面 401−1、401−2、・・・・401−10 レー
ン 500 入力波形 502 カーソル 504 出力波形 700、701、702、703、705 スポット分
割線 861 ノイズ除去部 862 極小値検出部 863 極大値検出部 864 しきい値決定部 865 スポット決定部 868 スポット境界座標決定部
───────────────────────────────────────────────────── フロントページの続き (72)発明者 佐藤 真 秋田県秋田市中通り2−2−32 住友生命 ビル4F 株式会社富士フイルムソフト開 発センター秋田内 (72)発明者 染 真人 神奈川県足柄上郡開成町宮台798番地 富 士写真フイルム株式会社内

Claims (19)

    【特許請求の範囲】
  1. 【請求項1】 少なくとも一つの試料を一次元的に展開
    して得られた試料中の標識物質の位置情報を、電気信号
    に変換して得られたデータに基づき、生成された生化学
    画像の画像データを記憶する画像データ記憶手段と、前
    記画像データ記憶手段に記憶された画像データの中から
    画像データを選択して、所定の処理を施すことにより、
    平面座標系に展開された表示画像データを生成する表示
    画像データ生成手段と、前記表示画像データに基づき、
    画像を表示する表示手段と、前記表示画像データを、前
    記表示画像データが展開される前記平面座標系の一方の
    座標軸方向に分割して、前記一方の座標軸方向の定量す
    べき領域を画定するレーンを形成するための第1の分割
    座標データを生成する第1の画像データ分割手段とを備
    えたことを特徴とする生化学画像解析装置。
  2. 【請求項2】 さらに、前記一方の座標軸方向に分割さ
    れた前記表示画像データを、前記平面座標系の他方の座
    標軸方向に分割し、前記レーンごとに、前記他方の座標
    軸方向の定量すべき領域を画定するスポットを形成する
    ための第2の分割座標データを生成する第2の画像デー
    タ分割手段を備えたことを特徴とする請求項1に記載の
    生化学画像解析装置。
  3. 【請求項3】 さらに、前記表示画像データの中から、
    前記表示手段に表示されるべき画像のノイズに対応する
    データを除去して、ノイズ除去画像データを生成する第
    1のノイズ除去手段を備え、前記第1の画像データ分割
    手段が、前記第1のノイズ除去手段によってノイズに対
    応するデータが除去されたノイズ除去データを、前記表
    示画像データが展開される前記平面座標系の一方の座標
    軸方向に分割して、前記一方の座標軸方向の定量すべき
    領域を画定するレーンを形成するための第1の分割座標
    データを生成するように構成されたことを特徴とする請
    求項1または2に記載の生化学画像解析装置。
  4. 【請求項4】 さらに、前記表示画像データの中から、
    前記表示手段に表示されるべき画像のノイズに対応する
    データを除去して、ノイズ除去画像データを生成する第
    2のノイズ除去手段を備え、前記第2の画像データ分割
    手段が、前記第2のノイズ除去手段によりノイズに対応
    するデータが除去されたノイズ除去データを、前記平面
    座標系の他方の座標軸方向に分割し、前記レーンごと
    に、前記他方の座標軸方向の定量すべき領域を画定する
    スポットを形成するための第2の分割座標データを生成
    するように構成されたことを特徴とする請求項2または
    3に記載の生化学画像解析装置。
  5. 【請求項5】 前記第1の画像データ分割手段が、前記
    一方の座標の値を共通にする画素データの有する濃度デ
    ータの値を累算し、前記一方の座標方向の周辺分布に対
    応するデータを算出する第1の周辺分布データ算出手段
    と、前記第1の周辺分布データ算出手段により算出され
    た周辺分布に対応するデータに基づいて、前記レーンの
    開始点の座標を算出するレーン開始点座標算出手段と、
    前記第1の周辺分布データ算出手段により算出された周
    辺分布に対応するデータに基づいて、前記レーン終了点
    の座標を算出するレーン終了点座標算出手段と、前記レ
    ーン開始点座標算出手段および前記レーン終了点座標算
    出手段により算出された前記開始点および前記終了点の
    座標に基づいて、前記レーンを画定する分割座標データ
    を生成するレーン分割座標データ生成手段とを備えたこ
    とを特徴とする請求項1ないし4のいずれか一項に記載
    の生化学画像解析装置。
  6. 【請求項6】 前記第2の画像データ分割手段が、前記
    第1の画像データ分割手段により生成された第1の分割
    座標データに基づいて画定されたレーンごとに、前記他
    方の座標の値を共通にする画素データの有する濃度デー
    タの値を累算し、前記他方の座標方向の周辺分布データ
    を算出する第2の周辺分布データ算出手段と、前記第2
    の周辺分布算出手段により算出された周辺分布に対応す
    るデータに基づいて、前記レーンごとのスポットを画定
    する第2の分割座標データを生成するスポット分割座標
    データ生成手段とを有することを特徴とする請求項2な
    いし5の何れか一項に記載の生化学画像解析装置。
  7. 【請求項7】 前記第1のノイズ除去手段が、前記第1
    の周辺分布データ算出手段により算出された周辺分布に
    対応するデータに対して、ノイズ除去処理をおこない、
    前記第1のノイズ除去手段によりノイズに対応するデー
    タが除去された周辺分布に対応するデータに基づいて、
    前記レーン開始点座標算出手段が、前記レーン終了点の
    座標を算出し、前記レーン終了点座標算出手段が、前記
    レーン終了点の座標を算出するように構成されたことを
    特徴とする請求項5または6に記載の生化学画像解析装
    置。
  8. 【請求項8】 前記第2のノイズ除去手段が、前記第2
    の周辺分布データ算出手段により算出された周辺分布に
    対応するデータに対して、ノイズ除去処理をおこない、
    前記第2のノイズ除去手段によりノイズに対応するデー
    タが除去された周辺分布に対応するデータに基づいて、
    前記スポット分割座標データ生成手段が、第2の分割座
    標データを生成するように構成されたことを特徴とする
    請求項6または7に記載の生化学画像解析装置。
  9. 【請求項9】 前記第1の画像分割手段が、前記第1の
    ノイズ除去手段によりノイズに対応するデータが除去さ
    れたノイズ除去データに一次元フーリエ変換を施して、
    濃度データのピークの座標を推定するピーク推定手段
    と、該ピーク推定手段により推定されたピークの座標を
    補正するピーク座標補正手段と、該ピーク座標補正手段
    により得られたピークの座標にしたがって、第1の分割
    座標データを生成するレーン分割座標データ生成手段を
    備えたことを特徴とする請求項1ないし8のいずれか1
    項に記載の生化学画像解析装置。
  10. 【請求項10】 さらに、前記第1の分割座標データと
    前記第2の分割座標データとに基づいて、前記各レーン
    に形成されたスポットが属すべきグループを示すグルー
    プデータを生成するグループデータ生成手段を備えたこ
    とを特徴とする請求項2ないし9のいずれか1項に記載
    の生化学画像解析装置。
  11. 【請求項11】 前記グループデータ生成手段が、各レ
    ーンの前記第2の分割座標データを対比して、前記各レ
    ーンに形成されたスポットをグループ化するスポットグ
    ループ化手段と、前記各レーンの前記第2の分割座標デ
    ータを比較して、各レーンをグループ化するレーングル
    ープ化手段と、該レーングループ化手段による前記各レ
    ーンのグループ化の結果に基づいて、前記スポットグル
    ープ化手段によりグループ化されたスポットのグループ
    を補正するスポットグループ化補正手段を備えたことを
    特徴とする請求項10に記載の生化学画像解析装置。
  12. 【請求項12】 前記スポット分割座標データ生成手段
    が、前記第2の周辺分布データ算出手段により算出され
    た周辺分布に対応するデータの値を所定値と比較して、
    領域開始点および領域終了点を検出する領域開始点/終
    了点検出手段と、前記領域開始点と前記領域終了点との
    間に存在する前記第2の周辺分布に対応するデータの微
    分値の符号が変化するピーク点を検出するピーク点検出
    手段と、前記領域開始点および領域終了点ならびに前記
    ピーク点に基づいて、前記スポットを形成するための第
    2の分割座標データを算出する分割座標データ算出手段
    とを有することを特徴とする請求項6ないし12のいず
    れか1項に記載の生化学画像解析装置。
  13. 【請求項13】 さらに、前記表示画像データの中か
    ら、前記表示手段に表示されるべき画像の背景ノイズに
    対応するデータを除去して、ノイズ除去画像データを生
    成する背景ノイズ除去手段を備え、前記第1の画像デー
    タ分割手段が、前記ノイズ除去画像データを、前記表示
    画像データが展開される前記平面座標系の一方の座標軸
    方向に分割して、前記一方の座標軸方向の定量すべき領
    域を画定するレーンを形成するための第1の分割座標デ
    ータを生成するように構成されたことを特徴とする請求
    項1、2、4、5、6、8、10、11、12のいずれ
    か1項に記載の生化学画像解析装置。
  14. 【請求項14】 前記背景ノイズ除去手段が、前記表示
    画像データを構成する複数の画素を含む領域を画定し、
    前記領域の所定の位置にある画素データの有する濃度デ
    ータが、前記領域の他の位置にある画素データの有する
    濃度データと異なる場合に、前記所定の位置にある画素
    データの有する濃度データの値を、所定の値に変換し、
    前記領域中の孤立した濃度データを有する画素データを
    除去する孤立点除去手段を備えたことを特徴とする請求
    項13に記載の生化学画像解析装置。
  15. 【請求項15】 前記背景ノイズ除去手段が、さらに、
    前記孤立点除去手段から出力されたデータを二値化する
    二値化手段と、前記二値化手段から出力されたデータか
    ら、複数の画素データを含む領域を画定し、前記領域の
    所定の位置にある画素データの有する濃度データが、前
    記領域の他の位置にある画素データの有する濃度データ
    と異なる場合に、前記所定の位置にある画素データの有
    する濃度データの値を、所定の値に変換し、該領域中の
    孤立した濃度データを有する画素データを除去する第2
    の孤立点除去手段を備えたことを特徴とする請求項14
    に記載の生化学画像解析装置。
  16. 【請求項16】 前記濃度データの値が、大きくなるに
    したがって、前記濃度データを有する画素データに対応
    する画像が、前記表示手段上に白く表示されるように構
    成され、前記背景ノイズ除去手段が、さらに、前記表示
    画像データを構成する画素データの有する濃度データ
    を、所定の範囲に拡張して、コントラストを強調するコ
    ントラスト強調手段と、前記コントラストの強調された
    表示画像データの濃度ヒストグラムを生成し、前記濃度
    ヒストグラムの面積の所定の割合を占める点より大きな
    値を有する濃度データに、濃度データの最大値を付与す
    る濃度変換手段とを有し、前記濃度変換手段により濃度
    データの変換されたデータが、前記孤立点除去手段に与
    えられることを特徴とする請求項13ないし15のいず
    れか1項に記載の生化学画像解析装置。
  17. 【請求項17】 前記濃度データの値が、大きくなるに
    したがって、前記濃度データを有する画素データに対応
    する画像が、前記表示手段上に白く表示されるように構
    成され、前記背景ノイズ除去手段が、さらに、前記表示
    画像データを構成する画素データの有する濃度データ
    を、所定の範囲に拡張して、コントラストを強調するコ
    ントラスト強調手段と、前記表示画像データから選択さ
    れた所定の領域に含まれる画素データの有する濃度デー
    タの平均値を算出し、前記平均値より大きな値を有する
    濃度データに、濃度データの最大値を付与する濃度変換
    手段を有し、前記濃度変換手段により濃度データの変換
    されたデータが、前記孤立点除去手段に与えられること
    を特徴とする請求項13ないし15のいずれか1項に記
    載の生化学画像解析装置。
  18. 【請求項18】 前記画像データが、蓄積性蛍光体シー
    トを用いて生成されたものであることを特徴とする請求
    項1ないし17のいずれか1項に記載の生化学画像解析
    装置。
  19. 【請求項19】 前記画像データが、オートラジオグラ
    フィ画像データおよび化学発光画像データよりなる群か
    ら選ばれた画像データにより構成されたことを特徴とす
    る請求項1ないし18のいずれか一項に記載の生化学画
    像解析装置。
JP22220795A 1995-01-09 1995-08-30 生化学画像解析装置 Expired - Fee Related JP3571819B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP22220795A JP3571819B2 (ja) 1995-01-09 1995-08-30 生化学画像解析装置
US08/581,261 US6144758A (en) 1995-01-09 1995-12-29 Biochemical image analyzing apparatus

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP1752195 1995-01-09
JP7-17521 1995-01-09
JP22220795A JP3571819B2 (ja) 1995-01-09 1995-08-30 生化学画像解析装置

Publications (2)

Publication Number Publication Date
JPH08247931A true JPH08247931A (ja) 1996-09-27
JP3571819B2 JP3571819B2 (ja) 2004-09-29

Family

ID=26354059

Family Applications (1)

Application Number Title Priority Date Filing Date
JP22220795A Expired - Fee Related JP3571819B2 (ja) 1995-01-09 1995-08-30 生化学画像解析装置

Country Status (2)

Country Link
US (1) US6144758A (ja)
JP (1) JP3571819B2 (ja)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6687394B1 (en) * 1999-04-08 2004-02-03 Fuji Photo Film Co. Ltd. Method and apparatus for quantifying image
JP2005181312A (ja) * 2003-12-19 2005-07-07 Palo Alto Research Center Inc サンプル媒体上のレチクルマークを用いて、希少細胞のスキャナ画像座標を顕微鏡座標に変換する方法
JP2010271873A (ja) * 2009-05-20 2010-12-02 Hypergear:Kk 画像の傾き角検出方法
JP2013224901A (ja) * 2012-04-23 2013-10-31 Tsubakimoto Chain Co 薄層クロマトグラフィ撮影装置
CN113177548A (zh) * 2021-05-08 2021-07-27 四川大学 一种针对免疫固定电泳的重点区域识别方法

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6498863B1 (en) 2000-09-20 2002-12-24 Media Cybernetics Inc. Method, system, and product for analyzing a digitized image of an array to create an image of a grid overlay
US6787761B2 (en) * 2000-11-27 2004-09-07 Surromed, Inc. Median filter for liquid chromatography-mass spectrometry data
US7583828B2 (en) * 2004-08-18 2009-09-01 General Electric Company Accumulation imaging
JP4238902B2 (ja) * 2006-09-04 2009-03-18 日本電気株式会社 文字ノイズ除去装置、文字ノイズ除去方法、文字ノイズ除去プログラム
JP2014220704A (ja) * 2013-05-09 2014-11-20 ソニー株式会社 画像処理装置および画像処理方法
CN108647272B (zh) * 2018-04-28 2020-12-29 江南大学 一种基于数据分布的小样本扩充对脱丁烷塔底丁烷浓度进行预测的方法
CN113674367B (zh) * 2021-08-20 2024-03-26 上海宝藤生物医药科技股份有限公司 一种电泳后脂蛋白胆固醇试剂扫描图的预处理方法

Family Cites Families (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3999047A (en) * 1972-09-05 1976-12-21 Green James E Method and apparatus utilizing color algebra for analyzing scene regions
US5028793A (en) * 1985-10-10 1991-07-02 Quantex Corporation Imaging screen for electrophoresis applications
US5281517A (en) * 1985-11-04 1994-01-25 Cell Analysis Systems, Inc. Methods for immunoploidy analysis
US5121436A (en) * 1987-08-14 1992-06-09 International Remote Imaging Systems, Inc. Method and apparatus for generating a plurality of parameters of an object in a field of view
JPS6460782A (en) * 1987-08-28 1989-03-07 Matsushita Refrigeration Silencer for enclosed motor compressor
JPS6460784A (en) * 1987-08-28 1989-03-07 Matsushita Refrigeration Enclosed motor compressor
JP3217343B2 (ja) * 1989-03-23 2001-10-09 オリンパス光学工業株式会社 画像処理装置
US5424186A (en) * 1989-06-07 1995-06-13 Affymax Technologies N.V. Very large scale immobilized polymer synthesis
US5133020A (en) * 1989-07-21 1992-07-21 Arch Development Corporation Automated method and system for the detection and classification of abnormal lesions and parenchymal distortions in digital medical images
JPH043952A (ja) * 1990-04-20 1992-01-08 Matsushita Electric Ind Co Ltd 表面解析方法および表面解析装置
JPH04232864A (ja) * 1990-07-12 1992-08-21 Bio Rad Lab Inc リン光体スクリーンを用いての生化学的アッセイにおける検出及びイメージング
US5539838A (en) * 1990-09-14 1996-07-23 Fuji Photo Film Co., Ltd. Abnormal pattern detecting apparatus pattern finding aparatus and liner pattern width calculating apparatus
US5185574A (en) * 1991-08-09 1993-02-09 Mayo Foundation For Medical Education And Research NMR measurements using recursive RF excitation
ES2152946T3 (es) * 1992-02-18 2001-02-16 Neopath Inc Metodo de identificacion de objetos mediante tecnicas de procesamiento de datos.
US5837832A (en) * 1993-06-25 1998-11-17 Affymetrix, Inc. Arrays of nucleic acid probes on biological chips
US5542003A (en) * 1993-09-13 1996-07-30 Eastman Kodak Method for maximizing fidelity and dynamic range for a region of interest within digitized medical image display
US5583973A (en) * 1993-09-17 1996-12-10 Trustees Of Boston University Molecular modeling method and system
JPH07299053A (ja) * 1994-04-29 1995-11-14 Arch Dev Corp コンピュータ診断支援方法
US5627907A (en) * 1994-12-01 1997-05-06 University Of Pittsburgh Computerized detection of masses and microcalcifications in digital mammograms

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6687394B1 (en) * 1999-04-08 2004-02-03 Fuji Photo Film Co. Ltd. Method and apparatus for quantifying image
JP2005181312A (ja) * 2003-12-19 2005-07-07 Palo Alto Research Center Inc サンプル媒体上のレチクルマークを用いて、希少細胞のスキャナ画像座標を顕微鏡座標に変換する方法
JP2012177697A (ja) * 2003-12-19 2012-09-13 Sri Internatl サンプル媒体上のレチクルマークを用いて、希少細胞のスキャナ画像座標を顕微鏡座標に変換する方法
JP2010271873A (ja) * 2009-05-20 2010-12-02 Hypergear:Kk 画像の傾き角検出方法
JP2013224901A (ja) * 2012-04-23 2013-10-31 Tsubakimoto Chain Co 薄層クロマトグラフィ撮影装置
CN113177548A (zh) * 2021-05-08 2021-07-27 四川大学 一种针对免疫固定电泳的重点区域识别方法
CN113177548B (zh) * 2021-05-08 2022-07-08 四川大学 一种针对免疫固定电泳的重点区域识别方法

Also Published As

Publication number Publication date
US6144758A (en) 2000-11-07
JP3571819B2 (ja) 2004-09-29

Similar Documents

Publication Publication Date Title
US6064755A (en) Image analyzing apparatus for producing density profile data of an image
JP3571819B2 (ja) 生化学画像解析装置
US6714669B2 (en) Image processing apparatus
US5946413A (en) Image analyzing apparatus for producing a profile of image data
US6343142B1 (en) Image analyzing apparatus
US6782140B1 (en) Image analyzing apparatus defining regions of interest
JP3851699B2 (ja) 画像解析装置
JPH0843538A (ja) 画像解析装置
US6256405B1 (en) Image forming apparatus
EP0677823B1 (en) Image analyzing apparatus
US6415038B1 (en) Image analyzing apparatus
JP2000292542A (ja) 画像定量方法および装置
US5831275A (en) Image forming apparatus
JP3453205B2 (ja) 画像表示装置および画像中の領域画定方法
JP3571801B2 (ja) 画像解析装置
JPH0814862A (ja) 画像解析装置
JPH11283019A (ja) 画像解析装置
JPH09212625A (ja) 画像処理装置
JP2000163591A (ja) 画像解析装置
JPH1051639A (ja) 画像処理装置
JP2001216497A (ja) 画像解析装置
JPH09190540A (ja) 画像形成装置
JPH09212671A (ja) 画像処理装置
JPH1049653A (ja) 画像表示装置
JPH08184676A (ja) 画像解析装置

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20040405

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20040428

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: 20040621

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20040625

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080702

Year of fee payment: 4

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080702

Year of fee payment: 4

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080702

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090702

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090702

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100702

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110702

Year of fee payment: 7

LAPS Cancellation because of no payment of annual fees