(57)【要約】【課題】各音響信号の周波数変動に柔軟に対応しながら精度よく、かつ高速に検索することを可能とする音響信号検索装置を提供する。【解決手段】音響信号検索装置１において、基準周波数検出部１１は、参照信号及び蓄積信号それぞれの時間周波数スペクトログラムから、時間間隔ごとの所定の基準特徴部を抽出し、抽出した所定の基準特徴部の周波数を時間間隔ごとの参照信号及び蓄積信号それぞれの基準周波数とする。特徴抽出部１２は、当該基準周波数を基準として一意に定められる周波数のパワースペクトル値に基づいて参照信号の信号全体における参照特徴量を抽出する。類似度計算部１３は、抽出した参照特徴量に基づいて蓄積信号の時間間隔ごとの特徴量を抽出し、抽出した蓄積信号の区間ごとの区間特徴量を算出し、算出した区間特徴量と、参照特徴量とに基づいて類似度を算出する。【選択図】図５
2008015002.tif 000002
【特許請求の範囲】
【請求項１】
入力される音響信号の時間周波数スペクトログラムに基づいて、当該音響信号の特徴量を抽出する音響信号特徴抽出装置における音響信号特徴抽出方法であって、
前記時間周波数スペクトログラムの時間間隔ごとに、前記時間周波数スペクトログラムの所定の基準特徴部を抽出し、
抽出した前記所定の基準特徴部の周波数を前記時間間隔ごとの基準周波数とし、
当該基準周波数を基準として一意に定められる周波数のパワースペクトル値に基づいて前記時間間隔ごとの音響信号の特徴量を抽出する
ことを特徴とする音響信号特徴抽出方法。
【請求項２】
前記所定の基準特徴部とは、パワースペクトル値が各時間間隔において最大値となる部分であることを特徴とする請求項１に記載の音響信号特徴抽出方法。
【請求項３】
前記所定の基準特徴部とは、時間軸で連続する幾つかの時間周波数スペクトログラムの中で、パワースペクトル値が最大値となる部分
であることを特徴とする請求項１に記載の音響信号特徴抽出方法。
【請求項４】
前記所定の基準特徴部とは、時間軸で連続する幾つかの時間周波数スペクトログラムのパワースペクトル値の重心となる部分
であることを特徴とする請求項１に記載の音響信号特徴抽出方法。
【請求項５】
参照音響信号に基づいて、当該参照音響信号に類似する音を含む区間を、入力される検索対象音響信号から検索する音響信号検索装置における音響信号検索方法であって、
前記参照音響信号及び前記検索対象音響信号それぞれの時間周波数スペクトログラムから、前記時間周波数スペクトログラムの時間間隔ごとの所定の基準特徴部を抽出し、
抽出した前記所定の基準特徴部の周波数を前記時間間隔ごとの前記参照音響信号及び前記検索対象音響信号それぞれの基準周波数とし、
前記参照音響信号の基準周波数を基準として一意に定められる周波数のパワースペクトル値に基づいて前記参照音響信号の前記時間間隔ごとの特徴量を抽出し、抽出した特徴量に基づいて前記参照音響信号の信号全体の時間長を一区間とする参照特徴量を算出し、
前記検索対象音響信号の基準周波数を基準として一意に定められる周波数のパワースペクトル値に基づいて前記検索対象音響信号の前記時間間隔ごとの特徴量を抽出し、
抽出した前記検索対象音響信号の前記特徴量に基づいて前記区間ごとの区間特徴量を算出し、
算出した区間特徴量と、前記参照特徴量とに基づいて類似度を算出し、
算出した類似度に基づいて、前記参照音響信号の音に類似する前記検索対象音響信号の区間を検索し、検索により検出した前記参照音響信号の音に類似する前記検索対象音響信号の区間を出力する
ことを特徴とする音響信号検索方法。
【請求項６】
入力される音響信号の時間周波数スペクトログラムに基づいて、当該音響信号の特徴量を抽出する音響信号特徴抽出装置であって、
前記時間周波数スペクトログラムの時間間隔ごとに、前記時間周波数スペクトログラムの所定の基準特徴部を抽出し、抽出した前記所定の基準特徴部の周波数を前記時間間隔ごとの基準周波数とする基準周波数検出手段と、
当該基準周波数を基準として一意に定められる周波数のパワースペクトル値に基づいて前記時間間隔ごとの音響信号の特徴量を抽出する特徴量抽出手段と、
を備えたことを特徴とする音響信号特徴抽出装置。
【請求項７】
参照音響信号に基づいて、当該参照音響信号に類似する音を含む区間を、入力される検索対象音響信号から検索する音響信号検索装置であって、
前記参照音響信号及び前記検索対象音響信号それぞれの時間周波数スペクトログラムから、前記時間周波数スペクトログラムの時間間隔ごとの所定の基準特徴部を抽出し、抽出した前記所定の基準特徴部の周波数を前記時間間隔ごとの前記参照音響信号及び前記検索対象音響信号それぞれの基準周波数とする基準周波数検出手段と、
前記参照音響信号の基準周波数を基準として一意に定められる周波数のパワースペクトル値に基づいて前記参照音響信号の前記時間間隔ごとの特徴量を抽出し、抽出した特徴量に基づいて前記参照音響信号の信号全体の時間長を一区間とする参照特徴量を算出し、前記検索対象音響信号の基準周波数を基準として一意に周波数のパワースペクトル値に基づいて前記検索対象音響信号の前記時間間隔ごとの特徴量を抽出する特徴抽出手段と、
前記特徴抽出手段が抽出した前記検索対象音響信号の前記特徴量に基づいて前記区間ごとの区間特徴量を算出し、算出した区間特徴量と、前記参照特徴量とに基づいて類似度を算出し、算出した類似度に基づいて、前記参照音響信号の音に類似する前記検索対象音響信号の区間を検索し、検索により検出した前記参照音響信号の音に類似する前記検索対象音響信号の区間を出力する類似度計算手段と、
を備えたことを特徴とする音響信号検索装置。
【請求項８】
入力される音響信号の時間周波数スペクトログラムに基づいて、当該音響信号の特徴量を抽出する音響信号特徴抽出装置のコンピュータに、
前記時間周波数スペクトログラムの時間間隔ごとに、前記時間周波数スペクトログラムの所定の基準特徴部を抽出するステップと、
抽出した前記所定の基準特徴部の周波数を前記時間間隔ごとの基準周波数とするステップと、
当該基準周波数を基準として一意に定められる周波数のパワースペクトル値に基づいて前記時間間隔ごとの音響信号の特徴量を抽出するステップと、
を実行させるための音響信号特徴抽出プログラム。
【請求項９】
参照音響信号に基づいて、当該参照音響信号に類似する音を含む区間を、入力される検索対象音響信号から検索する音響信号検索装置のコンピュータに、
参照音響信号及び検索対象音響信号の時間周波数スペクトログラムを入力するステップと、
前記参照音響信号及び前記検索対象音響信号それぞれの時間周波数スペクトログラムから、前記時間周波数スペクトログラムの時間間隔ごとの所定の基準特徴部を抽出するステップと、
抽出した前記所定の基準特徴部の周波数を前記時間間隔ごとの前記参照音響信号及び前記検索対象音響信号それぞれの基準周波数とするステップと、
前記参照音響信号の基準周波数を基準として一意に定められる周波数のパワースペクトル値に基づいて前記参照音響信号の前記時間間隔ごとの特徴量を抽出し、抽出した特徴量に基づいて前記参照音響信号の信号全体の時間長を一区間とする参照特徴量を算出するステップと、
前記検索対象音響信号の基準周波数を基準として一意に定められる周波数のパワースペクトル値に基づいて前記検索対象音響信号の前記時間間隔ごとの特徴量を抽出するステップと、
抽出した前記検索対象音響信号の前記特徴量に基づいて前記区間ごとの区間特徴量を算出するステップと、
前記類似度計算手段が、算出した区間特徴量と、前記参照特徴量とに基づいて類似度を算出するステップと、
算出した類似度に基づいて、前記参照音響信号の音に類似する前記検索対象音響信号の区間を検索し、検索により検出した前記参照音響信号の音に類似する前記検索対象音響信号の区間を出力するステップと、
を実行させるための音響信号検索プログラム。
【請求項１０】
請求項８に記載の音響信号特徴抽出プログラムを記録したコンピュータ読み取り可能な記録媒体。
【請求項１１】
請求項９に記載の音響信号検索プログラムを記録したコンピュータ読み取り可能な記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、音響信号である参照信号を検索キーとし、この参照信号よりも長い音響信号である蓄積信号の中から参照信号に類似した信号の位置を検索する音響信号検索に関するものであり、特に各音響信号の周波数変動（移調またはピッチシフト）に対応が可能な音響信号特徴抽出方法、抽出装置、抽出プログラム、該プログラムを記録した記録媒体、および該特徴を利用した音響信号検索方法、検索装置、検索プログラム、並びに該プログラムを記録した記録媒体に関する。
【背景技術】
【０００２】
近年、特徴抽出法と特徴照合法を要素とした音響信号検索技術が、放送のモニタリングやプレイリストの配信、或いは音楽認識サービスなどにおいて実用化されてきている。これは、音楽（またはその一部）を入力すれば自動的にその特徴を抽出し、データベースに格納された楽曲の特徴と照合することにより所望の情報を瞬時に検出できるというものである。
【０００３】
例えば、音響信号の高速な一致探索手法として、「時系列アクティブ探索（ＴＡＳ）法」（特許文献１参照）が知られている。これは時系列信号から特徴データとしてヒストグラムを生成し、生成したヒストグラムの重なり率から類似度を算出するものであり、また、特徴照合の際、時間軸方向で照合が不要な区間を算出して探索をスキップできることから、無駄な照合を省略し、実用上十分な探索精度を保ったまま探索速度を大幅に向上させることが可能になるというものである。
【０００４】
また、テレビドラマのＢＧＭやラジオのＤＪつき音楽番組のように、ナレーションなどがかぶさった音響信号から背景音楽を探索する技術として特許文献２に示される「分割一致探索（ＤＡＬ）法」がある。ＤＡＬ法の原理は、検索キーとなる音の信号を時間・周波数方向において小領域に分割し、各小領域について一致する小領域を検索対象の信号の中から探索するというものであり、この分割によりナレーションの重畳がない領域を見つけ出し、流れている音楽が検索キーと同じである区間を検出することができる。
【０００５】
しかしながら、これらの手法では、例えば、もとは同じ音響信号であっても、音程を変えて再生された音響信号については類似性を検出することができなかった。音程が異なるが同じ音楽を探索する手法としてハミング検索の手法がある。特にGhiasらは単旋律の音高系列を“Ｓ（同じ）”、“Ｕ（上がる）”、“Ｄ（下がる）”の相対音高の文字列で表すmelodic contour を用いた単旋律の照合手法、及びハミングを検索キーとし、類似する単旋律を検索する類似音楽の検索法を提案している（非特許文献１参照）。
【０００６】
melodic contourは旋律の表現としては楽譜などに比べ単純であるが、単旋律の識別、照合には有効であり、その他の研究においても単旋律の探索における特徴として用いられている。そして、これまでの類似音楽検索の多くは、ハミングなどで入力された単旋律の音高系列を検索キーとし、同じく単旋律の音高系列からなるデータベースから検索に類似するものを、単旋律同士の照合により探索するものであった。ここで、多くの手法において、データベースの音高系列は、ＭＩＤＩや楽譜情報、もしくは検索キーと同じくハミングによる主旋律の入力から得られる、既に主旋律などが抽出された単旋律の音高系列であった。そのため、通常、多重奏である通常の音楽ＣＤや放送中の音を用いた音響信号の探索に、これらのハミング検索の手法を用いるには、多重音からの主旋律抽出、もしくは多重奏と多重奏の間での主旋律の照合が必要になるなど困難な問題が多い。
【特許文献１】特許第３０６５３１４号明細書
【特許文献２】特開２００４−１０２０２３号公報
【非特許文献１】A. Ghias, J. Logan, D. Chamberlin and B. C. Smith “Query by humming:Musical information retrieval in an audio database,” Proc. ACM Multimedia '95, pp.231-236, San Francisco, USA, Nov. 1995.
【発明の開示】
【発明が解決しようとする課題】
【０００７】
上述したように、これまでの音響信号の一致探索では、同じ音響信号であっても音程を変えて再生された音響信号の探索は困難であるという問題があった。
【０００８】
また、入力された音響信号に対して周波数変動を付加するという方法では、照合に必要となる情報量が大幅に増大してしまうという問題があった。
【０００９】
本発明は、上記問題を解決すべくなされたもので、検索キーとして参照信号を入力した場合に、一致する信号が含まれる蓄積信号中の区間を、各音響信号の周波数変動に柔軟に対応しながら精度よく、かつ高速に検索することが可能な音響信号特徴抽出方法、抽出装置、抽出プログラム、該プログラムを記録した記録媒体、および該特徴を利用した音響信号検索方法、検索装置、検索プログラム、並びに該プログラムを記録した記録媒体を提供することにある。
【課題を解決するための手段】
【００１０】
上記問題を解決するために、本発明は、入力される音響信号の時間周波数スペクトログラムに基づいて、当該音響信号の特徴量を抽出する音響信号特徴抽出装置における音響信号特徴抽出方法であって、前記時間周波数スペクトログラムの時間間隔ごとに、前記時間周波数スペクトログラムの所定の基準特徴部を抽出し、抽出した前記所定の基準特徴部の周波数を前記時間間隔ごとの基準周波数とし、当該基準周波数を基準として一意に定められる周波数のパワースペクトル値に基づいて前記時間間隔ごとの音響信号の特徴量を抽出することを特徴とする音響信号特徴抽出方法である。
【００１１】
本発明は、上記に記載の発明において、前記所定の基準特徴部とは、パワースペクトル値が各時間間隔において最大値となる部分であることを特徴とする。
【００１２】
本発明は、上記に記載の発明において、前記所定の基準特徴部とは、時間軸で連続する幾つかの時間周波数スペクトログラムの中で、パワースペクトル値が最大値となる部分であることを特徴とする。
【００１３】
本発明は、上記に記載の発明において、前記所定の基準特徴部とは、時間軸で連続する幾つかの時間周波数スペクトログラムのパワースペクトル値の重心となる部分であることを特徴とする。
【００１４】
本発明は、参照音響信号に基づいて、当該参照音響信号に類似する音を含む区間を、入力される検索対象音響信号から検索する音響信号検索装置における音響信号検索方法であって、前記参照音響信号及び前記検索対象音響信号それぞれの時間周波数スペクトログラムから、前記時間周波数スペクトログラムの時間間隔ごとの所定の基準特徴部を抽出し、抽出した前記所定の基準特徴部の周波数を前記時間間隔ごとの前記参照音響信号及び前記検索対象音響信号それぞれの基準周波数とし、前記参照音響信号の基準周波数を基準として一意に定められる周波数のパワースペクトル値に基づいて前記参照音響信号の前記時間間隔ごとの特徴量を抽出し、抽出した特徴量に基づいて前記参照音響信号の信号全体の時間長を一区間とする参照特徴量を算出し、前記検索対象音響信号の基準周波数を基準として一意に定められる周波数のパワースペクトル値に基づいて前記検索対象音響信号の前記時間間隔ごとの特徴量を抽出し、抽出した前記検索対象音響信号の前記特徴量に基づいて前記区間ごとの区間特徴量を算出し、算出した区間特徴量と、前記参照特徴量とに基づいて類似度を算出し、算出した類似度に基づいて、前記参照音響信号の音に類似する前記検索対象音響信号の区間を検索し、検索により検出した前記参照音響信号の音に類似する前記検索対象音響信号の区間を出力することを特徴とする音響信号検索方法である。
【００１５】
本発明は、入力される音響信号の時間周波数スペクトログラムに基づいて、当該音響信号の特徴量を抽出する音響信号特徴抽出装置であって、前記時間周波数スペクトログラムの時間間隔ごとに、前記時間周波数スペクトログラムの所定の基準特徴部を抽出し、抽出した前記所定の基準特徴部の周波数を前記時間間隔ごとの基準周波数とする基準周波数検出手段と、当該基準周波数を基準として一意に定められる周波数のパワースペクトル値に基づいて前記時間間隔ごとの音響信号の特徴量を抽出する特徴量抽出手段と、を備えたことを特徴とする音響信号特徴抽出装置である。
【００１６】
本発明は、参照音響信号に基づいて、当該参照音響信号に類似する音を含む区間を、入力される検索対象音響信号から検索する音響信号検索装置であって、前記参照音響信号及び前記検索対象音響信号それぞれの時間周波数スペクトログラムから、前記時間周波数スペクトログラムの時間間隔ごとの所定の基準特徴部を抽出し、抽出した前記所定の基準特徴部の周波数を前記時間間隔ごとの前記参照音響信号及び前記検索対象音響信号それぞれの基準周波数とする基準周波数検出手段と、前記参照音響信号の基準周波数を基準として一意に定められる周波数のパワースペクトル値に基づいて前記参照音響信号の前記時間間隔ごとの特徴量を抽出し、抽出した特徴量に基づいて前記参照音響信号の信号全体の時間長を一区間とする参照特徴量を算出し、前記検索対象音響信号の基準周波数を基準として一意に周波数のパワースペクトル値に基づいて前記検索対象音響信号の前記時間間隔ごとの特徴量を抽出する特徴抽出手段と、前記特徴抽出手段が抽出した前記検索対象音響信号の前記特徴量に基づいて前記区間ごとの区間特徴量を算出し、算出した区間特徴量と、前記参照特徴量とに基づいて類似度を算出し、算出した類似度に基づいて、前記参照音響信号の音に類似する前記検索対象音響信号の区間を検索し、検索により検出した前記参照音響信号の音に類似する前記検索対象音響信号の区間を出力する類似度計算手段と、を備えたことを特徴とする音響信号検索装置である。
【００１７】
本発明は、入力される音響信号の時間周波数スペクトログラムに基づいて、当該音響信号の特徴量を抽出する音響信号特徴抽出装置のコンピュータに、前記時間周波数スペクトログラムの時間間隔ごとに、前記時間周波数スペクトログラムの所定の基準特徴部を抽出するステップと、抽出した前記所定の基準特徴部の周波数を前記時間間隔ごとの基準周波数とするステップと、当該基準周波数を基準として一意に定められる周波数のパワースペクトル値に基づいて前記時間間隔ごとの音響信号の特徴量を抽出するステップと、を実行させるための音響信号特徴抽出プログラムである。
【００１８】
本発明は、参照音響信号に基づいて、当該参照音響信号に類似する音を含む区間を、入力される検索対象音響信号から検索する音響信号検索装置のコンピュータに、参照音響信号及び検索対象音響信号の時間周波数スペクトログラムを入力するステップと、前記参照音響信号及び前記検索対象音響信号それぞれの時間周波数スペクトログラムから、前記時間周波数スペクトログラムの時間間隔ごとの所定の基準特徴部を抽出するステップと、抽出した前記所定の基準特徴部の周波数を前記時間間隔ごとの前記参照音響信号及び前記検索対象音響信号それぞれの基準周波数とするステップと、前記参照音響信号の基準周波数を基準として一意に定められる周波数のパワースペクトル値に基づいて前記参照音響信号の前記時間間隔ごとの特徴量を抽出し、抽出した特徴量に基づいて前記参照音響信号の信号全体の時間長を一区間とする参照特徴量を算出するステップと、前記検索対象音響信号の基準周波数を基準として一意に定められる周波数のパワースペクトル値に基づいて前記検索対象音響信号の前記時間間隔ごとの特徴量を抽出するステップと、抽出した前記検索対象音響信号の前記特徴量に基づいて前記区間ごとの区間特徴量を算出するステップと、前記類似度計算手段が、算出した区間特徴量と、前記参照特徴量とに基づいて類似度を算出するステップと、算出した類似度に基づいて、前記参照音響信号の音に類似する前記検索対象音響信号の区間を検索し、検索により検出した前記参照音響信号の音に類似する前記検索対象音響信号の区間を出力するステップと、を実行させるための音響信号検索プログラムである。
【００１９】
本発明は、上記に記載の音響信号特徴抽出プログラムを記録したコンピュータ読み取り可能な記録媒体である。
【００２０】
本発明は、上記に記載の音響信号検索プログラムを記録したコンピュータ読み取り可能な記録媒体である。
【発明の効果】
【００２１】
この発明によれば、時間周波数スペクトログラムの時間間隔ごとに、時間周波数スペクトログラムの所定の基準特徴部を抽出し、抽出した前記所定の基準特徴部の周波数を時間間隔ごとの基準周波数とし、当該基準周波数を基準として一意に定められる周波数のパワースペクトル値に基づいて前記時間間隔ごとの音響信号の特徴量を抽出する構成とした。これにより、音響信号の各時間周波数スペクトログラムから、パワースペクトルの所定の基準特徴部を検出して各時間間隔における基準周波数として設定し、当該基準周波数に基づいてパワースペクトルの特徴を抽出することにより、周波数による変動が発生していたとしても、その変動に影響を受けない特徴抽出を行うことができる。
【００２２】
また、本発明によれば、所定の基準特徴部とは、パワースペクトル値が各時間間隔において最大値となる部分である構成とした。これにより、パワースペクトル値が最大値となる周波数値を基準周波数とした特徴抽出を行うことが可能となる。
【００２３】
また、本発明によれば、所定の基準特徴部とは、時間軸で連続する幾つかの時間周波数スペクトログラムの中で、パワースペクトル値が最大値となる部分である構成とした。これにより、連続する時間の中で特徴抽出に最も適切な基準周波数を設定することが可能となる。
【００２４】
また、本発明によれば、所定の基準特徴部とは、時間軸で連続する幾つかの時間周波数スペクトログラムのパワースペクトル値の重心となる部分である構成とした。これにより、音響信号にノイズ等が含まれていたとしても特徴抽出に最も適切な基準周波数を設定することが可能となる。
【００２５】
また、本発明によれば、参照音響信号及び前記検索対象音響信号それぞれの時間周波数スペクトログラムから、時間周波数スペクトログラムの時間間隔ごとの所定の基準特徴部を抽出し、抽出した所定の基準特徴部の周波数を時間間隔ごとの参照音響信号及び検索対象音響信号それぞれの基準周波数とし、参照音響信号の基準周波数を基準として一意に定められる周波数のパワースペクトル値に基づいて参照音響信号の前記時間間隔ごとの特徴量を抽出し、抽出した特徴量に基づいて前記参照音響信号の信号全体の時間長を一区間とする参照特徴量を算出し、検索対象音響信号の基準周波数を基準として一意に定められる周波数のパワースペクトル値に基づいて検索対象音響信号の前記時間間隔ごとの特徴量を抽出し、抽出した前記検索対象音響信号の前記特徴量に基づいて区間ごとの区間特徴量を算出し、算出した区間特徴量と、参照特徴量とに基づいて類似度を算出し、算出した類似度に基づいて、参照音響信号の音に類似する検索対象音響信号の区間を検索し、検索により検出した参照音響信号の音に類似する前記検索対象音響信号の区間を出力する構成とした。これにより、参照音響信号および検索対象音響信号の各時間周波数スペクトログラムから、パワースペクトルの所定の基準特徴部を検出して各時間間隔における基準周波数として設定、すなわち、周波数による変動が発生していたとしても、その周波数変動に依存しない特徴を抽出可能な基準周波数を設定することが可能となる。そして、該基準周波数に基づき各信号の周波数スペクトル特徴を抽出し、抽出した周波数スペクトル特徴に基づいて周波数の上下の変動に対応した音響信号検出が可能となる。
【発明を実施するための最良の形態】
【００２６】
（本発明の原理）
最初に、図１から図４を参照して本発明の原理について説明する。本発明は、音響特徴の抽出の際、音響特徴を抽出する各時点で、音響信号の周波数の変動に応じて、音響特徴を抽出する周波数をずらす構成を有することを特徴としている。例えば、音響信号のパワースペクトログラムが、周波数軸を対数で表した場合に、図１のように表されるとする。このとき、音響特徴は、図２のように、ｆ１，ｆ２，…，ｆｃ，…，ｆｋの各周波数のパワーにとして抽出される。このとき、この音響信号の周波数をａ倍にして再生すると、図３のように、そのパワースペクトログラムは、パワースペクトログラム値の変化の特徴を維持したまま、周波数軸方向にａ倍だけずれたものになる。すなわち、このａ倍を移調やピッチシフトとすれば、ａの値さえ分かれば、移調やピッチシフトによる周波数変動の影響を吸収することが可能となる。
【００２７】
そこで、本発明では、音響特徴を抽出する基準となる周波数をｆｃを、例えば、その時点のスペクトルのピークにあわせ、当該ｆｃを基準として音響特徴を示す周波数としてｆ１，ｆ２，…，ｆｃ，…ｆｋを抽出し、抽出した周波数ごとのパワースペクトル値の平均値を算出し、算出した平均値で周波数ごとのパワースペクトル値を正規化する。そして、図４に示すように、ａ倍に相当する移調やピッチシフトが行われた音響信号についても、同じように基準周波数をその時点のスペクトルのピークにあわせることで、ａｆｃが基準周波数として得られ、当該ａｆｃを基準として音響特徴を示す周波数としてａｆ１，ａｆ２，…，ａｆｃ，…ａｆｋを抽出し、上記と同じように抽出した周波数ごとのパワースペクトル値を正規化する。このようにして得られた２つの音響特徴は、周波数値の違いはあるものの、ｆｃとａｆｃを基準として観測した場合、対応する周波数における正規化されたパワースペクトル値は同じ値となり、周波数変動による影響が吸収された音響特徴の抽出を行うことが可能となる。以下、当該原理を実現する本発明の実施形態を図面を参照して説明する。
【００２８】
（本発明の第１実施形態）
次に、本発明の第１実施形態に係る音響信号検索装置１について説明する。音響信号検索装置１は、検索したい音響信号、すなわち参照信号を入力し、これより長い音響信号である蓄積信号中から参照信号と類似する信号を含む区間を検出するものである。また、音響信号検索装置１は、上述した原理を適用して音響信号の周波数変動があっても、当該変動に対して頑健、つまり、その変動を吸収して検索を可能にするものである。なお、本明細書において蓄積信号中で参照信号に類似する信号を含む区間を検出するとは、この区間が始まる区間の先頭の時点を検出することと定義する。
【００２９】
図５は、第１実施形態に係る音響信号検索装置１の構成を示すブロック図である。音響信号検索装置１は、基準周波数検出部１１と、特徴抽出部１２と、類似度計算部１３とを備えている。音響信号検索装置１では、参照信号並びに蓄積信号が入力されると、まず最初に図示しない変換部により変換されて、それぞれの信号の時間周波数スペクトログラム（以下、スペクトログラムと記載）が出力される。
【００３０】
基準周波数検出部１１は、上記の変換部から参照信号と蓄積信号のそれぞれの信号のスペクトログラムが入力され、入力される参照信号及び蓄積信号のスペクトログラムからパワースペクトルの特徴的な部分を検出し、各時点における基準周波数ｆｃを設定する。特徴抽出部１２は、基準周波数検出部によって検出された各時点における基準周波数ｆｃを基準に、参照信号ならびに蓄積信号の特徴を算出する。類似度計算部１３は、上記参照信号と同じ長さの注目窓を蓄積信号に設定し、参照信号の特徴と注目窓内における蓄積信号の特徴との類似度を算出し、類似度を所定の閾値と比較することで類似性を判定する。
【００３１】
次に、上述した音響信号検索装置１の動作説明を行う。図６は、図５の音響信号検索装置１の動作例を示すフローチャートである。音響信号検索装置１は、検索キーである参照信号を読み込むことで入力を行い（ステップＳ２０１）、前述した変換部により、スペクトログラムが出力される(ステップＳ２０２)。ここで、変換部によるスペクトログラムへの変換は、例えば、高速フーリエ変換を行うことにより、あるいは帯域通過フィルタを複数用いることにより行われる。
【００３２】
変換部により出力された参照信号のスペクトログラムは、基準周波数検出部１１に入力され、基準周波数検出部１１は、図７（ａ）に示すようにスペクトログラムを構成する周波数スペクトルの各時点においてパワースペクトルが最も高くなる周波数を特徴的な部分として検出し、検出した周波数をこの時点における基準周波数ｆｃとして出力する（ステップＳ２０３）。
【００３３】
なお、基準周波数検出部１１は、図７（ａ）、（ｂ）に示すようにサンプリングされた全ての周波数スペクトルについて基準周波数を各々決定する代わりに、図７（ｃ）に示すように時間軸に連続する複数の周波数スペクトルのうち強度が最大となる周波数を、これらの周波数スペクトルの基準周波数として設定するようにしてもよい。他にも、時間軸方向において連続する複数の周波数スペクトルの強度の重心を求め、これに対応する周波数を、この時間内における周波数スペクトルの基準周波数として設定するようにしてもよい。
【００３４】
また、基準周波数検出部１１における基準周波数の決定方法は上述の方法に限られず、例えば、各サンプリング時点においてパワースペクトルのピークを複数検出し、強度がＸ番目の周波数を選択したり、あるいは、強度がＹ番目になる周波数から数Ｈｚずらした点を選択したりすることによって基準周波数を決定することができる。また、これら複数検出されるパワースペクトルの、ピーク間に現れる谷となる点を基準周波数としてもよい。
【００３５】
特徴抽出部１２は、参照信号について、基準周波数検出部１１によって検出された基準周波数ｆｃを基準に、周波数軸上所定の間隔でパワースペクトル値を複数選択する。例えば、対数軸上で等間隔にパワースペクトル値を複数選択するようにしてもよい。そして、特徴抽出部１２は、選択したパワースペクトル値を要素とする特徴ベクトルをこの時点における周波数スペクトル特徴（参照信号の特徴量）として抽出する。次に、特徴抽出部１２は、各時点における周波数スペクトル特徴となる特徴ベクトルをベクトル量子化により符号化する（ステップＳ２０４）。そして、特徴抽出部１２は、更にこれらの符号の所定時間内における出現回数を符号毎に計数したヒストグラム（以下、当該ヒストグラムをヒストグラム特徴と記載）を生成して、生成したヒストグラム特徴を音響信号の固有の情報、すなわち音響指紋として類似度計算部１３に送出する（ステップＳ２０５）。なお、ヒストグラムを生成する際の前記の所定時間は、参照信号の時間長に等しいものとする。
【００３６】
一方、音響信号検索装置１は、検索対象となる蓄積信号の音響波形信号を読み込んで入力し（ステップＳ２０６）、上記参照信号と同様の手法で時間周波数スペクトログラムを抽出する（ステップＳ２０７）。そして、基準周波数検出部１１と特徴抽出部１２は、参照信号と同様に基準周波数ｆｃを検出し（ステップＳ２０８）、各時点における周波数スペクトル特徴となる特徴ベクトルを求め、ベクトル量子化により符号化する（ステップＳ２０９）。特徴抽出部１２は、蓄積信号の場合には、ヒストグラム特徴を出力するのではなく、特徴ベクトルを符号化した符号系列を類似度計算部１３に出力する。
【００３７】
なお、一般的には、蓄積信号の符号系列は、蓄積信号データベースに記憶され、検索処理が開始されると、逐次蓄積信号データベースにアクセスして蓄積信号の符号系列を読み出す構成をとることになるが、この構成には限られない。例えば、参照信号と蓄積信号を同時に入力し、参照信号のヒストグラム特徴を類似度計算部１３に記憶しておき、リアルタイムに入力される蓄積信号の符号系列に対して照合を行うようにすることもできる。
【００３８】
類似度計算部１３は、参照信号のヒストグラム特徴が特徴抽出部１２から入力されると、入力される蓄積信号の符号系列について、その先頭より参照信号と同じ時間長の注目窓を設定する（ステップＳ２１０)。そして、注目窓内における符号の出現回数を符号毎に計数することによりヒストグラムを生成し（ステップＳ２１１）、生成したヒストグラムと参照信号のヒストグラム特徴との類似度を算出する（ステップＳ２１２）。算出する手法としては、例えば、時系列アクティブ探索（ＴＡＳ）法（特許文献１）を適用することが可能である。
【００３９】
ここで、図８を参照しつつＴＡＳ法について説明する。ＴＡＳ法は、参照信号ヒストグラム特徴とのヒストグラム重なり率による類似度ＳＩＲが、探索閾値θを超えるヒストグラム特徴を有する蓄積信号上の区間を全て探索するものである。具体的には、類似度計算部１３は、前述したように、始めに、蓄積信号の特徴である符号系列に先頭より参照信号と同じ時間長の注目窓を設定する。次に、参照信号のヒストグラム特徴ＨＲと蓄積信号符号系列の注目窓におけるヒストグラム特徴ＨＩとの類似度ＳＩＲを算出する。類似度ＳＩＲの算出は、次式（１）によって行うことができる。
【００４０】
【数１】
2008015002.tif 000003
【００４１】
ただし、式（１）において、Ｌは、ヒストグラムのビン数であり、Ｄはヒストグラムの総度数であり、ｈｌγ、ｈＲγは、それぞれＨＩ、ＨＲのγ番目のビンに含まれる度数を表している。なお、本実施形態では、図８の入力にあたるスペクトログラムとして、参照信号の周波数スペクトル特徴を時系列に並べたものと、蓄積信号の周波数スペクトル特徴を時系列に並べたものを適用し、それぞれの周波数スペクトル特徴に対してベクトル量子化を行うことになる。
【００４２】
類似度計算部１３は、算出された類似度の値が閾値θを越えた場合に、蓄積信号中の該区間に類似性があると判断して検出結果を出力する。類似度計算部１３は、設定された注目窓における類似度計算が終了すると、蓄積信号の符号系列に対し新たな注目窓を設定する。そのとき、すでに算出した類似度ＳＩＲに基づき、類似する可能性のないことが保証される蓄積信号上の区間との照合を省略し、次の注目窓をどこに設定できるかを表す値であるスキップ可能幅（移動可能量）ｋを算出する(ステップＳ２１３)。スキップ可能幅ｋは、次式（２）に基づいて算出することができる。
【００４３】
【数２】
2008015002.tif 000004
【００４４】
すなわち、注目窓の設定は、ヒストグラム重なり率による類似度から、注目窓の現在位置近傍における類似度の上限値を求め、当該上限値に基づいて蓄積信号における次の注目窓の位置を求めて設定することになる。
【００４５】
類似度計算部１３は、設定した注目窓の位置が蓄積信号符号系列の終端に達したか否かを判定し（ステップＳ２１４）、蓄積信号符号系列の終端に達するまでステップＳ２１０からステップＳ２１３の処理を繰り返し、蓄積信号上で適宜注目窓を移動させながら類似度計算を行うことになる。そして、注目窓が蓄積信号の終端に達し、該蓄積信号について照合処理が終わると、音響信号検索装置１は検出結果を出力して全ての検索処理を終了する（ステップＳ２１５）。
【００４６】
なお、注目窓の設定を行う順序は、ステップＳ２１０に示した蓄積信号の特徴量である符号系列を対象として行う構成に限定されることはなく、例えば、図９に示すフローチャートのように蓄積信号の信号波形に対して行ってもよい（図９：ステップＳ３１１）。図９の場合には、注目窓内の情報についてスペクトログラム抽出、基準周波数検出、特徴量抽出、ヒストグラム生成を行い（図９：ステップＳ３１２からステップＳ３１５）、参照信号のヒストグラム特徴との類似度計算を行う（ステップＳ２１２）こととなる。
【００４７】
また、図１０に示すフローチャートのように、まず、図６のフローチャートと同様に、蓄積信号から特徴抽出を行い（ステップＳ２０６からＳ２０９）、さらに各時点におけるヒストグラム特徴を予め算出しておくことにより（図１０：ステップＳ３２１）、注目窓を蓄積信号のヒストグラム系列に対して設定（図１０：ステップＳ３２２）するようにすることもできる。
【００４８】
また、上述した第１実施形態における音響信号検索装置１の別の形態として、図１１に示す音響信号検索装置１ａがある。音響信号検索装置１ａは、図５に示した音響信号検索装置１の構成に加えて、特徴抽出部１２より抽出されたヒストグラム特徴を符号化する符号化部１４を備えた構成となっている。音響信号検索装置１ａでは、ヒストグラム特徴をベクトル量子化などにより符号化し、各音響信号の音響指紋を符号系列として表現することで、類似度計算部１３ａにおいて、参照信号と蓄積信号中の区間との類似度を、当該符号系列間の類似性から算出することができ、照合に必要となる計算量を大幅に削減するものである。
【００４９】
なお、蓄積信号検索装置１ａの符号化部１４から出力される各音響信号の音響指紋である符号系列について、ヒストグラム特徴を符号化する代わりに、特徴抽出部１２において算出される周波数スペクトル特徴（あるいは複数のスペクトル特徴を連続した一連の特徴ベクトルとしたもの）を、直接ベクトル量子化して算出するようにしてもよい。
【００５０】
また、さらに、上述した第１実施形態における音響信号検索装置１の別の形態として図１２に示す音響信号検索装置１ｂがある。音響信号検索装置１ｂでは、図１１に示した音響信号検索装置１ａの構成に加えて、索引生成部１５を備えた構成となっている。索引生成部１５は、符号化部１４により出力された蓄積信号の各符号について、予め蓄積信号上の所在を示す索引情報を生成する。このような構成とすることで、参照信号の符号に類似する蓄積信号上の符号の位置は、索引情報により容易に検出でき、これにより類似度計算部１３ｂは、類似性のある符号を蓄積信号上で探索する処理を省略し、蓄積信号上の区間における類似度をより高速に算出できるようになる。なお、図１２において、実線は、参照信号のスペクトログラムを対象とした処理の流れを示したものであり、破線は、蓄積信号のスペクトログラムを対象とした処理の流れを示したものである。
【００５１】
図１３は、図１２に示した音響信号検索装置１ｂの動作例を示したフローチャートである。なお、当該フローチャートにおいて、図６で説明した動作と同じ動作については、同じステップ番号を付し、以下、異なる動作について説明する。
【００５２】
まず、音響信号検索装置１ｂにおいて符号化部１４は、特徴抽出部１２により算出された各音響信号のヒストグラム特徴を符号化する（ステップＳ６０１、ステップＳ６０２）。索引生成部１５は蓄積信号の各符号について索引情報を生成する（ステップＳ８０１）。類似度計算部１３ｂは、参照信号の符号に類似する蓄積信号の符号を索引検索により検出し（ステップＳ８０３：Ｙｅｓ）、検出した蓄積信号の符号から、類似する可能性のある蓄積信号中の照合区間を出力する（ステップＳ８０４）。次に、類似度計算部１３ｂは、出力した照合区間がなくなるまで、照合区間の選択を繰り返し、照合区間ごとに参照信号との類似度、すなわち区間類似度を算出する（ステップＳ８０５からＳ８０７）。そして、類似度計算部１３ｂは、区間類似度を所定の閾値である探索閾値と比較し、区間類似度が探索閾値を超えているものを、類似しているものとし、区間類似度が探索閾値を超えていないものを類似していないものとして検出結果を出力する（ステップＳ８０８、ステップＳ８０９）。なお、ステップＳ８０３において、類似度計算部１３ｂは、参照信号の符号に類似する蓄積信号の符号を索引検索により検出できない場合（ステップＳ８０３：Ｎｏ）、その結果を検出結果として出力する（ステップＳ８０９）。
【００５３】
（第１実施形態の実装例）
次に、上述した第１実施形態の実装例について説明する。
まず、音響信号としては、音楽ＣＤなどで用いられる音響信号（サンプリング周波数４４．１ｋＨｚ、１６ｂｉｔ量子化、ステレオ）をダウンサンプリング、モノラル化、量子化ビット数を変更した音響信号（サンプリング周波数１１．０２５ｋＨｚ、８ビット量子化、モノラル）を用いる。そして、この音響信号に対し、まず、１０２４点のＦＦＴ（Fast Fourier Transform）を１０ｍｓごとに行い、１０ｍｓごとにその時点のパワースペクトルを抽出する。そして、この１０ｍｓごとにスペクトルを時間方向に並べたものとして、スペクトログラムを得る。また、（（１／ｘ３）×ｆｃ，（１／ｘ２）×ｆｃ，（１／ｘ）×ｆｃ，ｘ×ｆｃ，ｘ２×ｆｃ，ｘ３×ｆｃ）の各周波数においてパワーを抽出し、抽出した各周波数のパワーの平均値を算出し、算出した平均値で各周波数のパワーを正規化したものを特徴ベクトルとする。ここで、ｘは７１／６である。
【００５４】
そして、ｆｃは先に抽出したスペクトルのうち１３００Ｈｚから上下１０％の周波数のうちでパワーが最大となる周波数とする。このｆｃに基づき、各時点で、（（１／ｘ３）×ｆｃ，（１／ｘ２）×ｆｃ，（１／ｘ）×ｆｃ，ｘ×ｆｃ，ｘ２×ｆｃ，ｘ３×ｆｃ）の各周波数におけるパワーを、ＦＦＴの結果から得られる各周波数のうち、これらの周波数に最も近い周波数のパワーを（（１／ｘ３）×ｆｃ，（１／ｘ２）×ｆｃ，（１／ｘ）×ｆｃ，ｘ×ｆｃ，ｘ２×ｆｃ，ｘ３×ｆｃ）でのパワーとし、各周波数のパワーの平均値を算出し、算出した平均値で各周波数のパワーを正規化したものを特徴ベクトルとして抽出する。
【００５５】
また、特徴ベクトルの各成分の対応する周波数について、その周波数を中心周波数とするバンドパスフィルタを用いて、再度パワースペクトルを抽出してもよい。
【００５６】
このようにして、特徴ベクトルを抽出することで、上下１０％までピッチが変動した音響信号でも探索を行うことができる。
【００５７】
なお、これらの特徴ベクトルの次元（上記実装例では７）やｘの値（上記実装例では７１／６）、ｆｃの設定範囲（上記実装例では１０％）は、探索の目的とする信号やピッチ変動の幅（上記実装例では上下１０％）に応じて、変更可能なものであり、上記した具体的な数値に限られないものとする。
【００５８】
また、特徴ベクトルの各成分が対応する周波数について、上記実装例では、第１成分から等比数列としたが、各成分間の比が決まっていれば、それらが等比である必要はない。例えば、第２成分の対応する周波数は第１成分が対応する周波数の１．２倍の周波数、第３成分の対応する周波数は第１成分の対応する周波数の２倍の周波数などとしてもよい。
【００５９】
（第２実施形態）
次に、本発明の第２実施形態に係る音響信号検索装置２について説明する。音響信号検索装置２は、検索したい音響信号、すなわち参照信号を入力し、これより長い信号である蓄積信号中から参照信号と類似する信号を含む区間を検出するものである。また、音響信号検索装置２は、音響信号に大きな妨害音（加法性雑音）が重畳し、周波数変動が発生した場合であっても、当該変動に対して頑健に、つまり、その変動を吸収して検索を可能にするものである。
【００６０】
図１４は、第２実施形態に係る音響信号検索装置２の構成を示すブロック図である。
図９における音響信号検索装置２は、スペクトログラム分割部２０と、基準周波数検出部２１と、特徴抽出部２２と、符号化部２４と、索引生成部２５と、類似度計算部２３とを備えている。なお、図１４において、実線は、参照信号のスペクトログラムを対象とした処理の流れを示したものであり、破線は、蓄積信号のスペクトログラムを対象とした処理の流れを示したものである。
【００６１】
スペクトログラム分割部２０は、入力される参照信号ならびに蓄積信号のスペクトログラムを、小領域スペクトログラムと呼ぶ小領域のスペクトログラムに分割し、それぞれ小領域参照信号スペクトログラムおよび小領域蓄積信号スペクトログラムとして出力する。基準周波数検出部２１は、参照信号ならびに蓄積信号の各小領域信号スペクトログラムからパワースペクトルの特徴的な部分を検出し、各時点における基準周波数ｆｃを設定する。特徴抽出部２２は、基準周波数検出部２１によって検出された各時点における基準周波数ｆｃを基準に、参照信号ならびに蓄積信号の小領域スペクトログラム特徴を算出する。
【００６２】
符号化部２４は、参照信号ならびに蓄積信号の各小領域スペクトログラム特徴を符号化し、それぞれ参照信号小領域符号および蓄積信号小領域符号として出力する。索引生成部は、蓄積信号小領域符号について、蓄積信号中の所在を示す索引情報を生成する。類似度計算部２３は、参照信号小領域符号に類似する符号を、後述する予め生成した小領域符号間類似度表を参照することにより検出し、検出した類似する符号について索引検索を行う。また、類似度計算部２３は、蓄積信号中に、類似する符号、すなわち小領域参照信号スペクトログラムに類似する小領域蓄積信号スペクトログラムが存在することを検出すると、該小領域蓄積信号スペクトログラムを含む蓄積信号中の区間と参照信号との類似度すなわち区間類似度を算出する。
【００６３】
次に、上述した音響信号検索装置２の動作について説明する。図１５は、音響信号検索装置２の動作例を示したフローチャートである。音響信号検索装置２は、参照信号ならびに蓄積信号の音響信号を読み込むことにより入力すると（ステップＳ９０１、ステップＳ９１１）、第１実施形態に係る音響信号検索装置１と同じく、変換部によりこれらの波形信号からパワースペクトルを得て、得られたパワースペクトルに基づいて時間周波数スペクトログラムが抽出される（ステップＳ９１１、ステップＳ９１２)。
【００６４】
次に、スペクトログラム分割部２０は、参照信号ならびに蓄積信号の各時間周波数スペクトログラムから一定の大きさ（狭帯域・短時間）のスペクトログラム、すなわち小領域スペクトログラム(small-region spectrogram)を周波数軸方向及び時間軸方向において所定の間隔で切り出し、それぞれ小領域参照信号スペクトログラム、小領域蓄積信号スペクトログラムとして出力する（ステップＳ９０３、ステップＳ９１３)。小領域スペクトログラムは、互いに重複するようにして出力するようにしてもよいし、重複させなくともよい。ここで、先頭の時点がｔｎであり、周波数帯がωｍである参照信号中の小領域参照信号スペクトログラムをＦｔｎ，ωｍと表す。同様に蓄積信号のスペクトログラムにおいて、先頭の時点がｔであり、周波数帯がωｍの上記Ｆｔｎ，ωｍと同じサイズの小領域蓄積信号スペクトログラムをＧｔ，ωｍと表す。
【００６５】
また、以下の説明において、参照信号のスペクトログラムにおいて、小領域スペクトログラムＦｔｎ，ωｍを切り出した全ての時点ｔｎの集合をＴＲ（ＴＲ＝｛ｔ１，ｔ２，…｝）とし、全ての周波数帯ωｍの集合をＷ（Ｗ＝｛ω１，ω２，…｝）として記載する。
【００６６】
基準周波数検出部２１は、スペクトログラム分割部２０により出力された参照信号ならびに蓄積信号の各小領域スペクトログラムについて、周波数帯ωｍにおけるパワースペクトルのピークを検出する（ステップＳ９０４、ステップＳ９１４）。
【００６７】
なお、小領域スペクトログラムの基準周波数ｆｃωｍの設定方法は、前述した第１実施形態における基準周波数検出部１１と同様に、小領域スペクトログラム中のサンプリングされた全ての周波数スペクトルについて基準周波数を各々決定する方法、または、ある一定時間においてパワースペクトル値が最大となる周波数もしくはパワースペクトル値の重心に対応する周波数を、この時間内の各時点における周波数スペクトルの基準周波数として設定する方法が可能である。また、第１実施形態と同様に、各サンプリング時点においてパワースペクトルのピークを複数検出し、強度がＸ番目の周波数を選択したり、あるいは強度がＹ番目になる周波数から数Ｈｚずらした点を選択したりすることによって基準周波数を決定する方法や、あるいは、パワースペクトルのピークに注目する代わりに、複数検出されるパワースペクトルの、ピーク間に現れる谷となる点を基準周波数とするような方法を用いるようにしてもよい。
【００６８】
次に、特徴抽出部２２は、参照信号ならびに蓄積信号の各小領域スペクトログラムについて、上記基準周波数検出部２１により検出された基準周波数ｆｃωｍを基準に、小領域スペクトログラム中の各時点における周波数スペクトル特徴を算出する（ステップＳ９０５、Ｓ９１５）。第２実施形態では、音量の変動に対応するため、周波数スペクトル特徴は、抽出した各時点において、その平均値で正規化する。そして、この正規化されたパワースペクトル値を要素とする周波数スペクトル特徴の特徴ベクトルをベクトル量子化により符号化し、これらの符号の小領域スペクトログラム内における出現回数を符号毎に計数したヒストグラム（ヒストグラム特徴）を生成する（ステップＳ９０６、ステップＳ９１６）。
【００６９】
符号化部２４は、特徴抽出部２２により算出された参照信号ならびに蓄積信号の各小領域スペクトログラムのヒストグラム特徴を、周波数帯ωｍ毎にベクトル量子化により符号化する（ステップＳ９０７、ステップＳ９１７）。具体的には、ヒストグラムのビンの値、すなわち各スペクトル符号の小領域スペクトログラム内での出現回数を要素とする特徴ベクトルをベクトル量子化して、あらかじめ生成されたベクトル量子化符号帳により、符号を割り当てる方法を用いる。なお、符号化部２４は、小領域スペクトログラムを符号化する際には、同じ周波数帯ωｍに属する小領域参照信号スペクトログラムならびに小領域蓄積信号スペクトログラムに対し同じベクトル符号帳を用いる。
【００７０】
符号化部２４によるベクトル量子化は、ヒストグラム特徴とベクトル量子化符号帳の各代表ベクトル間の距離を算出し、距離値が最小となる代表特徴に対応する符号を出力することによって行われる。以下、ヒストグラム特徴が符号化されて得られた符号を、参照信号小領域符号ならびに蓄積信号小領域符号と記載し、小領域参照信号スペクトログラムＦｔｎ，ωｍの参照信号小領域符号をｃ（Ｆｔｎ，ωｍ）、小領域蓄積信号スペクトログラムＧｔ，ωｍの蓄積信号小領域符号をｃ（Ｇｔ，ωｍ）と定義することとする。そして、帯域ωｍにおける各ｃ（Ｆｔｎ，ωｍ）、ｃ（Ｇｔ，ωｍ）には、符号ｑ（ωｍ，１）、ｑ（ωｍ，２），…，ｑ（ωｍ，ｚ）のいずれかが割振られることになる。
【００７１】
なお、これらの小領域スペクトログラムの符号化は、上記特徴抽出部２２において抽出された小領域参照信号スペクトログラム及び小領域蓄積信号スペクトログラム内の周波数スペクトル特徴(あるいは複数のスペクトル特徴を連続した一連の特徴ベクトルとしたもの)を、ヒストグラムを生成せずに直接ベクトル量子化し、各々参照信号小領域符号ならびに蓄積信号小領域符号を算出することも可能である。
【００７２】
索引生成部２５は、符号化部２４により出力された小領域蓄積信号スペクトログラムの各符号について、蓄積信号のどの時点に出現するか一覧化した索引情報を生成する（ステップＳ９２０）。索引情報は、図１６のように構成されており、符号ｑ（ωｍ，ｊ）が示すリスト（出現時点の配列；横方向の列)には、符号ｑ（ωｍ，ｊ）に分類された小領域蓄積信号スペクトログラムの蓄積信号上に出現する時点が、時系列に並べられて格納されている。
【００７３】
類似度計算部２３は、初めに小領域参照信号スペクトログラムに類似する小領域蓄積信号スペクトログラムを、参照信号小領域符号と蓄積信号小領域符号の符号間の類似度に基づいて索引情報から検出する。帯域ωｍにおける全ての小領域符号間の類似度は、符号の代表ベクトル間の距離から予め算出し、図１７に示す小領域符号間類似度表のように、小領域符号対毎にその類似度を表の上で定義、すなわち内部の記憶領域に記憶しておく。つまり、参照信号小領域符号ｑ（ωｍ,ｉ）と蓄積信号小領域符号ｑ（ωｍ,ｊ）との小領域符号間類似度υ（ωｍ,ｉ,ｊ）は、この表から一覧できるものである。
【００７４】
小領域符号間類似度υ（ωｍ,ｉ,ｊ）は、小領域符号ｑ（ωｍ,ｉ）とｑ（ωｍ,ｊ）に対応する代表ベクトル間の距離に基づいて算出され、類似度計算部２３により、距離値が小さい場合に、類似性が高いと判定され、距離値が大きい場合に類似性が低いものと判定される。本実施例においては、υ（ωｍ,ｉ,ｊ）を、０から１までの実数値として定義している。すなわち、帯域ωｍにおいて距離値が最大の場合には、υ（ωｍ,ｉ,ｊ）＝０であり、最小の場合には、υ（ωｍ,ｉ,ｊ）＝１となる。
【００７５】
類似度計算部２３は、参照信号小領域符号の集合から参照信号小領域符号を順次読み出して選択し（ステップＳ９３０）、図１７に示す小領域符号間類似度表を参照し（ステップＳ９３１）、当該参照信号小領域符号に類似する符号、すなわち小領域符号間類似度υ（ωｍ,ｉ,ｊ）が、所定の閾値である小領域用探索閾値ｓＰｔｈを超える符号を検出する（ステップＳ９３２）。ここでは、ある参照信号小領域符号について、同一の符号の値を有する符号を含む、少なくとも１つ以上の類似する符号を検出することになる。
【００７６】
なお、小領域用探索閾値ｓＰｔｈは、例えば、予め多数の音響信号において検出処理を繰り返し、検出漏れがない値を統計的に求めることにより設定することができる。また、小領域用探索閾値ｓＰｔｈは、全帯域で一定の値としているが、帯域ごとに異なる値を設定してもよい。
【００７７】
そして、類似度計算部２３は、ある一つの参照信号小領域スペクトログラムについて、その符号である参照信号小領域符号に類似する上記の少なくとも一つの符号の値を持つ小領域スペクトログラムが、蓄積信号中に存在するか否かを、図１６に示す索引情報により検索する（ステップＳ９３３）。蓄積信号中に、小領域参照信号スペクトログラムＦｔｎ，ωｍに類似する小領域蓄積信号スペクトログラムＧｔ＋ｔｎ,ωｍの存在が確認された場合、この小領域蓄積信号スペクトログラムＧｔ＋ｔｎ,ωｍを含む蓄積信号中の区間に、参照信号に対する類似性があると判定し、類似小領域蓄積信号スペクトログラムの蓄積信号中の出現位置ｔ＋ｔｎと、小領域参照信号スペクトル間の類似度ＳＰ（Ｆｔｎ，ωｍ, Ｇｔ＋ｔｎ,ωｍ）を出力する（ステップＳ９３４）。ここで、小領域スペクトル間の類似度すなわち小領域類似度は、次式（３）により定義される。
【００７８】
【数３】
2008015002.tif 000005
【００７９】
続いて、類似度計算部２３は、検出した類似蓄積信号小領域スペクトログラムの蓄積信号中の出現位置ｔ＋ｔｎから、類似性のある蓄積信号中の区間出現時点ｔを算出し（ステップＳ９３５）、小領域参照信号スペクトログラムとこれに類似する蓄積信号中の小領域スペクトログラムとの小領域類似度を、対応する区間出現時点tに、投票方式で積算する（ステップＳ９３６）。
【００８０】
類似度計算部２３は、上記ステップＳ９３０からステップＳ９３６の処理を繰り返し（ステップＳ９３７）、参照信号の全ての小領域スペクトログラムについて順次索引検出を行い、類似する小領域蓄積信号スペクトログラムが検出された場合に、この小領域蓄積信号スペクトログラムを含む蓄積信号中の区間出現時点ｔに対して小領域類似度を積算する。
【００８１】
類似度計算部２３は、各小領域参照信号スペクトログラムに類似する小領域蓄積信号スペクトログラムについて小領域類似度の積算処理が全て終了すると、次式（４）により、積算した小領域類似度を参照信号の全小領域参照信号スペクトログラム数で除算して正規化し、時点ｔに始まる蓄積信号中の区間と参照信号との区間類似度Ｓ（ｔ）を算出する（ステップＳ９４０）。
【００８２】
【数４】
2008015002.tif 000006
【００８３】
ここで、｜ＴＲ｜は時点ＴＲ（サンプリングされた時系列な離散的時点）の要素数を表し、｜ＴＷ｜は周波数帯域ωｍの要素数を表している。そして、類似度計算部２３は、以上のようにして算出された区間類似度Ｓ（ｔ）が、所定の値である探索閾値Ｓｔｈを超える蓄積信号上の時点ｔに始まる区間を、参照信号に類似する区間として出力する（ステップＳ９４１、ステップＳ９４２）。また、別の形態として、探索閾値を超えた複数の区間
を出力するかわりに、探索閾値を超え、もっとも大きな区間類似度を有する区間のみを出
力するようにしてもよい。
【００８４】
なお、探索閾値Ｓｔｈについては、実験的に任意で設定してもよいし、また別に、類似度計算部２３において算出された区間類似度Ｓ（ｔ）の統計分布から標準偏差を算出し、得られた区間類似度Ｓ（ｔ）の最大値に対して、−ｂσの値を探索閾値Ｓｔｈとして設定することも可能である。ここで、前記のｂは、実験を通じて適切な値が設定される。
【００８５】
（第２実施形態の実装例）
次に、上述した第２実施形態の実装例について説明する。
音響信号は、音楽ＣＤなどで用いられる音響信号（サンプリング周波数４４．１ｋＨｚ、１６ｂｉｔ量子化、ステレオ）をダウンサンプリング、モノラル化、量子化ビット数を変更した音響信号（サンプリング周波数１１．０２５ｋＨｚ、８ビット量子化、モノラル）を用いる。そして、この音響信号に対し、まず、１０２４点のＦＦＴ（Fast Fourier Transform）を１０ｍｓごとに行い、１０ｍｓごとにその時点のパワースペクトルを抽出する。そして、この１０ｍｓごとにスペクトルを時間方向に並べたものとして、スペクトログラムを得る。そして、周波数軸方向でのスペクトログラムの分割数は本実施形態に記載の装置で利用される４とし、各帯域において、（（１／ｘ３）×ｆｃ，（１／ｘ２）×ｆｃ，（１／ｘ）×ｆｃ，ｘ×ｆｃ，ｘ２×ｆｃ，ｘ３×ｆｃ）の各周波数においてパワーを抽出し、抽出した各周波数のパワーの平均値を算出し、算出した平均値で各周波数のパワーを正規化したものを特徴ベクトルとする。ここで、ｘは７１／２７としている。
【００８６】
そして、ｆｃは各帯域で各々６２０Ｈｚ、１０３０Ｈｚ、１７１０Ｈｚ、２８４０Ｈｚを中心（各帯域の浮動中心周波数とここでは記載）にして浮動的に設定するとし、各帯域で先に抽出したスペクトルのうち、浮動中心周波数から上下１０％の周波数のうちでパワーが最大となる周波数をｆｃとする。このｆｃに基づき、各時点で、（（１／ｘ３）×ｆｃ，（１／ｘ２）×ｆｃ，（１／ｘ）×ｆｃ，ｘ×ｆｃ，ｘ２×ｆｃ，ｘ３×ｆｃ）の各周波数におけるパワーを、ＦＦＴの結果から得られる各周波数のうち、これらの周波数に最も近い周波数のパワーを（（１／ｘ３）×ｆｃ，（１／ｘ２）×ｆｃ，（１／ｘ）×ｆｃ，ｘ×ｆｃ，ｘ２×ｆｃ，ｘ３×ｆｃ）でのパワーとし、各周波数のパワーの平均値を算出し、算出した平均値で各周波数のパワーを正規化したものを特徴ベクトルとして抽出する。
【００８７】
また、特徴ベクトルの各成分の対応する周波数について、その周波数を中心周波数とするバンドパスフィルタをもちいて、再度パワースペクトルを抽出してもよい。
【００８８】
このようにして、特徴ベクトルを抽出することで、上下１０％までピッチが変動した音響信号でも探索を行うことができる。
【００８９】
なお、これらのスペクトログラム分割数（上記実装例では４）、特徴ベクトルの次元（上記実装例では７）やｘの値（上記実装例では７１／２７）、ｆｃの設定範囲（上記実装例では１０％）は、探索の目的とする信号やピッチ変動の幅（上記実装例では上下１０％）に応じて、変更可能なものであり、具体的な数値に限られないものとする。
【００９０】
また、特徴ベクトルの各成分が対応する周波数について、上記実装例では、第１成分から等比数列としたが、各成分間の比が決まっていれば、それらが等比である必要はない。例えば、第２成分の対応する周波数は第１成分が対応する周波数の１．２倍の周波数、第３成分の対応する周波数は第１成分の対応する周波数の２倍の周波数などとしてもよい。
【００９１】
上記の第１及び第２実施形態の構成により、参照信号および蓄積信号の各時間周波数スペクトログラムから、パワースペクトルの特徴的な部分を検出して各時点における基準周波数を浮動的に設定、すなわち、周波数による変動が発生していたとしても、スペクトログラム全体の特徴を抽出可能な基準周波数を設定することが可能となる。そして、該基準周波数に基づき各信号の周波数スペクトル特徴を抽出し、抽出した周波数スペクトル特徴に基づいて周波数の上下の変動に対応した音響信号検出が可能となる。
【００９２】
また、第２の実施形態の構成により、各信号の時間周波数スペクトログラムを時間・周波数方向において小領域に分割し、小領域におけるピークを検出して特徴抽出を行い、参照信号と蓄積信号中のある区間との類似度を各小領域の類似性から判断するように構成することで、周波数変動がなされた音響信号に対し更に大きな妨害音が重畳していても精度よく、かつ高速に信号検出処理を行うことを可能としている。
【００９３】
また、第１及び第２実施形態の構成により、各信号の特徴を符号の集合として表現することにより、照合に掛かる情報量を大幅に低減することが可能となり、類似度の計算処理を高速化することが可能となる。更に、蓄積信号の符号を索引情報として予め用意しておくことで、類似性のある符号の位置を容易に検出できるようになり、参照信号に類似した信号を含有する蓄積信号中の区間検出をさらに高速に行うことを可能としている。
【００９４】
また、上記の実施形態によれば、テレビ放送におけるドラマの挿入曲や背景音楽のように、ストーリーの雰囲気に合わせて加工された音響信号についても検索が可能となり、マルチメディア放送のモニタリングが高速かつ高精度に行えるようになる。また、テレビ、ラジオで放送された音楽の楽曲情報やＣＭの情報を、多チャンネル同時にリアルタイムで検索しプレイリストを作成することが可能である。プレイリストは、例えば、放送局が視聴者からの問い合わせに対応するための広報資料、ある広告主が競合会社のＣＭ露出状況などを調査することのできるＣＭ確認用データ、レコード会社のマーケティングデータ、著作権料徴収のための基礎データとして幅広く活用できるものである。
【００９５】
また、上記の実施形態の構成により、著作権侵害物の検査、未知の音響情報の識別、インターネット著作権防止アプリケーションにも応用可能であり、メタデータ解析・配信サービスといった市場への展開が期待される。
【００９６】
また、本発明に記載の音響信号特徴抽出装置とは、音響信号検索装置１、１ａ、１ｂにおける基準周波数検出部１１、特徴抽出部１２を備えた構成であり、また、音響信号検索装置２における基準周波数検出部２１、特徴抽出部２２を備えた構成に該当する装置である。また、本発明に記載の参照音響信号は、参照信号に対応し、検索対象音響信号は、蓄積信号に対応する。
【００９７】
なお、上記の実施形態において、検出とは検索・認識を含み、特徴抽出は音響信号の分析を含む。また、上記の実施形態に記載の「各時点」とは、時間周波数スペクトログラムの時間間隔ごとの各時点の意味であり、時間軸において、時間周波数スペクトログラムの算出元となる一定の時間範囲を示すものとする。
【００９８】
また、実際に上述した音響信号検索装置１、１ａ、１ｂ、２を構築する際には、音響信号検索装置の機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータで実行することにより、蓄積信号のデータベースへの登録処理および入力信号検出処理を行っても良い。ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとずる。また、「コンピュータシステム」は、ホームページ提供環境（あるいは表示環境）を備えたＷＷＷシステムも含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ-ＲＯＭ等の可搬媒体、コンピュータに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムが送信された場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリ（ＲＡＭ）のように、一定時間プログラムを保持しているものも含むものとする。
【００９９】
また、上述したプログラムは、このプログラムを記憶装置等に格納したコンピュータシステムから、伝送媒体を介して、あるいは、伝送媒体中の伝送波により他のコンピュータシステムに伝送されてもよい。ここで、プログラムを伝送する「伝送媒体」は、インターネット等のネットワーク（通信網）や電話回線等の通信回線（通信線）のように情報を伝送する機能を有する媒体のことをいう。また、上記プログラムは、前述した機能の一部を実現するためのものであってもよい。さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル（差分プログラム）であってもよい。
【０１００】
また、本発明は、上述の実施形態によって限定されず、本発明の技術思想の範囲内で様々な変形が可能である。
【図面の簡単な説明】
【０１０１】
【図１】本発明の原理を説明するためのスペクトログラムの概略図（その１）である。
【図２】本発明の原理を説明するためのスペクトログラムの概略図（その２）である。
【図３】本発明の原理を説明するためのスペクトログラムの概略図（その３）である。
【図４】本発明の原理を説明するためのスペクトログラムの概略図（その４）である。
【図５】本発明の第１実施形態に係る音響信号検索装置の概略ブロック図である。
【図６】同実施形態に係る音響信号検索装置の動作を示すフローチャート（その１）である。
【図７】同実施形態に係る基準周波数検出の概要を説明する概念図である。
【図８】同実施形態に係るＴＡＳ法を説明する概念図である。
【図９】同実施形態に係る音響信号検索装置の動作を示すフローチャート（その２）である。
【図１０】同実施形態に係る音響信号検索装置の動作を示すフローチャート（その３）である。
【図１１】同実施形態に係る音響信号検索装置の他の実施形態（その１）である。
【図１２】同実施形態に係る音響信号検索装置の他の実施形態（その２）である。
【図１３】第１実施形態の他の実施形態（その２）の音響信号検索装置の動作を示すフローチャートである。
【図１４】本発明の第２実施形態に係る音響信号検索装置の概略ブロック図である。
【図１５】同実施形態に係る音響信号検索装置の動作を示すフローチャートである。
【図１６】同実施形態に係る索引情報の構成を説明するための概念図である。
【図１７】同実施形態に係る小領域符号間類似度表の構成を説明するための概念図である。
【符号の説明】
【０１０２】
１    音響信号検索装置
１１  基準周波数検出部
１２  特徴抽出部
１３  類似度計算部
【図１】
2008015002.tif 000007
【図２】
2008015002.tif 000008
【図３】
2008015002.tif 000009
【図４】
2008015002.tif 000010
【図５】
2008015002.tif 000011
【図６】
2008015002.tif 000012
【図７】
2008015002.tif 000013
【図８】
2008015002.tif 000014
【図９】
2008015002.tif 000015
【図１０】
2008015002.tif 000016
【図１１】
2008015002.tif 000017
【図１２】
2008015002.tif 000018
【図１３】
2008015002.tif 000019
【図１４】
2008015002.tif 000020
【図１５】
2008015002.tif 000021
【図１６】
2008015002.tif 000022
【図１７】
2008015002.tif 000023