JP2009112332A - データ変換装置およびデータ変換プログラム - Google Patents

データ変換装置およびデータ変換プログラム Download PDF

Info

Publication number
JP2009112332A
JP2009112332A JP2007285282A JP2007285282A JP2009112332A JP 2009112332 A JP2009112332 A JP 2009112332A JP 2007285282 A JP2007285282 A JP 2007285282A JP 2007285282 A JP2007285282 A JP 2007285282A JP 2009112332 A JP2009112332 A JP 2009112332A
Authority
JP
Japan
Prior art keywords
data
time
sensitivity
frequency components
wavelet
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
JP2007285282A
Other languages
English (en)
Other versions
JP5062627B2 (ja
Inventor
Yoshiki Nakamura
芳樹 中村
Hiroshi Furuya
寛 古谷
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.)
Tokyo Institute of Technology NUC
Original Assignee
Tokyo Institute of Technology NUC
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 Tokyo Institute of Technology NUC filed Critical Tokyo Institute of Technology NUC
Priority to JP2007285282A priority Critical patent/JP5062627B2/ja
Publication of JP2009112332A publication Critical patent/JP2009112332A/ja
Application granted granted Critical
Publication of JP5062627B2 publication Critical patent/JP5062627B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Eye Examination Apparatus (AREA)

Abstract

【課題】本発明の目的は、刺激量の時間変化から感度の時間変化を定量的に予測するのに適したデータ変換装置及びデータ変換プログラムを提供することにある。
【解決手段】 本発明のデータ変換装置は、刺激量の時間変化データをウェーブレット分解することにより、前記時間変化データに含まれる複数の時間周波数成分を抽出する分解手段(S3)と、刺激量の時間変化とヒトの感度の時間変化との関係から予め定められた重み係数により、前記抽出された複数の時間周波数成分を重み付けする重み付け手段(S5)と、前記重み付けされた複数の時間周波数成分をウェーブレット合成することにより、前記時間変化データ通りに前記刺激量を時間変化させたときにヒトにもたらされる感度を予測する合成手段(S6)と、を備えたことを特徴とする。
【選択図】 図2

Description

