"蝸牛遅延特性を模擬する蝸牛遅延フィルタを用いて、デジタルデータである音響信号に位相変調を施し、前記位相変調が施された音響信号に電子透かしデータを埋め込む電子透かしデータ埋め込み装置によって、デジタルデータである音響信号に電子透かしデータが埋め込まれた場合に、前記蝸牛遅延フィルタが模擬する蝸牛遅延特性を推定する蝸牛遅延特性推定手段と、前記蝸牛遅延特性推定手段により推定された蝸牛遅延特性に基づいて、音響信号に埋め込まれた前記電子透かしデータを検出する電子透かし検出手段とを備える、電子透かし検出装置。前記電子透かしデータ埋め込み装置が、複数の異なる蝸牛遅延フィルタを用いて音響信号に位相変調を施すことにより、複数の異なる位相変調された音響信号を生成し、電子透かしデータに応じて、前記複数の異なる位相変調された音響信号の中から一の音響信号を選択し、選択した音響信号同士を接合することにより、電子透かしデータを埋め込むように構成されており、前記蝸牛遅延特性推定手段が、前記複数の異なる蝸牛遅延フィルタがそれぞれ模擬する複数の異なる蝸牛遅延特性を推定するように構成され、前記電子透かし検出手段が、前記蝸牛遅延特性推定手段により推定された前記複数の異なる蝸牛遅延特性に基づいて、電子透かしデータが埋め込まれた音響信号が、前記複数の異なる蝸牛遅延フィルタのうちの何れの蝸牛遅延フィルタが適用されて位相変調が施されたかを判定することにより、電子透かしデータを検出するように構成されている、請求項1に記載の電子透かし検出装置。前記蝸牛遅延特性推定手段が、前記蝸牛遅延フィルタの零点を推定することにより、蝸牛遅延特性を推定するように構成されている、請求項1又は2に記載の電子透かし検出装置。前記蝸牛遅延特性推定手段が、チャープz変換を用いて、前記蝸牛遅延フィルタの零点を推定するように構成されている、請求項3に記載の電子透かし検出装置。前記蝸牛遅延特性手段により推定された蝸牛遅延特性の逆特性を有するフィルタを電子透かしデータが埋め込まれた音響信号に施すことにより、電子透かしデータが埋め込まれる前の音響信号を取得する原信号取得手段をさらに備える、請求項4に記載の電子透かし検出装置。電子透かしデータが埋め込まれた音響信号の位相変調に適用されたと前記電子透かし検出手段により判定された蝸牛遅延フィルタの逆フィルタを当該音響信号に施すことにより、電子透かしデータが埋め込まれる前の音響信号を取得する原信号取得手段をさらに備える、請求項2に記載の電子透かし検出装置。蝸牛遅延特性を模擬する蝸牛遅延フィルタを用いて、デジタルデータである音響信号に位相変調を施し、前記位相変調が施された音響信号に電子透かしデータを埋め込む電子透かしデータ埋め込み装置によって、デジタルデータである音響信号に電子透かしデータが埋め込まれた場合に、前記蝸牛遅延フィルタが模擬する蝸牛遅延特性を推定するステップaと、推定された蝸牛遅延特性に基づいて、音響信号に埋め込まれた前記電子透かしデータを検出するステップbとを有する、電子透かし検出方法。前記電子透かしデータ埋め込み装置が、複数の異なる蝸牛遅延フィルタを用いて音響信号に位相変調を施すことにより、複数の異なる位相変調された音響信号を生成し、電子透かしデータに応じて、前記複数の異なる位相変調された音響信号の中から一の音響信号を選択し、選択した音響信号同士を接合することにより、電子透かしデータを埋め込むように構成されており、前記ステップaにおいて、前記複数の異なる蝸牛遅延フィルタがそれぞれ模擬する複数の異なる蝸牛遅延特性を推定し、前記ステップbにおいて、前記ステップaにより推定された前記複数の異なる蝸牛遅延特性に基づいて、電子透かしデータが埋め込まれた音響信号が、前記複数の異なる蝸牛遅延フィルタのうちの何れの蝸牛遅延フィルタが適用されて位相変調が施されたかを判定することにより、電子透かしデータを検出する、請求項7に記載の電子透かし検出方法。前記ステップaにおいて、前記蝸牛遅延フィルタの零点を推定することにより、蝸牛遅延特性を推定する、請求項8又は9に記載の電子透かし検出方法。前記ステップaにおいてチャープz変換を用いて、前記蝸牛遅延フィルタの零点を推定する、請求項9に記載の電子透かし検出方法。一般に、電子音響透かし技術では、知覚不可能性埋め込み情報が利用者に知覚されず、埋め込みによる原信号の知覚可能な歪みが生じないこと、頑健性通常の信号変換処理及び埋め込み情報を削除するといった悪意のある攻撃に対して影響を受けないこと、及び秘匿性情報が埋め込まれていることに気付かせないこと、気付かれてもその情報を容易に検出されないことが要求されている。上記1のLSB法は、振幅情報に大きく影響を与えない下位ビットに情報を埋め込むため、知覚不可能性を満たすが、ビット変化に敏感なため頑健性に問題がある。また、上記2のDSS法の場合、スペクトル全体に情報を埋め込むため、信号変形処理には頑健であるが、埋め込まれた情報を容易に知覚できるため知覚不可能性に問題がある。上記3のECHO法は、エコー時間及び1次反射音の振幅を調整することで歪みがなく、知覚不可能な埋め込みを実現できるが、自己相関法及びケプストラム処理を利用することで透かし情報を容易に検出除去できるため、上記の従来の方法の中でもっとも頑健性秘匿性に欠ける。また、上記4のPPM法は、周期的な位相変調が比較的知覚され難いという聴覚特性に基づいているが、位相変調が高い周波数成分の位相スペクトルをランダムに歪ませるため、知覚不可能性に問題がある。他方、上記のCD法の場合、知覚不可能性、秘匿性、及び頑健性を十分に満足するものの、埋め込まれた情報を検出するために原信号を参照することが必要となるため、応用範囲が限定されるという問題がある。本発明は斯かる事情に鑑みてなされたものであり、その主たる目的は、CD法により埋め込まれた情報を、原信号を参照することなく検出することができる電子透かし検出装置及び電子透かし検出方法を提供することにある。上述した課題を解決するために、本発明の一の態様の電子透かし検出装置は、蝸牛遅延特性を模擬する蝸牛遅延フィルタを用いて、デジタルデータである音響信号に位相変調を施し、前記位相変調が施された音響信号に電子透かしデータを埋め込む電子透かしデータ埋め込み装置によって、デジタルデータである音響信号に電子透かしデータが埋め込まれた場合に、前記蝸牛遅延フィルタが模擬する蝸牛遅延特性を推定する蝸牛遅延特性推定手段と、前記蝸牛遅延特性推定手段により推定された蝸牛遅延特性に基づいて、音響信号に埋め込まれた前記電子透かしデータを検出する電子透かし検出手段とを備える。この態様において、前記電子透かしデータ埋め込み装置が、複数の異なる蝸牛遅延フィルタを用いて音響信号に位相変調を施すことにより、複数の異なる位相変調された音響信号を生成し、電子透かしデータに応じて、前記複数の異なる位相変調された音響信号の中から一の音響信号を選択し、選択した音響信号同士を接合することにより、電子透かしデータを埋め込むように構成されており、前記蝸牛遅延特性推定手段が、前記複数の異なる蝸牛遅延フィルタがそれぞれ模擬する複数の異なる蝸牛遅延特性を推定するように構成され、前記電子透かし検出手段が、前記蝸牛遅延特性推定手段により推定された前記複数の異なる蝸牛遅延特性に基づいて、電子透かしデータが埋め込まれた音響信号が、前記複数の異なる蝸牛遅延フィルタのうちの何れの蝸牛遅延フィルタが適用されて位相変調が施されたかを判定することにより、電子透かしデータを検出するように構成されていてもよい。また、前記態様において、前記蝸牛遅延特性推定手段が、前記蝸牛遅延フィルタの零点を推定することにより、蝸牛遅延特性を推定するように構成されていてもよい。また、前記態様において、前記蝸牛遅延特性推定手段が、チャープz変換を用いて、前記蝸牛遅延フィルタの零点を推定するように構成されていてもよい。また、前記態様において、前記蝸牛遅延特性手段により推定された蝸牛遅延特性の逆特性を有するフィルタを電子透かしデータが埋め込まれた音響信号に施すことにより、電子透かしデータが埋め込まれる前の音響信号を取得する原信号取得手段をさらに備えていてもよい。また、前記態様において、電子透かしデータが埋め込まれた音響信号の位相変調に適用されたと前記電子透かし検出手段により判定された蝸牛遅延フィルタの逆フィルタを当該音響信号に施すことにより、電子透かしデータが埋め込まれる前の音響信号を取得する原信号取得手段をさらに備えていてもよい。本発明の一の態様の電子透かし検出方法は、蝸牛遅延特性を模擬する蝸牛遅延フィルタを用いて、デジタルデータである音響信号に位相変調を施し、前記位相変調が施された音響信号に電子透かしデータを埋め込む電子透かしデータ埋め込み装置によって、デジタルデータである音響信号に電子透かしデータが埋め込まれた場合に、前記蝸牛遅延フィルタが模擬する蝸牛遅延特性を推定するステップaと、推定された蝸牛遅延特性に基づいて、音響信号に埋め込まれた前記電子透かしデータを検出するステップbとを有する。この態様において、前記電子透かしデータ埋め込み装置が、複数の異なる蝸牛遅延フィルタを用いて音響信号に位相変調を施すことにより、複数の異なる位相変調された音響信号を生成し、電子透かしデータに応じて、前記複数の異なる位相変調された音響信号の中から一の音響信号を選択し、選択した音響信号同士を接合することにより、電子透かしデータを埋め込むように構成されており、前記ステップaにおいて、前記複数の異なる蝸牛遅延フィルタがそれぞれ模擬する複数の異なる蝸牛遅延特性を推定し、前記ステップbにおいて、前記ステップaにより推定された前記複数の異なる蝸牛遅延特性に基づいて、電子透かしデータが埋め込まれた音響信号が、前記複数の異なる蝸牛遅延フィルタのうちの何れの蝸牛遅延フィルタが適用されて位相変調が施されたかを判定することにより、電子透かしデータを検出するようにしてもよい。また、前記態様において、前記ステップaにおいて、前記蝸牛遅延フィルタの零点を推定することにより、蝸牛遅延特性を推定するようにしてもよい。また、前記態様において、前記ステップaにおいてチャープz変換を用いて、前記蝸牛遅延フィルタの零点を推定するようにしてもよい。本発明に係る電子透かし検出装置及び電子透かし検出方法によれば、原信号を参照することなく、CD法により埋め込まれた電子透かしデータを検出することができる。以下、本発明の好ましい実施の形態を、図面を参照しながら説明する。なお、以下に示す各実施の形態は、本発明の技術的思想を具体化するための方法及び装置を例示するものであって、本発明の技術的思想は下記のものに限定されるわけではない。本発明の技術的思想は、特許請求の範囲に記載された技術的範囲内において種々の変更を加えることができる。本実施の形態に係る電子透かし検出装置は、原信号に埋め込まれた電子透かしデータをその原信号を参照することなく検出することができる装置である。このように原信号を参照することなく電子透かしデータを検出することを、本明細書ではブラインド検出と称する。以下、この電子透かし検出装置と、電子透かしデータを埋め込む電子透かし埋込装置について説明する。[電子透かし埋込装置の構成]図1は、本発明の実施の形態に係る電子透かし埋込装置の構成を示すブロック図である。図1に示すように、電子透かし埋込装置1は、CPU11、ROM12、RAM13、信号入力部14、信号出力部15及びハードディスク16を備えており、これらのCPU11、ROM12、RAM13、信号入力部14、信号出力部15及びハードディスク16は、バス17によって接続されている。CPU11は、ROM12及びハードディスク16に記憶されているコンピュータプログラムを実行する。これにより、電子透かし埋込装置1は、後述するような動作を実行し、音響信号に対する電子透かしデータの埋め込みを実現する。ROM12は、マスクROM、PROM、EPROM、又はEEPROM等によって構成されており、CPU11によって実行されるコンピュータプログラム及びこれに用いられるデータ等を記憶している。RAM13は、SRAMまたはDRAM等によって構成されており、ハードディスク16に記憶されているプログラムの読み出しに用いられる。また、RAM13は、CPU11がコンピュータプログラムを実行するときに、CPU11の作業領域としても利用される。信号入力部14は、処理対象となる原信号である音響信号及びその音響信号に埋め込まれる電子透かしデータの入力を、外部の装置から受け付ける。また、信号出力部15は、電子透かしデータが埋め込まれた音響信号以下、透かし入り音響信号というを外部の装置へ出力する。なお、本実施の形態においては、原信号である音響信号はデジタルデータである。しかし、当該音響信号がアナログデータであってもよく、その場合は、A/D変換機能を備えた信号入力部14が、入力された音響信号をA/D変換することによりデジタルデータに変換した上で、その後の処理を行うようにすればよい。ハードディスク16には、オペレーティングシステム及びアプリケーションプログラム等、並びにCPU11に実行させるための種々のコンピュータプログラムおよび当該コンピュータプログラムの実行に用いられるデータ等がインストールされている。このコンピュータプログラムには、電子透かしデータの埋め込みを行うための電子透かし埋込プログラム16Aが含まれる。ハードディスク16にインストールされる電子透かし埋込プログラム16Aは、フレキシブルディスクドライブ、CD-ROMドライブ、またはDVD-ROMドライブ等の外部記憶装置図示せずを介して、可搬型記録媒体から読み出される。なお、このように可搬型記録媒体によって提供されるのみならず、電気通信回線有線、無線を問わないを介して電子透かし埋込装置1と通信可能に接続された外部の装置から電子透かし埋込プログラム16Aを提供することも可能である。例えば、電子透かし埋込プログラム16Aがインターネット上のサーバコンピュータのハードディスク内に格納されている場合において、このサーバコンピュータに電子透かし埋込装置1がアクセスして、当該コンピュータプログラムをダウンロードし、これをハードディスク16にインストールすることも可能である。ハードディスク16には、例えば米マイクロソフト社が製造販売するWindows登録商標等のマルチタスクオペレーティングシステムがインストールされている。以下の説明において、本実施の形態に係る電子透かし埋込プログラム16Aは当該オペレーティングシステム上で動作するものとする。次に、上記の電子透かし埋込装置1の構成を、図2に示す機能ブロック図を参照しながら説明する。なお、以下において、nはサンプリング番号を、kは音響信号のフレーム番号をそれぞれ示している。図2に示すように、電子透かし埋込装置1は、音響信号xnをフレーム化するフレーム処理部101と、2つの蝸牛遅延フィルタ102a及び102bと、電子透かしデータskの値に応じて第1蝸牛遅延フィルタ102a及び第2蝸牛遅延フィルタ102bの何れかを選択するフィルタ選択部103とを備えている。フィルタ選択部103は、電子透かしデータのビット値が0である場合に第1蝸牛遅延フィルタ102aを選択し、同じく1である場合に第2蝸牛遅延フィルタ102bを選択する。第1蝸牛遅延フィルタ102a及び第2蝸牛遅延フィルタ102bでは、後述するようにして音響信号に群遅延を与える。このようにして群遅延が付与された音響信号が統合され、電子透かしデータが埋め込まれた音響信号である透かし入り音響信号ynが生成される。なお、本実施の形態において、これらのフレーム処理部101、第1蝸牛遅延フィルタ102a及び第2蝸牛遅延フィルタ102b、並びにフィルタ選択部103は、CPU11が電子透かし埋込プログラム16Aを実行することによって実現される。[蝸牛遅延フィルタ]以下、第1蝸牛遅延フィルタ102a及び第2蝸牛遅延フィルタ102bの詳細について説明する。これらの第1蝸牛遅延フィルタ102a及び第2蝸牛遅延フィルタ102bは、人間の聴覚の蝸牛遅延特性を模擬したデジタルフィルタであり、具体的には、振幅成分にはまったく影響を与えず、位相特性のみを変化させる全域通過フィルタで構成される。本実施の形態において、蝸牛遅延フィルタ102a及び102bは、以下の式1の伝達関数Hzにより定義される1次の無限インパルス応答型全域通過フィルタで構成される。ここで、bmはHmzのフィルタ係数を表している。このように、第1蝸牛遅延フィルタ102a及び第2蝸牛遅延フィルタ102bを1次の無限インパルス応答型全域通過フィルタで構成することにより、高速な処理が可能になる。なお、無限インパルス応答型全域通過フィルタの群遅延特性が蝸牛遅延特性をより正確に表していれば、フィルタ次数は1次以上であってもよく、また、フィルタのカスケード段数は1段以上であってもよい。第1蝸牛遅延フィルタ102a及び第2蝸牛遅延フィルタ102bにより与えられる群遅延γmωは以下の式2により算出される。図3は、本発明の実施の形態1における電子透かし埋込装置1が備える第1蝸牛遅延フィルタ102a及び第2蝸牛遅延フィルタ102bの特性を示すグラフである。図3において、縦軸は群遅延を、横軸は音響信号の周波数をそれぞれ示している。図3において、細い実線は、人間の聴覚における蝸牛遅延を1/10倍に縮小した蝸牛遅延特性を示している。また、太い実線は、フィルタ係数b=0795の場合に上記式1により定義される第1蝸牛遅延フィルタ102aの特性を示し、破線は、フィルタ係数b=0865の場合に同じく定義される第2蝸牛遅延フィルタ102bの特性を示している。なお、図3において細い実線で示されている蝸牛遅延特性は、TDauOWegnerVMellertandBKollmeierAuditorybrainstemresponseswithoptimizedchirpsignalscompensatingbasilarmembranedispersionJAcoustSocAm107153015402000を参考にして定めたものである。以上より、第1蝸牛遅延フィルタ102a及び第2蝸牛遅延フィルタ102bを音響信号にかけると、実際の蝸牛遅延の1/10倍の蝸牛遅延を当該音響信号に付与することになる。したがって、人間の実際の蝸牛遅延特性を近似するためには、このような蝸牛遅延フィルタを10段カスケード接続する必要がある。しかし、実際と同様の蝸牛遅延量を音響信号に与えることにすると、その音響信号を知覚する際の群遅延量は実際の蝸牛遅延量の2倍になってしまうため、遅延が大きすぎると考えられる。そこで、本実施の形態においては、上記のように実際の蝸牛遅延の1/10倍の蝸牛遅延を音響信号に与えることにしている。本実施の形態において、第1蝸牛遅延フィルタ102a及び第2蝸牛遅延フィルタ102bはそれぞれ、下記の式3及び式4にしたがって、原信号である音響信号xnに蝸牛遅延パターンを付与し、中間信号w0n及びw1nを得る。そして、フィルタ選択部103が、電子透かしデータのビット値に応じて、フレーム毎に中間信号w0n及びw1nを選択統合することにより、下記の式5に示す透かし入り音響信号ynを取得する。ただし、k-1ΔW<n≦kΔWを満足する。ここで、ΔW=fs/Nbitはフレーム長であり、fsは原信号のサンプリング周波数を、Nbitは1秒あたりの情報埋込ビットレートをそれぞれ表している。[電子透かし検出装置の構成]図4は、本発明の実施の形態に係る電子透かし検出装置の構成を示すブロック図である。図4に示すように、電子透かし検出装置2は、上記の電子透かし埋込装置1と同様に、CPU21、ROM22、RAM23、信号入力部24、及びハードディスク25を備えており、これらのCPU21、ROM22、RAM23、信号入力部24及びハードディスク25は、バス26によって接続されている。CPU21、ROM22及びRAM23のそれぞれについては、電子透かし埋込装置1が備えるCPU11、ROM12及びRAM13と同様であるので、説明を省略する。信号入力部24は、透かし入り音響信号の入力を外部の装置から受け付ける。この透かし入り音響信号は、信号入力部24に対して電子透かし埋込装置1から直接入力されてもよく、他の装置及び/又は通信ネットワーク等を介して入力されてもよい。ハードディスク25には、電子透かし埋込装置1の場合と同様に、オペレーティングシステム及びCPU21に実行させるための種々のコンピュータプログラム等がインストールされている。このコンピュータプログラムには、電子透かしデータの検出を行うための電子透かし検出プログラム25Aが含まれる。電子透かし埋込プログラム16Aの場合と同様に、ハードディスク25にインストールされる電子透かし検出プログラム25Aは、可搬型記録媒体によって提供されてもよく、電気通信回線を介して提供されてもよい。また、この電子透かし検出プログラム25Aは、電子透かし埋込プログラム16Aの場合と同様に、ハードディスク25にインストールされているオペレーティングシステム上で動作するものとする。次に、上記の電子透かし検出装置2の構成を、図5に示す機能ブロック図を参照しながら説明する。図5に示すように、電子透かし検出装置2は、電子透かし埋込装置1により生成された透かし入り音響信号ynをフレーム化するフレーム処理部201と、フレーム化された透かし入り音響信号ynに対して、チャープz変換を施す2つのチャープz変換部202a及び202bと、これらの第1チャープz変換部202a及び第2チャープz変換部202bによるチャープz変換の結果に基づいて、電子透かしデータのビット値を検出するビット値検出部203とを備えている。なお、本実施の形態において、これらのフレーム処理部201、第1チャープz変換部202a及び第2チャープz変換部202b、並びにビット値検出部203は、CPU21が電子透かし検出プログラム25Aを実行することによって実現される。[チャープz変換]第1チャープz変換部202a及び第2チャープz変換部202bが実行するチャープz変換CZTは、周波数スペクトルのフレキシブルな分析を可能とする手法として知られ例えば、WangTTThesegmentedchirpztransformanditsapplicationinspectrumanalysisIEEETransInstrumentationandmeasurement393183231990を参照、高速フーリエ変換FFTの実装にも活用されている。このチャープz変換は、離散フーリエ変換DFTと比較して、周波数分解能及び周波数応答のダイナミックレンジを自由に変えられるという特徴を有している。また、z平面上で任意のM点でのz変換を効率良く求めることができるという特徴も有している。一般に、チャープz変換は、z=rexpjωnでN点のDFTと結ばれる大きさr=1で正規化周波数ωn=2πn/Nのとき単位円周上のDFTと等価である関係にある。ここで、チャープz変換は下記の式6により表される。ただし、A=A0expj2πθ0、W=W0expj2πφ0である。ここで、θ0及びφ0は初期位相である。上述したように、A=1、M=N、W=exp-j2π/Nのとき、CZTはDFTに一致する。[ブラインド検出の原理]本実施の形態では、上記のチャープz変換を用いることにより、第1蝸牛遅延フィルタ102a及び第2蝸牛遅延フィルタ102bを用いて音響信号に埋め込まれた電子透かしデータのブラインド検出を実現する。以下、このブラインド検出の原理について説明する。第1蝸牛遅延フィルタ102a及び第2蝸牛遅延フィルタ102bの極及び零点は、図6に示すとおりに配置される。これらの蝸牛遅延フィルタ102a及び102bは、上述したように1次IIR型全域通過フィルタであり、その特徴として極図6中の×及び零点図6中の○は中心点から単位円に向かって垂線を描いたときに交差する半径及びその逆数bm及び1/bmの関係にある。一般に、bmの値が減少するにしたがい、極は中心点に近付き、零点は単位円から外側に向かって離れていく。反対に、bmの値が増加するにしたがい、極及び零点は互いに単位円に向かって近付いていく。この場合の群遅延量は、図3に示すように、bmの値の増加とともに増加する。なお、図6において、太字の○及び×は第1蝸牛遅延フィルタ102aの曲及び零点をそれぞれ示し、細字の○及び×は第2蝸牛遅延フィルタ102bの曲及び零点をそれぞれ示している。透かし入り音響信号ynは、上述したような遅延情報が埋め込まれた信号として観測されることになる。そのため、ynから遅延情報、すなわち遅延情報の付与に利用された蝸牛遅延フィルタの極及び零点の位置を推定することにより、ブラインド検出を実現することができる。なお、原信号xn自体も数列の特性として極及び零点を持つため音源が有界であるとして、その信号の減衰に関係する極など、観測信号ynから仮に極及び零点の位置を推定できたとしても、それはIIR型全域通過フィルタ蝸牛遅延フィルタによって付与されたものなのか、原信号そのものが持つものであるのかを見極める必要がある。チャープz変換を用いることにより、蝸牛遅延フィルタの極及び零点の位置を推定することができることを示すために、上記の式1の蝸牛遅延フィルタの零点r=1/bmを通るようにrを選択して、原信号xn及び遅延情報を埋め込んだ信号ynをチャープz変換A=r、M=N、W=exp-j2π/Nすることにより周波数分析を行う。以下、原信号である楽器音をxnとし、第1蝸牛遅延フィルタ102a及び第2蝸牛遅延フィルタ102bを利用してAISLabの電子透かしデータを埋め込んだ信号をynとする。ここでは、第1蝸牛遅延フィルタ102a及び第2蝸牛遅延フィルタ102bはいずれも直流成分のところに極及び零点を配置しており、r=1/b0又はr=1/b1としたチャープz変換の周波数分析を行う。なお、サンプリング周波数は44.1kHz、ビットレートはNbit=4bpsとして、1フレーム250msに1ビット相当の遅延情報を埋め込むものとする。図7は、その分析結果を示すグラフである。図7a乃至iは、左から右にフレーム1でのxn、フレーム1でのyn、フレーム2のynの周波数スペクトルを、上から下にr=1、r=1/b0、r=1/b1でのチャープz変換により分析した結果をそれぞれ示している。図7gに示すように、xnに関する分析結果では、極及び零点配置の周波数付近でのスペクトルには特段変化がみられない。他方、フレーム1のynではr=1/b1でのチャープz変換の結果図9hにおいて、フレーム2のynではr=1/b0でのチャープz変換の結果図9fにおいて、最も低い周波数領域直流成分から低周波数域までの範囲;例えば図3に示す遅延が見られる周波数帯のところでスペクトル成分が劇的に減少していることがわかる図中の矢印で示す箇所。これは、零点の影響によるディップくぼみに対応しているため、原理的にはその大きさは-∞dBになる。それ以外の分析r=1、r=1/b0フレーム1の場合、及びr=1/b1フレーム2の場合では、最も低い周波数のところでスペクトル成分の変化はほとんど見られないすなわち、-∞dB線形で0に近付かない。なお、この結果に関しては、他のフレーム及び他の対象信号でも同様のことが起こることが確認されている。以上より、対象信号に係わらず、蝸牛遅延フィルタの零点を交差するようにz平面上の軌跡に沿ってチャープz変換を行うことにより、ynから蝸牛遅延フィルタの零点の位置を推定することが可能であることが分かる。なお、原理的には、rを零点ではなく極の値にしてチャープz変換を行うことも可能である極の場合は∞dBのスペクトルピークを得ることになるが、計算機上でのダイナミックレンジ内のオーバーフローを検出しなければならないため、零点を用いる方が望ましい。零点を利用する場合は、ダイナミックレンジ内の0を探せばよいため、より容易な処理で足りることになる。本実施の形態では、第1チャープz変換部202aがr=1b0のz平面上の軌跡に沿ったチャープz変換を行い、第2チャープz変換部202bがr=1b1のz平面上の軌跡に沿ったチャープz変換を行う。これらのチャープz変換の結果を用いることにより、対象信号が、第1蝸牛遅延フィルタ102aフィルタ係数b0及び第2蝸牛遅延フィルタ102bフィルタ係数b1の何れにより群遅延が与えられたものであるのかを推定することが可能になる。[電子透かし埋込装置1及び電子透かし検出装置2の動作]次に、上述したように構成された本実施の形態の電子透かし埋込装置1及び電子透かし検出装置2の動作について、図8及び図9に示すフローチャートと図2及び図5とを参照しながら説明する。[電子透かし埋込処理]図8は、本発明の実施の形態における電子透かし埋込装置1が実行する電子透かし埋込処理の手順を示すフローチャートである。電子透かし埋込装置1は、フレーム処理部101において、外部から入力された音響信号原信号を各フレームに分割するS101。次に、電子透かし埋込装置1は、フィルタ選択部103において、電子透かしデータのビット値に応じて適用する蝸牛遅延フィルタの選択を行う。具体的には、外部から入力され、2進数表現のデータに変換された電子透かしデータのビット値が0及び1の何れであるかを判定しS102、その判定結果に応じて第1蝸牛遅延フィルタ102a及び第2蝸牛遅延フィルタ102bの何れかを選択する。なお、電子透かしデータとしては、例えば著作権者名等の著作権情報またはシリアルナンバー等が挙げられる。ステップS102において電子透かしデータのビット値が0であると判定した場合S102で0、電子透かし埋込装置1は、第1蝸牛遅延フィルタ102aを用いて、音響信号原信号に対して位相変調を施すS103。他方、電子透かしデータのビット値が1であると判定した場合S102で1、電子透かし埋込装置1は、第2蝸牛遅延フィルタ102bを用いて、音響信号原信号に対して位相変調を施すS104。これらのステップS103及びS104により、電子透かしデータが音響信号により埋め込まれることになる。次に、電子透かし埋込装置1は、当該フレームに埋め込む電子透かしデータのすべてのビットが処理されたか否かを判定するS105。ここでまだ処理されていないビットがあると判定した場合S105でNO、電子透かし埋込装置1は、ステップS102へ戻り、それ以降の処理を繰り返す。他方、すべてのビットが処理されたと判定した場合S105でYES、電子透かし埋込装置1は、ステップS103及びS104により電子透かしデータの各ビットが埋め込まれた音響信号を接合することにより、透かし入り音響信号を生成するS106。上記の電子透かし埋込処理をすべてのフレームについて行い、それらを接続することにより、透かし入り音響信号ynが生成される。なお、フレームの接続箇所に不連続点が生じることによりスペクトル拡散の原因でもある知覚不可能性に影響が出ることを防止するために、接続部前のフレームの後ろ数点1ms程度をスプラインSpline補間で滑らかにすることが望ましい。[電子透かし検出処理]次に、上記のようにして電子透かしデータが埋め込まれた透かし入り音響信号から、当該電子透かしデータを検出する電子透かし検出処理について説明する。本実施の形態では、上述したように、原信号を参照しないブラインド検出を行う。なお、電子透かし検出装置2は、電子透かし埋込装置1により電子透かしデータが埋め込まれた際のビットレートを示す情報を記憶しており、当該情報に基づいて下記のセグメントの設定を行うものとする。図9は、本発明の実施の形態における電子透かし検出装置2が実行する電子透かし検出処理の手順を示すフローチャートである。電子透かし検出装置2は、フレーム処理部201において、外部から入力された透かし入り音響信号を各フレームに分割するS201。次に、電子透かし検出装置2は、処理対象のセグメントを設定しS202、第1チャープz変換部202aにおいて、当該セグメントの音響信号に対してチャープz変換を行うS203。さらに、第2チャープz変換部202bにおいて、同じ音響信号に対してチャープz変換を行うS204。次に、電子透かし検出装置2は、ステップS203及びS204により得られた2つの周波数スペクトルのうちの何れが、最も低い周波数でのスペクトルの値が急激に減少しているか否かを判定し、その判定結果に基づき、当該音響信号に対して位相変調を施した蝸牛遅延フィルタの零点を推定するS205。本実施の形態の場合、上記のようにスペクトルの値が急激に減少しているのが第1チャープz変換部202aにより得られた周波数スペクトルである場合は当該零点が1/b0であると推定され、同じく第2チャープz変換部202bにより得られた周波数スペクトルである場合は当該零点が1/b1であると推定される。次に、電子透かし検出装置2は、ビット値検出部203において、ステップS205により推定された蝸牛遅延フィルタの零点が1/b0及び1/b1の何れであるかを判定すしS206、1/b0と判定した場合S206で1/b0はビット値0を検出するS207。他方、1/b1と判定した場合S206で1/b1はビット値1を検出するS208。その後、電子透かし検出装置2は、処理対象のフレームのすべてのセグメントについて処理を行ったか否かを判定するS209。ここで、まだ処理を行っていないセグメントがあると判定した場合S209でNO、電子透かし検出装置2は、ステップS202へ戻り、それ以降の処理を繰り返す。他方、すべてのセグメントについて処理を行ったと判定した場合S209でYES、電子透かし検出装置2は、ステップS207及びS208においてビット値検出部203により検出したビット値を接合することにより、電子透かしデータを復元するS210。以上のようにして、蝸牛遅延フィルタを用いて音響信号に埋め込まれた電子透かしデータをブラインド検出することができる。[他の手法との比較評価]次に、上述した本実施の形態の電子透かし埋込処理により埋め込まれた電子透かしデータの知覚不可能性と、同じく電子透かし検出処理によるビット検出の正確性とについて、他の手法と比較評価する。本発明者等は、RWC音楽データベース後藤、橋口、西村、岡、RWC研究用音楽データベース音楽ジャンルデータベースと楽器音データベース情処学研究報告、2002MUS45419262002の全102曲を評価用の原信号サンプリング周波数44.1kHz、16ビット量子化として利用して、客観評価実験を行った。ここでは、冒頭10秒間を元曲として、8文字の情報AISLabを透かし情報として各原信号に埋め込んだ。また、Nbit=4bpsをベースとし、12条件のNbitNbit4,8,16,32,64,128,256,512,1024,2048,4096,819bpsで、電子透かしデータを原信号の両チャンネルに埋め込み、その特性評価を行った。音質評価に関しては、YLinandWHAbdullaPerceptualevaluationofaudiowatermarkingusingobjectivequalitymeasureProcICASSP2008174517482008に基づき、オーディオ信号に対する知覚評価尺度PEAQPKabalAnexaminationandinterpretationofITURBS1387PerceptualevaluationofaudioqualityTSPLabTechnicalReportDeptElectricalComputerEngineeringMcGUniv2002及び対数スペクトル歪尺度LSDを利用した。比較対象の手法としては、代表的な電子音響透かし法であるLSB法、DSS法、ECHO法、PPM法を利用した。なお、これらの手法は、PPM法を除き、何れもブラインド検出法である。また、発明者等により非特許文献6及び7にて提案されているCD法も比較対象とした。以下、この比較対象のCD法をCDNonBlind法とし、本実施の形態の電子透かし検出方法をCDBlind法と表現する。図10は、上記の客観評価実験の結果を示すグラフであり、a乃至cはそれぞれPEAQ、LSD、ビット検出率についての実験結果を示している。なお、図10では、上記102曲についての平均値が示されている。まず、図10aに示す結果について検討する。PEAQのODGObjectiveDifferenceGrade値は0知覚不可能~-4非常に耳障りであるであるため、ここでは-1知覚される可能性があるが耳障りではないを知覚不可能性の閾値と定めた。図10aに示されるように、DSS法が最も悪く、ECHO法もビットレートが8bps以降から先で急激に悪くなっている。また、PPM法は全般的にODGが-2程度となっている。他方、LSB法は、今回の全てのビットレートにおいて良好な結果が得られている。また、CDNonBlind法では、ビットレートが4bpsでは全く問題がないが、128bpsあたりからODG値が減少し、1024bps程度以降で閾値-1を下回っている。これらに対し、本実施の形態のCDBlind法では、64bpsの時点で既に-1.0付近となり、bpsの増加とともに-3.0付近まで低下している。次に、図10bに示す結果について検討する。一般にLSDは1dB内の歪みであれば音質がよいといわれているため、ここでは、LSDの閾値を1dBに設定した。図10bに示すように、LSB法は、ビットレートを変えても埋め込みによる歪みの影響を受けておらず、良好な結果が得られている。他方、DSS法の場合、ビットレートの増加にかかわらず評価閾値の上にあり、音質評価では問題があることが分かる。ECHO法及びPPM法は、いずれも評価閾値内にあり、特に音質に関して問題があるとはいえない。また、CDNonBlind法は、すべてのビットレートで閾値内にあり、256bpsまでは0.5dB以内を維持するという良好な結果となっている。これらに対し、CDBlind法は、ビットレートの増加に対して単調増加しており、Nbit<1024bpsまでは閾値以下-1dBにあるものの、CDNonBlind法と比較すると若干大きな値となっている。しかし、4~64bpsの付近では、CDBlind法でのLSDがCDNonBlind法のものよりも若干小さい値になっている。なお、CDBlind法及びCDNonBlind法のLSDでの差は、図10aに示すPEAQの場合ほど大きくはなっていない。これは、単純なスペクトル歪みでみる場合と比べると、聴覚的な印象に基づく尺度では両者により違いが現れるためであると考えられる。最後に、図10cに示す結果について検討する。ここでは、ビット検出率の閾値を75とした。図10cに示すように、LSB法を除き、いずれの手法ともビットレートの増加に伴いビット検出率の低下が見られる。CDNonBlind法は、Nbit=1024bps程度で閾値を切るが、他の従来手法ではもっと低いビットレートで閾値を切っている。他方、本実施の形態のCDBlind法では、ビット検出率の低下はほとんど見られず、CDNonBlind法と比べても良好な結果が得られている。具体的には、Nbit<512まではほぼ100であり、1024bpsに至って98となっている。なお、上記の客観評価実験では、LSB法が最も良い結果を出しているが、LSB法の場合、埋め込みされた信号が少しでも改変された場合に検出できないため、頑健性に大きな問題があることが非特許文献6及び7等で指摘されている。これに対し、CDNonBlind法の場合、UnokiMImabeppuKHamadaDHaniuAandMiyauchiREmbeddinglimitationswithdigitalaudiowatermarkingmethodbasedoncochleardelaycharacteristicsJInformationHidingandMultimediaSignalProcessing21232011等に示されるように、十分な頑健性を備えている。しかしながら、CDNonBlind法ではブラインド検出ができないという問題があり、本実施の形態のCDBlind法では、この問題を解消しつつ、優れた知覚不可能性及び頑健性を得ることが可能である。[原信号取得処理]従来の多くの電子音響透かし技術では、電子透かしデータを原信号に埋め込んだ後、それを検出することのみが考慮され、検出後にその電子透かしデータを取り除くことまでは検討されていない。そのため、埋め込まれた電子透かしデータを取り除くための工夫はなされず、除去が困難な態様で電子透かしデータの埋め込みを行っている。このことから、従来の多くの技術は、非可逆的な電子音響透かし技術であるといえる。これに対し、本実施の形態では、原信号に対して蝸牛遅延フィルタにより位相変調を施すという比較的単純な処理で電子透かしデータの埋め込みを行っており、その後検出された電子透かしデータを利用することにより、簡易な方法で当該電子透かしデータを取り除いて原信号を取得することができる。このように、本実施の形態では、可逆型の電子音響透かし技術を実現することができる。以下、この原信号を取得するための処理について説明する。図11は、本発明の実施の形態における電子透かし検出装置2が実行する原信号取得処理の手順を示すフローチャートである。なお、以下では、電子透かし検出装置2が、電子透かし埋込装置1が備える第1蝸牛遅延フィルタ102a及び第2蝸牛遅延フィルタ102bの逆フィルタ、すなわち第1蝸牛遅延フィルタ102a及び第2蝸牛遅延フィルタ102bが模擬する蝸牛遅延特性の逆特性を有するフィルタを備えているものとする。電子透かし検出装置2は、フレーム処理部201において、外部から入力された透かし入り音響信号を各フレームに分割するS301。次に、電子透かし検出装置2は、上記の電子透かし検出処理により検出された電子透かしデータを参照しS302、その電子透かしデータのビット値が0及び1の何れであるかを判定するS303。ステップS303において電子透かしデータのビット値が0であると判定した場合S303で0、電子透かし検出装置2は、第1蝸牛遅延フィルタ102aの逆フィルタを用いて、透かし入り音響信号に対して位相変調を施すS304。他方、電子透かしデータのビット値が1であると判定した場合S303で1、電子透かし検出装置2は、第2蝸牛遅延フィルタ102bの逆フィルタを用いて、透かし入り音響信号に対して位相変調を施すS305。次に、電子透かし埋込装置1は、当該フレームに埋め込まれている電子透かしデータのすべてのビットについて処理がなされたか否かを判定するS306。ここでまだ処理がなされていないビットがあると判定した場合S306でNO、電子透かし検出装置2は、ステップS303へ戻り、それ以降の処理を繰り返す。他方、すべてのビットについて処理がなされたと判定した場合S306でYES、電子透かし検出装置2は、ステップS304及びS305により位相変調が施された音響信号を接合することにより、原信号を復元するS307。上記の原信号取得処理をすべてのフレームについて行い、それらを接続することにより、原信号が取得されることになる。なお、電子透かし埋込処理の場合と同様に、フレームの接続箇所に不連続点が生じることにより知覚不可能性に影響が出ることを防止するために、接続部前のフレームの後ろ数点1ms程度をスプライン補間で滑らかにすることが望ましい。[原信号取得処理の評価]上述した原信号取得処理により取得された信号が実際の原信号と一致しているか等の点について確認するために、上記の客観評価実験と同様の実験を行った。以下、この結果について検討する。図12は、CDNonBlind法及びCDBlind法における電子透かし埋込処理により生成された透かし入り音響信号についての上記客観評価実験の結果を示すグラフであり、a乃至cはそれぞれPEAQ、LSD、ビット検出率についての実験結果を示している。なお、図12には、上記102曲についての平均値が示されている。図12において、CDBlind法の結果は、上述したスプライン補間を行った場合BlindSplineありと行っていない場合BlindSplineなしとに分けて示されている。図12を参照すると、スプライン補間を行っている方が、PEAQ、LSD、ビット検出率の何れについても良い結果が出ていることが分かる。ただし、ビット検出率についてはほとんど違いがない。他方、図13は、本実施の形態の原信号取得処理により電子透かしデータを削除する前及び削除した後における上記の客観評価実験の結果を示すグラフであり、a乃至cはそれぞれPEAQ、LSD、SNRSignalNoiseRatioについての実験結果を示している。このSNRにおいて、Sは原信号、Nは原信号と回復信号上記の原信号取得処理により得られた信号との差を意味している。なお、ここでも、上記102曲についての平均値が示されている。図13を参照すると、電子透かしデータを削除する前よりも、削除した後の方が全般的に良好な結果となっていることが分かる。特に、図13cに示すSNRではそのことが顕著となっている。回復信号が原音に近付くほどSNRは高い値になるため、図13cに示す結果は、本実施の形態の原信号取得処理により取得された信号が原信号に近いこと、換言すると透かし入り音響信号から埋め込まれた電子透かしデータを効果的に削除することができたことを表しているといえる。このように、本実施の形態では、蝸牛遅延フィルタの逆フィルタを用いて位相変調を行うという簡易な処理で、透かし入り音響信号から電子透かしデータを除去して原信号を取得することができる。このように、原信号を取得することができるため、その原信号に新たな電子透かしデータを埋め込み、これを流通させるようなことも可能となる。これにより、埋め込み情報例えば、著作権情報、シリアルナンバー等の内容を更新することができる電子音響透かし技術を実現することができる。その他の実施の形態上記の各実施の形態においては、電子透かしデータの埋め込み処理及び検出処理がソフトウェアにより実現されているが、本発明はこれに限定されるわけではない。例えば、これらの処理の全部又は一部が、DSPDigitalSignalProcessor等の専用のハードウェア回路によって実現されてもよい。また、上記の各実施の形態においては、原信号であるモノラル音楽信号に対して電子透かしデータを埋め込んでいるが、本発明はこれに限られるわけではなく、ステレオ音楽信号の両チャンネルに対して電子透かしデータを埋め込むことも可能である。"