(57)【要約】【課題】時系列に変化する各音源の位置を一目瞭然にする音源位置表示装置を実現する。【解決手段】ＣＰＵ１０は、ＬチャンネルおよびＲチャンネルの音響波形データに周波数分析（高速フーリエ変換ＦＦＴ）を施して得たＬチャンネルおよびＲチャンネルのパワースペクトルから定位情報ｐａｎを取得し、取得した定位情報ｐａｎに応じて表示色を異ならせたパワースペクトルを曲頭から曲末尾までの時系列順に画面表示するので、時系列に変化する各音源の位置を一目瞭然にする。【選択図】図１
2014229965.tif 000002
【特許請求の範囲】
【請求項１】
音響波形データからパワースペクトルを取得するパワー取得手段と、
前記パワー取得手段により取得されたパワースペクトルから定位情報を取得する定位情報取得手段と、
前記定位情報取得手段により取得された定位情報に応じて表示色を異ならせたパワースペクトルを表示する表示手段と
を具備することを特徴とする音源位置表示装置。
【請求項２】
前記表示手段は、前記定位情報取得手段により取得された定位情報を色相に変換する色相変換手段を備え、当該色相変換手段により変換された色相で当該定位情報に対応するパワースペクトルを表示することを特徴とする請求項１記載の音源位置表示装置。
【請求項３】
前記色相変換手段は、前記定位情報取得手段により取得された定位情報に応じた三原色値（ｒ，ｇ，ｂ）を発生する色関数を備えることを特徴とする請求項１記載の音源位置表示装置。
【請求項４】
前記表示手段は、前記定位情報取得手段により取得された定位情報を明度に変換する明度変換手段を更に備えることを特徴とする請求項１記載の音源位置表示装置。
【請求項５】
前記表示手段は、前記定位情報取得手段により取得された定位情報を彩度に変換する彩度変換手段を更に備えることを特徴とする請求項１記載の音源位置表示装置。
【請求項６】
前記表示手段は、前記定位情報取得手段により取得された定位情報に応じて表示色を異ならせたパワースペクトルを時系列順に画面表示することを特徴とする請求項１記載の音源位置表示装置。
【請求項７】
前記表示手段は、前記定位情報取得手段により取得された定位情報が、所定範囲に含まれる場合に、前記色相変換手段により変換された色相を補正する色補正手段を更に備えることを特徴とする請求項１記載の音源位置表示装置。
【請求項８】
前記色補正手段は、前記定位情報取得手段により取得された定位情報が、特定の範囲に含まれる場合に、前記色相変換手段により変換された色相を強調することを特徴とする請求項７記載の音源位置表示装置。
【請求項９】
音響波形データからパワースペクトルを取得し、取得したパワースペクトルから定位情報を取得し、取得した定位情報に応じて表示色を異ならせたパワースペクトルを表示することを特徴とする音源位置表示方法。
【請求項１０】
コンピュータに、
音響波形データからパワースペクトルを取得するパワー取得ステップと、
前記パワー取得ステップで取得されたパワースペクトルから定位情報を取得する定位情報取得ステップと、
前記定位情報取得ステップで取得された定位情報に応じて表示色を異ならせたパワースペクトルを表示する表示ステップと
を実行させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ステレオ音場における音源の位置を一目瞭然にする音源位置表示装置、音源位置表示方法およびプログラムに関する。
【背景技術】
【０００２】
ステレオ音場における音源の位置を可視化表示する技術が知られている。例えば特許文献１には、入力楽音信号から抽出される各周波数帯域毎の定位情報で定まる２次元平面（定位軸−周波数軸平面）において、各周波数帯域毎のレベル分布を、当該２次元平面の高さ方向となるレベル軸で３次元座標で表現し、そのレベル軸（高さ方向）を区別可能な態様で画面表示させることによってボーカルや楽器などの音源の存在領域を視認させる技術が開示されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１１−１５８６７４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ところで、上記特許文献１に開示の技術では、２次元平面（定位軸−周波数軸平面）の高さ方向で各周波数帯域毎のレベル分布を表現するので、ボーカルや楽器単位の存在領域がグループ化した状態で画面表示される。この為、入力される音響信号に従って時系列に変化する各音源の位置を一目瞭然にすることが出来ないという問題が有る。
【０００５】
本発明は、このような事情に鑑みてなされたもので、時系列に変化する各音源の位置を一目瞭然にすることが出来る音源位置表示装置、音源位置表示方法およびプログラムを提供することを目的としている。
【課題を解決するための手段】
【０００６】
上記目的を達成するため、本発明の音源位置表示装置は、音響波形データからパワースペクトルを取得するパワー取得手段と、前記パワー取得手段により取得されたパワースペクトルから定位情報を取得する定位情報取得手段と、前記定位情報取得手段により取得された定位情報に応じて表示色を異ならせたパワースペクトルを表示する表示手段とを具備することを特徴とする。
【０００７】
また、本発明の音源位置表示方法では、音響波形データからパワースペクトルを取得し、取得したパワースペクトルから定位情報を取得し、取得した定位情報に応じて表示色を異ならせたパワースペクトルを表示することを特徴とする。
【０００８】
更に、本発明のプログラムでは、コンピュータに、音響波形データからパワースペクトルを取得するパワー取得ステップと、前記パワー取得ステップで取得されたパワースペクトルから定位情報を取得する定位情報取得ステップと、前記定位情報取得ステップで取得された定位情報に応じて表示色を異ならせたパワースペクトルを表示する表示ステップとを実行させることを特徴とする。
【発明の効果】
【０００９】
本発明では、時系列に変化する各音源の位置を一目瞭然にすることが出来る。
【図面の簡単な説明】
【００１０】
【図１】本発明の第１実施形態による音源位置表示装置１００の全体構成を示すブロック図である。
【図２】ＲＡＭ１２のワークエリアＷＥの構成を示すメモリマップである。
【図３】色相テーブルＣＴの一例を示す図である。
【図４】色相テーブルＣＴの他の一例を示す図である。
【図５】色相テーブルＣＴの他の一例を示す図である。
【図６】メインルーチンの動作を示すフローチャートである。
【図７】分析処理の動作を示すフローチャートである。
【図８】表示処理の動作を示すフローチャートである。
【図９】色値算出処理の動作を示すフローチャートである。
【図１０】色値算出処理の動作を示すフローチャートである。
【図１１】第１実施形態における音源位置表示の一例を示す図である。
【図１２】第２実施形態におけるＲＡＭ１２のワークエリアＷＥの構成を示すメモリマップである。
【図１３】第２実施形態の色補正テーブルＲＴの一例を示す図である。
【図１４】第２実施形態の色補正テーブルＲＴの他の一例を示す図である。
【図１５】第２実施形態による表示処理の動作を示すフローチャートである。
【図１６】第２実施形態による色補正処理の動作を示すフローチャートである。
【図１７】第２実施形態における音源位置表示の一例を示す図である。
【発明を実施するための形態】
【００１１】
以下、図面を参照して本発明の実施形態について説明する。
［第１実施形態］
Ａ．構成
図１は、本発明の第１実施形態による音源位置表示装置１００の全体構成を示すブロック図である。この図に示す音源位置表示装置１００は、ＣＰＵ１０、ＲＯＭ１１、ＲＡＭ１２、入力手段１３、表示手段１４、音響出力手段１５および音響入力手段１６を備える。
【００１２】
ＣＰＵ１０は、入力手段１３を介してユーザにより入力されるコマンドに応じて装置各部を制御する。ユーザにより入力されるコマンドの種類には、処理終了を指示する終了コマンド、再生機能の実行を指示する再生コマンド、音響入力機能の実行を指示する音響入力コマンド、演奏入力機能の実行を指示する演奏入力コマンド、本発明の要旨に係わる音源位置表示機能の実行を指示する音源位置表示コマンドがある。
【００１３】
再生機能とは、後述する音響入力機能により取り込んだ音響波形データや、後述する演奏入力機能により取り込んだ演奏データ（ＭＩＤＩデータ）を、音響出力手段１５から再生出力させる機能である。音響入力機能とは、音響入力手段１４によって、ステレオ音場として収音したステレオ成分の音響信号をサンプリングし、これにより得られる音響波形データ（Ｌチャンネル、Ｒチャンネル）をＲＡＭ１２の波形データエリアに保存する機能である。演奏入力機能とは、入力手段１３を介して外部ＭＩＤＩ楽器から入力される演奏データ（ＭＩＤＩデータ）を取り込んでＲＡＭ１２の演奏データエリアに保存する機能である。
【００１４】
音源位置表示機能とは、音響入力機能によってＲＡＭ１２に保存された音響波形データ（Ｌチャンネル、Ｒチャンネル）に周波数分析（高速フーリエ変換ＦＦＴ）を施してパワースペクトル（ＬチャンネルおよびＲチャンネル）を取得すると共に、ＬチャンネルおよびＲチャンネルのパワースペクトルから定位情報ｐａｎを取得し、取得した定位情報ｐａｎに応じてパワースペクトルの表示色を変化させ、それを時間軸（横軸）および周波数軸（縦軸）とした表示画面に表示する機能である。こうして、定位情報ｐａｎに応じて表示色が変化するパワースペクトルを時系列順に画面表示することで時系列に変化する各音源の位置を一目瞭然にすることが可能になる。
【００１５】
ＲＯＭ１１には、上述した再生機能、音響入力機能、演奏入力機能および音源位置表示機能を具現する各種プログラム（後述するメインルーチン、分析処理、表示処理および色値算出処理を含む）が記憶される。ＲＡＭ１２は、音響入力機能によって取り込んだ音響波形データ（Ｌチャンネル、Ｒチャンネル）を保存する音響波形データや、演奏入力機能によって取り込んだ演奏データ（ＭＩＤＩデータ）を保存するデータエリア（不図示）の他、ＣＰＵ１０の作業エリアとして各種レジスタ・フラグデータを一時記憶するワークエリアＷＥを備える。
【００１６】
ここで、図２を参照して、ＲＡＭ１２が備えるワークエリアＷＥに一時記憶される主要なレジスタデータ（変数）について説明する。図２は、主に音源位置表示機能を具現する際に使われるレジスタデータ（変数）である。レジスタｄｏＤａｔａ［ｎ］には、１フレーム分のパワースペクトル値がストアされる。但し、レジスタｄｏＤａｔａ［ｎ］の前半（０〜ｉＦｒａｍｅ２−１）には、Ｒチャンネルのパワースペクトル値がストアされ、レジスタｄｏＤａｔａ［ｎ］の後半（ｉＦｒａｍｅ２〜ｉＦｒａｍｅ２×２）には、Ｌチャンネルのパワースペクトル値がストアされる。
【００１７】
レジスタｄには、モノラル化したパワースペクトル値がストアされる。レジスタｄｆｒには、ポインタｊに対応する周波数（後述する）がストアされる。レジスタｉＦｒａｍｅ２には、フレーム長さ（分析窓幅）の半分の長さがストアされる。レジスタｄｏＦｓには、音響波形データ（Ｌチャンネル、Ｒチャンネル）のサンプリング周波数がストアされる。レジスタｔｉｍｅｎには、時間軸における表示幅（表示領域の幅）に相当するフレーム数がストアされる。レジスタｉＨｅｉｇｈｔには、表示領域の高さ（周波数を対数表示する縦軸の範囲）がストアされる。
【００１８】
レジスタｗｉｄｔｈには、時間軸の表示幅（表示ピクセル数）がストアされる。レジスタｒｅｓ１００ｃｅｎｔには、表示分解能に相当する半音分の表示ピクセル数がストアされる。レジスタｐａｎには、左右両チャンネルのパワースペクトルの差分を、左右両チャンネルのパワースペクトルの和で乗算して得られる定位情報がストアされる。以下、レジスタｐａｎの値を定位情報ｐａｎと称す。定位情報ｐａｎは、左端Ｌｅｆｔで「−１．０」、中央Ｃｅｎｔｅｒで「０．０」、右端Ｒｉｇｈｔで「１．０」となる。
【００１９】
レジスタｐａｎ、ｐｅａｋＰｏｓＳｈｆｉｔＲ、ｐｅａｋＰｏｓＳｈｉｆｔＬ、ｂｒｉｇｈｔｎｅｓｓ、ｂｒｉｇｈｔｎｅｓｓＯｆｆｓｅｔは、色相テーブルＣＴの色相特性を設定する。色相テーブルＣＴとは、パワースペクトルをカラー表示する際の色（色相）を、定位情報ｐａｎに応じて変換する。すなわち、定位情報ｐａｎを読み出しアドレスとして、対応する色値（ｒ（赤），ｇ（緑），ｂ（青））を読み出すデータテーブルである。
【００２０】
ここで、図３に図示する色相テーブルＣＴの一例を参照して、当該色相テーブルＣＴの色相特性を設定するレジスタｐｅａｋＰｏｓＳｈｆｉｔＲ、ｐｅａｋＰｏｓＳｈｉｆｔＬ、ｂｒｉｇｈｔｎｅｓｓ、ｂｒｉｇｈｔｎｅｓｓＯｆｆｓｅｔについて説明する。図３に図示する色相テーブルＣＴには、定位情報ｐａｎに応じて色値（ｒ，ｇ，ｂ）を決める色関数（Ｒ関数、Ｇ関数、Ｂ関数）が設けられる。なお、図３に図示する色値（ｒ，ｇ，ｂ）は、「１」を最大値として正規化している。
【００２１】
図３に図示する一例の場合、Ｒ関数は中央（定位情報ｐａｎ＝０）の色値ｒを最大値「１」とする２つの直線から構成される。すなわち左端（定位関数ｐａｎ＝−１．０）においてレジスタｂｒｉｇｈｔｎｅｓｓで指定される値（もしくはレジスタｂｒｉｇｈｔｎｅｓｓＯｆｆｓｅｔで指定されるオフセット値）から中央（定位情報ｐａｎ＝０）の色値ｒ「１」に至る傾きを有する第１の直線と、右端（定位情報ｐａｎ＝１．０）においてレジスタｂｒｉｇｈｔｎｅｓｓで指定される値（もしくはレジスタｂｒｉｇｈｔｎｅｓｓＯｆｆｓｅｔで指定されるオフセット値）から中央（定位情報ｐａｎ＝０）の色値ｒ「１」に至る傾きを有する第２の直線とから構成される。
【００２２】
Ｇ関数は、レジスタｐｅａｋＰｏｓＳｈｆｉｔＲの値で指定される定位位置の色値ｇを最大値「１」とし、Ｒ関数と同じ傾きを有する第１および第２の直線から構成される。なお、Ｇ関数を構成する第１の直線は、色値ｇ（最大値）「１」から「０」までの値域に対応する定義域（定位情報ｐａｎ）を有する。一方、Ｇ関数を構成する第２の直線は、レジスタｐｅａｋＰｏｓＳｈｆｉｔＲの値で指定される定位位置から右端（定位情報ｐａｎ＝１．０）までの定義域に対応する値域を有する。
【００２３】
Ｂ関数は、レジスタｐｅａｋＰｏｓＳｈｆｉｔＬの値で指定される定位位置の色値ｂを最大値「１」とし、Ｒ関数と同じ傾きを有する第１および第２の直線から構成される。なお、Ｂ関数を構成する第１の直線は、レジスタｐｅａｋＰｏｓＳｈｆｉｔＬの値で指定される定位位置から左端（定位情報ｐａｎ＝−１．０）までの定義域に対応する値域を有する。一方、Ｂ関数を構成する第２の直線は、色値ｂ（最大値）「１」から「０」までの値域に対応する定義域（定位情報ｐａｎ）を有する。
【００２４】
さて、このような色相特性を有する色相テーブルＣＴにおいて、例えば入力する定位情報ｐａｎの値が「０．５」の場合には、図３に図示する当該色相テーブルＣＴから色値（ｒ１，ｇ１，ｂ１）が出力される。色相テーブルＣＴの色相特性は、分析対象となる音楽（音響信号）の種類に応じて調整する。具体的には、Ｇ関数およびＢ関数の最大値位置と第１および第２の直線の傾きを異ならせる。Ｇ関数およびＢ関数の最大値位置を変えることでパワースペクトルの表示色を異ならせ、これにより音源位置を浮かび上がらせるように表示することが可能になる。
【００２５】
Ｇ関数およびＢ関数における第１および第２の直線の傾きは、音源の広がりに合わせて調整する。第１および第２の直線の傾きを急峻にすると、ピアノなど広がり幅が大きい音源を捉えきれなくなり、一方、なだらかな傾きにすれば、トランペットなど広がり幅の狭い音源を捉え難くなる為、ある程度の広がり幅の音源を捉えつつ広がり幅の狭い音源も捕捉可能な傾きを経験的に見つけることになる。
【００２６】
例えば、図３に図示した色相テーブルＣＴよりも、ＲＧＢ各関数を構成する第１および第２の直線の傾きを大きくした場合には、図４に図示する色相特性の色相テーブルＣＴとなる。さらに、図４に図示する色相テーブルＣＴにおいて、Ｇ関数の最大値を左端（定位情報ｐａｎ＝−１．０）、Ｂ関数の最大値を右端（定位情報ｐａｎ＝１．０）に設定した場合には、図５に図示する色相特性の色相テーブルＣＴとなる。
【００２７】
次に、再び図１を参照して実施形態の構成について説明を進める。図１において、入力手段１３は、前述した各種コマンド（終了コマンド、再生コマンド、音響入力コマンド、演奏入力コマンドおよび音源位置表示コマンド）を入力する各種操作子や、演奏入力用の鍵盤およびＭＩＤＩインタフェース（不図示）を備える。
【００２８】
表示手段１４は、ＣＰＵ１０から供給される表示制御信号に基づき、定位情報ｐａｎに応じて表示色が変化するパワースペクトルを時系列順に画面表示する。なお、表示画面は、後述するように、横軸を時間軸、縦軸を対数周波数軸とする。時間軸の長さは分析対象とした音楽一曲分に相当し、対数周波数軸の範囲は、８８鍵のピアノにおける半音分の周波数を分解能ｒｅｓ１００ｃｅｎｔとする為、都合８８×分解能ｒｅｓ１００ｃｅｎｔ以上の表示領域となる。
【００２９】
音響出力手段１５は、音源装置、Ｄ／Ａ変換器、ノイズ除去用のフィルタ、アンプおよびステレオ再生用スピーカから構成される。この音響出力手段１５が音響波形データを再生する場合には、ＣＰＵ１０の制御の下にＲＡＭ１２のデータエリアから読み出される音響波形データ（Ｌチャンネル、Ｒチャンネル）をＤ／Ａ変換器にてアナログ信号形式に変換した後、ノイズ除去用のフィルタおよびアンプを介してステレオ再生用スピーカから放音出力する。一方、音響出力手段１５が演奏データを再生する場合には、ＣＰＵ１０の制御の下にＲＡＭ１２のデータエリアから読み出される演奏データに基づき音源装置が楽音波形データを発生し、この楽音波形データをＤ／Ａ変換器にてアナログ信号形式に変換した後、ノイズ除去用のフィルタおよびアンプを介してステレオ再生用スピーカから放音出力する。
【００３０】
音響入力手段１６は、ステレオ録音用のマイク、アンプおよびＡ／Ｄ変換器から構成される。音響入力手段１６は、ＣＰＵ１０の制御の下に、ステレオ音場として収音したステレオ成分の音響信号をサンプリングし、これにより得られる音響波形データ（Ｌチャンネル、Ｒチャンネル）をＲＡＭ１２の波形データエリアに記憶する。
【００３１】
Ｂ．動作
次に、図６〜図１１を参照して、上記構成による音源位置表示装置１００のＣＰＵ１０が実行するメインルーチン、分析処理、表示処理および色値算出処理の各動作を説明する。
【００３２】
（１）メインルーチンの動作
図６は、メインルーチンの動作を示すフローチャートである。音源位置表示装置１００がパワーオン（電源投入）されると、ＣＰＵ１０は図６に図示するメインルーチンを実行し、ステップＳＡ１に進み、装置各部を初期化する。続いて、ステップＳＡ２〜ＳＡ６では、入力手段１３を介してユーザ操作により入力されるコマンドの種類（終了コマンド、再生コマンド、音響入力コマンド、演奏入力コマンドおよび音源位置表示コマンド）を判別する。コマンド入力されなければ、ステップＳＡ２〜ＳＡ６の各判断結果が何れも「ＮＯ」となり、コマンド入力待ちの状態で待機する。
【００３３】
そして、コマンド入力待ちの状態において、再生コマンドが入力された場合には、ステップＳＡ３の判断結果が「ＹＥＳ」になり、ステップＳＡ７に進み、後述のステップＳＡ８（音響入力処理）にて取り込んだ音響波形データや、後述のステップＳＡ９（演奏入力処理）にて取り込んだ演奏データ（ＭＩＤＩデータ）を、音響出力手段１５から再生出力する再生処理を実行した後、上記ステップＳＡ２に処理を戻す。
【００３４】
音響入力コマンドが入力された場合には、ステップＳＡ４の判断結果が「ＹＥＳ」になり、ステップＳＡ８に進み、音響入力手段１４を用い、ステレオ音場として収音したステレオ成分の音響信号をサンプリングし、これにより得られる音響波形データ（Ｌチャンネル、Ｒチャンネル）をＲＡＭ１２の波形データエリアに保存する音響入力処理を実行した後、上記ステップＳＡ２に処理を戻す。
【００３５】
演奏入力コマンドが入力された場合には、ステップＳＡ５の判断結果が「ＹＥＳ」になり、ステップＳＡ９に進み、入力手段１３を介して外部ＭＩＤＩ楽器から入力される演奏データ（ＭＩＤＩデータ）を取り込んでＲＡＭ１２の演奏データエリアに保存する演奏入力処理を実行した後、上記ステップＳＡ２に処理を戻す。
【００３６】
音源位置表示コマンドが入力された場合には、ステップＳＡ６の判断結果が「ＹＥＳ」になり、ステップＳＡ１０に進み、音響入力処理（ステップＳＡ６）においてＲＡＭ１２に保存された音響波形データ（Ｌチャンネル、Ｒチャンネル）を読み出す。
【００３７】
次いで、ステップＳＡ１１では、読み出した音響波形データ（Ｌチャンネル、Ｒチャンネル）に周波数分析（高速フーリエ変換ＦＦＴ）を施してパワースペクトル（ＬチャンネルおよびＲチャンネル）を取得し、取得したＬチャンネルおよびＲチャンネルのパワースペクトルから定位情報ｐａｎを取得する分析処理を実行する。
【００３８】
続いて、ステップＳＡ１２では、取得した定位情報ｐａｎに応じてパワースペクトルの表示色を変化させ、それを時間軸（横軸）および対数周波数軸（縦軸）とした表示画面に表示する表示処理を実行した後、上記ステップＳＡ２に処理を戻す。そして、終了コマンドが入力されると、ステップＳＡ２の判断結果が「ＹＥＳ」になり、本ルーチンを終了する。
【００３９】
（２）分析処理の動作
次に、図７を参照して分析処理の動作を説明する。上述したメインルーチンのステップＳＡ１１（図６参照）を介して本処理が実行されると、ＣＰＵ１０は、図７に図示するステップＳＢ１に進み、ポインタｉをゼロリセットする。なお、ポインタｉは、「０」の場合にＬチャンネルの音響波形データを指定し、「１」の場合にＲチャンネルの音響波形データを指定する。
【００４０】
続いて、ステップＳＢ２では、ポインタｉが「２」より小さいか否か、すなわち左右両チャンネルの音響波形データについて分析し終えたかどうかを判断する。ポインタｉが「２」より小さく、左右両チャンネルの音響波形データについて分析し終えていなければ、判断結果は「ＹＥＳ」になり、ステップＳＢ３に進み、ポインタｊをゼロリセットする。
【００４１】
なお、ポインタｊは、分析フレームの切り出し数（オーバーラップ数）を表す。本実施形態では、一つのチャンネルの音響波形データについて、２つのフレーム（ｊ＝０，１）で重複切り出しする態様としている。したがって、このステップＳＢ４では、一つのチャンネルの音響波形データについて分析し終えたかどうかを判断する。
【００４２】
そして、ポインタｉが「２」より小さく、、一つのチャンネルの音響波形データについて分析し終えていなければ、判断結果が「ＹＥＳ」になり、ステップＳＢ５に進む。ステップＳＢ５では、ポインタｉで指定されるチャンネルの音響波形データに対し、ポインタｊ番目のフレームの切り出しを行い、続くステップＳＢ６では、それに窓関数（例えばハニング窓関数）を乗算する。
【００４３】
次いで、ステップＳＢ７では、ポインタｊ番目のフレームにより切り出されて窓関数が乗算された音響波形データに対してＦＦＴ（高速フーリエ変換）を施し、続くステップＳＢ８では、ＦＦＴで得られた実数部と虚数部の二乗和平方根を算出してパワースペクトル値を取得する。そして、ステップＳＢ９では、取得したパワースペクトル値をＲＡＭ１２のレジスタｄｏＤａｔａにストアする。続いて、ステップＳＢ１０に進み、ポインタｊをインクリメントして歩進させた後、上述のステップＳＢ４に処理を戻す。
【００４４】
そして、歩進されたポインタｊの値が「２」より小さい場合、すなわち一つのチャンネルの音響波形データについて分析し終えていなければ、上記ステップＳＢ４の判断結果が「ＹＥＳ」となり、上述したステップＳＢ５〜ＳＢ１０を繰り返す。これにより、Ｌチャンネルの音響波形データについてパワースペクトル値を取得し終える。
【００４５】
さて、こうしてＬチャンネルの音響波形データのパワースペクトル値を取得し終え、これに応じて歩進されたポインタｊの値が「２」になったとする。そうすると、上記ステップＳＢ４の判断結果が「ＮＯ」になり、ステップＳＢ１１に進み、ポインタｉをインクリメントして歩進させた後、前述のステップＳＢ２に処理を戻す。歩進されたポインタｉの値が「２」より小さい場合、すなわち左右両チャンネルの音響波形データについて分析し終ていなければ、上記ステップＳＢ２の判断結果が「ＹＥＳ」になり、ステップＳＢ３に進み、ポインタｊをゼロリセットする。
【００４６】
以後、Ｒチャンネルの音響波形データについて、２つのフレーム（ｊ＝０，１）で重複切り出しするパワースペクトル分析を行うべく、上記ステップＳＢ４〜ＳＢ１０の処理を繰り返し実行してＲチャンネルの音響波形データから取得したパワースペクトル値をＲＡＭ１２のレジスタｄｏＤａｔａにストアする。そして、上記ステップＳＢ４〜ＳＢ１０の処理を繰り返して歩進されたポインタｊの値が「２」になると、上記ステップＳＢ２の判断結果が「ＮＯ」になり、本処理を終える。
【００４７】
このように、分析処理では、ＲＡＭ１２から音響波形データ（Ｌチャンネル、Ｒチャンネル）を読み出して周波数分析（高速フーリエ変換ＦＦＴ）を施し、ＦＦＴで得られた実数部と虚数部の二乗和平方根を算出してパワースペクトル値（ＬチャンネルおよびＲチャンネル）を取得し、レジスタｄｏＤａｔａにストアするようになっている。なお、レジスタｄｏＤａｔａの前半にＬチャンネルのパワースペクトル値が、後半にＲチャンネルのパワースペクトル値がそれぞれ格納される。
【００４８】
（３）表示処理の動作
次に、図８を参照して表示処理の動作を説明する。図８は、表示処理の動作を示すフローチャートである。前述したメインルーチンのステップＳＡ１２（図６参照）を介して本処理が実行されると、ＣＰＵ１０は図７に図示するステップＳＣ１に進み、カウンタｉをゼロリセットする。なお、カウンタｉは「０」から横軸の「表示領域幅ｔｉｍｅｎ」までのフレーム数を計数する。
【００４９】
続いて、ステップＳＣ２では、フレーム数が横軸の「表示領域幅ｔｉｍｅｎ」より小さいか否か、すなわち表示領域幅ｔｉｍｅｎまで表示し終えたかどうかを判断する。表示領域幅ｔｉｍｅｎまで表示し終えていなければ、判断結果は「ＹＥＳ」となり、ステップＳＣ３に進む。
【００５０】
ステップＳＣ３では、カウンタｊをゼロリセットする。なお、カウンタｊは、「０」から縦軸の表示領域の高さｉＨｅｉｇｈｔまで累算する。次いで、ステップＳＣ４では、カウンタｊの値が表示領域の高さｉＨｅｉｇｈｔより小さいか否か、つまりカウンタｉで指定されるフレーム（時間軸）において、縦軸の表示領域の高さｉＨｅｉｇｈｔまで表示し終えたかどうかを判断する。
【００５１】
縦軸の表示領域の高さｉＨｅｉｇｈｔまで表示し終えていなければ、上記ステップＳＣ４の判断結果は「ＹＥＳ」になり、ステップＳＣ５に進む。ステップＳＣ５では、次式（１）に基づき縦軸の目盛りに相当するカウンタｊの値を周波数に換算してレジスタｄｆｒにストアする。
ｄｒｆ ← ４４０×２＾（（ｊ／ｒｅｓ１００ｃｅｎｔ−６９＋２１）／１２）…（１）
なお、上記（１）式において、「２＾」は２のべき乗を表し、ｒｅｓ１００ｃｅｎｔは半音当たりの表示ピクセル数を表す。ｊ／ｒｅｓ１００ｃｅｎｔは、縦軸上の位置（カウンタｊの値）をＭＩＤＩノートナンバに換算する。したがって、上記（１）式は、その換算されたＭＩＤＩノートナンバ（カウンタｊの値）を周波数に換算する。
【００５２】
次に、ステップＳＣ６では、次式（２）に従ってフレームを指定するインデックスｉｎｄｅｘを算出する。
ｉｎｄｅｘ ← （ｉＦｒａｍｅ２×ｄｆｒ）／（ｄｏＦｓ／２）…（２）
なお、上記（２）式において、ｉＦｒａｍｅ２は１フレームの半分の長さ、ｄｆｒは上記ステップＳＣ５で算出した縦軸上の位置（カウンタｊの値）の周波数換算値、ｄｏＦｓはサンプリング周波数である。
【００５３】
続いて、ステップＳＣ７では、次式（３）に従い、インデックスｉｎｄｅｘで指定されるフレームにおけるＬチャンネルおよびＲチャンネルの各パワースペクトルの平均値、つまりモノラル化したパワーｄを算出する。
ｄ ← （ｄｏＤａｔｅ［ｉｎｄｅｘ］＋
ｄｏＤａｔｅ［ｉｎｄｅｘ＋ｉＦｒａｍｅ２］）／２…（３）
【００５４】
次いで、ステップＳＣ８では、次式（４）に基づき、ＬチャンネルのパワースペクトルとＲチャンネルのパワースペクトルとの差分値を、ＬチャンネルおよびＲチャンネルの各パワースペクトルの和（絶対値）で除算して定位情報ｐａｎを算出する。
ｐａｎ ← （ｄｏＤａｔａ［ｉｎｄｅｘ］−ｄｏＤａｔｅ［ｉｎｄｅｘ＋ｉＦｒａｍｅ２］）／ＡＢＳ（ｄｏＤａｔａ［ｉｎｄｅｘ］＋ｄｏＤａｔｅ［ｉｎｄｅｘ＋ｉＦｒａｍｅ２］）…（４）
【００５５】
そして、ステップＳＣ９に進むと、次式（５）に従い、モノラル化したパワーｄの対数値ｌｏｇ（ｄ）に応じた色値変換係数ｉｃｏｌｉｎｄｅｘを算出する。
ｉｃｏｌｉｎｄｅｘ ← −ｉＺｏｏｍ＋（ｉｃｏｌｎ＋ｉＺｏｏｍ）×ｌｏｇ（ｄ）／ｍａｘｌｏｇ…（５）
なお、上記（５）式において、ｉＺｏｏｍおよびｉｃｏｌｎは色値変換定数である。
【００５６】
こうして、上記ステップＳＣ８およびステップＳＣ９において、定位情報ｐａｎおよび色値変換係数ｉｃｏｌｉｎｄｅｘが得られると、ステップＳＣ１０を介して色値算出処理を実行する。色値算出処理では、後述するように、色相テーブルＣＴ（図３参照）を参照して定位情報ｐａｎを色値（ｒ，ｇ，ｂ）に変換してレジスタｉｃｏｌ（ｒ，ｇ，ｂ）にストアする。
【００５７】
次いで、ステップＳＣ１１では、上記ステップＳＣ１０で得られたレジスタｉｃｏｌ（ｒ，ｇ，ｂ）を、画面描画用のビットマップレジスタｂｉｔｍａｐ［ｉ，ｊ×ｗｉｄｔｈ］にセットして画面表示する。なお、ｗｉｄｔｈは縦軸の表示ピクセル数である。また、ステップＳＣ１１では、１オクターブ毎に目盛りとなる音程ガイドラインを描画する。
【００５８】
続いて、ステップＳＣ１２では、カウンタｊをインクリメントして歩進させた後、前述のステップＳＣ４に処理を戻す。以後、カウンタｊの値が表示領域の高さｉＨｅｉｇｈｔを超えるまで上記ステップＳＣ４〜ＳＣ１２を繰り返すことによって、カウンタｉで指定されるフレーム（時間軸）において、定位情報ｐａｎに応じて表示色が変化するパワースペクトルを縦軸（対数周波数軸）に表示する。
【００５９】
そして、縦軸の表示領域の高さｉＨｅｉｇｈｔまで表示し終えると、上記ステップＳＣ４の判断結果が「ＮＯ」になり、ステップＳＣ１３に進み、カウンタｉをインクリメントして歩進させた後、前述のステップＳＣ２に処理を戻す。以後、歩進されたカウンタｉの値が横軸（時間軸）の「表示領域幅ｔｉｍｅｎ」に達するまで上記ステップＳＣ３〜ＳＣ１３を繰り返すことによって、定位情報ｐａｎに応じて表示色が変化するパワースペクトルを、曲頭から曲末尾までの時系列順に画面表示する。そして、歩進されたカウンタｉの値が横軸（時間軸）の「表示領域幅ｔｉｍｅｎ」を超えると、上記ステップＳＣ２の判断結果が「ＮＯ」になり、本処理を終える。
【００６０】
（４）色値算出処理の動作
次に、図９を参照して色値算出処理の動作を説明する。図９は、色値処理の動作を示すフローチャートである。上述した表示処理のステップＳＣ１０（図８参照）を介して本処理が実行されると、ＣＰＵ１０は図９に図示するステップＳＤ１に進み、前述の表示処理で算出した色値変換係数ｉｃｏｌｉｎｄｅｘを、レジスタｉｎｄｅｘに置換する。以下、レジスタｉｎｄｅｘの値を色値変換係数ｉｎｄｅｘと称する。
【００６１】
続いて、ステップＳＤ２以降では、定位情報ｐａｎの値に基づき色相テーブルＣＴ（図３参照）の色値変換領域を判別し、判別した色値変換領域における色関数（Ｒ関数、Ｇ関数、Ｂ関数）に従って定位情報ｐａｎを色値（ｒ，ｇ，ｂ）に変換する。以下、色値変換領域が「ｐａｎ≧０」の場合と、「ｐａｎ＜０」の場合とに分けて動作説明を進める。
【００６２】
ａ．色値変換領域がｐａｎ≧０の場合
色値変換領域がｐａｎ≧０であると、ステップＳＤ２の判断結果が「ＹＥＳ」になり、ステップＳＤ３に進み、次式（６）に基づき定位情報ｐａｎをＧ関数に代入して色値ｇを算出する。
ｇ ← （（１−（ｐａｎ＋ｐｅａｋＰｏｓＳｈｉｆｔＲ））×ｂｒｉｇｈｔｎｅｓｓ＋ｂｒｉｇｈｔｎｅｓｓＯｆｆｓｅｔ−ｂｒｉｇｈｔｎｅｓｓ）×ｉｎｄｅｘ／ｉｃｏｌｎ…（６）
なお、上記（６）式において、色値変換係数ｉｎｄｅｘ／色値変換定数ｉｃｏｌｎは色値幅を決める乗算係数である。
【００６３】
次いで、ステップＳＤ４では、次式（７）に基づき定位情報ｐａｎをＲ関数に代入して色値ｒを算出する。
ｒ ← （（１−ｐａｎ）×ｂｒｉｇｈｔｎｅｓｓ＋ｂｒｉｇｈｔｎｅｓｓＯｆｆｓｅｔ）×ｉｎｄｅｘ／ｉｃｏｌｎ…（７）
【００６４】
続いて、ステップＳＤ５では、定位情報ｐａｎが「１−ｐｅａｋＰｏｓＳｈｉｆｔＬ」以下であるか否かを判断する。定位情報ｐａｎが「１−ｐｅａｋＰｏｓＳｈｉｆｔＬ」以下であると、判断結果は「ＹＥＳ」になり、ステップＳＤ６に進み、次式（８）に基づき定位情報ｐａｎをＢ関数に代入して色値ｂを算出してステップＳＤ８に進む。
ｂ ← （（ｐａｎ−ｐｅａｋＰｏｓＳｈｉｆｔＬ−０）×ｂｒｉｇｈｔｎｅｓｓ＋ｂｒｉｇｈｔｎｅｓｓＯｆｆｓｅｔ）×ｉｎｄｅｘ／ｉｃｏｌｎ…（８）
【００６５】
一方、定位情報ｐａｎが「１−ｐｅａｋＰｏｓＳｈｉｆｔＬ」より大きければ、上記ステップＳＤ５の判断結果が「ＮＯ」になり、ステップＳＤ７に進み、次式（９）に基づき定位情報ｐａｎをＢ関数に代入して色値ｂを算出してステップＳＤ８に進む。
ｂ ← （（１−（ｐａｎ−ｐｅａｋＰｏｓＳｈｉｆｔＬ））×ｂｒｉｇｈｔｎｅｓｓ＋ｂｒｉｇｈｔｎｅｓｓＯｆｆｓｅｔ）×ｉｎｄｅｘ／ｉｃｏｌｎ…（９）
【００６６】
そして、ステップＳＤ８では、上記ステップＳＤ２〜ＳＤ７において、定位情報ｐａｎに応じて色相テーブルＣＴ（図３参照）から得た色値（ｒ，ｇ，ｂ）をレジスタｉｃｏｌにストアして本処理を終える。
【００６７】
ｂ．色値変換領域がｐａｎ＜０の場合
色値変換領域がｐａｎ＜０であると、ステップＳＤ２の判断結果が「ＮＯ」になり、図１０に図示するステップＳＤ９に進み、次式（１０）に基づき定位情報ｐａｎをＢ関数に代入して色値ｂを算出する。
ｂ ← （（１＋ｐａｎ−ｐｅａｋＰｏｓＳｈｉｆｔＬ）×ｂｒｉｇｈｔｎｅｓｓ＋ｂｒｉｇｈｔｎｅｓｓＯｆｆｓｅｔ−ｂｒｉｇｈｔｎｅｓｓ）×ｉｎｄｅｘ／ｉｃｏｌｎ…（１０）
【００６８】
続いて、ステップＳＤ１１では、次式（１１）に基づき定位情報ｐａｎをＲ関数に代入して色値ｒを算出する。
ｒ ← （（１＋ｐａｎ）×ｂｒｉｇｈｔｎｅｓｓ＋ｂｒｉｇｈｔｎｅｓｓＯｆｆｓｅｔ）×ｉｎｄｅｘ／ｉｃｏｌｎ…（１１）
【００６９】
次いで、ステップＳＤ１１では、定位情報ｐａｎが「−１＋ｐｅａｋＰｏｓＳｈｉｆｔＲ」以上であるか否かを判断する。定位情報ｐａｎが「−１＋ｐｅａｋＰｏｓＳｈｉｆｔＲ」以上であると、判断結果は「ＹＥＳ」になり、ステップＳＤ１２に進み、次式（１２）に基づき定位情報ｐａｎをＧ関数に代入して色値ｇを算出してステップＳＤ８に進む。
ｇ ← （（０−（ｐａｎ＋ｐｅａｋＰｏｓＳｈｉｆｔＲ））×ｂｒｉｇｈｔｎｅｓｓ＋ｂｒｉｇｈｔｎｅｓｓＯｆｆｓｅｔ）×ｉｎｄｅｘ／ｉｃｏｌｎ…（１２）
【００７０】
一方、定位情報ｐａｎが「−１＋ｐｅａｋＰｏｓＳｈｉｆｔＲ」より小さければ、上記ステップＳＤ１１の判断結果が「ＮＯ」になり、ステップＳＤ１３に進み、次式（１３）に基づき定位情報ｐａｎをＧ関数に代入して色値ｇを算出してステップＳＤ８に進む。
ｇ ← （（１＋（ｐａｎ＋（１−ｐｅａｋＰｏｓＳｈｉｆｔＲ）））×ｂｒｉｇｈｔｎｅｓｓ＋ｂｒｉｇｈｔｎｅｓｓＯｆｆｓｅｔ）×ｉｎｄｅｘ／ｉｃｏｌｎ…（１３）
【００７１】
そして、ステップＳＤ８では、上記ステップＳＤ９〜ＳＤ１３において、定位情報ｐａｎに応じて色相テーブルＣＴ（図３参照）から得た色値（ｒ，ｇ，ｂ）をレジスタｉｃｏｌにストアして本処理を終える。
【００７２】
以上説明したように、第１実施形態によれば、ＬチャンネルおよびＲチャンネルの音響波形データに周波数分析（高速フーリエ変換ＦＦＴ）を施して得たＬチャンネルおよびＲチャンネルのパワースペクトルから定位情報ｐａｎを取得し、取得した定位情報ｐａｎに応じて表示色を異ならせたパワースペクトルを曲頭から曲末尾までの時系列順に画面表示するので、例えば図１１に図示する表示画面例のように、時系列に変化する各音源の位置を一目瞭然にすることが可能になる。
【００７３】
なお、図１１はカラー表示出来ない為にモノクロ表示となっているが、実際には前述した色相テーブルＣＴ（図３参照）に基づき定位情報ｐａｎに応じてパワースペクトルの表示色を異ならせている為、どの色（どの位置）で表示されるかで時系列に変化する各音源の位置が一目瞭然になる訳である。
【００７４】
なお、上述した第１実施形態では、ＲＡＭ１２の波形データエリアに保存しておいた音響波形データ（Ｌチャンネル、Ｒチャンネル）を処理対象としたが、これに限定されずステレオ音場として収音したステレオ成分の音響信号をサンプリングした音響波形データ（Ｌチャンネル、Ｒチャンネル）についてパワースペクトル分析および表示処理を施してリアルタイムに変化する音源位置を表示することも勿論可能である。
【００７５】
また、本実施形態では、色相テーブルＣＴ（図３参照）に基づき定位情報ｐａｎに応じてパワースペクトルの表示色を異ならせているが、これに限らず、更に定位情報ｐａｎを明度に変換するテーブルや、定位情報ｐａｎを彩度に変換するテーブルを設け、これらを組み合わせて定位情報ｐａｎに応じてパワースペクトルの表示色、明度および彩度を異ならせることでより一層音源の位置を明確に表現する態様も可能である。
【００７６】
［第２実施形態］
次に、第２実施形態を説明する。第２実施形態による音源位置表示装置１００は、図１に図示した第１実施形態と同一の構成なので、その説明については省略する。第２実施形態が第１実施形態と相違する点は、色補正テーブルＲＴを備えることにある。ここで、図１２〜図１４を参照して色補正テーブルＲＴの構成を説明する。
【００７７】
先ず、図１２は、第２実施形態による音源位置表示装置１００のＲＡＭ１２が備えるワークエリアＷＥに一時記憶される主要なレジスタデータ（変数）の構成を示す図である。この図において、第１実施形態と相違する点は、色補正テーブルＲＴと、当該色補正テーブルＲＴの色補正特性を指定する値を保持するレジスタｅｍｐｈａｓｉｚｅＬＲ、ｅＬＲＣｒａｔｉｏ、ｔｒａｐｅｚｏｉｄを具備することにある。
【００７８】
色補正テーブルＲＴとは、定位情報ｐａｎを色値（ｒ，ｇ，ｂ）の補正係数ｅｍｐｈａｓｉｚｅに変換する。補正係数ｅｍｐｈａｓｉｚｅは、前述の色相テーブルＣＴによって定位情報ｐａｎに応じて変換された色値（ｒ，ｇ，ｂ）に各々乗算される。
【００７９】
ここで、図１３に図示する色補正テーブルＲＴの一例を参照して、当該色補正テーブルＲＴの色補正特性を指定する値を保持するレジスタｅｍｐｈａｓｉｚｅＬＲ、ｅＬＲＣｒａｔｉｏ、ｔｒａｐｅｚｏｉｄについて説明する。色補正テーブルＲＴは、定位情報ｐａｎを補正係数ｅｍｐｈａｓｉｚｅに変換する台形の補正関数を備える。この台形の補正関数は、レジスタｔｒａｐｅｚｏｉｄにストアされる上底の半値幅、レジスタｅｍｐｈａｓｉｚｅＬＲにストアされる上底中心位置およびレジスタｅＬＲＣｒａｔｉｏにストアされる脚（斜辺）の傾きで指定される。
【００８０】
すなわち、台形の補正関数において、上底に含まれる範囲の定位情報ｐａｎに対応して「１」の補正係数ｅｍｐｈａｓｉｚｅを発生する。つまり上底に含まれる範囲の定位情報ｐａｎについては色値補正が為されず、上底から下底に向かって右脚の傾きを有する直線（又は上底から下底に向かって左脚の傾きを有する直線）に対応する範囲の定位情報ｐａｎでは色値をリニアに減衰させる補正係数ｅｍｐｈａｓｉｚｅを発生する。なお、左脚（左斜辺）の傾きは、右脚の傾きと逆符号となる。
【００８１】
図１３に図示した色補正テーブルＲＴは、右脚（右斜辺）の傾きを表すレジスタｅＬＲＣｒａｔｉｏの値が「０」以下の一例である。右脚の傾きが「０」より大きい場合には、図１４に図示するように、図１３に図示した台形を上下反転した形の補正関数となる。具体的には、下底範囲に含まれる範囲の定位情報ｐａｎについては「０」の補正係数ｅｍｐｈａｓｉｚｅを発生する。つまり下底に含まれる範囲の定位情報ｐａｎについては色値を「０」とする。下底から上底に向かって右脚の傾きを有する直線（又は下底から上底に向かって左脚の傾きを有する直線）に対応する範囲の定位情報ｐａｎでは色値をリニアに増加させる補正係数ｅｍｐｈａｓｉｚｅを発生する。
【００８２】
このような補正特性を備える色補正テーブルＲＴを用いることによって、特定の範囲に含まれる定位情報ｐａｎに対応した色値を強調する為、その特定の範囲に位置するパワースペクトルの表示色が強調され、これにより音源位置を浮かび上がらせるように表示することが可能になる。
【００８３】
以下、こうした色補正テーブルＣＴを用いる第２実施形態の表示処理の動作について説明する。図１５は、第２実施形態による表示処理の動作を示すフローチャートである。この図に示す第２実施形態の動作が、図８に図示した第１実施形態の表示処理と相違する点は、ステップＳＥ１０の色値算出処理の後に、ステップＳＥ１４の色補正処理を設けたことにある。第２実施形態の表示処理におけるステップＳＥ１〜ＳＥ１３は、前述した第１実施形態の表示処理におけるステップＳＣ１〜ＳＣ１３と同一であるので、その説明については省略する。
【００８４】
図１６は、第２実施形態による色補正処理の動作を示すフローチャートである。図１５に図示した第２実施形態による表示処理のステップＳＥ１４を介して色補正処理が実行されると、ＣＰＵ１０は、図１６に図示するステップＳＦ１に進み、定位情報ｐａｎからレジスタｅｍｐｈａｓｉｚｅＬＲＣにストアされる上底中心位置を減算した値が「０．０」以上であるか否かを判断する。以下、ｐａｎ−ｅｍｐｈａｓｉｚｅＬＲＣ≧０．０の場合と、ｐａｎ−ｅｍｐｈａｓｉｚｅＬＲＣ＜０の場合とに分けて動作説明を進める。
【００８５】
ａ．ｐａｎ−ｅｍｐｈａｓｉｚｅＬＲＣ≧０．０の場合
この場合、上記ステップＳＦ１の判断結果が「ＹＥＳ」になり、ステップＳＦ２に進み、レジスタｅＬＲＣｒａｔｉｏにストアされる脚（斜辺）の傾きが「０．０」以上であるか否かを判断する。脚（斜辺）の傾きが「０．０」以上であると、上記ステップＳＦ２の判断結果は「ＹＥＳ」になり、ステップＳＦ３に進む。ステップＳＦ３では、次式（１４）に基づき定位情報ｐａｎを色補正テーブルＲＴ（台形の補正関数）に代入して補正係数ｅｍｐｈａｓｉｚｅを算出した後、ステップＳＦ８に進む。
ｅｍｐｈａｓｉｚｅ ← １．０−ｅＬＲＣｒａｔｉｏ×
（ｐａｎ−ｔｒａｐｅｚｏｉｄ−ｅｍｐｈａｓｉｚｅＬＲＣ）…（１４）
【００８６】
一方、脚（斜辺）の傾きが「０．０」より小さいと、上記ステップＳＦ２の判断結果が「ＮＯ」になり、ステップＳＦ４に進む。ステップＳＦ４では、次式（１５）に基づき定位情報ｐａｎを色補正テーブルＲＴ（台形の補正関数）に代入して補正係数ｅｍｐｈａｓｉｚｅを算出した後、ステップＳＦ８に進む。
ｅｍｐｈａｓｉｚｅ ← ０．０−ｅＬＲＣｒａｔｉｏ×
（ｐａｎ−ｔｒａｐｅｚｏｉｄ−ｅｍｐｈａｓｉｚｅＬＲＣ）…（１５）
【００８７】
ｂ．ｐａｎ−ｅｍｐｈａｓｉｚｅＬＲＣ＜０の場合
この場合、上記ステップＳＦ１の判断結果が「ＮＯ」になり、ステップＳＦ５に進み、レジスタｅＬＲＣｒａｔｉｏにストアされる脚（斜辺）の傾きが「０．０」以上であるか否かを判断する。脚（斜辺）の傾きが「０．０」以上であると、上記ステップＳＦ５の判断結果は「ＹＥＳ」になり、ステップＳＦ６に進む。ステップＳＦ６では、次式（１６）に基づき定位情報ｐａｎを色補正テーブルＲＴ（台形の補正関数）に代入して補正係数ｅｍｐｈａｓｉｚｅを算出した後、ステップＳＦ８に進む。
ｅｍｐｈａｓｉｚｅ ← １．０−ｅＬＲＣｒａｔｉｏ×
（ｅｍｐｈａｓｉｚｅＬＲＣ−ｐａｎ−ｔｒａｐｅｚｏｉｄ）…（１６）
【００８８】
一方、脚（斜辺）の傾きが「０．０」より小さいと、上記ステップＳＦ５の判断結果が「ＮＯ」になり、ステップＳＦ７に進む。ステップＳＦ７では、次式（１７）に基づき定位情報ｐａｎを色補正テーブルＲＴ（台形の補正関数）に代入して補正係数ｅｍｐｈａｓｉｚｅを算出した後、ステップＳＦ８に進む。
ｅｍｐｈａｓｉｚｅ ← ０．０−ｅＬＲＣｒａｔｉｏ×
（ｅｍｐｈａｓｉｚｅＬＲＣ−ｐａｎ−ｔｒａｐｅｚｏｉｄ−）…（１７）
【００８９】
そして、ステップＳＦ８では、算出された補正係数ｅｍｐｈａｓｉｚｅが「０．０」〜「１．０」の範囲内であれば、その算出された補正係数ｅｍｐｈａｓｉｚｅをそのまま使うが、「０．０」より小さい場合には補正係数ｅｍｐｈａｓｉｚｅを「０．０」に下限値設定し、一方、「１．０」より大きい場合には補正係数ｅｍｐｈａｓｉｚｅを「１．０」に上限値設定する。次いで、ステップＳＦ９では、前述した色値算出処理で得られた色値（ｒ，ｇ，ｂ）にそれぞれ補正係数ｅｍｐｈａｓｉｚｅを乗算して補正色値（ｒ，ｇ，ｂ）を生成し、続くステップＳＦ１０では、生成された補正色値（ｒ，ｇ，ｂ）をレジスタｉｃｏｌにストアして本処理を終える。
【００９０】
以上説明したように、第２実施形態によれば、ＬチャンネルおよびＲチャンネルの音響波形データに周波数分析（高速フーリエ変換ＦＦＴ）を施して得たＬチャンネルおよびＲチャンネルのパワースペクトルから定位情報ｐａｎを取得し、取得した定位情報ｐａｎに応じてパワースペクトルの表示色を異ならせ、更に特定の範囲に含まれる定位情報ｐａｎに対応した色値のパワースペクトルを強調しつつ曲頭から曲末尾までの時系列順に画面表示するので、例えば図１７に図示する表示画面例のように、時系列に変化する各音源の位置を一目瞭然にすることが可能になる。
【００９１】
なお、図１７はカラー表示出来ない為にモノクロ表示となっているが、実際には前述した色相テーブルＣＴ（図３参照）に基づき定位情報ｐａｎに応じてパワースペクトルの表示色を異ならせ、更に色補正テーブルＲＴ（図１３参照）に基づき特定の範囲に含まれる定位情報ｐａｎに対応した色値のパワースペクトルを強調している為、どの色（どの位置）で表示されるかで時系列に変化する各音源の位置が一目瞭然になる訳である。
【００９２】
また、上述した第２実施形態では、色値算出処理で得られた色値（ｒ，ｇ，ｂ）に対して一律に補正係数ｅｍｐｈａｓｉｚｅを乗算して補正色値（ｒ，ｇ，ｂ）を生成する態様としたが、これに限らず、例えば色値ｒ、色値ｇおよび色値ｂ毎の補正テーブルを設けておき、定位情報ｐａｎに応じて各補正テーブルからそれぞれ得られる各色値毎の補正係数ｅｍｐｈａｓｉｚｅを、対応する色値（ｒ，ｇ，ｂ）に各々乗算して補正色値（ｒ，ｇ，ｂ）を得ることも可能である。このようにした場合、特定の範囲に含まれる定位情報ｐａｎに対応した特定色のパワースペクトルを強調表示したり抑制表示したりすることで所望の音源位置を明確に表現し得ることも可能になる。
【００９３】
以上、本発明の実施の一形態について説明したが、本発明はそれに限定されるものではなく、本願出願の特許請求の範囲に記載された発明とその均等の範囲に含まれる。以下では、本願出願当初の特許請求の範囲に記載された各発明について付記する。
【００９４】
（付記）
［請求項１］
音響波形データからパワースペクトルを取得するパワー取得手段と、
前記パワー取得手段により取得されたパワースペクトルから定位情報を取得する定位情報取得手段と、
前記定位情報取得手段により取得された定位情報に応じて表示色を異ならせたパワースペクトルを表示する表示手段と
を具備することを特徴とする音源位置表示装置。
【００９５】
［請求項２］
前記表示手段は、前記定位情報取得手段により取得された定位情報を色相に変換する色相変換手段を備え、当該色相変換手段により変換された色相で当該定位情報に対応するパワースペクトルを表示することを特徴とする請求項１記載の音源位置表示装置。
【００９６】
［請求項３］
前記色相変換手段は、前記定位情報取得手段により取得された定位情報に応じた三原色値（ｒ，ｇ，ｂ）を発生する色関数を備えることを特徴とする請求項１記載の音源位置表示装置。
【００９７】
［請求項４］
前記表示手段は、前記定位情報取得手段により取得された定位情報を明度に変換する明度変換手段を更に備えることを特徴とする請求項１記載の音源位置表示装置。
【００９８】
［請求項５］
前記表示手段は、前記定位情報取得手段により取得された定位情報を彩度に変換する彩度変換手段を更に備えることを特徴とする請求項１記載の音源位置表示装置。
【００９９】
［請求項６］
前記表示手段は、前記定位情報取得手段により取得された定位情報に応じて表示色を異ならせたパワースペクトルを時系列順に画面表示することを特徴とする請求項１記載の音源位置表示装置。
【０１００】
［請求項７］
前記表示手段は、前記定位情報取得手段により取得された定位情報が、所定範囲に含まれる場合に、前記色相変換手段により変換された色相を補正する色補正手段を更に備えることを特徴とする請求項１記載の音源位置表示装置。
【０１０１】
［請求項８］
前記色補正手段は、前記定位情報取得手段により取得された定位情報が、特定の範囲に含まれる場合に、前記色相変換手段により変換された色相を強調することを特徴とする請求項７記載の音源位置表示装置。
【０１０２】
［請求項９］
音響波形データからパワースペクトルを取得し、取得したパワースペクトルから定位情報を取得し、取得した定位情報に応じて表示色を異ならせたパワースペクトルを表示することを特徴とする音源位置表示方法。
【０１０３】
［請求項１０］
コンピュータに、
音響波形データからパワースペクトルを取得するパワー取得ステップと、
前記パワー取得ステップで取得されたパワースペクトルから定位情報を取得する定位情報取得ステップと、
前記定位情報取得ステップで取得された定位情報に応じて表示色を異ならせたパワースペクトルを表示する表示ステップと
を実行させることを特徴とするプログラム。
【符号の説明】
【０１０４】
１０  ＣＰＵ
１１  ＲＯＭ
１２  ＲＡＭ
１３  入力手段
１４  表示手段
１５  音響出力手段
１６  音響入力手段
１００  音源位置表示装置
【図１】
2014229965.tif 000003
【図２】
2014229965.tif 000004
【図３】
2014229965.tif 000005
【図４】
2014229965.tif 000006
【図５】
2014229965.tif 000007
【図６】
2014229965.tif 000008
【図７】
2014229965.tif 000009
【図８】
2014229965.tif 000010
【図９】
2014229965.tif 000011
【図１０】
2014229965.tif 000012
【図１１】
2014229965.tif 000013
【図１２】
2014229965.tif 000014
【図１３】
2014229965.tif 000015
【図１４】
2014229965.tif 000016
【図１５】
2014229965.tif 000017
【図１６】
2014229965.tif 000018
【図１７】
2014229965.tif 000019