本発明は、ヒトの感度の時間予測などに適用されるデータ変換装置及びデータ変換プログラムに関する。
また、本発明は、ヒトに対する刺激量の時間設計などに適用されるデータ変換装置及びデータ変換プログラムに関する。
ヒトの知覚には、与えられた刺激量の時間変化を検出するという特性がある。但し、或る時点におけるヒトの感度は、その時点で与えられた刺激量だけでなく、それ以前に与えられた刺激量にも影響される。
例えば、昼間にトンネルへ入った直後には周囲が見え難くいが、次第に見え易くなる(暗順応)。また、映画館から明るい外に出た直後には周囲が見え難いが、次第に見え易くなる(明順応)。
或いは、携帯型音楽プレーヤを大音量で使用し続けると音量が物足りなく感じる、塩味の強い料理を食べ続けてていると物足りなく感じる、帰宅直後には感じていた家の臭いが暫くすると感じ無くなる、テレビモニタの表示色が青味がかっていたとしても使用し始めるとその青味を感じ無くなる、といった感度が悪くなる方向への順応(慣れ)もある。
そこで近年、ヒトの感度の時間変化を定量的に予測することが望まれるようになってきた。例えば、ヒトの眼の感度の時間変化が予測できれば、施設内の照明輝度を時間設計することで、利用者に与える印象の時間変化を設計者が制御することも可能である。
特開2004−61150号公報
しかしながら現時点では、明るさ感の空間分布を定量的に予測する方法は提案されているものの(特許文献1などを参照)、明るさ感度の時間変化を定量的に予測する方法は提案されていない。
本発明の目的は、刺激量の時間変化から、その時間変化がヒトにもたらす感度を定量的に予測するのに適したデータ変換装置及びデータ変換プログラムを提供することにある。
また、本発明の目的は、感度の時間変化から、その時間変化をヒトに与えるための刺激量を定量的に予測するのに適したデータ変換装置及びデータ変換プログラムを提供することにある。
本発明のデータ変換装置は、刺激量の時間変化データをウェーブレット分解することにより、前記時間変化データに含まれる複数の時間周波数成分を抽出する分解手段と、刺激量の時間変化とヒトの感度の時間変化との関係から予め定められた重み係数により、前記抽出された複数の時間周波数成分を重み付けする重み付け手段と、前記重み付けされた複数の時間周波数成分をウェーブレット合成することにより、前記時間変化データ通りに前記刺激量を時間変化させたときにヒトにもたらされる感度を予測する合成手段と、を備えたことを特徴とする。
また、本発明の別のデータ変換装置は、ヒトの感度の時間変化データをウェーブレット分解することにより、前記時間変化データに含まれる複数の時間周波数成分を抽出する分解手段と、ヒトの感度の時間変化と刺激量の時間変化との関係から定められた重み係数により、前記抽出された複数の時間周波数成分を重み付けする重み付け手段と、前記重み付けされた複数の時間周波数成分をウェーブレット合成することにより、前記時間変化データ通りに前記感度を時間変化させるためにヒトに与えるべき刺激量を予測する合成手段と、を備えたことを特徴とする。
なお、本発明の何れかのデータ変換装置において、前記重み付け手段は、前記複数の時間周波数成分の成分量が正である場合と負である場合とで異なる値の重み係数を使用してもよい。
また、前記ウェーブレット分解及び前記ウェーブレット合成には、直交系のウェーブレットが用いられてもよい。
また、前記ウェーブレット分解及び前記ウェーブレット合成には、双直交系のウェーブレットが用いられてもよい。
また、前記刺激量は、ヒトの眼に呈示される光の輝度であり、前記感度は、ヒトの眼の輝度分解能であってもよい。
また、本発明のデータ変換プログラムは、刺激量の時間変化データをウェーブレット分解することにより、前記時間変化データに含まれる複数の時間周波数成分を抽出する分解手順と、刺激量の時間変化とヒトの感度の時間変化との関係から定められた重み係数により、前記抽出された複数の時間周波数成分を重み付けする重み付け手順と、前記重み付けされた複数の時間周波数成分をウェーブレット合成することにより、前記時間変化データ通りに前記刺激量を時間変化させたときにヒトにもたらされる感度を予測する合成手順と、をコンピュータに実行させることを特徴とする。
また、本発明の別のデータ変換プログラムは、ヒトの感度の時間変化データをウェーブレット分解することにより、前記時間変化データに含まれる複数の時間周波数成分を抽出する分解手順と、ヒトの感度の時間変化と刺激量の時間変化との関係から定められた重み係数により、前記抽出された複数の時間周波数成分を重み付けする重み付け手順と、前記重み付けされた複数の時間周波数成分をウェーブレット合成することにより、前記時間変化データ通りに前記感度を時間変化させるためにヒトに与えるべき刺激量を予測する合成手順とをコンピュータに実行させることを特徴とする。
なお、本発明の何れかのデータ変換プログラムにおいて、前記重み付け手順では、前記複数の時間周波数成分の成分量が正である場合と負である場合とで異なる値の重み係数を使用してもよい。
また、前記ウェーブレット分解及び前記ウェーブレット合成には、直交系のウェーブレットが用いられてもよい。
また、前記ウェーブレット分解及び前記ウェーブレット合成には、双直交系のウェーブレットが用いられてもよい。
また、前記刺激量は、ヒトの眼に呈示される光の輝度であり、前記感度は、ヒトの眼の輝度分解能であってもよい。
本発明によれば、刺激量の時間変化から、その時間変化がヒトにもたらす感度を定量的に予測するのに適したデータ変換装置及びデータ変換プログラムが実現する。
また、本発明によれば、感度の時間変化から、その時間変化をヒトに与えるための刺激量を定量的に予測するのに適したデータ変換装置及びデータ変換プログラムが実現する。
[第1実施形態]
以下、本発明の第1実施形態を説明する。本実施形態は、輝度の時間変化データを感度の時間変化データへ変換するデータ変換装置の実施形態である。
図1は、データ変換装置10の概略構成を示すブロック図である。図1に示すとおり、データ変換装置10は、データ変換プログラムがインストールされたコンピュータであり、処理対象である入力データや処理結果である出力データを記憶するメモリ10Aと、後述するフローチャート(図2)の手順に従ってデータ変換の演算処理を実行する演算部10Bとで構成される。なお、コンピュータにデータ変換プログラムをインストールするには、データ変換プログラムが記録された記録媒体(CD−ROMなど)を用いれば良い。又は、インターネットを介してダウンロード可能な搬送波(データ変換プログラムを含む)を用いても良い。
図2は、演算部10Bによる変換処理の手順を示すフローチャートである。以下、各ステップを順に説明する。
ステップS1:演算部10Bは、入力データをメモリ10Aに取り込む。この入力データは、輝度の時間変化データ(図3)であり、例えば、或る施設に装備された照明の輝度の時間変化データである。ここでは、この時間変化データの輝度軸の単位を輝度値の実数値とし、時間軸の単位を経過時間ではなく遡及時間とする。また、時間軸の刻みを0.025秒とし、データ数を113とする。
ステップS2:演算部10Bは、輝度の時間変化データ(図3)の輝度軸を、順応時における増分閾値の対数値に変換する。輝度の増分閾値とは、その輝度の下でヒトが知覚可能な輝度増加量の最小値のことであり、順応時における増分閾値とは、その輝度に順応しているヒトの増分閾値のことである。変換後の時間変化データは、図4に示すとおりである。
また、本ステップにおいて、演算部10Bは、図5に示すとおり、輝度の時間変化データの時間軸上に被変換点を定め、その被変換点の過去側の所定範囲のデータを参照する。以下、このデータを「参照データD(0)」とおく。なお、参照データD(0)のデータ範囲は、ヒトの眼が順応するのに必要な時間以上であり、例えば1.25秒である。ここでは時間軸の刻みが0.025秒なので、データ範囲を1.25秒とすると、参照データD(0)のデータ数は50となる。
ステップS3(1回目):演算部10Bは、図6(a)に示すとおり、参照データD(0)に対し分解レベル1のウェーブレット分解を施す。これによって、参照データD(0)に含まれる最も時間周波数の高い成分(分解レベル1のディテールデータDH(1))と、それ以外の成分(分解レベル1の近似データDL(1))とが得られる。このうちディテールデータDH(1)は、メモリ10Aに記憶される。なお、本ステップでは、ウェーブレット分解のマザーウェーブレットを、以下の式で表されるHaarとする。因みに、Haarは直交ウェーブレットの一種である。
Figure 2009112332
ステップS4:演算部10Bは、ウェーブレット分解の分解レベルが最終レベルに達したか否かを判別し、最終レベルに達していない場合はステップS3に戻り、最終レベルに達していた場合はステップS5へ移行する。なお、本実施形態では、参照データD(0)のデータ数が50なので、ウェーブレット分解の最終レベルを「6」とする。よって、ステップS3は6回繰り返される。
ステップS3(2回目):演算部10Bは、図6(b)に示すとおり、分解レベル1の近似データDL(1)に対し分解レベル2のウェーブレット分解を施す。これによって、近似データDL(1)に含まれる最も時間周波数の高い成分(分解レベル2のディテールデータDH(2))と、それ以外の成分(分解レベル2の近似データDL(2))とが得られる。このうちディテールデータDH(2)は、メモリ10Aに記憶される。
ステップS3(3回目):演算部10Bは、図6(c)に示すとおり、分解レベル2の近似データDL(2)に対し分解レベル3のウェーブレット分解を施す。これによって、近似データDL(2)に含まれる最も時間周波数の高い成分(分解レベル3のディテールデータDH(3))と、それ以外の成分(分解レベル3の近似データDL(3))とが得られる。このうちディテールデータDH(3)は、メモリ10Aに記憶される。
ステップS3(4回目):演算部10Bは、図6(d)に示すとおり、分解レベル3の近似データDL(3)に対し分解レベル4のウェーブレット分解を施す。これによって、近似データDL(3)に含まれる最も時間周波数の高い成分(分解レベル4のディテールデータDH(4))と、それ以外の成分(分解レベル4の近似データDL(4))とが得られる。このうちディテールデータDH(4)は、メモリ10Aに記憶される。
ステップS3(5回目):演算部10Bは、図6(e)に示すとおり、分解レベル4の近似データDL(4)に対し分解レベル5のウェーブレット分解を施す。これによって、近似データDL(4)に含まれる最も時間周波数の高い成分(分解レベル5のディテールデータDH(5))と、それ以外の成分(分解レベル5の近似データDL(5))とが得られる。このうちディテールデータDH(5)は、メモリ10Aに記憶される。
ステップS3(6回目):演算部10Bは、図6(f)に示すとおり、分解レベル5の近似データDL(5)に対し分解レベル6のウェーブレット分解を施す。これによって、近似データDL(5)に含まれる最も時間周波数の高い成分(分解レベル6のディテールデータDH(6))と、それ以外の成分(分解レベル6の近似データDL(6))とが得られる。これらのディテールデータDH(6)及び近似データDL(6)は、メモリ10Aに記憶される。
ステップS5:演算部10Bは、ディテールデータDH(1)〜DH(6)及び近似データDL(6)に対し係数処理を施すことによりディテールデータDH’(1)〜DH’(6)及び近似データDL’(6)を取得し、それらのディテールデータDH’(1)〜DH’(6)及び近似データDL’(6)をメモリ10Aに格納する。本ステップにおける係数処理は、以下の式で表される。
Figure 2009112332
この式中で各ディテールデータDH(1)〜DH(6)に乗算される重み係数を表にすると、図7のとおりとなる。図7に示すとおり、重み係数の値は、ディテールデータ毎(分解レベル毎)に用意されている。各分解レベルの重み係数の値は、輝度の時間変化データに含まれる各時間周波数成分がヒトの感度に与える影響度に相当する。
また、図7に示すとおり、重み係数の値は、ディテールデータの成分量が正であるときと負であるときとで異なるが、これは、同じ時間周波数成分であっても輝度が低から高へと変化するときと、高から低へと変化するときとでは、ヒトの眼の感度に与える影響度が異なることに起因する。
なお、図7に示す重み係数は、輝度の時間変化(図15(a))と、その時間変化が感度に及ぼす時間変化(図15(b))との関係に基づき予め決定されたものである。その決定方法については後に説明する。
ステップS6(1回目):演算部10Bは、図8(a)に示すとおり、分解レベル6のディテールデータDH’(6)及び近似データDL’(6)に対し合成レベル1のウェーブレット合成を施し、分解レベル5の近似データDL’(5)を取得する。なお、本ステップで使用されるマザーウェーブレットは、ステップS3で使用したのと同じウェーブレット(Haar)である。
ステップS7:演算部10Bは、合成レベルが最終レベルに達したか否かを判別し、最終レベルに達していない場合はステップS6に戻り、最終レベルに達していた場合はステップS8へ移行する。なお、上述したステップS3では分解レベルの最終レベルを「6」としたので、本ステップでも合成レベルの最終レベルを「6」とする。
ステップS6(2回目):演算部10Bは、図8(b)に示すとおり、分解レベル5のディテールデータDH’(5)及び近似データDL’(5)に対し合成レベル2のウェーブレット合成を施し、分解レベル4の近似データDL’(4)を取得する。
ステップS6(3回目):演算部10Bは、図8(c)に示すとおり、分解レベル4のディテールデータDH’(4)及び近似データDL’(4)に対し合成レベル3のウェーブレット合成を施し、分解レベル3の近似データDL’(3)を求める。
ステップS6(4回目):演算部10Bは、図8(d)に示すとおり、分解レベル3のディテールデータDH’(3)及び近似データDL’(3)に対し合成レベル4のウェーブレット合成を施し、分解レベル2の近似データDL’(2)を求める。
ステップS6(5回目):演算部10Bは、図8(e)に示すとおり、分解レベル2のディテールデータDH’(2)及び近似データDL’(2)に対し合成レベル5のウェーブレット合成を施し、分解レベル1の近似データDL’(1)を求める。
ステップS6(6回目):演算部10Bは、図8(f)に示すとおり、分解レベル1のディテールデータDH’(1)及び近似データDL’(1)に対し合成レベル6のウェーブレット合成を施し、分解レベル0の近似データDL’(0)を求める。そして、演算部10Bは、その近似データDL’(0)のうち被変換点に対応するデータを、その被変換点の感度データとしてメモリ10Aに記憶する。
ステップS8:演算部10Bは、全ての被変換点に関する変換処理(ステップS3〜S7)が終了したか否かを判別し、終了していなかった場合は被変換点を変更してからステップS3に戻り、終了していた場合はフローを終了する。
以上の結果、全ての被変換点の感度データがメモリ10Aに記憶される。全ての被変換点の感度データを時間順に並べると、図9に示すような感度の時間変化データが得られる。
この時間変化データの感度軸は、増分閾値(=ヒトが知覚可能な輝度増加量の最小値)によって表されているので、図9に示すカーブのうち落ち込んでいる箇所ほど感度が高く、膨らんでいる箇所ほど感度が低いとみなせばよい。
この感度の時間変化データの描くカーブ(図9)を、輝度の時間変化データの描くカーブ(図3)と比較すると明らかなとおり、照明が点灯された直後(t=50)には感度が急激に低下し、それから徐々に感度が回復していることがわかる。また、照明が再び消灯されると(t=25)、緩やかに感度が高まっていくことがわかる。
したがって、本実施形態のデータ変換装置10によれば、施設内の照明の輝度の時間変化から、その時間変化がヒトにもたらす感度の時間変化を、定量的に予測することができる。
なお、本実施形態のデータ変換装置10は、ウェーブレット分解とウェーブレット合成とに直交ウェーブレットを用いるので、輝度の時間変化データから感度の時間変化データへの変換誤差は非常に小さい。
したがって、本実施形態で取得した感度の時間変化データ(図9)を、同じ直交ウェーブレットを用いて逆変換すれば、元の輝度の時間変化データ(図3)に戻すこともできる。つまり、本実施形態の変換方法を応用すれば、輝度の時間変化データと感度の時間変化データとの双方向変換も可能である(逆変換の詳細は第2実施形態で説明する)。
[第2実施形態]
以下、本発明の第2実施形態を説明する。本実施形態は、感度の時間変化データを輝度の時間変化データへ変換するデータ変換装置の実施形態である。本実施形態のデータ変換装置の構成は、第1実施形態のそれ(図1)と同じであり、動作のみが異なる。よって、ここでは動作のみを説明する。
図10は、演算部10Bによる変換処理の手順を示すフローチャートである。以下、各ステップを順に説明する。
ステップS21:演算部10Bは、入力データをメモリ10Aに取り込む。この入力データは、感度の時間変化データ(図11)であり、例えば、或る施設の照明設計者が利用者に与えようとしている感度の時間変化データである。ここでは、この時間変化データの感度軸の単位を輝度値の増分閾値の対数値とし、時間軸の単位を経過時間ではなく遡及時間とする。また、時間軸の刻みを0.025秒とし、データ数を113とする。
また、本ステップにおいて、演算部10Bは、感度の時間変化データ(図11)の時間軸上に被変換点を定め、その被変換点の過去側の所定範囲のデータを参照する。以下、このデータを「参照データD(0)」とおく。なお、参照データD(0)のデータ範囲は、ヒトの眼の感度が順応するのに必要な時間以上であり、例えば1.25秒である。ここでは時間軸の刻みが0.025秒なので、データ範囲を1.25秒とすると、参照データD(0)のデータ数は50となる。
ステップS22:演算部10Bは、第1実施形態のステップS3と同様のウェーブレット分解(図6参照)を行う。なお、本ステップでも、ウェーブレット分解のマザーウェーブレットを、Haarとする。
ステップS23:演算部10Bは、ウェーブレット分解の分解レベルが最終レベルに達したか否かを判別し、最終レベルに達していない場合はステップS22に戻り、最終レベルに達していた場合はステップS24へ移行する。なお、本実施形態でも、参照データD(0)のデータ数が50なので、ウェーブレット分解の最終レベルを「6」とする。よって、ステップS22は6回繰り返される。繰り返しの結果、分解レベル1〜6のディテールデータDH(1)〜DH(6)及び分解レベル6の近似データDL(6)が取得され、メモリ10Aに記憶される。
ステップS24:演算部10Bは、第1実施形態のステップS5と同様、ディテールデータDH(1)〜DH(6)及び近似データDL(6)に対し係数処理を施すことによりディテールデータDH’(1)〜DH’(6)及び近似データDL’(6)を取得し、それらのディテールデータDH’(1)〜DH’(6)及び近似データDL’(6)をメモリ10Aに格納する。但し、本ステップにおける係数処理は、第1実施形態の係数処理とは異なり、以下の式で表される。
Figure 2009112332
この式中で各ディテールデータDH(1)〜DH(6)に乗算される重み係数を表にすると、図12のとおりとなる。図12に示すとおり、重み係数の値は、ディテールデータ毎(分解レベル毎)に用意されている。本実施形態の変換処理は、第1実施形態の変換処理の逆変換に相当するので、図12を前述した図7と比較すると明らかなとおり、本実施形態の重み係数の値は、第1実施形態の重み係数の逆数となっている。
なお、図12に示す重み係数は、輝度の時間変化(図5(a))と、その時間変化が感度に及ぼす時間変化(図15(b))との関係に基づき予め決定されたものである。その決定方法については後に説明する。
ステップS25:演算部10Bは、第1実施形態のステップS6と同様のウェーブレット合成(図8参照)を行う。なお、本ステップで使用されるマザーウェーブレットは、ステップS22で使用されたのと同じウェーブレット(Haar)である。
ステップS26:演算部10Bは、合成レベルが最終レベルに達したか否かを判別し、最終レベルに達していない場合はステップS25に戻り、最終レベルに達していた場合はステップS27へ移行する。なお、上述したステップS22では分解レベルの最終レベルを「6」としたので、本ステップでも合成レベルの最終レベルを「6」とする。これによって、分解レベル0の近似データDL’(0)が求まり、その近似データDL’(0)のうち被変換点に対応するデータが、その被変換点の輝度データとしてメモリ10Aに記憶される。
ステップS27:演算部10Bは、全ての被変換点に関する変換処理(ステップS22〜S26)が終了したか否かを判別し、終了していなかった場合は被変換点を変更してからステップS22に戻り、終了していた場合はステップS28へ移行する。その結果、全ての被変換点の輝度データがメモリ10Aに記憶される。全ての被変換点の輝度データを時間順に並べると、例えば図13に示すような輝度の時間変化データが得られる(図13は概念図)。但し、この時間変化データの輝度軸は、順応時における増分閾値の対数値(=ヒトが知覚可能な輝度増加量の最小値)によって表されている。
ステップS28:演算部10Bは、輝度の時間変化データ(図13)の輝度軸の単位を輝度値に変換し、変換後の時間変化データ(図14)をメモリ10Aへ記憶してからフローを終了する(図14は概念図)。なお、この変換は、第1実施形態のステップS2における変換の逆変換に相当する。
得られた時間変化データの描くカーブ(図14)を、感度の時間変化データの描くカーブ(図11)と比較すると明らかなとおり、感度を高感度aから低感度bへと直線的に変化させようとしたならば、輝度を緩やかに上昇させつつ、必要以上の輝度にまで上昇させてから(図13(A))、必要な輝度に落とし込めばよい(図13(B))ことがわかる。
したがって、本実施形態のデータ変換装置10によれば、眼の感度の時間変化から、その時間変化を利用者にもたらすのに必要な照明の輝度の時間変化を、定量的に予測することができる。
[重み係数の決定方法]
以下、第1実施形態の重み係数(図7)の決定方法を説明する。
第1実施形態の重み係数(図7)の決定は、実測データによる重回帰分析によって行われる。使用される実測データは、輝度の時間変化の実測データ(図15(a))と、その時間変化がヒトの感度に及ぼす時間変化の実測データ(図15(b))とである。これらの実測データには、なるべく多くの時間周波数成分が含まれていることが望ましい。よって、実測データ(図15(a))のカーブは、例えばパルス状とされる。
また、実測データ(図15(a)、(b))の時間軸の刻み(時間分解能)は、第1実施形態の入力データのそれと同じである。また、重回帰分析で想定されるマザーウェーブレットの種類は、第1実施形態で使用されるマザーウェーブレットの種類と同じ(ここではHaar)である。
但し、重回帰分析で想定される参照データのデータ範囲は、重み係数の決定精度を上げるため、第1実施形態で定めた参照データのデータ範囲よりも十分に広いことが望ましい。完全な暗順応は数十分かかるとも言われているため、そのデータ範囲は、数十分又はそれ以上に設定されることが望ましい。また、参照データのデータ範囲が広いので、重回帰分析で想定される最終分解レベル(及び最終合成レベル)は、その分だけ第1実施形態の最終分解レベル(及び最終合成レベル)よりも高く設定される。
以上の条件で重回帰分析を行えば、第1実施形態の重み係数(図7)は、高精度に算出される。そして、第2実施形態の重み係数(図12)は、第1実施形態の重み係数(図7)の逆数をとるだけで簡単に算出することができる。
[変形例]
なお、上記した何れかの実施形態では、ウェーブレット分解及びウェーブレット合成のマザーウェーブレットにHaarを用いたが、他の直交系ウェーブレットや、双直交系ウェーブレットを使用してもよい。但し、マザーウェーブレットの種類が異なると、前述した重み係数(図7、図12)の値も異なる。
また、上記した何れかの実施形態では、参照データのデータ数を50としたので、最終分解レベル(最終合成レベル)を「6」としたが、参照データのデータ数と最終分解レベル(最終合成レベル)との組み合わせを他の組み合わせにしてもよい。因みに、データ数が増える、つまり最終分解レベル(最終合成レベル)が高くなるほど、使用すべき重み係数の個数も増える。
また、上記した第2実施形態では、参照データのデータ数及び最終分解レベル(最終合成レベル)の組み合わせが第1実施形態のそれと同じであったが、異なってもよい。
また、上記した第2実施形態では、入力データの時間分解能が第1実施形態のそれと同じであったが、異なってもよい。但し、その場合は、第2実施形態の重み係数が、第1実施形態の重み係数の逆数にはならないので、第1実施形態の重み係数とは別に第2実施形態の重み係数も重回帰分析で決定される必要がある。
また、上記した第2実施形態では、マザーウェーブレットの種類が第1実施形態のそれと同じであったが、異なってもよい。但し、その場合は、第2実施形態の重み係数が、第1実施形態の重み係数の逆数にはならないので、第1実施形態の重み係数とは別に第2実施形態の重み係数も重回帰分析で決定される必要がある。
また、上記した何れかの実施形態では、ウェーブレット合成に係数処理後のディテールデータの全部を使用したが、係数処理後のディテールデータの一部(例えば、成分量の小さいもの)を無視し、他の一部のみを使用してもよい。これによって、演算が簡略化される。但し、係数処理後のディテールデータの全てをウェーブレット合成に使用する方が、正確な変換を行うことができる。
また、上記した第1実施形態では、輝度の時間変化データを感度の時間変化データへ変換する前に、輝度軸を対数値に変換したが、対数変換の他に例えば1/3乗などのPower Functionを用いてもよい。Power Functionは、均等知覚空間の表現に合わせて設定すればよい。
また、上記した第2実施形態では、感度の時間変化データの縦軸を対数値とし、感度の時間変化データを輝度の時間変化データへ変換した後に、その輝度軸を実数値に変換したが、対数変換の他に例えば1/3乗などのPower Functionを用いてもよい。Power Functionは、均等知覚空間の表現に合わせて設定すればよい。
また、上記した第1実施形態では、輝度の時間変化データを感度の時間変化データへ変換するデータ変換プログラム(第1のデータ変換プログラム)を説明し、第2実施形態では、感度の時間変化データを輝度の時間変化データへ変換するデータ変換プログラム(第2のデータ変換プログラム)を説明したが、輝度の時間変化データと感度の時間変化データとの間の双方向変換が可能なデータ変換プログラムを実現してもよい。
また、上述した変換処理を高速化すれば、データ変換装置10を次のように動作させることも可能である。すなわち、データ変換装置10は、輝度の時間変化データを示すカーブと、感度の時間変化データを示すカーブとを同時にモニタへ表示させる。この状態で、一方のカーブをユーザが任意に変形させると、データ変換装置10は、他方のカーブをリアルタイムで変形させる。このような動作プログラムは、照明設計に使用されるCG作成ソフトウエアなどに組み込まれるとよい。
このCG作成ソフトウエアを使用すれば、設計者は、照明器具の出力の時間変化を、自分の思い描いた通りに設計することができる。
また、上記した何れかの実施形態では、データ変換プログラムがインストールされたコンピュータをデータ変換装置10としたが、データ変換装置10を専用のハードウェア(LSI)によりチップ化してもよい。チップ化することにより、上述した変換処理を高速に行うことが容易となる。
また、上記した何れかの本実施形態で説明したデータ変換は、視覚と輝度との間に適用されるものであったが、特定の色感度と特定の色度との間や、特定の味覚と特定の味成分量との間や、温熱感覚と温度との間や、特定の嗅覚と特定の臭い成分量と間にも適用することができる。以下、幾つかの応用例を示す。
(1)視覚
・美術館のアプローチの各位置に設置された照明の輝度を、展示室内に入った直後に必要な感度が利用者に与えられるよう設計する。
・美術館で照明の明るい美術品を先に見せつつ、美術品と美術品との間に暗い部分を設け、その暗い部分の頻度や輝度などを適切に設計することにより、照明の暗い美術品に利用者が到達するまでの間に順応を図る。
・違和感のない色光の時間変化パターンを設計する。或いは、注意を引きやすい色光の時間変化パターンを設計する。
・自動車の運転時、夜間のスポーツ時などに、事故の生じ易い状況を予測する。
・自動車の運転時、夜間のスポーツ時などに、事故を防止するために照明を時間設計する。
(2)聴覚
・携帯オーディオの音量を、使用者に気づかれない程度の時間変化パターンで下げることにより、使用感を低下させずに事故や騒音を防ぐ。
(3)味覚
・塩味の濃いものばかりを食べると薄味のものに満足できなくなる様子を定量的に予測する。
・食事内容に不満感を与えずに生活習慣病を防止できるような食事計画をする。
(4)温熱感覚
・冷暖房機の自動温度制御において、夏は違和感を与えにくいような時間変化パターンで温度を上げ、冬は違和感を与えにくいような時間変化パターンで温度を下げることにより、省エネルギーを図る。
(5)嗅覚
・家の臭いが気にならなくなる様子を定量的に予測する。
・注意を引きやすい臭いの変化パターンを予測する。
・麻薬捜査犬の嗅覚をベストに保つための環境整備をする。
データ変換装置10の概略構成を示すブロック図である。 第1実施形態の変換処理の手順を示すフローチャートである。 第1実施形態の入力データ(輝度の時間変化データ)を示す図である。 輝度軸の変換後の時間変化データを示す図である。 参照データを説明する図である。 ウェーブレット分解の説明図である。 第1実施形態の重み係数を示す図である。 ウェーブレット合成を説明する図である。 第1実施形態の出力データ(感度の時間変化データ)を示す図である。 第2実施形態の変換処理の手順を示すフローチャートである。 第2実施形態の入力データ(感度の時間変化データ)の概念図である。 第2実施形態の重み係数を示す図である。 輝度の時間変化データの概念図である。 第2実施形態の出力データ(輝度の時間変化データ)の概念図である。 重回帰分析で使用される実測データの一例を示す図である。

