JP3905117B1 - 動状態解析システム及び方法、並びにコンピュータプログラム - Google Patents

動状態解析システム及び方法、並びにコンピュータプログラム Download PDF

Info

Publication number
JP3905117B1
JP3905117B1 JP2006122555A JP2006122555A JP3905117B1 JP 3905117 B1 JP3905117 B1 JP 3905117B1 JP 2006122555 A JP2006122555 A JP 2006122555A JP 2006122555 A JP2006122555 A JP 2006122555A JP 3905117 B1 JP3905117 B1 JP 3905117B1
Authority
JP
Japan
Prior art keywords
difference
image
analysis
value
axis
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP2006122555A
Other languages
English (en)
Other versions
JP2007293720A (ja
Inventor
西川  孝
Original Assignee
株式会社菱化システム
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 株式会社菱化システム filed Critical 株式会社菱化システム
Priority to JP2006122555A priority Critical patent/JP3905117B1/ja
Application granted granted Critical
Publication of JP3905117B1 publication Critical patent/JP3905117B1/ja
Publication of JP2007293720A publication Critical patent/JP2007293720A/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Analysis (AREA)

Abstract

【課題】動状態解析システムにおいて、画像の動状態或いは変化状態を高精度で解析できるようにする。
【解決手段】動状態解析システム(100)は、時間軸又は時間と相関関係にある数量軸上で所定間隔をおいて並ぶ少なくとも二枚の画像の差分を、画素単位で演算する差分演算部(111)と、画像上における少なくとも一方向に交わる配列をなす複数の画素についての、該演算された差分の平均値や最大値などの差分統計値を、前記一方向に沿った座標軸上の各点に対応付けて演算する差分統計値演算手段(112)と、該演算された差分統計値を、座標軸と時間軸又は数量軸とで張られた面に置くことで、解析用画像を生成する解析用画像生成手段(113)とを備える。
【選択図】図1

Description

本発明は、例えば時間に対して動く画像や、時間と相関関係にある数量軸に対して変化する画像における動状態や変化状態を解析する動状態解析システム及び方法、並びにコンピュータプログラムに関する。
この種の動状態解析システムは、伝統的な動状態数値化法の一つである「時空間画像法」を用いる。時空間画像法を用いる場合、水平座標軸及び垂直座標軸で張られる二次元平面の濃度値画像を第三軸としての時間軸上に並べることにより得られる時空間画像の水平又は垂直方向で時間軸に沿った方向に切断して得られる時空間断面画像をパターン解析して得られる数値で、動状態を記述している。このような方法の一種である「時系列差分画像法」では特に、時系列上の二枚の画像の差分をとり、画像内の変化のある部分のみを抽出し、その差分画像をパターン解析することにより動状態を数値化している。
従来この種の動状態解析システムとしては、例えば車両における自動操縦や操縦支援を行う技術分野において、車載カメラで捕らえた画像中における動状態数値化の対象物を、パターン認識技術によって特定した上で、その対象物の動きを上記数値化方法により数値化するシステムが開発されている(特許文献1参照)。
加えて、このような動状態数値化システムは、例えば、プラズマディスプレイにおける蛍光体の層厚や層間ギャップの測定、積層フィルムから構成される偏光板の温度による層厚変化の測定など、積層構造の層断面計測システムの主要部としても機能する。ここでは、計測対象物の表面に垂直な方向(言い換えれば、深度方向)を走査方向として、一定走査速度で撮像される複数の表面画像或いは光学界面画像を、上述の如き動画を構成する連続画像に見立ててその動状態を解析すれば、実際には動状態ではなく、複数の表面画像或いは光学界面画像における各層の変化状態を解析できる訳である。
特開2001―175845号公報
しかしながら、上述した背景技術によれば、画像或いは画像データにおけるS/N(Signal to Noise)比が必ずしも十分高くないという技術的問題点がある。例えば、前述の時系列差分画像法によれば、差分或いは差分値における雑音成分が相対的に大きいため、この雑音成分が大なり小なり差分画像をパターン解析する際の障害となっている。また特許文献1の如く、対象物をパターン認識した上で、その対象物の動きを追うことで動状態を数値化するシステムの場合には、そもそも対象物を特定しない限り、そのような数値化はできないという技術的問題点もある。
加えて、前述した層断面計測システムによれば、画像を取得する際に干渉プロファイルを用いるのが通常であり、その感度が一般に悪い或いは信号が微弱であることに起因して、数値化する際の雑音が顕著である。この結果、層厚や層間ギャップを精度良く計測することは困難或いは実践上不可能であるという技術的問題点がある。
本発明は上述の問題点に鑑み成されたものであり、例えば画像の動状態或いは変化状態を高精度で解析可能である動状態解析システム及び方法、並びにコンピュータをそのようなシステムとして機能させるコンピュータプログラムを提供することを課題とする。
本発明の動状態解析システムは上記課題を解決するために、時間軸上で所定間隔をおいて並ぶ少なくとも二枚の画像の差分を、画素単位で演算する差分演算手段と、前記画像上における少なくとも一方向に交わる配列をなす複数の画素についての、前記演算された差分の平均値、前記演算された差分の最大値及び前記演算された差分が閾値以上となる画素の個数のうち少なくとも一つである、差分統計値を、前記一方向に沿った座標軸上の各点に対応付けて演算する差分統計値演算手段と、該演算された差分統計値を、前記座標軸と前記時間軸とで張られた面に置くことで、解析用画像を生成する解析用画像生成手段とを備える。
本発明の動状態解析システムによれば、時間軸上で所定間隔をおいて並ぶ少なくとも二枚の画像を含む、複数の画像が動状態解析の対象となる。このような少なくとも二枚の画像は、例えば細胞の活性度評価用、車両の自動操縦用、手話入力用の動画を構成しており、時間軸上にて一定間隔で撮像或いは取得、入力、生成等される複数のフレーム画像のうち、相前後する少なくとも二枚のフレーム画像である
解析時には先ず、例えばプロセッサ、メモリ等を備えてなる、差分演算手段によって、時間軸上で所定間隔をおいて並ぶ少なくとも二枚の画像の差分が、画素単位で演算される。ここに「少なくとも二枚の画像の差分」とは、少なくとも二枚の画像間で、同一画素(即ち、同一位置或いは同一アドレスの画素)における、単色の又はRGB系やCMY系等のカラー各色の、濃度情報、輝度情報、彩度情報、色彩情報、色調情報等の、画像を規定するパラメータの差分或いは差分値を意味する。この際、「二枚の画像の差分」或いは「二枚のみの画像の差分」とは、一の画像における特定画素について、前又は後にある他の画像における同一の特定画素との間での、演算の対象となるパラメータについての差分或いは差分値である。そして「三枚以上の画像の差分」とは、一の画像における特定画素について、前又は後にある他の複数の画像における同一の特定画素との間での、演算の対象となるパラメータについての平均値等の基準値(即ち、3枚以上の画像における特定画素のパラメータに基づく基準値)からのズレとしての、差分或いは差分値である。単純な場合には、二枚の画像間での白黒表示における濃度の差分が、演算される。また、このような差分が演算される画像は、デジタル画像でもよいし、アナログ画像でもよい。更に、本発明に係る「差分」或いは「差分値」とは、時間軸に対する微小変化である「微分」或いは「微分値」を含む概念である。このように本発明に係る「差分」とは、狭義には、二枚の画像間における同一画素についての演算の対象となるパラメータの単純な差分を意味し、広義には、三枚以上の画像間における同一画素についての演算の対象となるパラメータの基準値からの差分或いはズレも意味し得る。
続いて、例えばプロセッサ、メモリ等を備えてなる、差分統計値演算手段によって、画像上における一方向に沿って設定された座標軸の各点において、該一方向に交わる(例えば、直交する)配列をなす複数の画素の差分統計値が演算される。ここに「一方向」は、例えば一枚の画像であるビットマップ画像やラスター画像における、該一枚の画像の水平方向、垂直方向、水平垂直の両方向、斜めなどの任意の方向或いは任意の二方向などであり、好ましくは当該動状態解析システムの解析対象となる動画の性質に応じて適宜に設定される性質のものである。また、このような画素の配列は、典型的には、同一画像、即ち同一時刻の画像或いは同一機会に取得された画像を構成する、画素の配列であり、該一方向に交わる一次元配列(即ち、一列の画素列)でもよいし、このような一次元配列を含む又は同一画像内で平面的に広がる2次元配列(即ち、複数列の画素列又は所定平面エリア内の画素群)であってもよい。言い換えれば、係る画素の配列は、一つの画像上において上記一方向に交わる方向に延びる線分上に少なくとも位置する。
但し、このような画素の配列は、複数の画像を時間軸に沿って配列した場合における、該複数の画像に跨るものであってもよい。即ち、このような画素の配列は、典型的には、同一画像上における一次元的な配列又は二次元的な配列であるが、複数の画像に跨る一次元的な配列、二次元的な配列又は三次元的な配列であってもよい。本発明に係る「差分統計値」は、差分の平均値、差分の最大値及び演算された差分が閾値以上となる画素の個数のうち、少なくとも一つである。差分に係る、これら平均値等以外の統計値として知られている既存の統計値を、ここでの差分統計値として用いることも可能である。
続いて、例えばプロセッサ、メモリ等を備えてなる、解析用画像生成手段によって、該演算された差分統計値が、座標軸と時間軸とで張られた面に置かれることで、解析用画像が生成される。ここでは、例えば座標軸と時間軸との二次元平面上の各点に、差分統計値が値により異なる色濃度やカラーなどで示される、解析用画像が生成される。数学的には、係る座標軸と時間軸とを二つの直交軸とする二次元平面上の各点について対応する差分統計値が割り当てられ、言い換えれば、係る二次元平面に対して更にこれに垂直な第3の軸として差分統計値を示す数量軸が設定され、係る二次元平面の各点における差分統計値の大小に応じて変化する色や濃度等が各点に割り当てられることとなり、これにより、色や濃度等のパターンが顕在化される。
このような一連の処理の中で、差分統計値演算手段によって差分統計値を演算する際には特に、画像或いは画像データにおけるランダムな又はランダム性が相対的に高いノイズ成分を、画像或いは画像データにおけるランダム性が相対的に低い信号成分に対して、相対的に弱めることができる。即ち、平均値、最大値又は差分が閾値以上となる画素の個数である、本発明に係る「差分統計値」は、単なる差分値或いは単なる微分値と比べて、S/N比率が高い。言い換えれば、画像或いは画像データにおける信号成分が微弱であっても、差分統計値を、演算する段階で、信号成分を相対的に強めている。逆に言えば、本発明に係る「差分統計値」は、これら平均値、最大値又は差分が閾値以上となる画素の個数に限らず、多少なりとも差分そのものと比べて、S/N比率が高くなるように或いは信号成分が相対的に強められるように該差分が統計処理された値であればよいことになる。
特に本発明では、差分統計値の一例として差分の「平均」を用いることで、前述の通りS/N比を向上させ、視野内で均一であるが微弱な動状態を顕著な信号として抽出できる。
差分統計値の他の例として差分の「最大値」を用いることで、視野内で局所的に発生する動状態を信号として抽出できる。尚、この場合に、差分の「平均」を用いると、動状態が局所的であるが故に、抽出される信号は弱められてしまう。
更に、差分統計値の他の例として差分に係る「閾値以上となる画素の個数」を用いることで、差分が極端に大きい動状態の信号と閾値以上で相対的に小さい動状態の信号とを、同じ重みとして扱うことができるため、相対的に小さい動状態の信号が、極端に大きい動状態の信号に隠れてしまうことを防止できる。尚、この場合に、前述の差分の「平均」や「最大値」を用いると、相対的に弱い動状態の信号が極端に大きい動状態の信号に隠されてしまう。
以上の結果、解析用画像生成手段により生成される解析用画像上では、元の画像上で各画素の濃度等そのものや、特に二枚の画像間の画素単位での差分そのものを直接解析する場合に比較して、画像或いは画像データにおけるノイズ成分が大きくても、より高精度で、画像の動状態或いは変化状態を解析可能となる。
従って、このような動状態解析システムを例えば、細胞の活性度評価システムに適用すれば、高精度で細胞の活性度を評価でき、自動操縦用システムに適用すれば、高精度で動く物体の状態を評価でき、手話入力用システムに適用すれば、高精度で手話入力或いは手話解読を行うことができる
特に本発明によれば、動く物体の状態を解析するために、先ず動く物体を特定する必要がない。言い換えれば、動く物体を特定すること自体も解析の一環として可能である。このように、先ず動く物体を特定する必要がある、従来の自動操縦用システムにおける時空間断面画像のパターン解析技術と比べて、本発明は格段に優れている。
本発明の動状態解析システムの一態様では、前記生成された解析用画像を、パターン解析するパターン解析手段を更に備える。
この態様によれば、例えばプロセッサ、メモリ等を備えてなる、パターン解析手段によって、解析用画像生成手段により生成された解析用画像を、自動的に解析可能となり、実践上有利である。例えば、パターン解析手段によって、細胞の活性度、動く物体の特定やその動状態、手話入力或いは手話解読等を、解析用画像上で自動的に実行可能となる。
本発明の動状態解析システムの他の態様では、動画を撮像して複数の画像を生成する撮像手段を更に備え、前記差分演算手段は、前記複数の画像に含まれる前記少なくとも二枚の画像の差分を演算する。
この態様によれば、例えば、光学顕微鏡、車載カメラ、テレビカメラ等を含んでなる撮像手段によって、動画が撮像され、時間軸上で所定間隔或いは一定間隔で並ぶ複数の画像が生成される。すると、差分演算手段によって、係る複数の画像に含まれる少なくとも二枚の画像の差分が各画素について演算され、更に、差分統計値演算手段による演算及び解析用画像生成手段による解析用画像の生成が行われる。従って高精度で、例えば、細胞の活性度の評価、動く物体の状態の評価、手話入力或いは手話解読等を、実行可能となる。
本発明の動状態解析方法は上記課題を解決するために、時間軸上で所定間隔をおいて並ぶ少なくとも二枚の画像の差分を、画素単位で演算する差分演算工程と、前記画像上における少なくとも一方向に交わる配列をなす複数の画素についての、前記演算された差分の平均値、前記演算された差分の最大値及び前記演算された差分が閾値以上となる画素の個数のうち少なくとも一つである、差分統計値を、前記一方向に沿った座標軸上の各点に対応付けて演算する差分統計値演算工程と、該演算された差分統計値を、前記座標軸と前記時間軸とで張られた面に置くことで、解析用画像を生成する解析用画像生成工程とを備える。
本発明の動状態解析方法によれば、上述した本発明の動状態解析システムの場合と同様に、より高精度で画像の動状態を解析可能となる。
尚、本発明の動状態解析方法においても、上述した本発明の動状態解析システムにおける各種態様と同様の各種態様を採ることが可能である。
本発明のコンピュータプログラムは上記課題を解決するために、コンピュータを上述した本発明に係る動状態解析システム(但し、その各種態様を含む)として機能させる。
本実施形態のコンピュータプログラムによれば、当該コンピュータプログラムを格納するCD−ROM、DVD−ROM等の記録媒体から、当該コンピュータプログラムを、動状態解析システムに備えられたコンピュータに読み込んで実行させれば、或いは、当該コンピュータプログラムを通信手段を介してダウンロードさせた後に実行させれば、上述した本発明に係る動状態解析システムを比較的簡単に構築できる。これにより、上述した本発明に係る動状態解析システムの場合と同様に、より高精度で画像の動状態を解析可能となる。
本発明の作用及び他の利得は次に説明する実施するための最良の形態から明らかにされよう。
以下、本発明の実施形態を図面に基づいて説明する。
<第1実施形態>
本発明の動状態解析システム及び方法、或いは、コンピュータを動状態解析システムとして機能させることが可能なコンピュータプログラム(以下、適宜、単にプログラムと称することもある)に係る第1実施形態を、図1から図7を参照して説明する。
先ず、図1を参照して、本実施形態に係る動状態解析システムの構成について説明する。図1は、本実施形態に係る動状態解析システムのブロック図である。
図1において、動状態解析システム100は、制御部100c及び撮像装置101を備える。制御部100cは、フレームバッファ102、差分演算部111、差分統計値演算部112、解析用画像生成部113、及びパターン解析部114を備える。動状態解析システム100は更に、表示装置116、プリンタ117、及び記憶装置118を備えて構成されている。
撮像装置101は、本発明に係る「撮像手段」の一例であり、例えばCCDカメラ等からなる車載カメラ、倒立型光学顕微鏡、テレビカメラ、ビデオカメラ等から構成されている。撮像装置101は、動画や、連続写真或いは連続顕微鏡写真を撮像し、時間軸上で所定間隔或いは一定間隔で並ぶ複数の画像を生成するように構成されている。
フレームバッファ102は、このように生成された動画等を、例えばフレーム画像の単位で逐次蓄える。或いは、連続写真を一枚ずつ逐次蓄える。より具体的には、フレームバッファ102は、時間軸又は時間と相関関係にある数量軸上で所定間隔をおいて並ぶ二枚の画像を含む、複数の画像を、少なくとも一枚ずつ或いは複数枚ずつ蓄える。
差分演算部111は、本発明に係る「差分演算手段」の一例であり、フレームバッファ102に蓄えられており、時間軸又は時間と相関関係にある数量軸上で所定間隔をおいて並ぶ二枚の画像の差分を、画素単位で演算するように構成されている。
差分統計値演算部112は、本発明に係る「差分統計値演算手段」の一例であり、画像上における少なくとも一方向に交わる配列をなす複数の画素についての、差分演算部111により演算された差分の差分統計値を、係る一方向に沿った座標軸上の各点に対応付けて演算するように構成されている。
ここで、このような座標軸について、図2を参照して説明を加える。ここに図2は、時間軸上で所定間隔をおいて並ぶ複数のフレーム画像における、一方向或いは座標軸、及び画素の配列を図式的に示す概念図である。
図2において、動画(又は連続写真或いは連続顕微鏡写真)DMは、時間軸t上における各時刻t1, t2,…, ti,…, tnに取得された複数のフレーム画像DM1, DM2,…, DMi,…, DMnから構成されている。ここで、フレーム画像DM1に着目すると、そのy軸に沿う方向を本発明に係る「一方向」の一例として、座標軸501をフレーム画像DM1上に定義すると共にこれに直交する座標軸503をフレーム画像DM1上に定義する。ここでは特に、フレーム画像DM1におけるx方向の中心線の位置に、座標軸501を定義する。すると、このフレーム画像DM1上にある各画素P(t1, x, y)についての差分は、この画素P(t1, x, y)の値(例えば、色別の濃度値)と、フレーム画像DM2上にある対応する画素P(t2, x, y)の値との差となる。即ち、差分演算部111(図1参照)により、このような差分が、フレーム画像DM1上にある座標軸503上の各画素Pについて逐次演算されることとなり、更に、差分統計値演算部112により、座標軸501上にある画素P(t1, xc, y)についての差分統計値が、座標軸503上にある複数の画素Pの差分を用いての差分統計値として演算される。更に、このような座標軸501上にある画素P(t1, xc, y)の差分統計値の演算を各yに対して行うことで、即ち図2において座標軸503を座標軸501に沿って上下に移動させつつ座標軸501上にある画素P(t1, xc, y)についての差分統計値の演算を、該座標軸503上にある複数の画素Pの差分を用いて行うことで、一枚のフレーム画像DM1におけるx方向の中心線である座標軸501上の各画素P(t1, xc, y)についての差分統計値が得られる。更に、このような演算を、複数のフレーム画像DM1, DM2,…, DMi,…, DMnについて同様に行うことで、最終的には、各フレーム画像DMiにおけるx方向の中心線である座標軸501上の各画素P(t1, xc, y)についての差分統計値が得られる。
再び図1に戻り、差分統計値演算部112は、画素P(t1, xc, y)に対する差分統計値として、座標軸503上にある複数の画素Pの差分の平均値又は最大値を演算するように構成されている。或いは、座標軸503上にある複数の画素Pの中で、差分が閾値以上となる画素の個数を、差分統計値として演算してもよい。
解析用画像生成部113は、本発明に係る「解析用画像生成手段」の一例であり、差分統計値演算部112により演算された差分統計値を、座標軸501と時間軸tとで張られた面に置くことで、解析用画像を生成するように構成されている。より具体的には、座標軸501上にある各画素P(t1, xc, y)の差分統計値を、座標軸501(言い換えれば「方向軸」)と、時間軸又は時間と相関関係にある数量軸とで張られる解析用画像の画素の数値とする。
ここで、このような解析用画像について、図3を参照して説明を加える。ここに図3は、解析用画像を図式的に示す概念図である。
図3において、解析用画像DTMは、図2に示したx方向の中心線である座標軸501のところにxの値が固定されており、yの値及びtの値をパラメータとする、一枚の解析用画像DTMである。ここでは、図2においてx方向に沿って配列された複数の画素についての差分が、座標軸501上の画素P(t, xc, y)に対する差分統計値を演算するために専ら用いられているので、当該解析用画像DTMでの画素P(t, xc, y)における信号成分が相対的に強められる結果となっている。尚、図3中、動状態パターンPATを粗いドット状で示しているが、実際の動状態パターンPATは、解析用画像DTMにおいて、ドットが識別できない程度に細いか、又はドット状ではなく連続的或いはアナログ的なパターンであってもよい。
再び図1に戻り、パターン解析部114は、本発明に係る「パターン解析手段」の一例であり、解析用画像生成部113により生成された解析用画像DTMを、パターン解析し、その解析結果を所定フォーマットの電子データとして生成するように構成されている。
表示装置116は、例えばLCD(液晶表示装置)、CRT(カソードレイチューブ)等からなり、パターン解析部114によるパターン解析結果たる電子データを、所定フォーマットで表示する。プリンタ117は、例えばレーザプリンタ、インクジェットプリンタ等からなり、パターン解析部114によるパターン解析結果たる電子データを、所定フォーマットで印刷する。記憶装置118は、例えばハードディスク装置、光ディスク装置等からなり、パターン解析部114によるパターン解析結果たる電子データを、所定フォーマットで記憶するように構成されている。
制御部100cは、以上の如き機能を有する、差分演算部111、差分統計値演算部112、解析用画像生成部113、及びパターン解析部114をハードウエア的に備えていてもよい。或いは、好ましくはCPU、メモリ等を備えたコンピュータをこれら各部として機能させるコンピュータプログラムが格納されたCD−ROM、DVD−ROM等の記録媒体からの該コンピュータプログラムのインストールや、或いはインターネット等経由のダウンロードによるインストールなどによって、これら各部がコンピュータ内に論理的に構築されてもよい。
尚、図2及び図3に示したように、差分統計値が演算される画素Pを、x方向の中心線上に固定された座標軸501上の画素Pにのみ限定する(その結果として、座標軸501に対応する一枚の解析用画像DTMのみが最終的に得られる)必要はない。即ち、差分統計値が演算される画素Pを、x軸に沿って所定間隔で配列された複数の座標軸501上の画素Pとしてもよい(その結果として、複数の座標軸501に対応する複数枚の解析用画像DTM(x)が最終的に得られることになる)。
例えば図4に示すように、図2の場合と同様に差分を演算した後に、フレーム画像DM1上にある座標軸503上の各画素Pについて、その左右において座標軸503上で幅504(但し、フレーム画像のx方向の長さより短い幅)内にある複数の画素Pの差分を用いて、夫々演算してもよい。即ち、所定間隔で或いは画素単位で、座標軸503上に並ぶ複数の座標Pについての差分統計値を、自らを中心に幅504内に位置する他の複数の画素の差分を用いて演算してもよい。このような演算を、座標軸501に沿って配列された複数の座標軸503について行い、更に、これを複数のフレーム画像DMについて行えば、例えば、図5に示したように、動状態パターンPATを画像として夫々有する複数の解析用画像DTM(x)を生成できる。図2及び図3の場合と異なり図4及び図5の場合には、x方向の中心線付近についてのみならず、所望のx座標における複数の解析用画像DTM(x)が得られる。このように本実施形態の骨子であるS/N比の改善と迅速に動きのある場所の特定を行うという特徴を発揮させるためには、幅504は、座標軸503に対して2つ以上であって、動きの探索領域を全て包含する個数に設定することも可能である。尚、図4及び図5は、図2及び図3と夫々同趣旨の概念図であり、図2及び図3と同様の構成要素には同様の参照符号が付与されている。
ここで図6を参照して、解析用画像生成部113により生成される解析用画像について、具体例を挙げて説明を加える。ここに、図6(a)は、車載カメラにより撮像される動画DMを示す平面図であり、図6(b)及び図6(c)は、その動画について差分演算部111による演算及び差分統計値演算部112による演算を経て、解析用画像生成部113により生成された、x方向の解析用画像DTMX及びy方向の解析用画像DTMYを夫々、概念的に示す概念図である。
図6(a)に示すように、動画DM中には、白い大きな矢印で示されている不図示の自車の進行方向の前方に、動くものとして左から右へ進行中の車両D101があり、更に、自車の進行に伴って相対的に手前側に動く街灯D102とビルD103とがあり、加えて、直進している自車の進行に伴って殆ど動かない太陽D104がある。
図6(b)及び図6(c)に示すように、x方向の解析用画像DTMX及びy方向の解析用画像DTMYでは、車両D101、街灯D102、ビルD103及び太陽D104の動状態が、帯状のパターンとして顕在化されている。先ず、解析用画像DTMX及びDTMY上で設定された閾値以上の画素を探すことで、動く対象物を特定することが可能となる。更に、解析用画像DTMX及びDTMYの帯状のパターンの方向と傾きと解析することにより、不図示の自車の走行による相対移動であるか、又は他の対象物の移動であるかを特定することが可能となる。より具体的には、自車の走行により相対移動する対象物(例えば、街灯D102やビルD103)は、x方向及びy方向夫々の消失点座標を指向する帯状のパターンとして解析することができ、自車の振動により相対移動する無限遠に位置すると考えられる対象物(例えば、太陽D104)は、x座標及びy座標を変えない縦の帯状のパターンとして解析することができ、他の移動する対象物(例えば、車両D101)は、これら以外の方向の帯状のパターンとして解析することができる。
次に、以上のように構成された動状態解析システム100の解析動作について、図1等に加えて、図7のフローチャートを参照して説明を加える。
図7において、先ず撮像装置101によって、車載カメラや倒立型光学顕微鏡による動画や連続写真など、時間軸又は時間と相関関係にある数量軸上で所定間隔をおいて並ぶ二枚の画像を含む、複数の画像が撮像される(ステップS11)。次に、差分演算部111によって、時間軸又は数量軸上で所定間隔をおいて並ぶ二枚の画像の差分が、各画素Pについて演算される(ステップS12)。尚、ここでの差分は、例えば単色の又はカラー各色の濃度値の差分である。続いて、全ての画素Pについて差分演算が終了したかが判定され(ステップS13)、終了していなければ(ステップS13:NO)、差分の演算が続行される(ステップS12)。
このステップS13における判定の結果、全ての画素Pについて差分演算が終了していれば(ステップS13:Yes)、差分統計値演算部112によって、画像上における一方向に沿って設定された座標軸の各点において、該一方向に直交する配列をなす複数の画素P(図2参照)を用いての差分統計値が演算される(ステップS14)。尚、ここでの差分統計値は、例えば、前述したように平均値又は最大値である。或いは前述したように閾値以上となる画素の個数である。続いて、差分統計値を演算すべき座標軸上の全ての画素Pについて差分統計値演算が終了したかが判定され(ステップS15)、終了していなければ(ステップS15:NO)、差分統計値の演算が続行される(ステップS14)。
このステップS15における判定の結果、差分統計値を演算すべき座標軸上の全ての画素Pについて差分統計値演算が終了していれば(ステップS15:Yes)、解析用画像生成部113によって、該演算された差分統計値が、座標軸と時間軸又はこれと相関関係にある数量軸とで張られた面に置かれることで、例えば解析用画像DTM(図3参照)や解析用画像DTM(x)(図5参照)が生成される。このような解析用画像上では、差分統計値の値に応じた色や濃度等の割り当てにより、解析用画像DTMやDTM(x)上でパターンPATが顕在化される(図3、図5、並びに図6(b)及び(c)参照)。
続いて、このように生成された解析用画像を、パターン解析部114によって自動でパターン解析することで、元となる画像における動状態或いは変化状態が解析され(ステップS17)、その解析結果が、表示装置116或いはプリンタ117によって表示或いは出力され、或いは、記憶装置118に記憶される(ステップS18)。
以上のように第1実施形態によれば、ステップS14における差分統計値演算部112によって演算した差分統計値を用いることで、信号成分をノイズ成分に対して強められる。この際、例えば差分統計値として「最大値」を用いれば、視野内で局所的に発生する動く対象物を感度良く捉えられ、画像の動状態が解析可能となる。特に本発明によれば、例えば図6に例示した車両D101の如き、対象物がどこにあるかを予め特定しなくても、この対象物を特定する処理も含めて、解析可能であるので、実用上大変有利である。
<第2実施形態>
本発明の動状態解析システム及び方法、或いは、コンピュータを動状態解析システムとして機能させることが可能なコンピュータプログラムに係る第2実施形態を、図8から図13を参照して説明する。第2実施形態では、第1実施形態と比較して、画像の取得に関する構成、及び画像データの属性が異なっており、更にこれに伴い処理の一部も異なっている。よって、第2実施形態について、第1実施形態と重複する説明を省略すると共に、図面上における共通箇所には同一符号を付して示し、基本的に異なる点についてのみ、図8から図13を参照して説明する。
先ず、図8を参照して、本実施形態に係る動状態解析システムの構成について説明する。図8は、図1と同趣旨の、本実施形態に係る動状態解析システムのブロック図である。
図8において、動状態解析システム200は、制御部200c、光学顕微鏡装置201、及び駆動装置204を備える。制御部200cは、光学顕微鏡装置201からの画像が取り込まれる画像メモリ202と、駆動装置204における走査を制御する走査制御装置203とを備え、その他、第1実施形態の場合と同様に、差分演算部111、差分統計値演算部112、解析用画像生成部113、及びパターン解析部114を備える。動状態解析システム100は更に、第1実施形態の場合と同様に、表示装置116、プリンタ117、及び記憶装置118を備えて構成されている。
光学顕微鏡装置201は、本発明に係る「撮像手段」の一例であり、多層構造体800に光源光を照射し、その反射光と参照光との干渉を利用して、光学的に多層構造体800内に存在する各反射界面からの微小な反射光を受光することで、各画像を取得する。より具体的には、多層構造体800の表面に垂直な、即ち多層構造体800の深度方向を走査方向として、その焦点位置或いは撮像位置が走査方向に沿って走査されつつ、各深度における表面画像或いは光学界面画像を撮像可能に構成されている。
画像メモリ202は、このように取得された複数の画像を、少なくとも一枚ずつ或いは複数枚ずつ、走査線に対応付けて、即ち各深度別に逐次に蓄える。
駆動装置204は、このような走査方向に沿って一定速度で連続的或いは間欠的な走査を行わせるべく、光学顕微鏡装置201の焦点位置或いは撮像位置を変更させる駆動を行い、走査制御装置203は、駆動装置204による駆動動作を、光学顕微鏡装置201における撮像動作と同期制御するように構成されている。
尚、駆動装置204は、光学顕微鏡装置201を多層構造体800の表面に沿って(即ち、X方向やY方向に)移動可能に構成されてもよく、このように構成すれば、より大きな多層構造体800の所望の平面位置における断面図を、最終的に得られる解析用画像として、生成可能となる。但し、例えば数mm角や数百μm角など、多層構造体800がその平面サイズにおいて光学顕微鏡装置201の撮像可能な平面範囲と比べて十分に小さければ、このように光学顕微鏡装置201を多層構造体800の表面に沿って移動させる必要はない。
このため、走査制御装置203による制御下で、駆動装置204により光学顕微鏡装置201が一定方向(即ち、深度方向)に連続的又は断続的或いは間欠的に走査されつつ、光学顕微鏡装置201による画像取得が順次行われ、時間軸と相関関係にある数量軸上で所定間隔或いは一定間隔で並ぶ複数の画像が取得される。即ち、各深度の光学界面図が取得される。そして、このように取得された画像或いは画像データが、画像メモリ202に走査線に対応付けられて、即ち各深度に対応付けられて逐次に蓄えられ、以下、第1実施形態の場合と基本的に同様に、差分演算部111及び差分統計値演算部112による演算、更に解析用画像生成部113やパターン解析部114による処理が行われることになる。
ここで第2実施形態における差分や差分統計値の演算について、図9を参照して説明を加える。ここに図9は、時間と相関関係にある数量軸上で所定間隔をおいて並ぶ複数の断面画像における、一方向或いは座標軸、及び画素の配列を図式的に示す概念図である。
図9において、光学顕微鏡装置201の対物レンズ600は、太い矢印ARの如き光源光を照射し、その反射光を受光し、参照光との干渉を利用しての光学界面図SVを取得する。この際、駆動装置204により駆動されることで、その対物レンズ600の焦点位置或いは撮像位置は、走査方向(即ち、深度方向或いはZ軸方向)に、図中細い矢印SCNで示したように順次走査され、各深度の光学界面図SV1、SV2、…、SVi、…、SVmを取得することになる。ここで、光学界面図SV1に着目すると、そのx軸に沿う方向を本発明に係る「一方向」の一例として、座標軸601を光学界面図SV1上に定義すると共にこれに直交する座標軸603を光学界面図SV1上に定義する。ここでは特に、光学界面図SV1におけるy方向の中心線の位置に、座標軸601を定義する。また、本発明に係る「時間軸と相関関係にある数量軸」の一例を、ここでのz軸と定義する。すると、この光学界面図SV上にある各画素Q(x、y、z1)についての差分は、この画素Q(x、y、z1)の値(例えば、色別の濃度)と一定時間ごとに矢印SCN方向にΔZだけずらして得られる画素(x、y、z2)(但し、z2=z1+Δz)の値との差となる。即ち、差分演算部111(図8参照)により、このような差分が光学界面図SV1上にある座標軸603上の各画素Pについて逐次演算されることとなり、更に、差分統計値演算部112により、座標軸603上にあり且つ座標軸601上にある画素Q(x、yc,z1)についての差分統計値が、座標軸603上にある複数の画素Pの差分を用いての差分統計値として演算される。更に、このような座標軸601上にある画素Q(x、yc,z1)の差分統計値の演算を各xに対して行うことで、一枚の光学界面図SV1におけるy方向の中心線である座標軸601上の各画素Q(x、yc,z1)についての差分統計値が得られる。更に、このような演算を、複数の光学界面図SV1、SV2、…、SVi、…、SVmについて同様に行うことで、最終的には、各光学界面図SViにおけるy方向の中心線である座標軸601上の各画素Q(x、yc,z1)についての差分統計値が得られる。
更に第2実施形態における解析用画像生成部113により生成される解析用画像について、図10を参照して説明を加える。ここに図10は、解析用画像を図式的に示す概念図である。
図10において、解析用画像STMは、図9に示したy方向の中心線である座標軸601のところにyの値が固定されており、xの値及びzの値をパラメータとする、一枚の解析用画像STM、即ち多層構造体800のy方向の中央断面である。ここでは、図9においてy方向に沿って配列された複数の画素についての差分が、座標軸601上の画素Q(x、yc,z)に対する差分統計値を演算するために専ら用いられているので、画素Q(x、yc,z)における信号成分が相対的に強められる結果となっている。尚、図10中、動状態パターンPATを粗いドット状で示しているが、実際の動状態パターンPATは、解析用画像STMの全体に対して、ドットが識別できない程度に細いか、又は或いはドット状ではなく連続的或いはアナログ的なパターンであってもよい。
尚、図9及び図10に示したように、差分統計値が演算される画素Qを、y方向の中心線上に固定された座標軸601上の画素Qにのみ限定する(その結果として、座標軸601に対応する一枚の解析用画像STMのみが最終的に得られる)必要はない。即ち、差分統計値が演算される画素Qを、y軸に沿って所定間隔で配列された複数の座標軸601上の画素Qとしてもよい(その結果として、複数の座標軸601に対応する複数枚の解析用画像STM(y)が最終的に得られることになる)。
例えば図11に示すように、図9の場合と同様に差分を演算した後に、光学界面図SV1上にある座標軸603上の各画素Qについて、その上下において座標軸603上で幅604(但し、光学界面図のy方向の長さより短い幅)内にある複数の画素Qの差分を用いて、夫々演算してもよい。即ち、所定間隔で或いは画素単位で、座標軸603上に並ぶ複数の座標Qについての差分統計値を、自らを中心に幅604内に位置する他の複数の画素の差分を用いて演算してもよい。このような演算を、座標軸601に沿って配列された複数の座標軸603について行い、更に、これを複数の光学界面図SVについて行えば、例えば、図12に示したように、動状態パターンPATを画像として夫々有する複数の解析用画像STM(y)を生成できる。このように図9及び図10の場合と異なり、図11及び図12の場合には、中心線付近についてのみならず、所望のy座標における複数の解析用画像STM(y)が得られる。尚、図11及び図12は、図9及び図10と夫々同趣旨の概念図であり、図2及び図3と同様の構成要素には同様の参照符号が付与されている。
ここで図13を参照して、解析用画像生成部113により生成される解析用画像について、具体例を挙げて説明を加える。図13(a)は、本実施形態により生成された解析用画像STMを示す平面図であり、図13(b)は、比較例における同趣旨の平面図である。
図13(a)に示すように、解析用画像STM(y)中では、層701〜706について、各層の配置や、特に層704が薄くなっている層703と層705との間の層間ギャップGPについても、高い精度で、その層厚が計測可能とされている。これは主に、差分統計値演算部112によって差分統計値を演算する際に画像データにおけるランダム性が高いノイズ成分を弱めることにより、達成されるものと考察される。
図13(b)に示すように、本実施形態の構成から差分統計値を演算する処理を除去して、差分をそのまま使って解析用画像を生成した比較例によれば、層701〜706について、各層の配置が明確ではなく、特に層704が薄くなっている領域では、不明確さが顕著であり、層間ギャップGPについても、高精度での層厚計測が不可能である。これは主に、差分統計値演算部112によって差分統計値を演算する際に画像データにおけるランダム性が高いノイズ成分を弱めることができていないためと考察される。
以上のように第2実施形態によれば、差分統計値演算部112によって差分統計値を演算する際に、画像データにおけるランダム性が高いノイズ成分を弱めることができるので、解析用画像上では、より高精度で、画像の動状態或いは変化状態を解析可能となる。また、微小ボイドや微小異物が分布していると想定される透明積層フィルムに対して本実施形態を適用し、この際、差分統計値として前述の「最大値」又は「閾値以上の画素の個数」を用いることにより、その深さ方向の分布状態が解析可能となる。
<各種実施例>
次に以上のように構成された実施形態を適用した各種実施例について説明する。
<<第1実施例>>
前述の第1実施形態を、細胞の活性度評価システムに適用した第1実施例について、図14から図16を参照して説明する。ここに図14は、本実施例のシステムのブロック図であり、図15は、本実施例のシステムにより撮像される、シャーレ内における細胞サンプルが動く様子を概念的に示すx−y平面図であり、図16(a)及び(b)は、本実施例により生成される、x軸及びy軸に夫々沿った方向の時間的変化を示す解析用画像である。
図14において、細胞の活性度評価システム10は、シャーレ内の細胞サンプル11を評価対象とするものであり、倒立型光学顕微鏡の照明装置12、対物レンズ13及び撮像装置14と、計測及び制御用のコンピュータ15とを備えて構成される。
このように構成されているので、シャーレ内の底部に設置された透明ガラス上の細胞サンプル11の透過顕微鏡像が、撮像装置14で観測可能となる。この際、シャーレ内の環境を調整することで細胞サンプル11を生きた状態に保つことができ、図15に例示する如く動く細胞サンプル11の画像が、動画或いは連続顕微鏡写真として、撮像装置14により取り込まれる。更に、前述の差分演算部、差分統計値演算部及び解析用画像生成部として機能するコンピュータ15によって、設定された時間毎に取り込まれた画像の差分の絶対値を算出し、更に差分統計値の一例として差分の最大値を算出することで、図16(a)に例示する如き、横軸が原画像のX方向の座標であり且つ縦軸が時間軸の解析用画像と、図16(b)に例示する如き、横軸が原画像のY方向の座標であり且つ縦軸が時間軸の解析用画像を作成する。これらの解析用図面上では、動きがない箇所は、基本的に何らのパターンも描かれることの無い平素な箇所或いは領域となり、動きがある箇所は、該動きに応じて何らのかのパターンが描かれることになる。
本実施例によれば、差分統計値の一例として差分の最大値を用いて生成された解析用画像の細線化処理により、細胞サンプル11の動きの総量、動きのパターン(例えば、直線運動か回転運動か)、動きの速度等を、高い精度で算出でき、細胞の活性度を数値化できる。
<<第2実施例>>
前述の第1実施形態を、自動操縦用動状態解釈システムに適用した第2実施例について、図17及び図18を参照して説明する。ここに図17は、本実施例のシステムのブロック図であり、図18は、本実施例により生成される、x軸及びy軸に夫々沿った方向の時間的変化を示す解析用画像である。
図17において、車両25に搭載された自動操縦用動状態解釈システム20は、車載カメラ21と、計測及び制御用のコンピュータ22とを備える。
このように構成されているので、前述の差分演算部、差分統計値演算部及び解析用画像生成部として機能するコンピュータ22によって、取り込まれた画像の差分の絶対値を算出し、更に差分統計値の一例として差分の最大値を算出することで、図18に例示する如き、横軸が原画像のX方向の座標であり且つ縦軸が時間軸である解析用画像と、横軸が原画像のY方向の座標であり且つ縦軸が時間軸である解析用画像とを作成する。これらの解析用図面上では、動きがない箇所は、基本的に何らのパターンも描かれることの無い平素な箇所或いは領域となり、動きがある箇所は、該動きに応じて何らのかのパターンが描かれることになる。
本実施例によれば、差分統計値の一例として差分の最大値を用いて生成された解析用画像上での流れのパターンから、自車である車両25の進行方向、速度、消失点の座標などを知ることが可能となる。更に、例えば図6に示したように、自車以外の車両25から見て視野を横切る対象物についても、その存在と方向、速度などを知ることも可能となる。
<<第3実施例>>
前述の第1実施形態を、手話入力システムに適用した第3実施例について、図19及び図20を参照して説明する。ここに図19は、本実施例のシステムのブロック図であり、図20(a)及び(b)は、本実施例により生成される、x軸及びy軸に夫々沿った方向の時間的変化を示す解析用画像である。
図19において、手話入力システム30は、話者35の手36の動き37を動画又は連続写真として撮像するビデオカメラ31と、計測及び制御用のコンピュータ32とを備える。
このように構成されているので、前述の差分演算部、差分統計値演算部及び解析用画像生成部として機能するコンピュータ32によって、取り込まれた画像の差分の絶対値を算出し、更に差分統計値の一例として差分の最大値を算出することで、図20(a)に例示する如き、横軸が原画像のX方向の座標であり且つ縦軸が時間軸である解析用画像と、図20(b)に例示する如き、横軸が原画像のY方向の座標であり且つ縦軸が時間軸である解析用画像を作成する。これらの解析用図面上では、動きがない箇所は、基本的に何らのパターンも描かれることの無い平素な箇所或いは領域となり、動きがある箇所は、該動きに応じて何らのかのパターンが描かれることになり、手話の区切りは、動きの一時的に停止するパターンとして顕在化している。ここで特に、予め、手話における文字、単語、句、節、定型文章などと、図20(a)及び(b)に出現するパターンとの対応関係を特定しておき、例えば電子データとしてテーブル化しておけば、その後、実際に手話入力を行う際に、このテーブルを参照してのパターンマッチングによって、図20(a)及び(b)の如き解析用画像に基づいて、対応する手話の内容を解釈できる。
本実施例によれば、差分統計値の一例として差分の最大値を用いて生成された解析用画像上での流れのパターンから、手話者35の手36の動き37を撮像することで、手話入力を、自動的に行うことが可能となる。
<<第4実施例>>
前述の第2実施形態を、層断面計測システムに適用した第4実施例について、図21及び図22を参照して説明する。ここに図21は、本実施例のシステムのブロック図であり、図22は、本実施例のシステムにより生成される解析用画像である、x−z層断面及びy−z層断面である。
図21において、層断面計測システム40は、透明多層構造体サンプル43を評価対象とするものであり、光学顕微鏡装置41と、光学顕微鏡装置41の焦点位置或いは撮像位置を透明多層構造体サンプル43の表面に垂直な走査方向に、相対的に順次移動させる駆動装置44と、計測及び制御用のコンピュータ45とを備えて構成される。
このように構成されているので、光学顕微鏡装置41の撮像位置を、透明多層構造体サンプル43に対して、その表面の法線方向或いは深さ方向に順次に走査させることにより、或いは必要に応じて、更にその表面に沿って順次に走査させることにより、透明多層構造体サンプル43の連続的な顕微鏡像が、光学顕微鏡装置41で観測可能となる。ここで撮像される連続的な画像は、時系列でZ軸を変化させて取得されるので、透明多層構造体サンプル43の各深度における光学界面画像を表すことになる。この際、前述の差分演算部、差分統計値演算部、解析用画像生成部、及び走査制御部として機能するコンピュータ45によって、取り込まれた画像の差分の絶対値を算出し、更に差分統計値を算出することで、図22(a)に例示する如き、横軸が原画像のX方向の座標であり且つ縦軸が原画像のZ方向の座標である、解析用画像と、図22(b)に例示する如き、横軸が原画像のY方向の座標であり且つ縦軸が原画像のZ方向の座標である、解析用画像を作成する。これらの解析用図面上では、積層状態に変化がない箇所、即ち反射界面が存在しない箇所は、基本的に何らのパターンも描かれることの無い平素な箇所或いは領域となり、積層状態に変化がある箇所、即ち反射界面が存在する箇所には、積層パターンが現れることになる。
本実施例によれば、差分統計値を用いて生成された解析用画像により、ミクロンオーダといった非常に薄い層の位置や層厚、層間ギャップ或いは空隙の位置や厚みなどを、高い精度で計測できる。
尚、本発明は、上述した実施形態に限られるものではなく、請求の範囲及び明細書全体から読み取れる発明の要旨、或いは思想に反しない範囲で適宜変更可能であり、そのような変更を伴う動状態数値化システム及び方法、並びにコンピュータプログラムもまた、本発明の技術的範囲に含まれるものである。
本発明の第1実施形態に係る動状態解析システムのブロック図である。 第1実施形態の、時間軸上で並ぶ複数のフレーム画像における、座標軸及び画素の配列の一例を図式的に示す概念図である。 第1実施形態における解析用画像の一例を図式的に示す概念図である。 第1実施形態の、時間軸上で並ぶ複数のフレーム画像における、座標軸及び画素の配列の他の例を図式的に示す概念図である。 第1実施形態における解析用画像の他の例を図式的に示す概念図である。 図6(a)は、車載カメラにより撮像される動画を示す平面図であり、図6(b)及び(c)は夫々、その動画について解析用画像生成部により生成された解析用画像を概念的に示す概念図である。 第1実施形態の動状態解析システムの動作を示すフローチャートである。 本発明の第2実施形態に係る動状態解析システムのブロック図である。 第2実施形態の、数量軸上で並ぶ複数の断面画像における、座標軸及び画素の配列の一例を図式的に示す概念図である。 第2実施形態における解析用画像の一例を図式的に示す概念図である。 第2実施形態の、数量軸上で並ぶ複数の断面画像における、座標軸及び画素の配列の他の例を図式的に示す概念図である。 第2実施形態における解析用画像の他の例を図式的に示す概念図である。 図13(a)は、本実施形態により生成された解析用画像を示す平面図であり、図13(b)は、比較例における同趣旨の平面図である。 本発明の第1実施例のシステムのブロック図である。 第1実施例のシステムにより撮像される細胞サンプルが動く様子を概念的に示すx−y平面図である。 第1実施例により生成される、x軸及びy軸に夫々沿った方向の時間的変化を示す解析用画像である。 本発明の第2実施例のシステムのブロック図である。 第2本実施例により生成される、x軸及びy軸に夫々沿った方向の時間的変化を示す解析用画像である。 本発明の第3実施例のシステムのブロック図である。 第3実施例により生成される、x軸及びy軸に夫々沿った方向の時間的変化を示す解析用画像である。 本発明の第4実施例のシステムのブロック図である。 第4実施例のシステムにより生成される解析用画像である、x−z層断面及びy−z層断面である。
符号の説明
10…細胞の活性度評価システム、20…自動操縦用動状態解釈システム、30…手話入力システム、40…層断面計測システム、100…動状態解析システム、101…撮像装置、100c…制御部、102…フレームバッファ、111…差分演算部、112…差分統計値演算部、113…解析用画像生成部、114…パターン解析部

Claims (5)

  1. 時間軸上で所定間隔をおいて並ぶ少なくとも二枚の画像の差分を、画素単位で演算する差分演算手段と、
    前記画像上における少なくとも一方向に交わる配列をなす複数の画素についての、前記演算された差分の平均値、前記演算された差分の最大値及び前記演算された差分が閾値以上となる画素の個数のうち少なくとも一つである、差分統計値を、前記一方向に沿った座標軸上の各点に対応付けて演算する差分統計値演算手段と、
    該演算された差分統計値を、前記座標軸と前記時間軸とで張られた面に置くことで、解析用画像を生成する解析用画像生成手段と
    を備えることを特徴とする動状態解析システム。
  2. 前記生成された解析用画像を、パターン解析するパターン解析手段を更に備えることを特徴とする請求項1に記載の動状態解析システム。
  3. 動画を撮像して複数の画像を生成する撮像手段を更に備え、
    前記差分演算手段は、前記複数の画像に含まれる前記少なくとも二枚の画像の差分を演算する
    ことを特徴とする請求項1又は2に記載の動状態解析システム。
  4. 時間軸上で所定間隔をおいて並ぶ少なくとも二枚の画像の差分を、画素単位で演算する差分演算工程と、
    前記画像上における少なくとも一方向に交わる配列をなす複数の画素についての、前記演算された差分の平均値、前記演算された差分の最大値及び前記演算された差分が閾値以上となる画素の個数のうち少なくとも一つである、差分統計値を、前記一方向に沿った座標軸上の各点に対応付けて演算する差分統計値演算工程と、
    該演算された差分統計値を、前記座標軸と前記時間軸とで張られた面に置くことで、解析用画像を生成する解析用画像生成工程と
    を備えることを特徴とする動状態解析方法。
  5. コンピュータを、請求項1から3のいずれか一項に記載の動状態解析システムとして機能させることを特徴とするコンピュータプログラム。
JP2006122555A 2006-04-26 2006-04-26 動状態解析システム及び方法、並びにコンピュータプログラム Expired - Fee Related JP3905117B1 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2006122555A JP3905117B1 (ja) 2006-04-26 2006-04-26 動状態解析システム及び方法、並びにコンピュータプログラム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2006122555A JP3905117B1 (ja) 2006-04-26 2006-04-26 動状態解析システム及び方法、並びにコンピュータプログラム

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2006222396A Division JP3900368B1 (ja) 2006-08-17 2006-08-17 層断面計測システム及び方法、並びにコンピュータプログラム

Publications (2)

Publication Number Publication Date
JP3905117B1 true JP3905117B1 (ja) 2007-04-18
JP2007293720A JP2007293720A (ja) 2007-11-08

Family

ID=38036293

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2006122555A Expired - Fee Related JP3905117B1 (ja) 2006-04-26 2006-04-26 動状態解析システム及び方法、並びにコンピュータプログラム

Country Status (1)

Country Link
JP (1) JP3905117B1 (ja)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4807354B2 (ja) * 2007-12-25 2011-11-02 住友電気工業株式会社 車両検出装置、車両検出システム及び車両検出方法
JP6750183B2 (ja) 2016-09-01 2020-09-02 公立大学法人会津大学 画像距離算出装置、画像距離算出方法および画像距離算出用プログラム
WO2018043437A1 (ja) * 2016-09-01 2018-03-08 公立大学法人会津大学 画像距離算出装置および画像距離算出用プログラムを記録したコンピュータ読み取り可能な非一過性の記録媒体
KR20230096994A (ko) * 2020-10-26 2023-06-30 스미도모쥬기가이고교 가부시키가이샤 표시방법, 제어장치, 표시장치, 및 프로그램

Also Published As

Publication number Publication date
JP2007293720A (ja) 2007-11-08

Similar Documents

Publication Publication Date Title
JP5997039B2 (ja) 欠陥検査方法および欠陥検査装置
US5892854A (en) Automatic image registration using binary moments
JP2004117078A (ja) 障害物検出装置及び方法
JPH11331609A (ja) ラスタ化された画像デ―タ内のエイリアスが除去された画素の検知方法
JP4872769B2 (ja) 路面判別装置および路面判別方法
JP3905117B1 (ja) 動状態解析システム及び方法、並びにコンピュータプログラム
US20180275414A1 (en) Display device and display method
JP2006138814A (ja) 液面検出方法
CN1677096A (zh) 印刷电路板检查装置
JP2000155839A (ja) 標識領域画像抽出方法、装置、標識領域画像抽出プログラムを記録した記録媒体
JP2007292715A (ja) 層断面計測システム及び方法、並びにコンピュータプログラム
JPH11248431A (ja) 3次元モデル生成装置および3次元モデル生成プログラムを記録したコンピュータ読取可能な記録媒体
JPH09223218A (ja) 走路検出方法及び装置
CN115409861A (zh) 基于图像处理的激光雷达地面点云分割方法、装置、设备及介质
JP2001056438A (ja) 共焦点走査型顕微鏡、共焦点走査型顕微鏡における画像表示方法および画像表示処理をするための処理プログラムを記録した記録媒体
JP2008293504A (ja) 対象物認識装置
JPH08210811A (ja) 円ないし楕円の検出方法およびその装置
JPH0520593A (ja) 走行レーン認識装置と先行車認識装置
JP4116597B2 (ja) バーコード認識装置、バーコード認識方法およびバーコード認識プログラム。
JP2009270984A (ja) 位置検出装置、位置検出プログラム、位置検出方法
JP7185839B2 (ja) 画像処理装置、画像処理方法及び画像処理プログラム
JP2589837B2 (ja) 合焦点条件検出のための処理領域決定方法
JP2006107018A (ja) 画像解析方法及び装置、画像処理方法及びシステム、これらの動作プログラム
CN114972005A (zh) 投影效果判断方法、设备、存储介质及程序产品
CN118037812A (zh) 一种中厚板板头不规则形状自适应测量方法和系统

Legal Events

Date Code Title Description
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: 20061226

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20070110

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

Year of fee payment: 4

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

Free format text: PAYMENT UNTIL: 20110119

Year of fee payment: 4

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

Free format text: PAYMENT UNTIL: 20120119

Year of fee payment: 5

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

Free format text: PAYMENT UNTIL: 20130119

Year of fee payment: 6

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

Free format text: PAYMENT UNTIL: 20130119

Year of fee payment: 6

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

Free format text: PAYMENT UNTIL: 20140119

Year of fee payment: 7

LAPS Cancellation because of no payment of annual fees