(57)【特許請求の範囲】
【請求項１】
周波数成分値を変更することでコンテンツに透かし情報の埋め込み処理を行う電子透かし埋め込み装置において，
透かし情報を埋め込むコンテンツの空間領域の値を入力し，空間領域から周波数領域への変換を行うにあたって，埋め込み処理で必要となる成分位置の周波数成分値のみを算出する手段と，
上記成分位置における透かし情報の埋め込み処理前の周波数成分値と埋め込み処理後の周波数成分値との差分値を算出する手段と，
上記周波数成分の差分値を空間領域に戻したときの空間領域での差分値と元の空間領域の値とを加算する手段とを備える
ことを特徴とする電子透かし埋め込み装置。
【請求項２】
周波数成分値を変更することでコンテンツに透かし情報の埋め込み処理を行う電子透かし埋め込み方法において，
透かし情報を埋め込むコンテンツの空間領域の値を入力し，空間領域から周波数領域への変換を行うにあたって，埋め込み処理で必要となる成分位置の周波数成分値のみを算出する過程と，
上記成分位置における透かし情報の埋め込み処理前の周波数成分値と埋め込み処理後の周波数成分値との差分値を算出する過程と，
上記周波数成分の差分値を空間領域に戻したときの空間領域での差分値と元の空間領域の値とを加算する過程とを有する
ことを特徴とする電子透かし埋め込み方法。
【請求項３】
周波数成分値を変更することでコンテンツに透かし情報を埋め込む処理を行うためのプログラムであって，
透かし情報を埋め込むコンテンツの空間領域の値を入力し，空間領域から周波数領域への変換を行うにあたって，埋め込み処理で必要となる成分位置の周波数成分値のみを算出する処理と，
上記成分位置における透かし情報の埋め込み処理前の周波数成分値と埋め込み処理後の周波数成分値との差分値を算出する処理と，
上記周波数成分の差分値を空間領域に戻したときの空間領域での差分値と元の空間領域の値とを加算する処理とを，
計算機に実行させるための電子透かし埋め込み用プログラム。
【請求項４】
周波数成分値を変更することでコンテンツに透かし情報を埋め込む処理を行うためのプログラムを記録した媒体であって，
透かし情報を埋め込むコンテンツの空間領域の値を入力し，空間領域から周波数領域への変換を行うにあたって，埋め込み処理で必要となる成分位置の周波数成分値のみを算出する処理と，
上記成分位置における透かし情報の埋め込み処理前の周波数成分値と埋め込み処理後の周波数成分値との差分値を算出する処理と，
上記周波数成分の差分値を空間領域に戻したときの空間領域での差分値と元の空間領域の値とを加算する処理とを，
計算機に実行させるためのプログラムを記録した
ことを特徴とする電子透かし埋め込み用プログラムの記録媒体。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は，電子透かし埋め込み装置に関わり，特に，画像の空間領域から周波数領域への変換およびその逆変換の演算量を削減することで，装置の小型化・低消費電力を実現する電子透かし埋め込み装置を提供するものである。
【０００２】
【従来の技術】
電子透かしとは，動画像などのコンテンツを知覚できない程度の微少量だけ変更することで，透かし情報と呼ばれるデータをコンテンツに埋め込む技術である。画像に対する電子透かしの手法は大きく分けて二つに大別される。一つは，空間領域において画素値を直接変更することで透かし情報を埋め込む手法であり，もう一つは，直交変換によって得られる周波数領域において周波数成分値を変更することで透かし情報を埋め込む手法である。
【０００３】
周波数成分値を変更する手法では，埋め込み処理，読み取り処理ともに離散コサイン変換（ＤＣＴ）などの直交変換を行う。さらに埋め込み処理では，逆変換も行う。
【０００４】
このような周波数成分値を変更する手法の一つとして，次の文献１で提案されている手法がある。
文献１：Hiroshi Ogawa, Takao Nakamura, Atsuki Tomioka, and Youichi Takashima. Digital watermarking technique for motion pictures based on quantization. IEICE Trans. Fundamentals, Vol.E83-A, No.1, pp.77-89, January 2000.
上記文献１では，直交変換として離散コサイン変換（ＤＣＴ）が用いられている。この文献１の手法を図４，図５および図６を用いて説明する。
【０００５】
図４は，電子透かし埋め込み処理を説明する図である。まず，動画像の入力画像４０１をＮ×Ｎ画素のブロック（Ｂk ）４０４に分割する。あらかじめ，すべてのブロックにおいて，ブロック（Ｂk ）４０４に埋め込む透かし情報（ｗ）４０３のデータ位置ｉk を決定しておく。同様に，透かし情報（ｗ）４０３の位置ｉk のデータｗikを埋め込む周波数領域（Ｄk ）４０５の中での成分の位置（ｕk ，ｖk ）を決定しておく。
【０００６】
各ブロック（Ｂk ）４０４に対してＮ×Ｎ・ＤＣＴを施す。Ｎ×Ｎ・ＤＣＴによって得られた周波数領域（Ｄk ）４０５の成分位置（ｕk ，ｖk ）の値ｄuk,vk を，透かし情報（ｗ）４０３のデータ位置ｉk のデータｗikの値に応じて，ｄ′uk,vk に変更する。ｄuk,vk をｄ′uk,vk に変更した周波数領域（Ｄ′k ）４０６に対して，Ｎ×Ｎ逆ＤＣＴ（ＩＤＣＴ）をかけて空間領域のブロック（Ｂ′k ）４０７に戻す。
【０００７】
すべてのブロックに対して上記の処理を施すことで，透かし埋め込み後の画像４０２を得る。
【０００８】
図５は読み取り処理を説明する図である。まず，動画像の入力画像５０１をＮ×Ｎ画素のブロック（Ｂk ）５０４に分割する。あらかじめ，すべてのブロックにおいて，ブロック（Ｂk ）５０４から読み取る透かし情報（ｗ）５０３のデータ位置ｉk を決定しておく。同様に，透かし情報（ｗ）５０３の位置ｉk のデータｗikを読み取る周波数領域（Ｄk ）５０５の中での成分の位置（ｕk ，ｖk ）を決定しておく。
【０００９】
各ブロック（Ｂk ）５０４に対してＮ×Ｎ・ＤＣＴを施す。Ｎ×Ｎ・ＤＣＴによって得られた周波数領域（Ｄk ）５０５の成分位置（ｕk ，ｖk ）の値ｄuk,vk を解析し，透かし情報（ｗ）５０３のデータ位置ｉk のデータｗikを読み取る。すべてのブロックに対して上記の処理を施すことで透かし情報（ｗ）５０３を得る。
【００１０】
図６に，従来の電子透かし埋め込み／読み取り装置の構成図を示す。埋め込み処理における処理の流れを，図６を用いて説明する。
【００１１】
入力画像の各ブロックＢk は，Ｎ×Ｎ・ＤＣＴ器６１に入力される。Ｎ×Ｎ・ＤＣＴ器６１で，入力されたブロックＢk は周波数領域Ｄk に変換され，成分値変更／解析器６２へと出力される。
【００１２】
成分値変更／解析器６２では，周波数領域Ｄk と成分位置（ｕk ，ｖk ），および透かしデータｗの位置ｉk のデータｗikを入力し，周波数成分値ｄuk,vk をｄ′uk,vk に変更する。ｄuk,vk をｄ′uk,vk に変更したＤ′k をＮ×Ｎ・ＩＤＣＴ器６３へと出力する。Ｎ×Ｎ・ＩＤＣＴ器６３では，Ｄ′k を空間領域のブロックＢ′k に変換し，出力する。
【００１３】
次に読み取り処理における処理の流れを，図６を用いて説明する。入力画像の各ブロックＢk は，Ｎ×Ｎ・ＤＣＴ器６１に入力される。Ｎ×Ｎ・ＤＣＴ器６１で，入力されたブロックＢk は周波数領域Ｄk に変換され，成分値変更／解析器６２へと出力される。
【００１４】
成分値変更／解析器６２では，周波数領域Ｄk と成分位置（ｕk ，ｖk ）とを入力し，周波数成分値ｄuk,vk を解析し，透かしデータｗの位置ｉk のデータｗikを読み取る。読み取り処理では，Ｎ×Ｎ・ＩＤＣＴ器６１は使用されない。
【００１５】
【発明が解決しようとする課題】
従来の技術では，電子透かし埋め込み／読み取り処理において，各ブロックＢを周波数領域Ｄに変換するために，Ｎ×Ｎの２次元ＤＣＴを行っていた。また，埋め込み処理において変更を加えた周波数領域Ｄ′を空間領域Ｂ′に戻すために，Ｎ×Ｎの２次元ＩＤＣＴを行っていた。
【００１６】
空間領域のブロックＢ＝｛ｂi,j ｝（０≦ｉ＜Ｎ，０≦ｊ＜Ｎ）を周波数領域Ｄ＝｛ｄu,v ｝（０≦ｕ＜Ｎ，０≦ｖ＜Ｎ）に変換するＮ×Ｎの２次元ＤＣＴは，次の式（１）で表される。
【００１７】
【数１】
0003737723.tif 000002
【００１８】
ここで，０≦ｕ＜Ｎ，０≦ｖ＜Ｎである。これは行列表現を使って表すと，
Ｄ＝ＡＢＡT                                     （２）
と書くことができる。ここでＡ，Ｂ，Ｄは，すべてＮ行Ｎ列の行列であり，ＡT は，Ａの転置行列である。Ａは，ＤＣＴ行列であり，第ｉ行第ｊ列の要素が，
【００１９】
【数２】
0003737723.tif 000003
【００２０】
である行列である。ただし，
【００２１】
【数３】
0003737723.tif 000004
【００２２】
である。
【００２３】
通常，Ｎ×Ｎの２次元ＤＣＴは，Ｎ回のＮ点１次元ＤＣＴを２回直列に行うことで実現する。Ｎ点１次元ＤＣＴを直接的に行うとＮ2 回の乗算が必要となる。Ｎ×Ｎの２次元ＤＣＴでは，約２Ｎ3 回の乗算が必要となる。ＩＤＣＴも同様に２Ｎ3 回の乗算が必要となる。
【００２４】
電子透かしの埋め込み／読み取り処理では，特定の周波数成分値ｄu,v （ｕ，ｖはあらかじめブロックごとに決定）のみを対象とする。それにもかかわらず，従来技術では，埋め込み処理においてＮ×Ｎの２次元ＤＣＴ／ＩＤＣＴを，また読み取り処理においてＮ×Ｎの２次元ＤＣＴを，それぞれ行っていた。そのため，演算量が膨大となり，その結果，回路規模を大きくしたり動作周波数を上げたりすることが必要となって，それによる消費電力の増大を招いていた。
【００２５】
本発明は，上記の点を鑑みなされたもので，透かし処理で必要となる周波数成分値のみを処理の対象とすることで大幅に演算量を削減し，その結果として，小型かつ低消費電力の装置を提供することを目的としている。
【００２６】
【課題を解決するための手段】
本発明の電子透かし埋め込み装置は，上記課題を解決するものであり，電子透かしデータを埋め込む周波数成分値のみＤＣＴ演算を行って計算し，電子透かしデータを埋め込んだ後，電子透かしデータを埋め込んだ周波数成分値と埋め込む前の周波数成分値との差分を求め，この差分をＩＤＣＴ演算して空間領域成分値の差分に変換し，電子透かしデータを埋め込む前の空間領域成分値と加算することにより，電子透かしデータを埋め込んだ空間領域成分値を求めることを特徴とする。
【００２７】
また，本発明に関連する電子透かし読み取り装置は，電子透かしデータが埋め込まれている周波数成分値のみＤＣＴ演算を行って計算し，その周波数成分値から電子透かしデータを読み取ることを特徴とする。
【００２８】
具体的には，以下のとおりである。本発明では，埋め込み／読み取り処理のＤＣＴにおいて，すべてのｕ，ｖに対して空間領域のブロックＢを周波数領域Ｄに変換することはせずに，成分値の変更／解析に必要となる特定の周波数成分値ｄu,v のみを求める。
【００２９】
特定の周波数成分値ｄu,v のみを求めるために，式（２）を，次の式（５）のように変形する。
【００３０】
ｄu,v ＝Ａv ＢＡT u                             （５）
ここで，Ａv はＡの第ｖ行からなる行ベクトルである。同様にＡT u はＡT の第ｕ列からなる列ベクトルである。
【００３１】
また，本発明では，埋め込み処理のＩＤＣＴにおいて，すべてのｕ，ｖに対して周波数領域Ｄ′を空間領域のブロックＢ′に変換することはせずに，埋め込み処理による成分値の差分Δｄu,v ＝ｄ′u,v −ｄu,v のみから空間領域のブロックの差分ΔＢを求める（式（６））。
【００３２】
ΔＢ＝ＡT u Δｄu,v Ａv                         （６）
ΔＢを元のブロックＢに加算することで埋め込み後のブロックＢ′を得る（式（７））。
【００３３】
Ｂ′＝Ｂ＋ΔＢ                                  （７）
次に，空間領域のブロックＢ＝｛ｂi,j ｝（０≦ｉ＜Ｎ，０≦ｊ＜Ｎ）から特定の周波数成分値ｄu,v を求める演算手法を説明する。
【００３４】
式（５）の演算は，
Ａv ＝（ａ(v) 0 ，ａ(v) 1 ，... ，ａ(v) N-1 ）
ＡT u ＝（ａ(u) 0 ，ａ(u) 1 ，... ，ａ(u) N-1 ）T 
とすると，以下のように書くことができる。
【００３５】
【数４】
0003737723.tif 000005
【００３６】
Ａv ＝（ａ(v) 0 ，ａ(v) 1 ，... ，ａ(v) N-1 ）には，式（９）のような対称性がある。
【００３７】
【数５】
0003737723.tif 000006
【００３８】
新たなＮ／２行Ｎ列の行列Ｃ＝｛ｃi,j ｝（０≦ｉ＜Ｎ／２，０≦ｊ＜Ｎ）を以下のように定義する。
【００３９】
【数６】
0003737723.tif 000007
【００４０】
このとき，式（８）は，以下のように変形できる。
【００４１】
【数７】
0003737723.tif 000008
【００４２】
ＡT u ＝（ａ(u) 0 ，ａ(u) 1 ，... ，ａ(u) N-1 ）T にも式（９）のような対称性がある。さらに新たなＮ／２行Ｎ／２列の行列Ｅ＝｛ｅi,j ｝（０≦ｉ＜Ｎ／２，０≦ｊ＜Ｎ／２）を以下のように定義する。
【００４３】
【数８】
0003737723.tif 000009
【００４４】
このとき，式（１１）は，さらに以下のように変形できる。
【００４５】
【数９】
0003737723.tif 000010
【００４６】
このように空間領域Ｂと周波数成分の位置（ｕ，ｖ）が与えられたとき，まずＥを求める。得られたＥに関して，式（１３）の右辺を計算することで，ｄu,v を求めることとする。
【００４７】
次に，周波数成分値の差分Δｄu,v から空間領域の差分ブロックΔＢ＝｛Δｂi,j ｝を求める演算手法を説明する。
【００４８】
式（６）の演算は，以下のように書くことができる。
【００４９】
Δｂi,j ＝ａ(v) j Δｄu,v ａ(u) i               （１４）
ただし，０≦ｉ＜Ｎ，０≦ｊ＜Ｎである。ａ(v) j およびａ(u) i の対称性から，まずΔｂi,j （０≦ｉ＜Ｎ／２，０≦ｊ＜Ｎ／２）を求める。
【００５０】
次に，Δｂi,j の正負をｕ，ｖの偶数／奇数に応じて適切に変えることで，ΔｂN-1-i,j ，Δｂi,N-1-j ，ΔｂN-1-i,N-1-j を得ることとする。
【００５１】
ΔｂN-1-i,j ＝（−１）u Δｂi,j                 （１５）
Δｂi,N-1-j ＝（−１）v Δｂi,j                 （１６）
ΔｂN-1-i,N-1-j ＝（−１）u （−１）v Δｂi,j   （１７）
このように周波数成分値の差分Δｄu,v と周波数成分の位置（ｕ，ｖ）とが与えられたとき，まずΔｂi,j （０≦ｉ＜Ｎ／２，０≦ｊ＜Ｎ／２）を求める。得られたΔｂi,j に関して，式（１５），式（１６），式（１７）を計算することで，ΔＢを求めることとする。
【００５２】
本発明の作用は，以下のとおりである。式（５）を直接的に求めた場合，乗算回数がＮ（Ｎ＋１）回となる。同様に式（６）を直接的に求めた場合，乗算回数がＮ（Ｎ＋１）回となる。
【００５３】
従来の手法では，Ｎ×Ｎの２次元ＤＣＴ／ＩＤＣＴを直接的に行うと，２Ｎ3 回の乗算が必要であるから，本発明の手法によって演算量を削減することができる。
【００５４】
また，式（５）を直接的に求めることはせず，式（１３）と変形してから演算する場合，乗算回数はＮ／２（Ｎ／２＋１）回となる。同様に式（６）を直接的に求めることはせず，Δｂi,j （０≦ｉ＜Ｎ／２，０≦ｊ＜Ｎ／２）を求めてから，式（１５），式（１６），式（１７）を用いてΔＢを求める場合，乗算回数はＮ／２（Ｎ／２＋１）回となる。
【００５５】
このようにＮ×Ｎの２次元ＤＣＴ／ＩＤＣＴをすべて行う従来の手法と比較し，演算量を削減することができる。
【００５６】
以上説明したように，電子透かし埋め込み／読み取り処理において，成分値の変更／解析に必要となる特定の周波数成分値のみを求め，埋め込み処理において，周波数成分の差分値を空間領域に戻し元のブロックと加算することによって，本発明の目的である演算量の制限を実現することができる。
【００５７】
以上の処理は，ハードウェア回路を用いて実現することができ，またソフトウェアプログラムを計算機に実行させることによって実現することもできる。計算機に実行させるためのプログラムは，計算機が読み取り可能な可搬媒体メモリ，半導体メモリ，ハードディスク等の適当な記録媒体に格納することができる。
【００５８】
【発明の実施の形態】
〔電子透かし埋め込み／読み取り装置の構成例〕
図１を用いて，電子透かし埋め込み／読み取り装置による処理を説明する。図１は，本発明の電子透かし埋め込み／読み取り装置の構成例を示す図である。
【００５９】
最初に，埋め込み処理における処理の流れを図１を用いて説明する。入力画像の各ブロックＢは，ＤＣＴ器１に入力される。成分位置ｕ，ｖもＤＣＴ器１に入力される。ＤＣＴ器１では，Ｎ×Ｎ・ＤＣＴをすべて行うことはせず，成分値の変更または解析において必要となる特定の周波数成分ｄu,v のみを求める。求めた周波数成分ｄu,v は，成分値変更／解析器２へと出力される。
【００６０】
成分値変更／解析器２では，周波数成分ｄu,v と透かしデータｗの位置ｉのデータｗi とを入力し，周波数成分値ｄu,v をｄ′u,v に変更する。成分値変更／解析器２は，ＩＤＣＴ器３へ周波数成分値の差分Δｄu,v ＝ｄ′u,v −ｄu,v を出力する。
【００６１】
ＩＤＣＴ器３には，周波数成分値の差分Δｄu,v と成分位置ｕ，ｖとが入力される。ＩＤＣＴ器３は，Δｄu,v を空間領域のブロックΔＢに変換し，差分加算器４へ出力する。
【００６２】
差分加算器４では，元のブロックＢと差分ブロックΔＢとを画素単位で加算して，埋め込み後のブロックＢ′として出力する。
【００６３】
次に，読み取り処理における処理の流れを図１を用いて説明する。入力画像の各ブロックＢは，ＤＣＴ器１に入力される。成分位置ｕ，ｖもＤＣＴ器１に入力される。ＤＣＴ器１でブロックＢから周波数成分ｄu,v が求められ，成分値変更／解析器２へと出力される。
【００６４】
成分値変更／解析器２では，入力した周波数成分値ｄu,v を解析し，透かしデータｗの位置ｉのデータｗi を読み取る。読み取り処理では，ＩＤＣＴ器３，差分加算器４は使用されない。
【００６５】
〔ＤＣＴ器の構成例〕
図２は，本発明で用いるＤＣＴ器の構成例を示す。図２を用いて，周波数成分ｄu,v のみを求める図１のＤＣＴ器１の構成例を説明する。
【００６６】
ＤＣＴ器１は，空間領域の値を入力する入力レジスタ１０１〜１０４と，入力値を成分位置に応じて加減算する加減算器１１０と，加減算器１１０の結果を保持する加減算レジスタ１０５と，加減算レジスタ１０５の出力とレジスタファイル１５０の出力のいずれかを選択する選択器１２０と，ＤＣＴ行列の要素を記憶する行列要素記憶部１３０と，選択器１２０の出力と行列要素記憶部１３０からの出力とを乗算する乗算器１４０と，乗算結果を一時的に蓄えるレジスタファイル１５０と，乗算結果を加算する加算器１６０と，出力レジスタ１０６と，出力レジスタ１０６またはレジスタファイル１５０の出力を選択する選択器１７０とから構成される。
【００６７】
動作は，以下のとおりである。空間領域のブロックＢ＝｛ｂi,j ｝（０≦ｉ＜Ｎ，０≦ｊ＜Ｎ）のうち，ｂi,j ，ｂi,N-1-j ，ｂN-1-i,j ，ｂN-1-i,N-1-j の４つの値が入力され，入力レジスタ１０１，１０２，１０３，１０４に保存される。
【００６８】
加減算器１１０では，上記入力レジスタ１０１，１０２，１０３，１０４から値を入力し，垂直の成分位置ｖおよび水平の成分位置ｕに応じて加減算を行う。演算結果は加減算レジスタ１０５に格納する。
【００６９】
選択器１２０では，入力の加減算結果が加減算レジスタ１０５から出力されている間は加減算レジスタ１０５からの入力を選択し，その後は，後述のレジスタファイル１５０からの出力を選択する。
【００７０】
行列要素記憶部１３０は，ＤＣＴ行列Ａの各要素の値が格納されたＲＯＭ等によって構成されている。行列要素記憶部１３０では，ＤＣＴ行列の要素を以下の順序で出力する。選択器１２０が加減算レジスタ１０５を選択しているときには，ａ(v) 0 ，ａ(v) 1 ，... ，ａ(v) N/2-1 の順序で出力し，選択器１２０がレジスタファイル１５０を選択しているときには，ａ(u) 0 ，ａ(u) 1 ，... ，ａ(u) N/2-1 の順序で出力する。
【００７１】
乗算器１４０では，乗算を行い，乗算結果を加算器１６０に出力する。選択器１７０では，選択器１２０が加減算レジスタ１０５を選択しているときには，レジスタファイル１５０からの出力を選択し，選択器１２０がレジスタファイル１５０を選択しているときには，出力レジスタ１０６からの出力を選択する。
【００７２】
加算器１６０は，乗算器１４０の乗算結果と，選択器１７０の出力とを入力し，加算結果をレジスタファイル１５０と出力レジスタ１０６とに出力する。
【００７３】
レジスタファイル１５０は，Ｎ／２個のレジスタを持ち，途中結果を格納する。Ｎ／２個の値の格納が終了すれば，格納順に選択器１２０へ値を出力する。最終的に加算器１６０による累算結果が，周波数成分値の差分Δｄu,v として出力レジスタ１０６から出力される。
【００７４】
〔ＩＤＣＴ器の構成例〕
図３は，本発明のＩＤＣＴ器の構成例を示す図である。図３を用いて，図１に示すＩＤＣＴ器３の構成例を説明する。
【００７５】
ＩＤＣＴ器３は，周波数成分の差分値Δｄu,v を入力する入力レジスタ３０１と，入力された値とレジスタファイル３４０の出力とを選択する選択器３１０と，ＤＣＴ行列の要素を記憶する行列要素記憶部３２０と，選択器３１０の出力と行列要素記憶部３２０からの出力とを乗算する乗算器３３０と，乗算結果を一時的に蓄えるレジスタファイル３４０と，成分位置に応じて乗算結果の符号を反転する符号反転部３５１〜３５４と，出力値を保持する出力レジスタ３０２〜３０５とから構成される。
【００７６】
動作は，以下のとおりである。周波数成分の差分値Δｄu,v が入力レジスタ３０１に人力される。選択器３１０では，後述するレジスタファイル３４０にＮ／２個のデータが格納されている間は，入力レジスタ３０１を選択し，その後は，レジスタファイル３４０からの出力を選択する。
【００７７】
行列要素記憶部３２０では，ＤＣＴ行列の要素を以下の順序で出力する。選択器３１０が入力レジスタ３０１を選択しているときには，ａ(v) 0 ，ａ(v) 1 ，... ，ａ(v) N/2-1 の順序で出力し，選択器３１０がレジスタファイル３４０を選択しているときには，ａ(u) 0 ，ａ(u) 1 ，... ，ａ(u) N/2-1 の順序で出力する。
【００７８】
乗算器３３０では，乗算を行う。選択器３１０が入力レジスタ３０１を選択しているときには，結果をレジスタファイル３４０に保存し，選択器３１０がレジスタファイル３４０を選択しているときには，結果を符号反転部３５１，３５２，３５３，３５４に出力する。
【００７９】
符号反転部３５１，３５２，３５３，３５４では，垂直の成分位置ｖおよび水平の成分位置ｕに応じて符号を反転またはそのままとし，出力レジスタ３０２，３０３，３０４，３０５から空間領域のブロックの差分値として出力する。
【００８０】
【発明の効果】
以上説明したように，本発明によれば，電子透かし埋め込み／読み取り処理において，画像の空間領域から周波数領域への変換およびその逆変換の演算量を削減することができるため，装置の小型化・低消費電力化が可能となる。
【図面の簡単な説明】
【図１】本発明の電子透かし埋め込み／読み取り装置の構成例を示す図である。
【図２】本実施の形態のＤＣＴ器の構成例を示す図である。
【図３】本実施の形態のＩＤＣＴ器の構成例を示す図である。
【図４】電子透かし埋め込み処理を説明する図である。
【図５】電子透かし読み取り処理を説明する図である。
【図６】従来の電子透かし埋め込み／読み取り装置の構成例を示す図である。
【符号の説明】
１      特定の成分値のみを求めるＤＣＴ器
２      成分値変更／解析器
３      特定の成分値から空間領域に戻すＩＤＣＴ器
４      差分加算器
６１    Ｎ×Ｎ・ＤＣＴ器
６２    成分値変更／解析器
６３    Ｎ×Ｎ・ＩＤＣＴ器
４０１  入力画像
４０２  埋め込み後の出力画像
４０３  埋め込まれる透かし情報
４０４  Ｎ×Ｎのブロック
４０５  ブロックを周波数領域に変換したもの
４０６  周波数領域の特定の成分値を変更したもの
４０７  埋め込み後のブロック
５０１  入力画像
５０３  読み出される透かし情報
５０４  Ｎ×Ｎのブロック
５０５  ブロックを周波数領域に変換したもの
【図１】
0003737723.tif 000011
【図２】
0003737723.tif 000012
【図３】
0003737723.tif 000013
【図４】
0003737723.tif 000014
【図５】
0003737723.tif 000015
【図６】
0003737723.tif 000016