Claims (12)

  1. 刺激量の時間変化データをウェーブレット分解することにより、前記時間変化データに含まれる複数の時間周波数成分を抽出する分解手段と、
    刺激量の時間変化とヒトの感度の時間変化との関係から予め定められた重み係数により、前記抽出された複数の時間周波数成分を重み付けする重み付け手段と、
    前記重み付けされた複数の時間周波数成分をウェーブレット合成することにより、前記時間変化データ通りに前記刺激量を時間変化させたときにヒトにもたらされる感度を予測する合成手段と、
    を備えたことを特徴とするデータ変換装置。
  2. ヒトの感度の時間変化データをウェーブレット分解することにより、前記時間変化データに含まれる複数の時間周波数成分を抽出する分解手段と、
    ヒトの感度の時間変化と刺激量の時間変化との関係から定められた重み係数により、前記抽出された複数の時間周波数成分を重み付けする重み付け手段と、
    前記重み付けされた複数の時間周波数成分をウェーブレット合成することにより、前記時間変化データ通りに前記感度を時間変化させるためにヒトに与えるべき刺激量を予測する合成手段と、
    を備えたことを特徴とするデータ変換装置。
  3. 請求項1又は請求項2に記載のデータ変換装置において、
    前記重み付け手段は、
    前記複数の時間周波数成分の成分量が正である場合と負である場合とで異なる値の重み係数を使用する
    ことを特徴とするデータ変換装置。
  4. 請求項1〜請求項3の何れか一項に記載のデータ変換装置において、
    前記ウェーブレット分解及び前記ウェーブレット合成には、
    直交系のウェーブレットが用いられる
    ことを特徴とするデータ変換装置。
  5. 請求項1〜請求項3の何れか一項に記載のデータ変換装置において、
    前記ウェーブレット分解及び前記ウェーブレット合成には、
    双直交系のウェーブレットが用いられる
    ことを特徴とするデータ変換装置。
  6. 請求項1〜請求項5の何れか一項に記載のデータ変換装置において、
    前記刺激量は、ヒトの眼に呈示される光の輝度であり、
    前記感度は、ヒトの眼の輝度分解能である
    ことを特徴とするデータ変換装置。
  7. 刺激量の時間変化データをウェーブレット分解することにより、前記時間変化データに含まれる複数の時間周波数成分を抽出する分解手順と、
    刺激量の時間変化とヒトの感度の時間変化との関係から定められた重み係数により、前記抽出された複数の時間周波数成分を重み付けする重み付け手順と、
    前記重み付けされた複数の時間周波数成分をウェーブレット合成することにより、前記時間変化データ通りに前記刺激量を時間変化させたときにヒトにもたらされる感度を予測する合成手順と、
    をコンピュータに実行させることを特徴とするデータ変換プログラム。
  8. ヒトの感度の時間変化データをウェーブレット分解することにより、前記時間変化データに含まれる複数の時間周波数成分を抽出する分解手順と、
    ヒトの感度の時間変化と刺激量の時間変化との関係から定められた重み係数により、前記抽出された複数の時間周波数成分を重み付けする重み付け手順と、
    前記重み付けされた複数の時間周波数成分をウェーブレット合成することにより、前記時間変化データ通りに前記感度を時間変化させるためにヒトに与えるべき刺激量を予測する合成手順と、
    をコンピュータに実行させることを特徴とするデータ変換プログラム。
  9. 請求項7又は請求項8に記載のデータ変換プログラムにおいて、
    前記重み付け手順では、
    前記複数の時間周波数成分の成分量が正である場合と負である場合とで異なる値の重み係数を使用する
    ことを特徴とするデータ変換プログラム。
  10. 請求項7〜請求項9の何れか一項に記載のデータ変換プログラムにおいて、
    前記ウェーブレット分解及び前記ウェーブレット合成には、
    直交系のウェーブレットが用いられる
    ことを特徴とするデータ変換プログラム。
  11. 請求項7〜請求項9の何れか一項に記載のデータ変換プログラムにおいて、
    前記ウェーブレット分解及び前記ウェーブレット合成には、
    双直交系のウェーブレットが用いられる
    ことを特徴とするデータ変換プログラム。
  12. 請求項7〜請求項11の何れか一項に記載のデータ変換プログラムにおいて、
    前記刺激量は、ヒトの眼に呈示される光の輝度であり、
    前記感度は、ヒトの眼の輝度分解能である
    ことを特徴とするデータ変換プログラム。
JP2007285282A 2007-11-01 2007-11-01 データ変換装置およびデータ変換プログラム Expired - Fee Related JP5062627B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2007285282A JP5062627B2 (ja) 2007-11-01 2007-11-01 データ変換装置およびデータ変換プログラム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2007285282A JP5062627B2 (ja) 2007-11-01 2007-11-01 データ変換装置およびデータ変換プログラム

Publications (2)

Publication Number Publication Date
JP2009112332A true JP2009112332A (ja) 2009-05-28
JP5062627B2 JP5062627B2 (ja) 2012-10-31

Family

ID=40780175

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2007285282A Expired - Fee Related JP5062627B2 (ja) 2007-11-01 2007-11-01 データ変換装置およびデータ変換プログラム

Country Status (1)

Country Link
JP (1) JP5062627B2 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5351438B2 (ja) * 2008-05-23 2013-11-27 日本放送協会 ディスプレイ制御装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1049791A (ja) * 1996-08-05 1998-02-20 Nippon Telegr & Teleph Corp <Ntt> 時系列信号処理方法および装置
JP2003330912A (ja) * 2002-05-10 2003-11-21 National Institute Of Advanced Industrial & Technology 時系列信号の解析方法およびプログラム
JP2003334250A (ja) * 2002-05-20 2003-11-25 Matsushita Electric Ind Co Ltd 照明装置および照明制御装置
WO2006132014A1 (ja) * 2005-06-06 2006-12-14 Tokyo Institute Of Technology 画像変換装置および画像変換プログラム
JP2007142670A (ja) * 2005-11-16 2007-06-07 Olympus Corp 画像処理システム、画像処理プログラム

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1049791A (ja) * 1996-08-05 1998-02-20 Nippon Telegr & Teleph Corp <Ntt> 時系列信号処理方法および装置
JP2003330912A (ja) * 2002-05-10 2003-11-21 National Institute Of Advanced Industrial & Technology 時系列信号の解析方法およびプログラム
JP2003334250A (ja) * 2002-05-20 2003-11-25 Matsushita Electric Ind Co Ltd 照明装置および照明制御装置
WO2006132014A1 (ja) * 2005-06-06 2006-12-14 Tokyo Institute Of Technology 画像変換装置および画像変換プログラム
JP2007142670A (ja) * 2005-11-16 2007-06-07 Olympus Corp 画像処理システム、画像処理プログラム

Also Published As

Publication number Publication date
JP5062627B2 (ja) 2012-10-31

Similar Documents

Publication Publication Date Title
JP4016072B2 (ja) 画像変換装置および画像変換プログラム
WO2010018477A3 (en) Model enhanced imaging
Shin et al. Sound quality evaluation of the booming sensation for passenger cars
US20170313237A1 (en) An Illumination Standard Calculation Method And System for A Tunnel Middle Section Based On Safe Visual Recognition
JP2011503780A (ja) 照明雰囲気を自動的にレンダリングするための光制御システム及び方法
JP2012504827A5 (ja)
CN106782544A (zh) 语音交互设备及其输出方法
JP2015518974A5 (ja)
RU2011118457A (ru) Определение и/или представление показателя риска для здоровья
JP2007292665A (ja) グレア評価装置およびグレア評価プログラム
CN110709903B (zh) 觉醒引导控制装置以及觉醒引导系统
JP7262002B2 (ja) 覚醒誘導システム
JP5062627B2 (ja) データ変換装置およびデータ変換プログラム
CN114828359A (zh) 基于音乐的氛围灯展示方法、装置、设备及存储介质
EP1843299A3 (en) Method for automatic generation of finite element mesh from integrated circuit layout data
EP2151714A3 (en) Original data generation program and original data generation method
Wan et al. The speed control effect of highway tunnel sidewall markings based on color and temporal frequency
WO2012086548A1 (ja) 照明装置
KR20110076745A (ko) 졸음운전 경고 장치 및 그 방법
JP2006139782A5 (ja)
CN111989917A (zh) 电子设备及其控制方法
RU2012133293A (ru) Способ для моделирования снижения дозы при сборе данных рентгеновской системой, компьютерная система и рентгеновская система
US20100262569A1 (en) Data converting apparatus and medium having data converting program
JP2008093131A (ja) 眩しさ評価装置
CN116610212A (zh) 一种多模态娱乐交互方法、装置、设备及介质

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20101022

TRDD Decision of grant or rejection written
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20120621

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20120703

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20120731

R150 Certificate of patent or registration of utility model

Ref document number: 5062627

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

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

Free format text: PAYMENT UNTIL: 20150817

Year of fee payment: 3

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees