JP2004040237A - 画像処理装置及び方法、並びにコンピュータプログラム及びコンピュータ可読記憶媒体 - Google Patents
画像処理装置及び方法、並びにコンピュータプログラム及びコンピュータ可読記憶媒体 Download PDFInfo
- Publication number
- JP2004040237A JP2004040237A JP2002191133A JP2002191133A JP2004040237A JP 2004040237 A JP2004040237 A JP 2004040237A JP 2002191133 A JP2002191133 A JP 2002191133A JP 2002191133 A JP2002191133 A JP 2002191133A JP 2004040237 A JP2004040237 A JP 2004040237A
- Authority
- JP
- Japan
- Prior art keywords
- noise
- visible
- multiplexed
- image
- pixel
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T1/00—General purpose image data processing
- G06T1/0021—Image watermarking
- G06T1/005—Robust watermarking, e.g. average attack or collusion attack resistant
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2201/00—General purpose image data processing
- G06T2201/005—Image watermarking
- G06T2201/0051—Embedding of the watermark in the spatial domain
Abstract
【課題】原画像にノイズを多重化させ、そのノイズを多重化した分布による目視可能な付加情報を可逆的に埋め込むと共に、ノイズ除去に対する高い耐性を有する画像を生成することを可能にする。
【解決手段】多値画像データにノイズを多重化し、そのノイズを多重化した分布による目視可能な付加情報を埋め込む場合、画素単位に第１のノイズ多重化処理を行うのか、第２のノイズ多重化処理を行うのかを示す情報を、前記付加情報として入力する（Ｓ１０２）。そして、その付加情報に基づき、多値画像データ中の注目画素が、いずれの位置にあるのか否かを判定する（Ｓ１０６）。そして、この判定結果に基づいて、第１のノイズ付加処理（Ｓ１０８）、或いは、第２のノイズ付加処理（Ｓ１１０）のいずれか一方を行う。
【選択図】 図１
【解決手段】多値画像データにノイズを多重化し、そのノイズを多重化した分布による目視可能な付加情報を埋め込む場合、画素単位に第１のノイズ多重化処理を行うのか、第２のノイズ多重化処理を行うのかを示す情報を、前記付加情報として入力する（Ｓ１０２）。そして、その付加情報に基づき、多値画像データ中の注目画素が、いずれの位置にあるのか否かを判定する（Ｓ１０６）。そして、この判定結果に基づいて、第１のノイズ付加処理（Ｓ１０８）、或いは、第２のノイズ付加処理（Ｓ１１０）のいずれか一方を行う。
【選択図】 図１
Description
【０００１】
【発明の属する技術分野】
本発明は、画像の著作権保護を実現するために、静止画像または動画像のデジタル画像データに対して、画質劣化を加える変換、画質劣化を除去する逆変換を行う画像処理装置及び方法、並びに、コンピュータプログラム及びコンピュータ可読記録媒体に関するものである。
【０００２】
【従来の技術】
画像をデジタルデータとして扱うデジタル画像には、従来のアナログ画像と比較して、コンピュータ等によって画質を劣化すること無く簡単にコピーでき、通信回線を通じて電送することができるという特徴がある。しかし、その一方で、このような特徴により、著作権等に絡むデジタル画像が安易に不正コピーされ再配布されることにもなる。
【０００３】
これを防ぐための方法の一つとして電子透かしと呼ばれる手法がある。電子透かしには、大きく分類して著作権情報、ユーザ情報等の透かし情報を目に見えない形で埋め込む不可視型と、画像上に著作権を保有する会社のロゴ等の透かし画像を目に見える形で形成する可視型とがある。
【０００４】
不可視型の電子透かしでは、透かし情報が埋め込まれているということが、埋め込み画像を一見しただけでは認識できない、或いは認識し難い。このため、透かし情報の削除は行われにくいものの、不正コピー、不正配布は可視型に比べ行われ易い。但し、例えデジタル画像が不正にコピー又は配布された場合であっても、そのデジタル画像中には透かし情報が残っているので、透かし情報として埋め込まれたユーザＩＤ等により、不正ユーザを特定することができる。
【０００５】
一方、可視型電子透かしでは、透かし情報は、デジタル画像上に目で見える形で書き込まれているので、そのままでは利用し辛く、不正コピーや不正配布を思いとどまらせる効果がある。可視型の電子透かしの埋め込み方法として、従来、著作権所有者のロゴ等の著作情報を表した画像の画素値を、原画像の画素値と置き換えることにより、著作情報を原画像に埋め込むという手法がよく用いられてきた。この手法の欠点としては、原画像の画素値が失われるため、差分情報なしには原画像を復元できないことが挙げられる。
【０００６】
特開２０００−１８４１７３では、可視型電子透かしを埋め込みたい位置において、画素の全てまたは一部のビットと埋め込みシリアル系列の間で演算処理（暗号化）を行い、可視の透かし画像形状情報を画像に埋め込む方法が提案されている。上記手法に基づき、完全に可逆な可視の電子透かしを実現可能である。
【０００７】
また、特開平８−２５６３２１では、ＪＰＥＧまたはＭＰＥＧ圧縮符号化方式で圧縮符号化された画像のビット列の一部分を抽出し、その抽出部分以外を参照せずに、独立に定義された変換方式によって、抽出されたビット列を直接変換し、書き戻す方法が提案されている。この技術は通称、「半開示」と呼ばれ、開示する画像情報の画質をコントロールした状態でユーザに提供できる。この方法においてビット列を変換する領域を指定することで可逆な可視の電子透かしを実現できる。
【０００８】
【発明が解決しようとする課題】
上記のような可視型電子透かしにおいて、鍵情報を持つ特定の利用者のみが除去を行え、その他の利用者には除去が難しいことが望ましい。
【０００９】
しかし、可視型電子透かしは以下のような攻撃でも除去を試みることが出来る。
【００１０】
一般に自然画像の特徴は近隣の画素値は相関性が高いという特徴をもつ。従って、可視型電子透かしが埋め込まれた領域の周囲を解析し、可視型電子透かしが埋め込まれた領域の画素値を、周囲の画素値を元に計算した値で置き換えることにより、可視型電子透かしの一部を除去することが可能であった。特に平坦で変化の少ない箇所の場合、上記の方法で埋め込まれた可視型電子透かしを簡単に除去することが可能である。つまり、従来までの可視型電子透かしの場合、上記の攻撃に対して適切な対処方法が存在しなかった。
【００１１】
本発明はかかる課題に鑑み成されたものであり、原画像にノイズを多重化させ、そのノイズを多重化した分布による目視可能な付加情報を可逆的に埋め込むと共に、ノイズ除去に対する高い耐性を有する画像を生成することが可能な画像処理装置及び方法、並びにコンピュータプログラム及びコンピュータ可読記憶媒体を提供しようとするものである。
【００１２】
また、他の発明は、上記のようにして可視付加情報が多重化された画像から付加情報を除去することでオリジナル画像、もしくはオリジナル画像に近い画像まで復元することが可能な画像処理装置及び方法、並びにコンピュータプログラム及びコンピュータ可読記憶媒体を提供しようとするものである。
【００１３】
【課題を解決するための手段】
かかる課題を解決するため、例えば本発明の画像処理装置は以下の構成を備える。すなわち、
多値画像データにノイズを多重化し、ノイズを多重化した分布による目視可能な付加情報を埋め込む画像処理装置であって、
第１の可視強度でもってノイズを多重化する第１のノイズ多重化手段と、
前記第１の可視強度とは異なる第２の可視強度でもってノイズを多重化する第２のノイズ多重化手段と、
所定の画素領域単位に、前記第１、第２の可視強度のいずれで多重化するかを示す情報を、前記付加情報として入力する入力手段と、
前記付加情報に基づき、多値画像データ中の注目画素領域が、前記第１、第２の可視強度のいずれで多重化すべきかを判定する判定手段と、
該判定手段の判定結果に基づいて、注目画素領域に対し、前記第１、第２のノイズ多重化手段のいずれか一方を選択してノイズ多重化を行わせる制御手段とを備える。
【００１４】
【本発明の実施の形態】
以下、添付図面に従って本発明に係る実施形態を詳細に説明する。
【００１５】
＜前提の説明＞
本実施形態においては、可視型電子透かしの埋め込みに用いる成分を、カラー画像を構成する輝度成分とする場合について説明する。輝度成分を操作する場合には、色成分の変化が無い為、人の目には明るさに変化が加わったように見える。従って、可視型電子透かしの埋め込みとしては好適な成分である。
【００１６】
しかしながら、可視型電子透かしの埋め込みに用いる成分は輝度成分だけに限らない。Ｒ（赤）、Ｇ（緑）、Ｂ（青）の各成分を、画像の特徴を保存しつつ、可視型電子透かしが人の目に好ましく見えるように、バランスよく操作することも可能である。その他の成分（例えばＣ（シアン）、Ｍ（マゼンダ）、Ｙ（イエロー）等）についても同様である。
【００１７】
また、説明を簡単にする為、入力される画像は８ビットグレースケール画像とするが、Ｒ（赤）、Ｇ（緑）、Ｂ（青）やまたは、Ｙ（輝度）、Ｕ、Ｖ（色差２成分）の色成分から構成される画像においても、本発明の実施の方法で処理可能である。
【００１８】
図２２は、本実施形態に於いて画像に埋め込む可視型電子透かしを表す透かし画像形状情報の一例である。図示の場合、単純な文字列「ＡＢＣ、ＤＥＦ」であるが、透かし画像形状情報としては著作権者のロゴ、画像撮影日時、個人名、企業名、ロゴ、印象的な効果を与える模様等のどのような画像情報になるであろう。また透かし画像形状情報は、画像の関心領域（例えば、医療用画像の患部等）であってもよい。
【００１９】
また本発明においては、透かし画像形状情報は、図２２のように第１、第２の可視型電子透かしが埋めこまれる位置を規定する各画素１ビット（２値）の情報を持つマスク画像（図２２では白（＝１）のアルファベットの領域が可視型電子透かしが埋めこまれる領域を示す）として説明を行う。
【００２０】
＜第１の実施形態＞
以下、本発明の第１の実施形態を図面と共に説明する。
【００２１】
図２５は実施形態における画像処理を行う情報処理装置のブロック構成図である。図示において、１は装置全体の制御を司るＣＰＵであり、２はブートプログラムやＢＩＯＳ等を記憶しているＲＯＭである。３はＣＰＵ１のワークエリアとして使用されるＲＡＭであり、ここにＯＳ、画像処理に係るプログラム等がロードされ、実行されることになる。４はＯＳ、画像処理プログラム、画像データファイル（処理前、処理後を含む）を格納するための外部記憶装置としてのハードディスク装置である。５が画像入力装置であって、例えばイメージスキャナ、デジタルカメラ、画像ファイルを格納している記憶媒体（メモリカードや、フレキシブルディスク、ＣＤＲＯＭ等）、或いは、ネットワーク上から画像をダウンロードするためのインタフェースであっても構わない。６は、画像を表示したり、各種操作を行うためのＧＵＩを提供する表示装置であり、７はキーボード、８は表示画面上の所望とする位置を指定したり、各種メニューを選択するために用いられるポインティングデバイスである。
【００２２】
上記構成における装置に電源が投入され、ＯＳがＲＡＭ３上にロードされ、ユーザの指示、或いは、予め自動起動するように設定しておくことで、本実施形態における画像処理プログラムがＲＡＭ３にロードされ、実行されることになる。
【００２３】
図１は、本発明の第１の実施形態による、画像処理である可逆ノイズ加算の処理を示すフローチャートである。
【００２４】
同図において、まず、ステップＳ１０２は初期設定のステップであり、画素位置と画素値で構成される複数の画素からなる原画像Ｉと、画像中にノイズを加算する領域の形状を示す画素位置からなる透かし画像形状情報Ｍと、２進数で表現される所定のシリアルビット系列を生成するための乱数キーＲ１及びＲ２と、加算するノイズの強度を定義する可視強度値Ｓ１及びＳ２と、入力原画像Ｉと等しい出力画像Ｗを設定する。
【００２５】
なお、原画像Ｉは画像入力装置５より直接入力したもので良いし、ＨＤＤ４に一旦保存した画像ファイルでも構わない。また、画像形状情報Ｍは、ＨＤＤ４に予め記憶されているものを利用するものとするが、ユーザが自由に作成しても良い。乱数キーＲ１、Ｒ２については、乱数発生するための関数（プログラム）を実行すればよいし、可視強度値Ｓ１、Ｓ２はキーボード等から入力しても良いし、予めこれらをファイルとしてＨＤＤに保存してあるものを採用しても構わない。また、出力画像Ｗの出力先は、ＨＤＤ４であるものとする。更に、シリアルビット系列は画像全体で固定としてもよいが、秘匿性を高めるため、乱数キーＲ１、Ｒ２に基づき画像の埋め込み位置に応じて変化するものとする。
【００２６】
また、このステップＳ１０２では、以下の処理に先立ち、原画像Ｉ中の特定の画素を順次選択していくよう設定するが、この選択順は、原画像Ｉの左上隅を開始位置とし、そこから右水平方向の１ラインをスキャンし、ラインエンドになった場合には、次のライン（２ライン目）の左から右方向へスキャンすることを繰り返し行うものとする。これは、ノイズ除去でも同じであるし、後述する各実施形態でも同様である。画像形状情報も同様である。ただし、画像形状情報のサイズが原画像よりも小さい場合には、画像形状情報を反復的に使用するものとする。
【００２７】
次にステップＳ１０４において、入力画像の未処理の画素を選択し、ステップＳ１０６において、その画素位置が、透かし画像形状情報のどの位置に相当するかを判定する。すなわち、画像形状情報中の“１”（図２２における白文字内）の位置にあるのか、“０”（図２２の黒領域）の位置にあるのかを判断する。形状情報内（白画素位置）にあると判断した場合にはステップＳ１０８に進んで第１のノイズ付加処理（ノイズ多重化処理）を行ない、そうでない場合にはステップＳ１１０に進んで第２のノイズ付加処理（ノイズ多重化処理）を行う。
【００２８】
ここで、ステップＳ１０８の第１のノイズ付加処理では、ステップＳ１０２で設定された乱数キーＲ１および可視強度値Ｓ１に基づき、第１のノイズ付加処理を実行する。ステップＳ１１０の第２のノイズ付加処理では、ステップ１０２で設定された乱数キーＲ２および可視強度値Ｓ２に基づき、第２のノイズ付加処理を実行する。
【００２９】
いずれの処理を行った場合にも、処理はステップＳ１１２に進み、入力画像の全画素についてのノイズの付加処理が終了したと判断するまで、ステップＳ１０８、Ｓ１１０を行う。
【００３０】
以上、本実施形態における可逆ノイズ付加処理の概要について簡単に述べたが、次に可逆ノイズが付加された画像からそのノイズを除去する可逆ノイズ除去の処理手順について図２を参照して説明する。なお、ノイズを除去する装置の構成は、図２５に示したものと実質的に同様であるので、その詳述は省略する。
【００３１】
同図において、まず、ステップＳ２０２において初期設定を行う。ここでは、画素位置と画素値で構成される複数の画素からなる原画像Ｉ、埋め込む画像の形状を示す画素位置からなる透かし画像形状情報Ｍ、２進数で表現される所定のシリアルビット系列を生成するための乱数キーＲ１及びＲ２、可視型電子透かしの強度を定義する異なる可視強度値Ｓ１及びＳ２、そして、入力原画像Ｉと等しい出力画像Ｗを設定する。なお、Ｒ１とＲ２は同じ鍵でもよいし、異なる鍵でもよい。また、原画像Ｉは主走査方向の画素数×副走査方向の画素数が定義された上で、主走査方向の１ライン毎の画素データが順次並んだ構造を成しているので、或る画素データを着目するとき、その画素位置は自ずと決まっている。従って、原画像が画素位置をその構成情報の１つとしているのは「画素位置」は便宜的なものである。
【００３２】
次にステップＳ２０４において、入力画像（ノイズが多重化された画像）の未処理の画素を選択する。
【００３３】
次にステップＳ２０６において、そして、その選択画素の位置が透かし画像形状情報のどの位置に相当するかを判定する。すなわち、画像形状情報中の“１”（図２２における白文字内）の位置にあるのか、“０”（図２２の黒領域）の位置にあるのかを判断する。形状情報内（白画素位置）にあると判断した場合にはステップＳ２０８に進んで第１のノイズ除去処理を行ない、そうでない場合にはステップＳ２１０に進んで第２のノイズ除去処理を行う。
【００３４】
そして、ステップＳ２１２では全ての画素が処理されたかを判定し、未処理の画素があるならステップＳ２０４に戻り、前述の処理を全ての画素が終了するまで続行する。
【００３５】
以上、本実施形態における可視型電子透かし除去処理の概要について簡単に述べた。
【００３６】
次に、本実施形態における、第１のノイズ付加処理（ステップＳ１０８）と第２のノイズ付加処理（ステップＳ１１０）の内部処理を詳しく述べる。
【００３７】
なお、第１のノイズ付加処理と第２のノイズ付加処理は、それぞれ同じ原理に基づき行うものとするが、異なる原理に基づきノイズ付加を行ってもよい。要は、結果として付加されるノイズが人の目に異なって知覚されればよい（実施形態では、可視強度値Ｓ１と可視強度値Ｓ２は、画像に異なる強度のノイズを付加するので、異なる知覚強度のノイズが多重化されることになる）。
【００３８】
本実施形態の特徴は、第１と第２のノイズ付加処理において付加されるノイズ量を異ならせ、画像の閲覧者に透かし画像形状を明示するものである。
【００３９】
図３は本実施の形態における、第１と第２のノイズ付加処理の内部処理の手順を示すフローチャートである。
【００４０】
ノイズ付加処理では、まず、演算ビット領域決定処理（ステップＳ３０４）において、可視強度値Ｓ（Ｓ１またはＳ２のいずれか）に従って、演算ビット領域を決定する。
【００４１】
次に演算処理３０６において、選択された演算ビット領域のビット値と乱数キーＲ（Ｒ１またはＲ２のいずれか）から生成されるシリアルビット系列の間で所定の演算処理を実行する。
【００４２】
最後に書き込み処理（ステップＳ３０８）において、演算処理（ステップＳ３０６）の演算結果を選択された画素の演算ビット領域に書き込む。
【００４３】
図４は可視強度値と演算ビット領域の関係を図的に示した図である。図中、４０２、４０４、４０６、４０８は８ビットで表現される各画素値（ここでは各画素の輝度成分の値とする）を２進数のビットで表現した図である。最上位ビット（ＭＳＢ， Ｍｏｓｔ Ｓｉｇｎｉｆｉｃａｎｔ Ｂｉｔ）をＢ７とし、最下位ビット（ＬＳＢ， Ｌｅａｓｔ Ｓｉｇｉｎｉｆｉｃａｎｔ Ｂｉｔ）をＢ０とし、中間のビットをそれぞれ大きいものから、Ｂ６、Ｂ５、Ｂ４、Ｂ３、Ｂ２、Ｂ１としている。
【００４４】
図４では、可視強度値が「１」のとき、４０２に示すようにＭＳＢを除くＢ６以下の７ビットが演算ビット領域に指定されている（ＭＳＢは不変であることを意味する）。同様に、可視強度値が「２」のときにはＢ５以下の６ビットが、可視強度値が「３」のときＢ４以下の５ビットが、可視強度値が「４」のときＢ３以下の４ビットが、それぞれ４０４、４０６、４０８に図示するように演算ビット領域に指定されている。なお、Ｓ＝２〜４において上位の「白」のビットも不変＝演算対象外）を意味する。
【００４５】
図示の関係において、演算ビット領域（サイズ）が大きいほど、その輝度の変化する割合が高くなり、オリジナル画素値に対する変更量が相対的に大きくなることがわかるであろう。
【００４６】
本実施形態における第１、第２のノイズ付加処理は、可視強度値Ｓに応じてこのように演算ビット領域が変化するように設計されている。例えば、第１のノイズ付加処理での可視強度値Ｓ１は例えば図４のＳ＝２を採用し、第２のノイズ付加処理での可視強度Ｓ２を、図４のＳ＝４を採用する。この結果、図２２に示す画像形状情報の白い部分では可視強度が大きなノイズが多重化されることで、目視できる程度の変化を与えることができるようになる。逆に、画像形状情報の黒い部分は可視強度が小さくなり、その変化量は相対的に小さなものとなり、その変化はあまり目立たないようにする。すなわち、原画像に付加する大部分は、画像形状情報における黒い部分であろうから、その部分についてはオリジナル画像の持つ雰囲気を損ねないようにするため、可視強度を小さなものした。
【００４７】
次に図５を用いながら演算処理（ステップＳ３０６）の内部処理の概要を説明する。図示の場合、可視強度Ｓ＝２の場合を示している。
【００４８】
５０２は入力画素値、５０４は乱数キーＲ（Ｒ１、Ｒ２のいずれか）から生成されるシリアルビット系列、５０６は入力画素とシリアルビット系列の対応するビット位置の排他的論理和（ＸＯＲ）、５０８は演算処理後の出力画素を示す。シリアルビット系列５０４はノイズ付加の鍵であると同時に、画素を復号化する際の鍵でもある。なお演算ビット領域以外は演算対象外であるので、対応するシリアルビット系列は必要ない（あっても無視する）。
【００４９】
５０６の太枠内は、入力画素の演算ビット領域と対応するシリアルビット系列（乱数）のビット領域を演算処理結果（ここでは排他的論理和）である。
【００５０】
５０８は５０６の演算処理結果を対応する入力画素の演算ビット領域の値に書き込んだ結果である。
【００５１】
図示の場合、演算結果の画素値と原画素値との差（以下でＢは２進数を意味する）は、
１００１１１０１（Ｂ）−１０１１０１０１（Ｂ）
＝１５７−１８１＝−２４
となり、着目画素値は「−２４」変化したことを意味する。
【００５２】
Ｂ５，Ｂ４，Ｂ３，Ｂ２，Ｂ１，Ｂ０の６ビットが演算ビット領域である場合、演算ビット領域が全て反転すると、最大で ２＾５ ＋ ２＾４ ＋ ２＾３ ＋ ２＾２ ＋ ２＾１ ＋ ２＾０ ＝ ３２ ＋ １６ ＋ ８ ＋４ ＋ ２ ＋ １ ＝ ６３ （ｘ＾ｙはｘのｙ乗を示す）の画素値の変化が実現する。従って、演算ビット領域を指定することで、演算処理で変化する画素値の最大量を決定することが出来る。可視強度値はノイズ付加処理により生成されるノイズの量を表す指標でもあるのは容易に理解できよう。
【００５３】
図６は、図２の可視型電子透かし除去装置における第１のノイズ除去処理（ステップＳ２０８）と第２のノイズ除去処理（ステップＳ２１０）の内部処理手順を示したフローチャートである。
【００５４】
図３とほぼ同様の手順でノイズの除去を行うことが出来るため、簡単に述べる。
【００５５】
ノイズ除去処理では、先ず、演算ビット領域決定処理において、可視強度Ｓ（Ｓ１またはＳ２）に従って、演算ビット領域を決定する（ステップＳ６０４）。
【００５６】
次に演算処理（ステップＳ６０６）において、選択された演算ビット領域のビット値と乱数キーＲ（Ｒ１またはＲ２）から生成されるシリアルビット系列の間で所定の演算処理（復元処理）を実行する。復元が可能であることは、図５の例の場合、処理済み画素値とシリアルビット列とを排他論理和を行って同図における「入力画素値」にまで完全に復元できることからも明らかであろう。
【００５７】
そして、最後に、書き込み処理（ステップＳ６０８）において、演算処理（ステップＳ６０６）の演算結果を選択された画素の演算ビット領域に書き込む。
【００５８】
以上、本実施形態におけるノイズ除去処理について述べた。
【００５９】
本実施形態では、第１の演算処理と第２の演算処理において、異なる可視強度値Ｓ１、Ｓ２を用いて、透かし画像形状内部／外部で演算処理（暗号化）を施す。
【００６０】
画像全体が演算処理（暗号化）されているにも関わらず、第１の領域（透かし画像形状内）と第２の領域（透かし画像形状外）のノイズ量が異なるように設定されている為、透かし画像形状情報を画像の利用者に提示することが可能である。
【００６１】
このとき可視強度値Ｓ１とＳ２の差が大きいほど、画像の利用者に透かし画像形状情報をはっきりと提示することが可能である。
【００６２】
図２２に示す画像形状情報は、簡単に説明するため、単純な文字列ＡＢＣとＤＥＦとを２行にして示したが、実際はオリジナル画像の作成者や撮影者名等が可視的に付加（多重化）されることになる。背景となる部分（図２２における黒色に相当する部分）は面積的に大きいから、それによってオリジナル画像の雰囲気が余り変更されないことが望ましい。従って、背景となる部分に適用される可視強度Ｓ２は、付加情報部分の可視強度Ｓ１よりも、小さいことが望ましい。この意味から、図４の例で示すのであれば、Ｓ１＝Ｓ＝１とし、Ｓ２＝Ｓ＝４とすることが望ましい。
【００６３】
以上の結果、本実施形態では可視型電子透かしの埋め込みによって画像全体が演算処理（暗号化）されている。従って、従来の周囲の画素から演算処理された可視型電子透かし領域の画素を計算して、可視型電子透かしを除去するという攻撃に耐性を持つ。
【００６４】
また本実施形態においては、透かし画像形状内部と外部でそれぞれ所定の演算処理を実行するとして説明したが、透かしの除去を難しくする目的では、透かし画像形状内部とその周囲でそれぞれ所定の演算処理を実行し、それ以外の領域では演算処理を施さないとしてもよい。
【００６５】
以上、実施形態では、
（１）透かし画像形状情報を画像の利用者に提示する、
（２）透かし画像形状情報の除去を困難にする、
という２つの目的を同時に実現することが可能である。
【００６６】
上記例の具体的な利用方法は、例えば、次のようなものである。
【００６７】
インターネット上で画像をサービスするサーバを設置し、図１の処理を施した付加情報を多重化した画像、及び、それぞれの画像を復元するための特有の情報（画像形状情報Ｍ、可視強度値Ｓ１、Ｓ２、乱数キーＲ１、Ｒ２）を記憶管理しておき、ユーザ（クライアント）は所望とする画像を選択し、ダウンロードできるようしておく。このダウンロードした画像には、上記のように付加情報（例えば撮影者や撮影日時等）が可視化して多重化され、尚且つ、その画像の全体の雰囲気を十分に伝えるものができるようになる。ユーザはダウンロードした画像について、オリジナルに復元したいことをサーバに伝えると（例えばブラウザ等で対応するボタンをクリックするとう操作で伝える）、サーバは、オリジナルそのものではなく、その画像に特有の復元に必要な鍵情報（画像形状情報Ｍ、乱数キーＲ１、Ｒ２、可視強度値Ｓ１、Ｓ２（もしくはそれぞれにおける演算ビット領域のサイズを指定する情報））を送信する（これにより、オリジナルそのものが漏洩することを防ぐことができる）。なお、送信する際に、これらの情報を秘密キーでもって暗号化しておけば更に安全なものとすることができる。ユーザ側のＰＣでは、これらの情報を受けとって、図２に示す処理を実行させる処理を行う。
【００６８】
なお、乱数キーＲ１、Ｒ２については、乱数そのものを送信するのでなく、乱数を発生する共通の関数をサーバ及びクライアント側で保有するようにし、その乱数を発生するため初期化パラメータのみを送信するようにする。勿論、全ての画像について同じパラメータとしてしまうと、サーバが保有している全画像を復元できてしまうので、個々の画像毎にこのパラメータを異なるようにしておく。
【００６９】
また、医療用画像（例えばＸ線透視画像）への適用した場合、予め患部部分を示す領域に対応する画像形状情報を用意しておけば良いであろう。この場合、患部については大きなノイズが多重化されるようにし、患部外が小さなノイズが多重化させるようにする。この結果、一個人の情報の漏洩を防ぐことができるようになるであろう。
【００７０】
＜第２の実施の形態＞
上記第１の実施形態では、それぞれ、透かし画像形状の内部と外部において付加されるノイズ量は、画素値に関わらず一定で可視強度値に従って変化した。
【００７１】
本第２の実施形態で述べるノイズ処理を用いても同様に、除去が困難な可逆の可視型電子透かしを実現するものであるが、第１の実施形態に対して異なる点は、ノイズ付加処理及びノイズ除去処理における演算ビット領域が画素値に従って変化する点である。
【００７２】
ところで、人間の視覚特性は、輝度値が小さいほど輝度値の変化に敏感で、輝度値が大きいほど輝度値の変化が目立たたないという特性がある。本第２の実施形態では、かかる特性をも利用するものである。
【００７３】
可視型電子透かしの埋め込みや除去に係る基本の処理手順は、第１の実施形態における図１、図２と同じであるので、その説明については省略し、第２の実施形態における第１、第２のノイズ付加処理を説明する。
【００７４】
図７は本第２の実施形態におけるノイズ付加処理の処理手順を示したフローチャートである。
【００７５】
まず、ステップＳ７０２で判定ビット領域の参照を行い、次にステップＳ７０４で判定ビット領域のビット値と可視強度値の値を元に演算ビット領域を決定する。そして、ステップＳ７０６では、演算ビット領域に演算処理を施し、最後にステップＳ７０８では、演算処理の結果を出力画像の画素の演算ビット領域に書き込む処理を実行する、という手順を踏む。
【００７６】
ここで、本第２の実施形態における判定ビット領域、演算ビット領域について説明する。
【００７７】
図８は８ビットで表現される各画素値（ここでは各画素の輝度成分の値とする）を２進数のビットで表現した図である。最上位ビット（ＭＳＢ， Ｍｏｓｔ ＳｉｇｎｉｆｉｃａｎｔＢｉｔ）をＢ７とし、最下位ビット（ＬＳＢ， Ｌｅａｓｔ Ｓｉｇｉｎｉｆｉｃａｎｔ Ｂｉｔ）をＢ０とし、中間のビットをそれぞれ大きいものから、Ｂ６、Ｂ５、Ｂ４、Ｂ３、Ｂ２、Ｂ１としている。
【００７８】
図示における判定ビット領域、演算ビット領域、非演算ビット領域について簡単に説明すると、判定ビット領域は演算ビット領域を決める為に参照するビット領域であり、演算ビット領域は、ノイズ付加のために所定の演算が行われるビット領域であり、非演算ビット領域は何の操作も行われないビット領域である。
【００７９】
本第２の実施形態では判定ビット領域（判定）は第１または第２のノイズ付加処理において、固定（ＭＳＢ（Ｂ７）及びＢ６の２ビット）であるとする。Ｂ７、Ｂ６とする理由は、原画像の成分の値はビット位置が大きいほど、その値に対して支配的であり、オリジナル画像の持つ雰囲気を保存するためである。
【００８０】
図７における、ステップＳ７０２の判定ビット領域参照処理では、この判定ビット領域のビット値を読み出す。これはこの８ビット表現の対象画素の概略値（おおよその輝度）を参照することに相当する。なお、ステップＳ７０２の時点では、演算ビット領域、非演算ビット領域は未決定である。
【００８１】
図９は、ステップ７０４の演算ビット領域決定処理で用いられる演算ビット領域決定テーブルＴの一例である。同図（ａ）が可視強度Ｓ１用のテーブルであり、同図（ｂ）が可視強度Ｓ２用のテーブルである。
【００８２】
演算ビット領域決定テーブル９０１では、判定ビット領域（Ｂ７、Ｂ６）のビット値は演算ビット領域と対応付けられて定義されている。９０１の演算ビット領域中のＹが記入されているビット位置は、演算ビット領域を示し、無記入のビット位置は、非演算ビット領域を示す。即ち演算ビット領域は、判定ビット領域のビット値の関数で与えられる。図示に示す如く、可視強度Ｓ２は、Ｓ１よりも演算ビット領域が総じて小さい（下位ビット側にシフトしている）。この意味は、可視強度Ｓ２の方が原画像に近い状態でノイズ多重化が行われることを意味する。
【００８３】
例えば、所定の可視強度値において、Ｂ７のビット値が１、 Ｂ６のビット値が０の場合、演算ビット領域決定テーブル９０１を参照した場合、同図（ａ）の場合（可視強度Ｓ１）は演算ビット領域はＢ４、Ｂ３、Ｂ２、Ｂ１、Ｂ０となる。また、同条件においては、同図（ｂ）によれば演算ビット領域はＢ２、Ｂ１、Ｂ０となることからも容易に理解できよう。そして、いずれの場合であっても、輝度が高い程、演算ビット領域は大きく、逆に、輝度が低いほど演算ビット領域は小さくなる。これは、先に説明したように、人間の視覚特性は、輝度値が小さいほど輝度値の変化に敏感で、輝度値が大きいほど輝度値の変化が目立たたないという特性に併せてのものである。
【００８４】
図１０はステップＳ７０４の演算ビット領域決定処理を実行する、演算ビット領域決定部のブロック図である。なお、以下では、可視強度Ｓ１における説明であるが、可視強度Ｓ２でも同様である。テーブル１００４の内容が、図９（ａ）であるか同図（ｂ）であるかの違いだけである。
【００８５】
演算ビット領域決定部１０００には、ステップＳ７０２の判定ビット領域参照処理手段から判定ビット領域のビット値１００２、初期設定で設定された可視強度値Ｓ（この場合はＳ１）、演算ビット領域決定テーブルＴがそれぞれ入力される。
【００８６】
演算ビット領域決定部１０００では、既に説明したように、判定ビット領域のビット値と演算ビット領域決定テーブルＴ、さらには可視強度値Ｓとに基づいて演算ビット領域を決定し、演算ビット領域情報として出力する。
【００８７】
図１１はステップＳ７０６の演算処理を実行する、演算部のブロック図である。
【００８８】
演算部１１００には、処理対象である入力画素値１１０２、乱数キーＲ（ここではＲ１）から生成される位置ごとに異なるシリアルビット系列１１０４、ステップ７０４で決定された演算ビット領域情報１１１０がそれぞれ入力される。演算手段１１００は可逆である所定の演算処理に基づき、演算処理（例えば排他的論理和）を施し、演算結果１１０６を出力する。
【００８９】
演算部１１００の動作の一例を以下、図５を用いて再び説明すると次の通りである。説明が重複するが、５０２は入力画素値、５０４は乱数キーＲから生成されるシリアルビット系列、５０６は入力画素とシリアルビット系列の対応するビット位置の排他的論理和（ＸＯＲ）、５０８は演算処理後の出力画素を示す。
【００９０】
まず入力画素５０２の判定ビット領域を参照し、演算ビット領域を決定する。５０２の太枠で囲まれたビット位置は演算ビット領域を示す。シリアルビット系列５０４は画素を復号化する際の鍵となる。なお演算ビット領域以外に対応するシリアルビット系列は必要ない（もしくは０）。
【００９１】
５０６の太枠内は、入力画素の演算ビット領域と対応するシリアルビット系列のビット領域を演算処理結果（ここでは排他的論理和）である。
【００９２】
５０８は５０６の演算処理結果を対応する入力画素の演算ビット領域の値に書き込んだ結果である。
【００９３】
このように本実施形態は、判定ビット領域から演算ビット領域を決定する処理を導入し、画像情報の一部を保存すると同時に、原画像の画素値の大きさに応じて画素値の最大変化量（Δｍａｘ）を決定することが出来る。
【００９４】
図９（ａ）の９０２は、演算ビット領域決定テーブル９０１に対応する画素値の最大変化量（Δｍａｘ）を示している。
【００９５】
一方、可視強度Ｓ２の場合の構成及び動作は、上記説明からすれば容易に理解できるであろうから、その説明は省略する。
【００９６】
図７のステップＳ７０６の演算処理により演算ビット領域のビットのランダマイズ的な効果が予測されるが、演算ビット領域を決定することで、ノイズを加算する透かし画像形状領域内部は、最大変化量Δｍａｘの範囲内でのランダム的な輝度変化を実現できる。
【００９７】
以上の通り、人間の視覚特性は、輝度値が小さいほど輝度値の変化に敏感で、輝度値が大きいほど輝度値の変化が目立たたないという特性がある。従って、演算ビット領域（最大変化量Δｍａｘ）は上記の人間の視覚特性を考慮して設計することが好ましい。
【００９８】
かかる点、図９（ａ）、（ｂ）によれば、人間の視覚特性を考慮して設計された演算ビット領域決定テーブルとなっているので、画質劣化を抑制するように作用する。
【００９９】
なお、空間上の色の変化が人間の知覚に変化と一致する均等知覚色空間ＣＩＥ１９７６Ｌ＊ｕ＊ｖ＊やＣＩＥ １９７６ Ｌ＊ａ＊ｂ＊（以降Ｌ＊ａ＊ｂ＊色空間と呼ぶ）が１９７６年にＣＩＥから推奨されている。
【０１００】
上記の均等知覚色空間は、上記の最大変化量Δｍａｘを決定する際にも有用である。
【０１０１】
次に本第２の実施形態におけるノイズ除去処理について述べる。
【０１０２】
図１２は第２の本実施形態における、ノイズ除去処理の処理を示すフローチャートである。本第２の実施形態におけるノイズ除去処理は、ほぼノイズ付加処理と同一の処理で実現できる。
【０１０３】
先ず、ステップＳ１２０２では、指定された原画像の画素値中の判定ビット領域のビット情報を読み出す（このとき、注目画素の位置が画像形状情報の白の領域にあるのか、黒の領域にあるのかは既に判定済みである点に注意されたい）。
【０１０４】
次にステップＳ１２０４では、ステップＳ１２０２で読み出した判定ビット領域の値に基づき、ステップＳ１２０６演算処理で処理対象となる演算ビット領域を決定する。
【０１０５】
次にステップＳ１２０６では、ステップＳ１２０４で決定された演算ビット領域のビットと初期設定で入力される乱数キーから生成されるシリアルビット系列の間で逆演算処理を行う。この演算は埋め込み時の演算と対応した逆演算処理（復号処理）である。
【０１０６】
次にステップＳ１２０８では、出力画像Ｅに対応する入力画素のビット領域の値を処理された演算ビット領域の値で書き込む、書き込み処理を実行する。
【０１０７】
以上、第２の本実施形態のノイズ除去処理の動作について述べた。
【０１０８】
本第２の実施形態によれば、第１の実施形態に加えて、原画像の透かし画像形状情報で示される領域において、原画像の特徴を維持しつつ、演算処理によって付加されるノイズ量が画像の階調に従って変化するように、画素値の変更量を設定することが出来る。なお、演算ビット領域を決定する際に演算ビット領域決定テーブルＴを用いたが、数式で表現される演算ビット領域決定関数Ｆを用いることも可能であり、この様な方法も本発明の範疇に含まれる。
【０１０９】
また、演算ビット領域を決定するテーブルを２つ示したが、そのテーブルが３以上あって、そのいずれを選択するかＳ１、Ｓ２で決定するようにしてもよい。
【０１１０】
＜第３の実施の形態＞
次に、本第３の実施形態について説明する。本第３の実施形態でも、付加されたノイズの除去が困難な可逆の可視型電子透かしを実現するものである。
【０１１１】
第１、第２の実施形態と異なる点は、注目画素に対してノイズを多重化する場合、その多重化する領域が可視強度Ｓ１、Ｓ２により決定され、最大でもその領域は７ビットであったが、本第３の実施形態では、演算処理により発生するノイズ量が可視強度値Ｓ１、Ｓ２だけでなく、注目画素の近傍領域の画素値に従って変化するようにすることで、この制限を撤廃できる点に在る。
【０１１２】
ノイズ付加や除去の全体の処理は、第１、第２の実施形態における図１、図２と同様であるものとし、その説明は省略し、第１、第２のノイズ付加処理について説明する。
【０１１３】
図１３は本第３の実施形態におけるノイズ付加処理（第１、第２のノイズ付加処理のいずれか）の処理手順を示したフローチャートである。本第３の実施形態のノイズ付加処理は、付加するノイズ量が可視強度値、及び近傍領域としての隣接画素値に従って変化する点を特徴とする。以下、本第３の実施形態のノイズ付加処理について説明する。
【０１１４】
ノイズ付加処理では、まず、ステップＳ１３０２において、初期設定で設定された近傍領域選択方法ＮＳに基づき、埋め込み対象画素の近傍領域を決定する。次に同じく初期設定で設定された近傍解析方法ＮＡに従い、近傍領域内の画素値を解析し、近傍領域解析値を生成する。ここで、近傍領域解析値は、近傍領域から得られる埋め込み対象画素の予測値である近傍領域画素値と、近傍領域の周波数的特性等を含む近傍領域特性値から成る。
【０１１５】
次にステップＳ１３０４では、ステップＳ１３０２で生成した近傍領域解析値や可視強度値Ｓ（Ｓ１、Ｓ２のいずれか）に基づき、ステップＳ１３０６の演算処理で処理対象となる演算ビット領域を決定する。
【０１１６】
次にステップＳ１３０６では、ステップＳ１３０４で決定された演算ビット領域のビットと初期設定で入力される乱数キーＲから生成されるシリアルビット系列の間で演算処理を行う。このとき演算は可逆の処理である必要がある。演算処理としては、排他的論理和やモジュロ加算、モジュロ乗算等の可逆な演算は全て含まれる。
【０１１７】
次にステップＳ１３０８では、出力画像Ｗに対応する入力画素のビット領域の値を処理された演算ビット領域の値で書き込む、書き込み処理を実行する。
【０１１８】
以上、本実施の形態におけるノイズ付加処理の概要を述べた。
【０１１９】
次に、近傍領域解析処理、近傍領域解析値について詳細に説明する。
【０１２０】
本第３の実施形態においては、埋め込み対象画素に可逆のノイズを埋め込むために隣接画素等から構成される近傍領域の解析に基づいて、演算処理する演算ビット領域を決定する。
【０１２１】
一般に自然画像では隣接画素の間では画素値に相関が大きい。即ち、隣接する画素にはよく似た画素値が配置されている場合が多い。そのため、自然画像では、近隣画素における人の目に知覚可能な変更量は、埋め込み対象画素における人の目に知覚可能な変更量として適切である可能性が高い。
【０１２２】
本第３の実施形態では、埋め込み対象画素と相関性の高い近隣画素を参照し、埋め込み対象画素の最大変化量Δｍａｘをきめ細かく設定することで、原画像のどの階調においてもほぼ同様に知覚されるノイズの加算を実現する。
【０１２３】
本実施形態において、埋め込み対象画素への可逆ノイズの埋め込みの為の演算ビット領域を決める近隣画素で構成される領域を近傍領域と呼ぶこととする。
【０１２４】
この近傍領域を構成する画素は実施形態では単一の画素とするが、複数の画素でもよい。また近傍領域は、埋め込み対象画素と相関が高いと予想される領域であればよく、必ずしも埋め込み対象画素に隣接している必要はないが、ある程度の距離以内にあることが望ましい。
【０１２５】
また、近傍領域内の画素を解析する場合、画素値だけでなく、近傍領域内の画素値の周波数特性や近傍領域内の画素値の分散等の統計的性質を利用してもよい。
【０１２６】
また高周波部分やテクスチャ領域において最大変化量Δｍａｘを大きく設定した演算領域決定テーブルＴ＿Ｎを設計すれば、高周波成分から成る領域やテクスチャ領域でも、人の目に均等に認識され易い可逆のノイズを付加することが出来る。
【０１２７】
次に本実施の形態における近傍領域選択解析処理について詳しく説明する。
【０１２８】
図１８は図１３のステップＳ１３０２の近傍領域選択解析処理を実行する近傍領域選択解析部の内部構成を示すブロック図であり、近傍領域選択部１８０２と近傍領域解析部１８０４で構成される。
【０１２９】
近傍領域選択部１８０２には、画像情報（画素位置、画素値）、埋め込み対象画素の位置情報、近傍領域選択方法ＮＳが入力される。近傍領域選択部１８０２は、前記の入力情報に基づき、近傍領域を決定するが、近傍領域は画像全体で固定でなく、画素位置や所定の鍵情報に従って変化させる形態にしてもよい。
【０１３０】
近傍領域選択部１８０２は後段の近傍領域解析部１８０４に近傍領域情報（画素の位置、画素値等）を出力する。
【０１３１】
後段の近傍領域解析部１８０４には、近傍領域情報（画素の位置、画素値等）と近傍領域解析方法ＮＡが入力され、前記の入力情報に基づき、近傍領域の画素値を解析する。そして、近傍領域解析部１８０４は近傍領域解析値（近傍領域画素値と近傍領域特性値）を出力する。
【０１３２】
次に近傍領域選択解析手段の処理について図を用いて具体的に説明する。
【０１３３】
図１５は８ビットグレースケールの入力画像１５０１と可視型電子透かしが埋め込まれた出力画像１５０２の関係を示している。
【０１３４】
図中の太枠で囲まれた画素（画素１３ａ、１４ａ、１５ａ、１８ａ、１９ａ、２０ａ）は、透かし画像形状内に相当し、可逆ノイズの埋め込み対象画素である。
【０１３５】
今、画素１３ａを着目画素とし、その演算ビット領域は、画素１３ａの近傍領域に基づいて選択され、その近傍領域はまず、近傍領域選択手段１８０２を用いて決定される。
【０１３６】
一例として、近傍領域選択部１８０２は、画素１３ａの近傍領域として左横の画素１２ａ（画素値１１２）を選択するとする（近傍領域として複数の画素領域を選択する場合については、後に説明する）。
【０１３７】
次に近傍領域解析部１８０４に画素１２ａ（画素値１１２）を入力する。図１５では説明を簡単にするため、近傍領域解析部１８０４は、入力画素値１１２をそのまま近傍領域解析値として出力するとする。
【０１３８】
次に、演算ビット領域決定処理では、前段の近傍領域選択解析処理から得られる近傍領域解析値を元に、埋め込み対象画素１３ａの演算ビット領域を決定する。
【０１３９】
次に、演算ビット領域決定処理について詳しく説明する。
【０１４０】
図１６はステップＳ１３０４の演算ビット領域決定処理を実行する、演算ビット領域決定部のブロック図である。
【０１４１】
演算ビット領域決定部１６００には、ステップＳ１３０２の近傍領域選択解析処理から入力される近傍領域解析値１６０２、初期設定で設定された可視強度値Ｓ１６０６、演算ビット領域決定テーブルＴ＿Ｎ１６０４がそれぞれ入力される。
【０１４２】
演算ビット領域決定部１６００では、近傍領域解析値１６０２と演算ビット領域決定テーブルＴ＿Ｎ１６０４、さらには可視強度値Ｓ１６０６とに基づいて、埋め込み対象画素の演算ビット領域を決定し、演算ビット領域情報１６０８として出力する。
【０１４３】
次に、演算ビット領域決定テーブルＴ＿Ｎについて説明する。
【０１４４】
演算ビット領域決定テーブルＴ＿Ｎは、ステップＳ１３０４の演算ビット領域決定処理で演算ビット領域を決定する際に用いられる参照テーブルである。
【０１４５】
図１７に近傍領域解析値（近傍領域内の画素の画素値）と対応付けられた演算ビット領域決定テーブルＴ＿Ｎの一例を示す。同図（ａ）は可視強度Ｓ１用の者であり、同図（ｂ）は可視強度Ｓ２用のものである。なお、図示では２種類示したが、３以上あって、その何れを使用するかを示す情報を可視強度Ｓ１、Ｓ２に割り当てても構わない。
【０１４６】
図１７の１７０１は近傍領域解析値を表し（図１７では簡単のために近傍領域画素値のみを表示している）、１７０２は近傍領域解析値に対応する埋め込み対象画素の演算ビット領域（図１７ではＹが記入されているビット位置が演算ビット領域）を表し、１７０３は演算ビット領域から計算される最大変化量Δｍａｘを表している。
【０１４７】
上記のように演算ビット領域決定テーブルＴ＿Ｎは、近傍領域解析値（図１７では近傍領域から計算される近傍領域画素値）及び可視強度値Ｓ１、Ｓ２が演算ビット領域と対応付けられた値を持っている。演算ビット領域決定部１６００では、上記の演算ビット領域決定テーブルＴ＿Ｎを参照し、入力される近傍領域解析値１６０２及び可視強度値Ｓ（１６０６）に対応する演算ビット領域を読み取り、出力する。
【０１４８】
最後にステップＳ１３０６の演算処理では、既に述べた方法で決定された演算ビット領域において、図５に図示したようなシリアルビット系列と演算ビット領域のビットの間でビット演算を実行する。ステップＳ１３０８の書き込み処理では、ステップＳ１３０６の演算結果を出力画像の対応する演算ビット領域に書き込む。
【０１４９】
次に、図１４を用いて、本実施の形態におけるノイズ除去処理の概要を簡単に述べる。
【０１５０】
ノイズ除去処理は、ステップＳ１４０２では、初期設定で設定された近傍領域選択方法ＮＳに基づき、埋め込み対象画素の近傍領域を決定する。次に同じく初期設定で設定された近傍解析方法ＮＡに従い、近傍領域内の画素値を解析し、近傍領域解析値を生成する。
【０１５１】
次にステップＳ１４０４では、ステップＳ１４０２で生成した近傍領域解析値や可視強度値Ｓ（Ｓ１、Ｓ２のいずれか）に基づき、ステップＳ１４０６の演算処理で処理対象となる演算ビット領域を決定する。
【０１５２】
次にステップＳ１４０６では、ステップＳ１４０４で決定された演算ビット領域のビットと初期設定で入力される乱数キーから生成されるシリアルビット系列の間で逆演算処理を行う。この演算は埋め込み時の演算と対応した逆演算処理（復号処理）である。
【０１５３】
ステップＳ１４０８では、入力画素のビット領域の値を処理された演算ビット領域の値を出力画像Ｅの対応画素に書き込む、書き込み処理を実行する。
【０１５４】
以上、本第３の実施形態のノイズ除去処理の動作について述べた。
【０１５５】
完全に可逆のノイズの付加を実現するためには、近傍領域選択方法に関して若干の条件が必要となる。即ち、逆演算時に、埋め込み時に演算処理された演算ビット領域を正しく知る必要がある。従って、近傍領域選択方法は、逆演算時（可逆ノイズ除去時）に演算ビット領域の決定に用いた近傍領域の画素値を参照できるように近傍領域を選択する必要がある。
【０１５６】
上記の条件を満たすような近傍領域の選択方法は、多数考えられるが、埋め込む処理及び除去処理で処理対象となる画素を順次選択する場合のスキャン方向が左から右へ進む場合には、埋め込み対象画素の左横を近傍領域とすることでる。以下、図１５を用いて埋め込みと原画像を必要としない除去方法の一例を説明する。
【０１５７】
まず埋め込みから簡単に説明する。なお、可視強度Ｓ１、Ｓ２の何れであっても処理内容は同じであるので、ここでは可視強度Ｓ１について説明することとする。
【０１５８】
今、埋め込み対象画素が画素１３ａであるとする。この状況において、その左隣の画素１２ａが近傍領域としているわけであるから、先ず、その左隣の画素１２ａ（画素値１１２）が選択、解析される。その結果、近傍領域解析値の近傍領域画素値が「１１２」として出力される。
【０１５９】
次に、演算ビット領域決定テーブルＴ＿Ｎ（図１７（ａ））の中で「１１２」を持つデータを探し出し、画素１３ａの演算ビット領域を決定する。そして、画素１３ａの演算ビット領域のビット値とシリアルビット系列の演算処理により、画素１３ａ（画素値１１６）に対して画素値１２６が計算され、出力画像１５０２の画素１３ａに画素値１２６が書き込まれる。
【０１６０】
次に画素１４ａに対する処理を行うことになる（注目画素位置を１つ右隣に進める）。入力画像１５０１の左横の画素１３ａ（変更前の画素１１６）を参照し、演算ビット領域決定テーブルＴ＿Ｎを用いて、画素１４ａの演算ビット領域を決定する。そして、演算ビット領域とシリアルビット系列の演算処理により、画素１４ａ（画素値１１４）に対して画素値９８が計算され、出力画像１５０２の画素１４ａに画素９８が書き込まれる。
【０１６１】
以上の処理を１ラインのエンドまで行なうことを各ラインについて行うことで、原画像に対するノイズの埋め込み（多重化）を行うことになる。
【０１６２】
次に上述のようにして付加されたノイズを除去するステップについて説明する。説明を簡単なものとするため、ここでも、可視強度Ｓ１の場合について説明する。
【０１６３】
まず、画素１３ａ（ノイズ加算後の画素値「１２６」の画素）が除去対象画素の場合、近傍領域として出力画像１５０２の左横の画素値１１２の画素１２ａ（既に復元処理が行われた画素である点に注意されたい）が選択、解析される。その結果、近傍領域解析値の近傍領域画素値が１１２として出力される。
【０１６４】
次に、演算ビット領域決定テーブルＴ＿Ｎ（図１７（ａ））を用いて、画素１３ａの演算ビット領域を決定する。そして演算ビット領域のビット値とシリアルビット系列の演算処理により、画素１３ａ（透かし埋め込み後の画素値１２６）に対して復元画素値１１６が計算され、出力画像１５０１（原画像）の画素１３ａに画素値１１６が書き込まれる。
【０１６５】
次に画素１４ａが除去対象画素の場合、入力画像１５０１の左横の画素１３ａではなく、既に復元された出力画像１５０２内の左横の画素１３ａ（復元画素値１１６）が選択、解析される。その結果、近傍領域解析値の近傍領域画素値が１１６として出力される。
【０１６６】
次に、演算ビット領域決定テーブルＴ＿Ｎを用いて、画素１４ａの演算ビット領域を決定する。そして演算ビット領域のビット値とシリアルビット系列の演算処理により、画素１４ａ（画素値９８）に対して画素値１１４が計算され、出力画像１５０１（原画像）の画素１４ａに画素値１１４が書き込まれる。
【０１６７】
このようにして、順次、復元処理済の近隣画素の画素値を選択、解析することで、埋め込み時と同じ演算ビット領域を決定することが出来、付加されたノイズを完全に除去することが出来る。
【０１６８】
簡単な説明の為に上記の説明では近傍領域として埋め込み対象画素の左横の画素を選択したが、直前のラインにおける、主走査方向の同じ位置の画素でも構わない。要するに、参照する近傍領域は復元処理が済んでいるものであれば良い。
【０１６９】
また、以下に説明するように、複数の画素にまたがる領域を近傍領域として選択、解析してもよい。
【０１７０】
例えば、近傍領域選択部は、画素１３ａの近傍領域として、画素７ａ、８ａ、１２ａを選択する。次に近傍領域解析部は、画素７ａ、８ａ、１２ａの画素値から、埋め込み対象画素１３ａの画素値を予測し（例えば平均値）、その予測値を近傍領域画素値としてもよい。
【０１７１】
それ以外にも、近傍領域選択部は、画素１３ａの近傍領域として、左上の画素１ａ、２ａ、６ａ、７ａの４画素を選択する方法なども考えられる。この場合、近傍領域解析部は、近傍領域内の分散値、周波数係数等も合わせて計算し、近傍領域特性値としてもよい。この場合、演算ビット領域決定テーブルＴ＿Ｎには、分散値、周波数係数等の近傍領域特性値にも対応する演算ビット領域が定義されているとする。
【０１７２】
以上は、スキャンする順序が１ラインの左から右、各ラインは上から下に向かうことを前提にしているが、スキャンする方向が異なればこれに限るものではない。
【０１７３】
図１７の演算ビット領域決定テーブルＴ＿Ｎの場合、近傍領域解析値（近傍領域画素値）である入力画素値１１２に対応する演算ビット領域は、Ｂ４， Ｂ３， Ｂ１となる。このとき、最大変化量（Δｍａｘ）を計算すると２＾４＋２＾３＋２＾１＝２６となる。例えば、埋め込み対象画素の画素値が１１２の場合には、１１２のＢ４，Ｂ３，Ｂ１はそれぞれ１，０，０であるから、変更量は＋方向には２＾３＋２＾１＝１０， −方向には２＾４＝１６ の幅を持つ。
【０１７４】
本第３の実施形態のノイズ付加方法は、埋め込み対象画素に近い値を持つ近傍領域画素値や埋め込み対象画素近傍の近傍領域特性値からなる近傍領域解析値に基づき、きめ細かな演算ビット領域の設定が可能であることが分かる。
【０１７５】
次に、本第３の実施形態のノイズ付加方法を用い、入力画像の全画素に可逆ノイズを付加する場合について述べる。
【０１７６】
以上の説明からもわかるように、注目画素に対し、その近傍領域が左画素とした場合、各ラインの先頭画素に着目したとき、近傍領域は存在しない。以下、近傍領域が存在しない場合の対処方法例を幾つか述べる。
【０１７７】
例えば、入力画像の端近くで近傍領域が存在しない画素が埋め込み対象画素の場合は、ノイズの付加を中止する方法をとってもよい。このとき、近傍領域が存在しない画素からノイズを除去する際には、ノイズは付加されていないことが分かっているので、近傍領域が得られるまでノイズの除去は実行しなくてよい。
【０１７８】
また、入力画像の端近くで近傍領域が存在しない画素が埋め込み対象画素の場合は、可視強度値Ｓに従ってのみ決められる固定の演算ビット領域に対して、演算ビット処理を行ってもよい。このとき、近傍領域が存在しない画素からノイズを除去する際には、可視強度値Ｓに従ってのみ決められる演算ビット領域に対して、逆演算処理を施し、ノイズの除去は実行するとよい。
【０１７９】
近傍領域選択方法ＮＳによって決められる近傍領域の大きさが得られるような入力画像の画素値に対しては、既に詳しく述べてきたように、近傍領域解析値に従って、演算ビット領域を決定する。
【０１８０】
人間の視覚特性は、輝度値が小さいほど輝度値の変化に敏感で、輝度値が大きいほど輝度値の変化が目立たたないという特性がある。従って、演算ビット領域（最大変化量Δｍａｘ）は上記の人間の視覚特性を考慮して設計することが好ましい。
【０１８１】
既に第１の実施形態でも述べた均等知覚色空間は、上記の演算ビット領域（最大変化量Δｍａｘ）を決定する際にも有用である。
【０１８２】
本実施の形態によれば、原画像の透かし画像形状情報で示される領域において、原画像の特徴を維持しつつ、演算処理によって付加されるノイズ量が画像の階調に従って変化するように、画素値の変更量を設定することが出来る。
【０１８３】
また本第３の実施形態におけるノイズ付加処理は、既に述べたようにノイズ除去処理を行う際に復元されている近傍領域が必要となる。
【０１８４】
従って、透かし画像形状内部に本第３の実施形態で述べたノイズ付加処理を用った後、透かし画像形状外部に第１または第２の実施の形態で述べたノイズ付加処理を用いてノイズを付加するとすれば、透かし画像形状外部のノイズが除去処理されることなしに、透かし画像形状内部のノイズ除去がより困難になる。従って、透かし画像形状内部のノイズ除去には、透かし画像形状内部のノイズ除去のための鍵情報と透かし画像形状外部のノイズ除去のための鍵情報が二重に必要になるという更なる効果も得られる。
【０１８５】
なお、演算ビット領域を決定する際に演算ビット領域決定テーブルＴ＿Ｎを用いたが、数式で表現される演算ビット領域決定関数Ｆを用いることも可能であり、この様な方法も本発明の範疇に含まれる。
【０１８６】
＜第４の実施の形態＞
以上説明した第１乃至第３の実施形態では、画素単位での埋め込みを想定したが、本第４の実施形態では、ＪＰＥＧやＪＰＥＧ２０００等の圧縮符号化画像に対して、可視型電子透かしを埋め込む場合について説明する。
【０１８７】
ＪＰＥＧやＪＰＥＧ２０００等の圧縮符号化方式では、入力色成分は規定されないが、一般にＲ（赤）、Ｇ（緑）、Ｂ（青）の色成分を、Ｙ（輝度）、Ｃｂ（色差）、Ｃｒ（色差）に変換後、離散コサイン変換や離散ウェーブレット変換を実行する場合が多い。
【０１８８】
従って、ＪＰＥＧやＪＰＥＧ２０００圧縮符号化されたカラー画像の輝度成分を表す周波数変換係数を可視型電子透かしの埋め込み成分として用いることで、特別な工夫なしに輝度値に埋め込むことが可能となる。
【０１８９】
ＪＰＥＧ圧縮符号化方式では、ブロック単位の圧縮符号化が行われている。例えば、ＪＰＥＧ圧縮符号化画像には、最小符号化単位（通常８ｘ８画素）が存在し、その単位毎に圧縮符号化の基本的な処理が行われている。従って、ＪＰＥＧ圧縮符号化画像に可視型電子透かしを埋め込む場合、透かし画像形状情報は、画素単位でなく、最小符号化単位とすると、第１乃至第３の実施の形態の手法を応用し易い。
【０１９０】
より詳しく説明すると、画像を８×８画素ブロック単位に周波数成分データに変換するために先ずＤＣＴ変換を行う。そして、この画素ブロックが、画像形状情報の０の領域にあるのか否かを判断し、「１」の領域にあると判断した場合には、ＤＣＴ変換した結果得られる直流成分ＤＣの値を構成するビットに対し、第１乃至第３の実施形態における第１のノイズ付加処理を行ない、「０」の領域にあると判断した場合には第１乃至第３の実施形態における第２のノイズ付加処理を行うことになる。
【０１９１】
図１９の１９０１は、ＪＰＥＧ圧縮符号化方式において、画像を最小符号化単位毎にブロック化して図示した図である。ＪＰＥＧ圧縮符号化画像の場合、最小符号化単位内（１９０１）で離散コサイン変換（ＤＣＴ）が実行される。１９０２は最小符号化単位ごとに得られるＤＣＴ係数を図示したものである。ＤＣＴ係数は直流（平均値）と交流（周波数成分）からなる。
【０１９２】
従って、最小符号化単位ごとに、ＤＣＴ係数のＤＣ成分（直流）に対して、第１の実施の形態で述べた演算ビット領域の演算処理を行えば、最小符号化単位内（１９０１）の平均値を変化させることが出来、ブロック単位の可視型電子透かし埋め込みが実現できる。
【０１９３】
従って透かし画像形状情報は、埋め込み対象最小符号化単位ブロックを指定する情報とすれば、第１乃至第３の実施形態を適用できる。
【０１９４】
なお、第３の実施形態に適用する場合、参照する近傍領域としては、注目画素ブロックの、例えば左隣の画素ブロックのＤＣ成分を参照すると良いであろう。
【０１９５】
以下、各画素ブロックのＤＣＴ変換後の直流成分にノイズを多重化した後は、公知の符号化処理を行ない、符号化することになる。ただし、ＪＰＥＧ圧縮符号化では、量子化処理により或る程度のデータを破棄することになる。通常、この量子化処理では、量子化ステップ（例えば、２とか４等）を決めているので、ＬＳＢもしくは、ＬＳＢから上位に向かう数ビットは除去されることになるので、ノイズ多重化する際には、可視強度Ｓ１、Ｓ２のいずれであっても、この除去されるビットは演算ビット領域から除外することになるであろう。
【０１９６】
一方、ノイズ除去する場合、逆量子化を行うことになるが、これによって得られる直流成分を、第１乃至第３の実施形態における注目画素と見たてて処理すれば良い。第３の実施形態に適用するのであれば、近傍領域は、例えば、復号済みの左隣の画素ブロックにおける直流成分を近傍領域として採用すれば良いであろう。ただし、第３の実施形態に適用する場合、処理単位が８×８画素単位となるので、その画素ブロックが左端にある場合には近傍領域は存在しないので、ノイズ埋め込み時、及び、ノイズ除去時においては、近傍領域を参照せず、自身の画素ブロックに対して固定の演算ビット領域で処理することになるであろう。
【０１９７】
一方、ＪＰＥＧ２０００圧縮符号化画像では、離散ウェーブレット変換（ＤＷＴ）を用いて、画像の形状情報を保持しつつ、画像を低周波から高周波成分に渡って段階的に帯域ごとに分割し、圧縮符号化する。
【０１９８】
図２０はＪＰＥＧ２０００圧縮符号化方式における離散ウェーブレット変換による帯域分割の様子を図的に示した図である。
【０１９９】
離散ウェーブレット変換では、画像に影響が大きい低周波画像成分はＬＬに集まり、ＬＬは原画像の画像特徴をよく保存している。従って、埋め込みに用いる要素を離散ウェーブレット変換（ＤＷＴ）の低周波成分（ＬＬ）とすると、第１乃至第３の実施形態と比較的近い形で可視型電子透かしを埋め込むことが可能である。
【０２００】
図２０の場合、ＬＬは２回ＤＷＴを行った結果の低周波成分のサブブロックであるが、そのサブブロックを構成している全ての画素に対して、第１乃至第３の実施形態と同様の埋め込みを行えばよい。
【０２０１】
なお、第１乃至第３の実施形態で説明した手法をＪＰＥＧ２０００圧縮符号化画像の離散ウェーブレット変換（ＤＷＴ）係数に適用する場合、ＤＷＴ係数は正負の値を取りうることを考慮し、演算ビット領域決定テーブルを設計する必要がある。
【０２０２】
また第１の実施形態をＪＰＥＧ２０００圧縮符号化方式に適応する場合、ノイズ付加処理対象となるＤＷＴ係数をツリー構造内から、可視強度値に応じて、選択してもよい。またサブバンド毎に量子化ステップが設定されている場合、量子化後のＤＷＴ係数を処理対象とすることで可逆のノイズ付加が可能となる。
【０２０３】
また、以下のようにしても良いであろう。
【０２０４】
図２１の２１０１は、図２０から所定の位置のＬＬとツリー構造の関係を満たすＤＷＴ係数を各サブバンドから抜き出したＤＷＴ係数列である。図中、２１０２は、可視強度値とノイズ付加処理対象となるＤＷＴ係数の関係を示す表である。２１０２では、可視強度値の強さに従って、ノイズ付加処理するＤＷＴ係数も変化している。ＬＬに対するノイズ付加処理は、第２の実施の形態のように判定ビット領域を設けて、演算ビット領域を決定した上でノイズ付加処理を施してもよい。その他のサブバンドのＤＷＴ係数は、演算ビット領域を所定の固定値とし、ノイズ付加処理を施してもよい。
【０２０５】
また、ＪＰＥＧ２０００圧縮符号化画像の場合、ＲＯＩ（Ｒｅｇｉｏｎ ｏｆＩｎｔｅｒｅｓｔ：関心領域）のために、画像サイズと同じサイズを持つ１ビットのビットプレーンが用意される（ＪＰＥＧ２０００の基本符号化システムでは、ＲＯＩ領域だけシフトアップして符号化される）。
【０２０６】
従って、関心領域が特に存在せず、透かし画像形状情報を透かし画像として画像閲覧者に示したい場合、透かし画像形状情報をＲＯＩに設定してもよい。
【０２０７】
例えば、著作権情報を示す可視のロゴ情報をＲＯＩに記述しておけば、コンテンツ配信で画像情報を伝送する際に、ロゴ情報から先に視聴者に示すことが出来、利用者にコンテンツの著作権者を明示する効果も実現できる。
【０２０８】
さらに、透かし画像形状情報は既にＲＯＩ情報として、画像と共に符号化されているため、可視型電子透かしを除去する際に必要な鍵情報が僅かですむというメリットがある。
【０２０９】
また、可視型電子透かしの除去に必要な透かし画像形状情報は、画像ファイルのヘッダ等の所定位置に添付することも可能である。このとき、可視型電子透かしが埋め込まれた画像から原画像に戻すときに、画像ファイル以外に必要な鍵情報だけですみ、配送する情報量を削減できる効果もある。
【０２１０】
また、可視型電子透かしの除去に必要な鍵（および透かし画像形状情報）は比較的情報量が少ない為、画像ファイルのヘッダ等の所定位置に添付することも可能である。このとき、特定の利用者のみ可視型電子透かしの除去を可能にするため、上記の鍵（および透かし画像形状情報）を所定の暗号方式（例えば、公開鍵暗号方式）で暗号化し、画像ファイルのヘッダ等の所定位置に添付してもよい。第１乃至第３の実施の形態では、暗号方式としては、非排他的論理和（ＸＯＲ演算）の場合のみ説明したが、所定の処理単位（例えば６４ビット）となるように複数の演算ビット領域を集め、ＤＥＳ等の共通鍵暗号方式や、公開鍵暗号方式を用いることも可能である。
【０２１１】
また、ＤＥＳ等の所定の処理単位毎に処理を行うブロック暗号方式に属する共通鍵暗号方式を用いる場合には、集めたビット数が所定の処理単位に満たない場合には、足りないビット数だけ０や１をパディングし所定単位を満たし、暗号化を行うとよい。このとき元の画素位置に格納できないビットは、ヘッダ等のファイルの所定位置に付加するとよい。
【０２１２】
あるいは、１〜数ビット単位で処理できる（共通鍵暗号方式に属する）ストリーム暗号方式に属する暗号方式を用いてもよい。
【０２１３】
このとき第１〜第３の実施の形態において、初期設定で入力されるのは乱数キーではなく、共通鍵暗号方式では共通鍵、公開鍵暗号方式では埋め込み時には公開鍵、抽出時に秘密鍵となる。
【０２１４】
本実施の形態では、暗号化方式としてＤＥＳを例に挙げたが、他の共通鍵暗号方式、例えば、ＡＥＳ、ＦＥＡＬ、ＩＤＥＡ、ＲＣ２、ＲＣ４、ＲＣ５、ＭＩＳＴＹ、シーザー型暗号、ビジネル暗号・ビューフォート暗号、プレイフェア暗号、ヒル暗号、バーナム暗号等を用いることもできる。
【０２１５】
また本実施の形態では、静止画像についてのみ説明を行ったが、動画像に対しても同様の原理が適用できることは容易に想像できる。例えば、ＭＰＥＧ圧縮符号化方式の場合、中間フレームを埋め込みの対象とすることで、比較的容易に可視型電子透かしを埋め込むことが出来る。ＭｏｔｉｏｎＪＰＥＧ２０００の場合は時間フレーム方向にＪＰＥＧ２０００圧縮符号化画像と同様の手法で可視型電子透かしの埋め込みを繰り返すことで埋め込み可能である。従って、動画像に対して可視型電子透かしを適用することも本発明の範疇であるとする。
【０２１６】
＜ノイズ多重化のサンプルの説明＞
最後に、第１の実施形態で述べたノイズ付加処理を用いて可視型電子透かしを埋め込んだサンプルを添付する。
【０２１７】
ただし、各画像は、本来は多値、すなわち、１画素が多ビットで構成される階調画像であるが、特許出願の明細書に添付する図面は２値画像のみであって多値画像が扱えない。従って、以下に説明する各画像はノイズ多重化した後の多値階調画像そのものではなく、それを更に誤差拡散処理でもって２値化した結果で示していることに留意されたい。
【０２１８】
図２３は横６４０画素×縦４８０画素の８ビットグレースケール画像（原画像）である。図２４は図２３の原画像に対し、図２２で示される透かし画像形状を埋め込んだ画像である。
【０２１９】
図２４では、透かし画像形状内部はＭＳＢ以外のビット位置が演算ビット領域であり、透かし画像形状外部はＭＳＢ， Ｂ６， Ｂ５以外のビット位置が演算ビット領域である。それぞれ鍵から生成されるシリアルビット列との間でＸＯＲ演算を実行し、演算結果で置き換え、可視型電子透かしが埋め込まれた画像（図２４）を生成している。
【０２２０】
図２４の透かし画像形状外部では、上位ビットを多く残していることから、階調性を維持しているが、透かし画像形状内部では、ＭＳＢしか残していないため、非常に階調性が乏しい画像となっている。しかし画像全体が演算処理されているにも関わらず、透かし画像形状内外でノイズ量が異なるため、画像の閲覧者は透かし画像形状情報ははっきりと見ることが出来る。
【０２２１】
上記の実施形態からも明らかなように、本実施形態を実現するその殆どは、ソフトウェアでもって実現できるものである。通常、パーソナルコンピュータ等の汎用情報処理装置にコンピュータプログラムを導入する場合、フロッピー（登録商標）ディスクやＣＤＲＯＭ、或いは、半導体メモリカード等のコンピュータ可読記憶媒体を装置にセットし、インストールプログラムを実行するか、或いは、システムにコピーすることになるので、かかるコンピュータ可読記憶媒体も本願発明の範疇に含まれるのは明らかである。
【０２２２】
また、コンピュータ上で稼働しているＯＳ等が処理の一部又は全部を行う場合、あるいは記憶媒体から読み出されたプログラムコードが、コンピュータに挿入された拡張機能ボードやコンピュータに接続された拡張機能ユニットに備わるメモリに書き込まれた後、そのプログラムコードの指示に基づいて、上記拡張機能ボードや拡張機能ユニットに備わるＣＰＵ等が処理の一部又は全部を行う場合にも、各実施の形態と同等の機能を実現できると共に、同等の効果を得ることができ、本発明の目的を達成することができる。
【０２２３】
以上説明したように、本実施形態によれば、画像及び鍵及び透かし画像形状情報を入力し、画像の構成要素の要素値を選択し、透かし画像形状情報内部であれば、第１のノイズ付加処理を構成要素の要素値に対し実行し、透かし画像形状情報外部であれば、第２のノイズ付加処理を構成要素の要素値に対し実行し、画像に可視の電子透かしを埋め込む。上記の埋め込みにより、可視型電子透かしは、可視型電子透かしが埋め込まれた領域の周囲を解析し、周囲の画素値を元に計算した値で可視型電子透かしの一部を除去することを困難にすると同時に、第１のノイズ付加処理と第２のノイズ付加処理で付加するノイズを視覚的に異ならせることにより、画像の閲覧者に可視型電子透かしを明示することが出来る。
【０２２４】
【発明の効果】
以上説明したように本発明によれば、原画像にノイズを多重化させ、そのノイズを多重化した分布による目視可能な付加情報を可逆的に埋め込むと共に、ノイズ除去に対する高い耐性を有する画像を生成することが可能になる。
【０２２５】
また、他の発明によれば、可視付加情報が多重化された画像から付加情報を除去することでオリジナル画像、もしくはオリジナル画像に近い画像まで復元することが可能になる。
【図面の簡単な説明】
【図１】実施形態における可視型電子透かし埋め込み処理を示すフローチャートである。
【図２】実施形態における可視型電子透かし除去処理を示すフローチャートである。
【図３】第１の実施形態におけるノイズ付加処理のフローチャートである。
【図４】第１の実施の形態における演算ビット領域を表す図である。
【図５】演算処理の具体例を示す図である。
【図６】第１の実施形態におけるノイズ除去処理のフローチャートである。
【図７】第２の実施形態におけるノイズ付加処理のフローチャートである。
【図８】第２の実施形態における判定ビット領域、演算ビット領域、非演算ビット領域の一例を示す図である。
【図９】第２の実施形態における演算ビット領域決定テーブルの一例を示す図である。
【図１０】第２の実施形態における演算ビット領域決定処理を実行する演算ビット領域決定部の内部構成を示すブロック図である。
【図１１】第２の実施形態における演算処理を実行する演算部の内部構成を示すブロック図である。
【図１２】第２の実施形態におけるノイズ除去処理のフローチャートである。
【図１３】第３の実施形態におけるノイズ付加処理のフローチャートである。
【図１４】第３の実施形態におけるノイズ除去処理のフローチャートである。
【図１５】第３の実施形態における原画像とノイズ多重化後の画像の関係を示す図である。
【図１６】第３の実施形態における演算ビット領域決定処理を実行する演算ビット領域決定部の内部構成を示すブロック図である。
【図１７】第３の実施形態における演算ビット領域決定テーブルの一例を示す図である。
【図１８】近傍領域選択解析処理を実行する近傍領域選択解析部の内部構成を示すブロック図である。
【図１９】ＪＰＥＧ圧縮符号化方式における最小符号化単位を図的に示した図である。
【図２０】ＪＰＥＧ２０００圧縮符号化方式における離散ウェーブレット変換による帯域分割の様子を図的に示した図である。
【図２１】所定のＬＬとツリー関係にある離散ウェーブレット変換係数と可視強度値との関係である。
【図２２】本実施形態における透かし画像形状情報の一例を示す図である。
【図２３】第１の実施形態を方法を適用する為の原画像サンプルである。
【図２４】第１の実施形態の方法を適用し、透かし画像形状を埋め込んだ画像サンプルである。
【図２５】実施形態における装置の具体的なブロック構成図である。
【発明の属する技術分野】
本発明は、画像の著作権保護を実現するために、静止画像または動画像のデジタル画像データに対して、画質劣化を加える変換、画質劣化を除去する逆変換を行う画像処理装置及び方法、並びに、コンピュータプログラム及びコンピュータ可読記録媒体に関するものである。
【０００２】
【従来の技術】
画像をデジタルデータとして扱うデジタル画像には、従来のアナログ画像と比較して、コンピュータ等によって画質を劣化すること無く簡単にコピーでき、通信回線を通じて電送することができるという特徴がある。しかし、その一方で、このような特徴により、著作権等に絡むデジタル画像が安易に不正コピーされ再配布されることにもなる。
【０００３】
これを防ぐための方法の一つとして電子透かしと呼ばれる手法がある。電子透かしには、大きく分類して著作権情報、ユーザ情報等の透かし情報を目に見えない形で埋め込む不可視型と、画像上に著作権を保有する会社のロゴ等の透かし画像を目に見える形で形成する可視型とがある。
【０００４】
不可視型の電子透かしでは、透かし情報が埋め込まれているということが、埋め込み画像を一見しただけでは認識できない、或いは認識し難い。このため、透かし情報の削除は行われにくいものの、不正コピー、不正配布は可視型に比べ行われ易い。但し、例えデジタル画像が不正にコピー又は配布された場合であっても、そのデジタル画像中には透かし情報が残っているので、透かし情報として埋め込まれたユーザＩＤ等により、不正ユーザを特定することができる。
【０００５】
一方、可視型電子透かしでは、透かし情報は、デジタル画像上に目で見える形で書き込まれているので、そのままでは利用し辛く、不正コピーや不正配布を思いとどまらせる効果がある。可視型の電子透かしの埋め込み方法として、従来、著作権所有者のロゴ等の著作情報を表した画像の画素値を、原画像の画素値と置き換えることにより、著作情報を原画像に埋め込むという手法がよく用いられてきた。この手法の欠点としては、原画像の画素値が失われるため、差分情報なしには原画像を復元できないことが挙げられる。
【０００６】
特開２０００−１８４１７３では、可視型電子透かしを埋め込みたい位置において、画素の全てまたは一部のビットと埋め込みシリアル系列の間で演算処理（暗号化）を行い、可視の透かし画像形状情報を画像に埋め込む方法が提案されている。上記手法に基づき、完全に可逆な可視の電子透かしを実現可能である。
【０００７】
また、特開平８−２５６３２１では、ＪＰＥＧまたはＭＰＥＧ圧縮符号化方式で圧縮符号化された画像のビット列の一部分を抽出し、その抽出部分以外を参照せずに、独立に定義された変換方式によって、抽出されたビット列を直接変換し、書き戻す方法が提案されている。この技術は通称、「半開示」と呼ばれ、開示する画像情報の画質をコントロールした状態でユーザに提供できる。この方法においてビット列を変換する領域を指定することで可逆な可視の電子透かしを実現できる。
【０００８】
【発明が解決しようとする課題】
上記のような可視型電子透かしにおいて、鍵情報を持つ特定の利用者のみが除去を行え、その他の利用者には除去が難しいことが望ましい。
【０００９】
しかし、可視型電子透かしは以下のような攻撃でも除去を試みることが出来る。
【００１０】
一般に自然画像の特徴は近隣の画素値は相関性が高いという特徴をもつ。従って、可視型電子透かしが埋め込まれた領域の周囲を解析し、可視型電子透かしが埋め込まれた領域の画素値を、周囲の画素値を元に計算した値で置き換えることにより、可視型電子透かしの一部を除去することが可能であった。特に平坦で変化の少ない箇所の場合、上記の方法で埋め込まれた可視型電子透かしを簡単に除去することが可能である。つまり、従来までの可視型電子透かしの場合、上記の攻撃に対して適切な対処方法が存在しなかった。
【００１１】
本発明はかかる課題に鑑み成されたものであり、原画像にノイズを多重化させ、そのノイズを多重化した分布による目視可能な付加情報を可逆的に埋め込むと共に、ノイズ除去に対する高い耐性を有する画像を生成することが可能な画像処理装置及び方法、並びにコンピュータプログラム及びコンピュータ可読記憶媒体を提供しようとするものである。
【００１２】
また、他の発明は、上記のようにして可視付加情報が多重化された画像から付加情報を除去することでオリジナル画像、もしくはオリジナル画像に近い画像まで復元することが可能な画像処理装置及び方法、並びにコンピュータプログラム及びコンピュータ可読記憶媒体を提供しようとするものである。
【００１３】
【課題を解決するための手段】
かかる課題を解決するため、例えば本発明の画像処理装置は以下の構成を備える。すなわち、
多値画像データにノイズを多重化し、ノイズを多重化した分布による目視可能な付加情報を埋め込む画像処理装置であって、
第１の可視強度でもってノイズを多重化する第１のノイズ多重化手段と、
前記第１の可視強度とは異なる第２の可視強度でもってノイズを多重化する第２のノイズ多重化手段と、
所定の画素領域単位に、前記第１、第２の可視強度のいずれで多重化するかを示す情報を、前記付加情報として入力する入力手段と、
前記付加情報に基づき、多値画像データ中の注目画素領域が、前記第１、第２の可視強度のいずれで多重化すべきかを判定する判定手段と、
該判定手段の判定結果に基づいて、注目画素領域に対し、前記第１、第２のノイズ多重化手段のいずれか一方を選択してノイズ多重化を行わせる制御手段とを備える。
【００１４】
【本発明の実施の形態】
以下、添付図面に従って本発明に係る実施形態を詳細に説明する。
【００１５】
＜前提の説明＞
本実施形態においては、可視型電子透かしの埋め込みに用いる成分を、カラー画像を構成する輝度成分とする場合について説明する。輝度成分を操作する場合には、色成分の変化が無い為、人の目には明るさに変化が加わったように見える。従って、可視型電子透かしの埋め込みとしては好適な成分である。
【００１６】
しかしながら、可視型電子透かしの埋め込みに用いる成分は輝度成分だけに限らない。Ｒ（赤）、Ｇ（緑）、Ｂ（青）の各成分を、画像の特徴を保存しつつ、可視型電子透かしが人の目に好ましく見えるように、バランスよく操作することも可能である。その他の成分（例えばＣ（シアン）、Ｍ（マゼンダ）、Ｙ（イエロー）等）についても同様である。
【００１７】
また、説明を簡単にする為、入力される画像は８ビットグレースケール画像とするが、Ｒ（赤）、Ｇ（緑）、Ｂ（青）やまたは、Ｙ（輝度）、Ｕ、Ｖ（色差２成分）の色成分から構成される画像においても、本発明の実施の方法で処理可能である。
【００１８】
図２２は、本実施形態に於いて画像に埋め込む可視型電子透かしを表す透かし画像形状情報の一例である。図示の場合、単純な文字列「ＡＢＣ、ＤＥＦ」であるが、透かし画像形状情報としては著作権者のロゴ、画像撮影日時、個人名、企業名、ロゴ、印象的な効果を与える模様等のどのような画像情報になるであろう。また透かし画像形状情報は、画像の関心領域（例えば、医療用画像の患部等）であってもよい。
【００１９】
また本発明においては、透かし画像形状情報は、図２２のように第１、第２の可視型電子透かしが埋めこまれる位置を規定する各画素１ビット（２値）の情報を持つマスク画像（図２２では白（＝１）のアルファベットの領域が可視型電子透かしが埋めこまれる領域を示す）として説明を行う。
【００２０】
＜第１の実施形態＞
以下、本発明の第１の実施形態を図面と共に説明する。
【００２１】
図２５は実施形態における画像処理を行う情報処理装置のブロック構成図である。図示において、１は装置全体の制御を司るＣＰＵであり、２はブートプログラムやＢＩＯＳ等を記憶しているＲＯＭである。３はＣＰＵ１のワークエリアとして使用されるＲＡＭであり、ここにＯＳ、画像処理に係るプログラム等がロードされ、実行されることになる。４はＯＳ、画像処理プログラム、画像データファイル（処理前、処理後を含む）を格納するための外部記憶装置としてのハードディスク装置である。５が画像入力装置であって、例えばイメージスキャナ、デジタルカメラ、画像ファイルを格納している記憶媒体（メモリカードや、フレキシブルディスク、ＣＤＲＯＭ等）、或いは、ネットワーク上から画像をダウンロードするためのインタフェースであっても構わない。６は、画像を表示したり、各種操作を行うためのＧＵＩを提供する表示装置であり、７はキーボード、８は表示画面上の所望とする位置を指定したり、各種メニューを選択するために用いられるポインティングデバイスである。
【００２２】
上記構成における装置に電源が投入され、ＯＳがＲＡＭ３上にロードされ、ユーザの指示、或いは、予め自動起動するように設定しておくことで、本実施形態における画像処理プログラムがＲＡＭ３にロードされ、実行されることになる。
【００２３】
図１は、本発明の第１の実施形態による、画像処理である可逆ノイズ加算の処理を示すフローチャートである。
【００２４】
同図において、まず、ステップＳ１０２は初期設定のステップであり、画素位置と画素値で構成される複数の画素からなる原画像Ｉと、画像中にノイズを加算する領域の形状を示す画素位置からなる透かし画像形状情報Ｍと、２進数で表現される所定のシリアルビット系列を生成するための乱数キーＲ１及びＲ２と、加算するノイズの強度を定義する可視強度値Ｓ１及びＳ２と、入力原画像Ｉと等しい出力画像Ｗを設定する。
【００２５】
なお、原画像Ｉは画像入力装置５より直接入力したもので良いし、ＨＤＤ４に一旦保存した画像ファイルでも構わない。また、画像形状情報Ｍは、ＨＤＤ４に予め記憶されているものを利用するものとするが、ユーザが自由に作成しても良い。乱数キーＲ１、Ｒ２については、乱数発生するための関数（プログラム）を実行すればよいし、可視強度値Ｓ１、Ｓ２はキーボード等から入力しても良いし、予めこれらをファイルとしてＨＤＤに保存してあるものを採用しても構わない。また、出力画像Ｗの出力先は、ＨＤＤ４であるものとする。更に、シリアルビット系列は画像全体で固定としてもよいが、秘匿性を高めるため、乱数キーＲ１、Ｒ２に基づき画像の埋め込み位置に応じて変化するものとする。
【００２６】
また、このステップＳ１０２では、以下の処理に先立ち、原画像Ｉ中の特定の画素を順次選択していくよう設定するが、この選択順は、原画像Ｉの左上隅を開始位置とし、そこから右水平方向の１ラインをスキャンし、ラインエンドになった場合には、次のライン（２ライン目）の左から右方向へスキャンすることを繰り返し行うものとする。これは、ノイズ除去でも同じであるし、後述する各実施形態でも同様である。画像形状情報も同様である。ただし、画像形状情報のサイズが原画像よりも小さい場合には、画像形状情報を反復的に使用するものとする。
【００２７】
次にステップＳ１０４において、入力画像の未処理の画素を選択し、ステップＳ１０６において、その画素位置が、透かし画像形状情報のどの位置に相当するかを判定する。すなわち、画像形状情報中の“１”（図２２における白文字内）の位置にあるのか、“０”（図２２の黒領域）の位置にあるのかを判断する。形状情報内（白画素位置）にあると判断した場合にはステップＳ１０８に進んで第１のノイズ付加処理（ノイズ多重化処理）を行ない、そうでない場合にはステップＳ１１０に進んで第２のノイズ付加処理（ノイズ多重化処理）を行う。
【００２８】
ここで、ステップＳ１０８の第１のノイズ付加処理では、ステップＳ１０２で設定された乱数キーＲ１および可視強度値Ｓ１に基づき、第１のノイズ付加処理を実行する。ステップＳ１１０の第２のノイズ付加処理では、ステップ１０２で設定された乱数キーＲ２および可視強度値Ｓ２に基づき、第２のノイズ付加処理を実行する。
【００２９】
いずれの処理を行った場合にも、処理はステップＳ１１２に進み、入力画像の全画素についてのノイズの付加処理が終了したと判断するまで、ステップＳ１０８、Ｓ１１０を行う。
【００３０】
以上、本実施形態における可逆ノイズ付加処理の概要について簡単に述べたが、次に可逆ノイズが付加された画像からそのノイズを除去する可逆ノイズ除去の処理手順について図２を参照して説明する。なお、ノイズを除去する装置の構成は、図２５に示したものと実質的に同様であるので、その詳述は省略する。
【００３１】
同図において、まず、ステップＳ２０２において初期設定を行う。ここでは、画素位置と画素値で構成される複数の画素からなる原画像Ｉ、埋め込む画像の形状を示す画素位置からなる透かし画像形状情報Ｍ、２進数で表現される所定のシリアルビット系列を生成するための乱数キーＲ１及びＲ２、可視型電子透かしの強度を定義する異なる可視強度値Ｓ１及びＳ２、そして、入力原画像Ｉと等しい出力画像Ｗを設定する。なお、Ｒ１とＲ２は同じ鍵でもよいし、異なる鍵でもよい。また、原画像Ｉは主走査方向の画素数×副走査方向の画素数が定義された上で、主走査方向の１ライン毎の画素データが順次並んだ構造を成しているので、或る画素データを着目するとき、その画素位置は自ずと決まっている。従って、原画像が画素位置をその構成情報の１つとしているのは「画素位置」は便宜的なものである。
【００３２】
次にステップＳ２０４において、入力画像（ノイズが多重化された画像）の未処理の画素を選択する。
【００３３】
次にステップＳ２０６において、そして、その選択画素の位置が透かし画像形状情報のどの位置に相当するかを判定する。すなわち、画像形状情報中の“１”（図２２における白文字内）の位置にあるのか、“０”（図２２の黒領域）の位置にあるのかを判断する。形状情報内（白画素位置）にあると判断した場合にはステップＳ２０８に進んで第１のノイズ除去処理を行ない、そうでない場合にはステップＳ２１０に進んで第２のノイズ除去処理を行う。
【００３４】
そして、ステップＳ２１２では全ての画素が処理されたかを判定し、未処理の画素があるならステップＳ２０４に戻り、前述の処理を全ての画素が終了するまで続行する。
【００３５】
以上、本実施形態における可視型電子透かし除去処理の概要について簡単に述べた。
【００３６】
次に、本実施形態における、第１のノイズ付加処理（ステップＳ１０８）と第２のノイズ付加処理（ステップＳ１１０）の内部処理を詳しく述べる。
【００３７】
なお、第１のノイズ付加処理と第２のノイズ付加処理は、それぞれ同じ原理に基づき行うものとするが、異なる原理に基づきノイズ付加を行ってもよい。要は、結果として付加されるノイズが人の目に異なって知覚されればよい（実施形態では、可視強度値Ｓ１と可視強度値Ｓ２は、画像に異なる強度のノイズを付加するので、異なる知覚強度のノイズが多重化されることになる）。
【００３８】
本実施形態の特徴は、第１と第２のノイズ付加処理において付加されるノイズ量を異ならせ、画像の閲覧者に透かし画像形状を明示するものである。
【００３９】
図３は本実施の形態における、第１と第２のノイズ付加処理の内部処理の手順を示すフローチャートである。
【００４０】
ノイズ付加処理では、まず、演算ビット領域決定処理（ステップＳ３０４）において、可視強度値Ｓ（Ｓ１またはＳ２のいずれか）に従って、演算ビット領域を決定する。
【００４１】
次に演算処理３０６において、選択された演算ビット領域のビット値と乱数キーＲ（Ｒ１またはＲ２のいずれか）から生成されるシリアルビット系列の間で所定の演算処理を実行する。
【００４２】
最後に書き込み処理（ステップＳ３０８）において、演算処理（ステップＳ３０６）の演算結果を選択された画素の演算ビット領域に書き込む。
【００４３】
図４は可視強度値と演算ビット領域の関係を図的に示した図である。図中、４０２、４０４、４０６、４０８は８ビットで表現される各画素値（ここでは各画素の輝度成分の値とする）を２進数のビットで表現した図である。最上位ビット（ＭＳＢ， Ｍｏｓｔ Ｓｉｇｎｉｆｉｃａｎｔ Ｂｉｔ）をＢ７とし、最下位ビット（ＬＳＢ， Ｌｅａｓｔ Ｓｉｇｉｎｉｆｉｃａｎｔ Ｂｉｔ）をＢ０とし、中間のビットをそれぞれ大きいものから、Ｂ６、Ｂ５、Ｂ４、Ｂ３、Ｂ２、Ｂ１としている。
【００４４】
図４では、可視強度値が「１」のとき、４０２に示すようにＭＳＢを除くＢ６以下の７ビットが演算ビット領域に指定されている（ＭＳＢは不変であることを意味する）。同様に、可視強度値が「２」のときにはＢ５以下の６ビットが、可視強度値が「３」のときＢ４以下の５ビットが、可視強度値が「４」のときＢ３以下の４ビットが、それぞれ４０４、４０６、４０８に図示するように演算ビット領域に指定されている。なお、Ｓ＝２〜４において上位の「白」のビットも不変＝演算対象外）を意味する。
【００４５】
図示の関係において、演算ビット領域（サイズ）が大きいほど、その輝度の変化する割合が高くなり、オリジナル画素値に対する変更量が相対的に大きくなることがわかるであろう。
【００４６】
本実施形態における第１、第２のノイズ付加処理は、可視強度値Ｓに応じてこのように演算ビット領域が変化するように設計されている。例えば、第１のノイズ付加処理での可視強度値Ｓ１は例えば図４のＳ＝２を採用し、第２のノイズ付加処理での可視強度Ｓ２を、図４のＳ＝４を採用する。この結果、図２２に示す画像形状情報の白い部分では可視強度が大きなノイズが多重化されることで、目視できる程度の変化を与えることができるようになる。逆に、画像形状情報の黒い部分は可視強度が小さくなり、その変化量は相対的に小さなものとなり、その変化はあまり目立たないようにする。すなわち、原画像に付加する大部分は、画像形状情報における黒い部分であろうから、その部分についてはオリジナル画像の持つ雰囲気を損ねないようにするため、可視強度を小さなものした。
【００４７】
次に図５を用いながら演算処理（ステップＳ３０６）の内部処理の概要を説明する。図示の場合、可視強度Ｓ＝２の場合を示している。
【００４８】
５０２は入力画素値、５０４は乱数キーＲ（Ｒ１、Ｒ２のいずれか）から生成されるシリアルビット系列、５０６は入力画素とシリアルビット系列の対応するビット位置の排他的論理和（ＸＯＲ）、５０８は演算処理後の出力画素を示す。シリアルビット系列５０４はノイズ付加の鍵であると同時に、画素を復号化する際の鍵でもある。なお演算ビット領域以外は演算対象外であるので、対応するシリアルビット系列は必要ない（あっても無視する）。
【００４９】
５０６の太枠内は、入力画素の演算ビット領域と対応するシリアルビット系列（乱数）のビット領域を演算処理結果（ここでは排他的論理和）である。
【００５０】
５０８は５０６の演算処理結果を対応する入力画素の演算ビット領域の値に書き込んだ結果である。
【００５１】
図示の場合、演算結果の画素値と原画素値との差（以下でＢは２進数を意味する）は、
１００１１１０１（Ｂ）−１０１１０１０１（Ｂ）
＝１５７−１８１＝−２４
となり、着目画素値は「−２４」変化したことを意味する。
【００５２】
Ｂ５，Ｂ４，Ｂ３，Ｂ２，Ｂ１，Ｂ０の６ビットが演算ビット領域である場合、演算ビット領域が全て反転すると、最大で ２＾５ ＋ ２＾４ ＋ ２＾３ ＋ ２＾２ ＋ ２＾１ ＋ ２＾０ ＝ ３２ ＋ １６ ＋ ８ ＋４ ＋ ２ ＋ １ ＝ ６３ （ｘ＾ｙはｘのｙ乗を示す）の画素値の変化が実現する。従って、演算ビット領域を指定することで、演算処理で変化する画素値の最大量を決定することが出来る。可視強度値はノイズ付加処理により生成されるノイズの量を表す指標でもあるのは容易に理解できよう。
【００５３】
図６は、図２の可視型電子透かし除去装置における第１のノイズ除去処理（ステップＳ２０８）と第２のノイズ除去処理（ステップＳ２１０）の内部処理手順を示したフローチャートである。
【００５４】
図３とほぼ同様の手順でノイズの除去を行うことが出来るため、簡単に述べる。
【００５５】
ノイズ除去処理では、先ず、演算ビット領域決定処理において、可視強度Ｓ（Ｓ１またはＳ２）に従って、演算ビット領域を決定する（ステップＳ６０４）。
【００５６】
次に演算処理（ステップＳ６０６）において、選択された演算ビット領域のビット値と乱数キーＲ（Ｒ１またはＲ２）から生成されるシリアルビット系列の間で所定の演算処理（復元処理）を実行する。復元が可能であることは、図５の例の場合、処理済み画素値とシリアルビット列とを排他論理和を行って同図における「入力画素値」にまで完全に復元できることからも明らかであろう。
【００５７】
そして、最後に、書き込み処理（ステップＳ６０８）において、演算処理（ステップＳ６０６）の演算結果を選択された画素の演算ビット領域に書き込む。
【００５８】
以上、本実施形態におけるノイズ除去処理について述べた。
【００５９】
本実施形態では、第１の演算処理と第２の演算処理において、異なる可視強度値Ｓ１、Ｓ２を用いて、透かし画像形状内部／外部で演算処理（暗号化）を施す。
【００６０】
画像全体が演算処理（暗号化）されているにも関わらず、第１の領域（透かし画像形状内）と第２の領域（透かし画像形状外）のノイズ量が異なるように設定されている為、透かし画像形状情報を画像の利用者に提示することが可能である。
【００６１】
このとき可視強度値Ｓ１とＳ２の差が大きいほど、画像の利用者に透かし画像形状情報をはっきりと提示することが可能である。
【００６２】
図２２に示す画像形状情報は、簡単に説明するため、単純な文字列ＡＢＣとＤＥＦとを２行にして示したが、実際はオリジナル画像の作成者や撮影者名等が可視的に付加（多重化）されることになる。背景となる部分（図２２における黒色に相当する部分）は面積的に大きいから、それによってオリジナル画像の雰囲気が余り変更されないことが望ましい。従って、背景となる部分に適用される可視強度Ｓ２は、付加情報部分の可視強度Ｓ１よりも、小さいことが望ましい。この意味から、図４の例で示すのであれば、Ｓ１＝Ｓ＝１とし、Ｓ２＝Ｓ＝４とすることが望ましい。
【００６３】
以上の結果、本実施形態では可視型電子透かしの埋め込みによって画像全体が演算処理（暗号化）されている。従って、従来の周囲の画素から演算処理された可視型電子透かし領域の画素を計算して、可視型電子透かしを除去するという攻撃に耐性を持つ。
【００６４】
また本実施形態においては、透かし画像形状内部と外部でそれぞれ所定の演算処理を実行するとして説明したが、透かしの除去を難しくする目的では、透かし画像形状内部とその周囲でそれぞれ所定の演算処理を実行し、それ以外の領域では演算処理を施さないとしてもよい。
【００６５】
以上、実施形態では、
（１）透かし画像形状情報を画像の利用者に提示する、
（２）透かし画像形状情報の除去を困難にする、
という２つの目的を同時に実現することが可能である。
【００６６】
上記例の具体的な利用方法は、例えば、次のようなものである。
【００６７】
インターネット上で画像をサービスするサーバを設置し、図１の処理を施した付加情報を多重化した画像、及び、それぞれの画像を復元するための特有の情報（画像形状情報Ｍ、可視強度値Ｓ１、Ｓ２、乱数キーＲ１、Ｒ２）を記憶管理しておき、ユーザ（クライアント）は所望とする画像を選択し、ダウンロードできるようしておく。このダウンロードした画像には、上記のように付加情報（例えば撮影者や撮影日時等）が可視化して多重化され、尚且つ、その画像の全体の雰囲気を十分に伝えるものができるようになる。ユーザはダウンロードした画像について、オリジナルに復元したいことをサーバに伝えると（例えばブラウザ等で対応するボタンをクリックするとう操作で伝える）、サーバは、オリジナルそのものではなく、その画像に特有の復元に必要な鍵情報（画像形状情報Ｍ、乱数キーＲ１、Ｒ２、可視強度値Ｓ１、Ｓ２（もしくはそれぞれにおける演算ビット領域のサイズを指定する情報））を送信する（これにより、オリジナルそのものが漏洩することを防ぐことができる）。なお、送信する際に、これらの情報を秘密キーでもって暗号化しておけば更に安全なものとすることができる。ユーザ側のＰＣでは、これらの情報を受けとって、図２に示す処理を実行させる処理を行う。
【００６８】
なお、乱数キーＲ１、Ｒ２については、乱数そのものを送信するのでなく、乱数を発生する共通の関数をサーバ及びクライアント側で保有するようにし、その乱数を発生するため初期化パラメータのみを送信するようにする。勿論、全ての画像について同じパラメータとしてしまうと、サーバが保有している全画像を復元できてしまうので、個々の画像毎にこのパラメータを異なるようにしておく。
【００６９】
また、医療用画像（例えばＸ線透視画像）への適用した場合、予め患部部分を示す領域に対応する画像形状情報を用意しておけば良いであろう。この場合、患部については大きなノイズが多重化されるようにし、患部外が小さなノイズが多重化させるようにする。この結果、一個人の情報の漏洩を防ぐことができるようになるであろう。
【００７０】
＜第２の実施の形態＞
上記第１の実施形態では、それぞれ、透かし画像形状の内部と外部において付加されるノイズ量は、画素値に関わらず一定で可視強度値に従って変化した。
【００７１】
本第２の実施形態で述べるノイズ処理を用いても同様に、除去が困難な可逆の可視型電子透かしを実現するものであるが、第１の実施形態に対して異なる点は、ノイズ付加処理及びノイズ除去処理における演算ビット領域が画素値に従って変化する点である。
【００７２】
ところで、人間の視覚特性は、輝度値が小さいほど輝度値の変化に敏感で、輝度値が大きいほど輝度値の変化が目立たたないという特性がある。本第２の実施形態では、かかる特性をも利用するものである。
【００７３】
可視型電子透かしの埋め込みや除去に係る基本の処理手順は、第１の実施形態における図１、図２と同じであるので、その説明については省略し、第２の実施形態における第１、第２のノイズ付加処理を説明する。
【００７４】
図７は本第２の実施形態におけるノイズ付加処理の処理手順を示したフローチャートである。
【００７５】
まず、ステップＳ７０２で判定ビット領域の参照を行い、次にステップＳ７０４で判定ビット領域のビット値と可視強度値の値を元に演算ビット領域を決定する。そして、ステップＳ７０６では、演算ビット領域に演算処理を施し、最後にステップＳ７０８では、演算処理の結果を出力画像の画素の演算ビット領域に書き込む処理を実行する、という手順を踏む。
【００７６】
ここで、本第２の実施形態における判定ビット領域、演算ビット領域について説明する。
【００７７】
図８は８ビットで表現される各画素値（ここでは各画素の輝度成分の値とする）を２進数のビットで表現した図である。最上位ビット（ＭＳＢ， Ｍｏｓｔ ＳｉｇｎｉｆｉｃａｎｔＢｉｔ）をＢ７とし、最下位ビット（ＬＳＢ， Ｌｅａｓｔ Ｓｉｇｉｎｉｆｉｃａｎｔ Ｂｉｔ）をＢ０とし、中間のビットをそれぞれ大きいものから、Ｂ６、Ｂ５、Ｂ４、Ｂ３、Ｂ２、Ｂ１としている。
【００７８】
図示における判定ビット領域、演算ビット領域、非演算ビット領域について簡単に説明すると、判定ビット領域は演算ビット領域を決める為に参照するビット領域であり、演算ビット領域は、ノイズ付加のために所定の演算が行われるビット領域であり、非演算ビット領域は何の操作も行われないビット領域である。
【００７９】
本第２の実施形態では判定ビット領域（判定）は第１または第２のノイズ付加処理において、固定（ＭＳＢ（Ｂ７）及びＢ６の２ビット）であるとする。Ｂ７、Ｂ６とする理由は、原画像の成分の値はビット位置が大きいほど、その値に対して支配的であり、オリジナル画像の持つ雰囲気を保存するためである。
【００８０】
図７における、ステップＳ７０２の判定ビット領域参照処理では、この判定ビット領域のビット値を読み出す。これはこの８ビット表現の対象画素の概略値（おおよその輝度）を参照することに相当する。なお、ステップＳ７０２の時点では、演算ビット領域、非演算ビット領域は未決定である。
【００８１】
図９は、ステップ７０４の演算ビット領域決定処理で用いられる演算ビット領域決定テーブルＴの一例である。同図（ａ）が可視強度Ｓ１用のテーブルであり、同図（ｂ）が可視強度Ｓ２用のテーブルである。
【００８２】
演算ビット領域決定テーブル９０１では、判定ビット領域（Ｂ７、Ｂ６）のビット値は演算ビット領域と対応付けられて定義されている。９０１の演算ビット領域中のＹが記入されているビット位置は、演算ビット領域を示し、無記入のビット位置は、非演算ビット領域を示す。即ち演算ビット領域は、判定ビット領域のビット値の関数で与えられる。図示に示す如く、可視強度Ｓ２は、Ｓ１よりも演算ビット領域が総じて小さい（下位ビット側にシフトしている）。この意味は、可視強度Ｓ２の方が原画像に近い状態でノイズ多重化が行われることを意味する。
【００８３】
例えば、所定の可視強度値において、Ｂ７のビット値が１、 Ｂ６のビット値が０の場合、演算ビット領域決定テーブル９０１を参照した場合、同図（ａ）の場合（可視強度Ｓ１）は演算ビット領域はＢ４、Ｂ３、Ｂ２、Ｂ１、Ｂ０となる。また、同条件においては、同図（ｂ）によれば演算ビット領域はＢ２、Ｂ１、Ｂ０となることからも容易に理解できよう。そして、いずれの場合であっても、輝度が高い程、演算ビット領域は大きく、逆に、輝度が低いほど演算ビット領域は小さくなる。これは、先に説明したように、人間の視覚特性は、輝度値が小さいほど輝度値の変化に敏感で、輝度値が大きいほど輝度値の変化が目立たたないという特性に併せてのものである。
【００８４】
図１０はステップＳ７０４の演算ビット領域決定処理を実行する、演算ビット領域決定部のブロック図である。なお、以下では、可視強度Ｓ１における説明であるが、可視強度Ｓ２でも同様である。テーブル１００４の内容が、図９（ａ）であるか同図（ｂ）であるかの違いだけである。
【００８５】
演算ビット領域決定部１０００には、ステップＳ７０２の判定ビット領域参照処理手段から判定ビット領域のビット値１００２、初期設定で設定された可視強度値Ｓ（この場合はＳ１）、演算ビット領域決定テーブルＴがそれぞれ入力される。
【００８６】
演算ビット領域決定部１０００では、既に説明したように、判定ビット領域のビット値と演算ビット領域決定テーブルＴ、さらには可視強度値Ｓとに基づいて演算ビット領域を決定し、演算ビット領域情報として出力する。
【００８７】
図１１はステップＳ７０６の演算処理を実行する、演算部のブロック図である。
【００８８】
演算部１１００には、処理対象である入力画素値１１０２、乱数キーＲ（ここではＲ１）から生成される位置ごとに異なるシリアルビット系列１１０４、ステップ７０４で決定された演算ビット領域情報１１１０がそれぞれ入力される。演算手段１１００は可逆である所定の演算処理に基づき、演算処理（例えば排他的論理和）を施し、演算結果１１０６を出力する。
【００８９】
演算部１１００の動作の一例を以下、図５を用いて再び説明すると次の通りである。説明が重複するが、５０２は入力画素値、５０４は乱数キーＲから生成されるシリアルビット系列、５０６は入力画素とシリアルビット系列の対応するビット位置の排他的論理和（ＸＯＲ）、５０８は演算処理後の出力画素を示す。
【００９０】
まず入力画素５０２の判定ビット領域を参照し、演算ビット領域を決定する。５０２の太枠で囲まれたビット位置は演算ビット領域を示す。シリアルビット系列５０４は画素を復号化する際の鍵となる。なお演算ビット領域以外に対応するシリアルビット系列は必要ない（もしくは０）。
【００９１】
５０６の太枠内は、入力画素の演算ビット領域と対応するシリアルビット系列のビット領域を演算処理結果（ここでは排他的論理和）である。
【００９２】
５０８は５０６の演算処理結果を対応する入力画素の演算ビット領域の値に書き込んだ結果である。
【００９３】
このように本実施形態は、判定ビット領域から演算ビット領域を決定する処理を導入し、画像情報の一部を保存すると同時に、原画像の画素値の大きさに応じて画素値の最大変化量（Δｍａｘ）を決定することが出来る。
【００９４】
図９（ａ）の９０２は、演算ビット領域決定テーブル９０１に対応する画素値の最大変化量（Δｍａｘ）を示している。
【００９５】
一方、可視強度Ｓ２の場合の構成及び動作は、上記説明からすれば容易に理解できるであろうから、その説明は省略する。
【００９６】
図７のステップＳ７０６の演算処理により演算ビット領域のビットのランダマイズ的な効果が予測されるが、演算ビット領域を決定することで、ノイズを加算する透かし画像形状領域内部は、最大変化量Δｍａｘの範囲内でのランダム的な輝度変化を実現できる。
【００９７】
以上の通り、人間の視覚特性は、輝度値が小さいほど輝度値の変化に敏感で、輝度値が大きいほど輝度値の変化が目立たたないという特性がある。従って、演算ビット領域（最大変化量Δｍａｘ）は上記の人間の視覚特性を考慮して設計することが好ましい。
【００９８】
かかる点、図９（ａ）、（ｂ）によれば、人間の視覚特性を考慮して設計された演算ビット領域決定テーブルとなっているので、画質劣化を抑制するように作用する。
【００９９】
なお、空間上の色の変化が人間の知覚に変化と一致する均等知覚色空間ＣＩＥ１９７６Ｌ＊ｕ＊ｖ＊やＣＩＥ １９７６ Ｌ＊ａ＊ｂ＊（以降Ｌ＊ａ＊ｂ＊色空間と呼ぶ）が１９７６年にＣＩＥから推奨されている。
【０１００】
上記の均等知覚色空間は、上記の最大変化量Δｍａｘを決定する際にも有用である。
【０１０１】
次に本第２の実施形態におけるノイズ除去処理について述べる。
【０１０２】
図１２は第２の本実施形態における、ノイズ除去処理の処理を示すフローチャートである。本第２の実施形態におけるノイズ除去処理は、ほぼノイズ付加処理と同一の処理で実現できる。
【０１０３】
先ず、ステップＳ１２０２では、指定された原画像の画素値中の判定ビット領域のビット情報を読み出す（このとき、注目画素の位置が画像形状情報の白の領域にあるのか、黒の領域にあるのかは既に判定済みである点に注意されたい）。
【０１０４】
次にステップＳ１２０４では、ステップＳ１２０２で読み出した判定ビット領域の値に基づき、ステップＳ１２０６演算処理で処理対象となる演算ビット領域を決定する。
【０１０５】
次にステップＳ１２０６では、ステップＳ１２０４で決定された演算ビット領域のビットと初期設定で入力される乱数キーから生成されるシリアルビット系列の間で逆演算処理を行う。この演算は埋め込み時の演算と対応した逆演算処理（復号処理）である。
【０１０６】
次にステップＳ１２０８では、出力画像Ｅに対応する入力画素のビット領域の値を処理された演算ビット領域の値で書き込む、書き込み処理を実行する。
【０１０７】
以上、第２の本実施形態のノイズ除去処理の動作について述べた。
【０１０８】
本第２の実施形態によれば、第１の実施形態に加えて、原画像の透かし画像形状情報で示される領域において、原画像の特徴を維持しつつ、演算処理によって付加されるノイズ量が画像の階調に従って変化するように、画素値の変更量を設定することが出来る。なお、演算ビット領域を決定する際に演算ビット領域決定テーブルＴを用いたが、数式で表現される演算ビット領域決定関数Ｆを用いることも可能であり、この様な方法も本発明の範疇に含まれる。
【０１０９】
また、演算ビット領域を決定するテーブルを２つ示したが、そのテーブルが３以上あって、そのいずれを選択するかＳ１、Ｓ２で決定するようにしてもよい。
【０１１０】
＜第３の実施の形態＞
次に、本第３の実施形態について説明する。本第３の実施形態でも、付加されたノイズの除去が困難な可逆の可視型電子透かしを実現するものである。
【０１１１】
第１、第２の実施形態と異なる点は、注目画素に対してノイズを多重化する場合、その多重化する領域が可視強度Ｓ１、Ｓ２により決定され、最大でもその領域は７ビットであったが、本第３の実施形態では、演算処理により発生するノイズ量が可視強度値Ｓ１、Ｓ２だけでなく、注目画素の近傍領域の画素値に従って変化するようにすることで、この制限を撤廃できる点に在る。
【０１１２】
ノイズ付加や除去の全体の処理は、第１、第２の実施形態における図１、図２と同様であるものとし、その説明は省略し、第１、第２のノイズ付加処理について説明する。
【０１１３】
図１３は本第３の実施形態におけるノイズ付加処理（第１、第２のノイズ付加処理のいずれか）の処理手順を示したフローチャートである。本第３の実施形態のノイズ付加処理は、付加するノイズ量が可視強度値、及び近傍領域としての隣接画素値に従って変化する点を特徴とする。以下、本第３の実施形態のノイズ付加処理について説明する。
【０１１４】
ノイズ付加処理では、まず、ステップＳ１３０２において、初期設定で設定された近傍領域選択方法ＮＳに基づき、埋め込み対象画素の近傍領域を決定する。次に同じく初期設定で設定された近傍解析方法ＮＡに従い、近傍領域内の画素値を解析し、近傍領域解析値を生成する。ここで、近傍領域解析値は、近傍領域から得られる埋め込み対象画素の予測値である近傍領域画素値と、近傍領域の周波数的特性等を含む近傍領域特性値から成る。
【０１１５】
次にステップＳ１３０４では、ステップＳ１３０２で生成した近傍領域解析値や可視強度値Ｓ（Ｓ１、Ｓ２のいずれか）に基づき、ステップＳ１３０６の演算処理で処理対象となる演算ビット領域を決定する。
【０１１６】
次にステップＳ１３０６では、ステップＳ１３０４で決定された演算ビット領域のビットと初期設定で入力される乱数キーＲから生成されるシリアルビット系列の間で演算処理を行う。このとき演算は可逆の処理である必要がある。演算処理としては、排他的論理和やモジュロ加算、モジュロ乗算等の可逆な演算は全て含まれる。
【０１１７】
次にステップＳ１３０８では、出力画像Ｗに対応する入力画素のビット領域の値を処理された演算ビット領域の値で書き込む、書き込み処理を実行する。
【０１１８】
以上、本実施の形態におけるノイズ付加処理の概要を述べた。
【０１１９】
次に、近傍領域解析処理、近傍領域解析値について詳細に説明する。
【０１２０】
本第３の実施形態においては、埋め込み対象画素に可逆のノイズを埋め込むために隣接画素等から構成される近傍領域の解析に基づいて、演算処理する演算ビット領域を決定する。
【０１２１】
一般に自然画像では隣接画素の間では画素値に相関が大きい。即ち、隣接する画素にはよく似た画素値が配置されている場合が多い。そのため、自然画像では、近隣画素における人の目に知覚可能な変更量は、埋め込み対象画素における人の目に知覚可能な変更量として適切である可能性が高い。
【０１２２】
本第３の実施形態では、埋め込み対象画素と相関性の高い近隣画素を参照し、埋め込み対象画素の最大変化量Δｍａｘをきめ細かく設定することで、原画像のどの階調においてもほぼ同様に知覚されるノイズの加算を実現する。
【０１２３】
本実施形態において、埋め込み対象画素への可逆ノイズの埋め込みの為の演算ビット領域を決める近隣画素で構成される領域を近傍領域と呼ぶこととする。
【０１２４】
この近傍領域を構成する画素は実施形態では単一の画素とするが、複数の画素でもよい。また近傍領域は、埋め込み対象画素と相関が高いと予想される領域であればよく、必ずしも埋め込み対象画素に隣接している必要はないが、ある程度の距離以内にあることが望ましい。
【０１２５】
また、近傍領域内の画素を解析する場合、画素値だけでなく、近傍領域内の画素値の周波数特性や近傍領域内の画素値の分散等の統計的性質を利用してもよい。
【０１２６】
また高周波部分やテクスチャ領域において最大変化量Δｍａｘを大きく設定した演算領域決定テーブルＴ＿Ｎを設計すれば、高周波成分から成る領域やテクスチャ領域でも、人の目に均等に認識され易い可逆のノイズを付加することが出来る。
【０１２７】
次に本実施の形態における近傍領域選択解析処理について詳しく説明する。
【０１２８】
図１８は図１３のステップＳ１３０２の近傍領域選択解析処理を実行する近傍領域選択解析部の内部構成を示すブロック図であり、近傍領域選択部１８０２と近傍領域解析部１８０４で構成される。
【０１２９】
近傍領域選択部１８０２には、画像情報（画素位置、画素値）、埋め込み対象画素の位置情報、近傍領域選択方法ＮＳが入力される。近傍領域選択部１８０２は、前記の入力情報に基づき、近傍領域を決定するが、近傍領域は画像全体で固定でなく、画素位置や所定の鍵情報に従って変化させる形態にしてもよい。
【０１３０】
近傍領域選択部１８０２は後段の近傍領域解析部１８０４に近傍領域情報（画素の位置、画素値等）を出力する。
【０１３１】
後段の近傍領域解析部１８０４には、近傍領域情報（画素の位置、画素値等）と近傍領域解析方法ＮＡが入力され、前記の入力情報に基づき、近傍領域の画素値を解析する。そして、近傍領域解析部１８０４は近傍領域解析値（近傍領域画素値と近傍領域特性値）を出力する。
【０１３２】
次に近傍領域選択解析手段の処理について図を用いて具体的に説明する。
【０１３３】
図１５は８ビットグレースケールの入力画像１５０１と可視型電子透かしが埋め込まれた出力画像１５０２の関係を示している。
【０１３４】
図中の太枠で囲まれた画素（画素１３ａ、１４ａ、１５ａ、１８ａ、１９ａ、２０ａ）は、透かし画像形状内に相当し、可逆ノイズの埋め込み対象画素である。
【０１３５】
今、画素１３ａを着目画素とし、その演算ビット領域は、画素１３ａの近傍領域に基づいて選択され、その近傍領域はまず、近傍領域選択手段１８０２を用いて決定される。
【０１３６】
一例として、近傍領域選択部１８０２は、画素１３ａの近傍領域として左横の画素１２ａ（画素値１１２）を選択するとする（近傍領域として複数の画素領域を選択する場合については、後に説明する）。
【０１３７】
次に近傍領域解析部１８０４に画素１２ａ（画素値１１２）を入力する。図１５では説明を簡単にするため、近傍領域解析部１８０４は、入力画素値１１２をそのまま近傍領域解析値として出力するとする。
【０１３８】
次に、演算ビット領域決定処理では、前段の近傍領域選択解析処理から得られる近傍領域解析値を元に、埋め込み対象画素１３ａの演算ビット領域を決定する。
【０１３９】
次に、演算ビット領域決定処理について詳しく説明する。
【０１４０】
図１６はステップＳ１３０４の演算ビット領域決定処理を実行する、演算ビット領域決定部のブロック図である。
【０１４１】
演算ビット領域決定部１６００には、ステップＳ１３０２の近傍領域選択解析処理から入力される近傍領域解析値１６０２、初期設定で設定された可視強度値Ｓ１６０６、演算ビット領域決定テーブルＴ＿Ｎ１６０４がそれぞれ入力される。
【０１４２】
演算ビット領域決定部１６００では、近傍領域解析値１６０２と演算ビット領域決定テーブルＴ＿Ｎ１６０４、さらには可視強度値Ｓ１６０６とに基づいて、埋め込み対象画素の演算ビット領域を決定し、演算ビット領域情報１６０８として出力する。
【０１４３】
次に、演算ビット領域決定テーブルＴ＿Ｎについて説明する。
【０１４４】
演算ビット領域決定テーブルＴ＿Ｎは、ステップＳ１３０４の演算ビット領域決定処理で演算ビット領域を決定する際に用いられる参照テーブルである。
【０１４５】
図１７に近傍領域解析値（近傍領域内の画素の画素値）と対応付けられた演算ビット領域決定テーブルＴ＿Ｎの一例を示す。同図（ａ）は可視強度Ｓ１用の者であり、同図（ｂ）は可視強度Ｓ２用のものである。なお、図示では２種類示したが、３以上あって、その何れを使用するかを示す情報を可視強度Ｓ１、Ｓ２に割り当てても構わない。
【０１４６】
図１７の１７０１は近傍領域解析値を表し（図１７では簡単のために近傍領域画素値のみを表示している）、１７０２は近傍領域解析値に対応する埋め込み対象画素の演算ビット領域（図１７ではＹが記入されているビット位置が演算ビット領域）を表し、１７０３は演算ビット領域から計算される最大変化量Δｍａｘを表している。
【０１４７】
上記のように演算ビット領域決定テーブルＴ＿Ｎは、近傍領域解析値（図１７では近傍領域から計算される近傍領域画素値）及び可視強度値Ｓ１、Ｓ２が演算ビット領域と対応付けられた値を持っている。演算ビット領域決定部１６００では、上記の演算ビット領域決定テーブルＴ＿Ｎを参照し、入力される近傍領域解析値１６０２及び可視強度値Ｓ（１６０６）に対応する演算ビット領域を読み取り、出力する。
【０１４８】
最後にステップＳ１３０６の演算処理では、既に述べた方法で決定された演算ビット領域において、図５に図示したようなシリアルビット系列と演算ビット領域のビットの間でビット演算を実行する。ステップＳ１３０８の書き込み処理では、ステップＳ１３０６の演算結果を出力画像の対応する演算ビット領域に書き込む。
【０１４９】
次に、図１４を用いて、本実施の形態におけるノイズ除去処理の概要を簡単に述べる。
【０１５０】
ノイズ除去処理は、ステップＳ１４０２では、初期設定で設定された近傍領域選択方法ＮＳに基づき、埋め込み対象画素の近傍領域を決定する。次に同じく初期設定で設定された近傍解析方法ＮＡに従い、近傍領域内の画素値を解析し、近傍領域解析値を生成する。
【０１５１】
次にステップＳ１４０４では、ステップＳ１４０２で生成した近傍領域解析値や可視強度値Ｓ（Ｓ１、Ｓ２のいずれか）に基づき、ステップＳ１４０６の演算処理で処理対象となる演算ビット領域を決定する。
【０１５２】
次にステップＳ１４０６では、ステップＳ１４０４で決定された演算ビット領域のビットと初期設定で入力される乱数キーから生成されるシリアルビット系列の間で逆演算処理を行う。この演算は埋め込み時の演算と対応した逆演算処理（復号処理）である。
【０１５３】
ステップＳ１４０８では、入力画素のビット領域の値を処理された演算ビット領域の値を出力画像Ｅの対応画素に書き込む、書き込み処理を実行する。
【０１５４】
以上、本第３の実施形態のノイズ除去処理の動作について述べた。
【０１５５】
完全に可逆のノイズの付加を実現するためには、近傍領域選択方法に関して若干の条件が必要となる。即ち、逆演算時に、埋め込み時に演算処理された演算ビット領域を正しく知る必要がある。従って、近傍領域選択方法は、逆演算時（可逆ノイズ除去時）に演算ビット領域の決定に用いた近傍領域の画素値を参照できるように近傍領域を選択する必要がある。
【０１５６】
上記の条件を満たすような近傍領域の選択方法は、多数考えられるが、埋め込む処理及び除去処理で処理対象となる画素を順次選択する場合のスキャン方向が左から右へ進む場合には、埋め込み対象画素の左横を近傍領域とすることでる。以下、図１５を用いて埋め込みと原画像を必要としない除去方法の一例を説明する。
【０１５７】
まず埋め込みから簡単に説明する。なお、可視強度Ｓ１、Ｓ２の何れであっても処理内容は同じであるので、ここでは可視強度Ｓ１について説明することとする。
【０１５８】
今、埋め込み対象画素が画素１３ａであるとする。この状況において、その左隣の画素１２ａが近傍領域としているわけであるから、先ず、その左隣の画素１２ａ（画素値１１２）が選択、解析される。その結果、近傍領域解析値の近傍領域画素値が「１１２」として出力される。
【０１５９】
次に、演算ビット領域決定テーブルＴ＿Ｎ（図１７（ａ））の中で「１１２」を持つデータを探し出し、画素１３ａの演算ビット領域を決定する。そして、画素１３ａの演算ビット領域のビット値とシリアルビット系列の演算処理により、画素１３ａ（画素値１１６）に対して画素値１２６が計算され、出力画像１５０２の画素１３ａに画素値１２６が書き込まれる。
【０１６０】
次に画素１４ａに対する処理を行うことになる（注目画素位置を１つ右隣に進める）。入力画像１５０１の左横の画素１３ａ（変更前の画素１１６）を参照し、演算ビット領域決定テーブルＴ＿Ｎを用いて、画素１４ａの演算ビット領域を決定する。そして、演算ビット領域とシリアルビット系列の演算処理により、画素１４ａ（画素値１１４）に対して画素値９８が計算され、出力画像１５０２の画素１４ａに画素９８が書き込まれる。
【０１６１】
以上の処理を１ラインのエンドまで行なうことを各ラインについて行うことで、原画像に対するノイズの埋め込み（多重化）を行うことになる。
【０１６２】
次に上述のようにして付加されたノイズを除去するステップについて説明する。説明を簡単なものとするため、ここでも、可視強度Ｓ１の場合について説明する。
【０１６３】
まず、画素１３ａ（ノイズ加算後の画素値「１２６」の画素）が除去対象画素の場合、近傍領域として出力画像１５０２の左横の画素値１１２の画素１２ａ（既に復元処理が行われた画素である点に注意されたい）が選択、解析される。その結果、近傍領域解析値の近傍領域画素値が１１２として出力される。
【０１６４】
次に、演算ビット領域決定テーブルＴ＿Ｎ（図１７（ａ））を用いて、画素１３ａの演算ビット領域を決定する。そして演算ビット領域のビット値とシリアルビット系列の演算処理により、画素１３ａ（透かし埋め込み後の画素値１２６）に対して復元画素値１１６が計算され、出力画像１５０１（原画像）の画素１３ａに画素値１１６が書き込まれる。
【０１６５】
次に画素１４ａが除去対象画素の場合、入力画像１５０１の左横の画素１３ａではなく、既に復元された出力画像１５０２内の左横の画素１３ａ（復元画素値１１６）が選択、解析される。その結果、近傍領域解析値の近傍領域画素値が１１６として出力される。
【０１６６】
次に、演算ビット領域決定テーブルＴ＿Ｎを用いて、画素１４ａの演算ビット領域を決定する。そして演算ビット領域のビット値とシリアルビット系列の演算処理により、画素１４ａ（画素値９８）に対して画素値１１４が計算され、出力画像１５０１（原画像）の画素１４ａに画素値１１４が書き込まれる。
【０１６７】
このようにして、順次、復元処理済の近隣画素の画素値を選択、解析することで、埋め込み時と同じ演算ビット領域を決定することが出来、付加されたノイズを完全に除去することが出来る。
【０１６８】
簡単な説明の為に上記の説明では近傍領域として埋め込み対象画素の左横の画素を選択したが、直前のラインにおける、主走査方向の同じ位置の画素でも構わない。要するに、参照する近傍領域は復元処理が済んでいるものであれば良い。
【０１６９】
また、以下に説明するように、複数の画素にまたがる領域を近傍領域として選択、解析してもよい。
【０１７０】
例えば、近傍領域選択部は、画素１３ａの近傍領域として、画素７ａ、８ａ、１２ａを選択する。次に近傍領域解析部は、画素７ａ、８ａ、１２ａの画素値から、埋め込み対象画素１３ａの画素値を予測し（例えば平均値）、その予測値を近傍領域画素値としてもよい。
【０１７１】
それ以外にも、近傍領域選択部は、画素１３ａの近傍領域として、左上の画素１ａ、２ａ、６ａ、７ａの４画素を選択する方法なども考えられる。この場合、近傍領域解析部は、近傍領域内の分散値、周波数係数等も合わせて計算し、近傍領域特性値としてもよい。この場合、演算ビット領域決定テーブルＴ＿Ｎには、分散値、周波数係数等の近傍領域特性値にも対応する演算ビット領域が定義されているとする。
【０１７２】
以上は、スキャンする順序が１ラインの左から右、各ラインは上から下に向かうことを前提にしているが、スキャンする方向が異なればこれに限るものではない。
【０１７３】
図１７の演算ビット領域決定テーブルＴ＿Ｎの場合、近傍領域解析値（近傍領域画素値）である入力画素値１１２に対応する演算ビット領域は、Ｂ４， Ｂ３， Ｂ１となる。このとき、最大変化量（Δｍａｘ）を計算すると２＾４＋２＾３＋２＾１＝２６となる。例えば、埋め込み対象画素の画素値が１１２の場合には、１１２のＢ４，Ｂ３，Ｂ１はそれぞれ１，０，０であるから、変更量は＋方向には２＾３＋２＾１＝１０， −方向には２＾４＝１６ の幅を持つ。
【０１７４】
本第３の実施形態のノイズ付加方法は、埋め込み対象画素に近い値を持つ近傍領域画素値や埋め込み対象画素近傍の近傍領域特性値からなる近傍領域解析値に基づき、きめ細かな演算ビット領域の設定が可能であることが分かる。
【０１７５】
次に、本第３の実施形態のノイズ付加方法を用い、入力画像の全画素に可逆ノイズを付加する場合について述べる。
【０１７６】
以上の説明からもわかるように、注目画素に対し、その近傍領域が左画素とした場合、各ラインの先頭画素に着目したとき、近傍領域は存在しない。以下、近傍領域が存在しない場合の対処方法例を幾つか述べる。
【０１７７】
例えば、入力画像の端近くで近傍領域が存在しない画素が埋め込み対象画素の場合は、ノイズの付加を中止する方法をとってもよい。このとき、近傍領域が存在しない画素からノイズを除去する際には、ノイズは付加されていないことが分かっているので、近傍領域が得られるまでノイズの除去は実行しなくてよい。
【０１７８】
また、入力画像の端近くで近傍領域が存在しない画素が埋め込み対象画素の場合は、可視強度値Ｓに従ってのみ決められる固定の演算ビット領域に対して、演算ビット処理を行ってもよい。このとき、近傍領域が存在しない画素からノイズを除去する際には、可視強度値Ｓに従ってのみ決められる演算ビット領域に対して、逆演算処理を施し、ノイズの除去は実行するとよい。
【０１７９】
近傍領域選択方法ＮＳによって決められる近傍領域の大きさが得られるような入力画像の画素値に対しては、既に詳しく述べてきたように、近傍領域解析値に従って、演算ビット領域を決定する。
【０１８０】
人間の視覚特性は、輝度値が小さいほど輝度値の変化に敏感で、輝度値が大きいほど輝度値の変化が目立たたないという特性がある。従って、演算ビット領域（最大変化量Δｍａｘ）は上記の人間の視覚特性を考慮して設計することが好ましい。
【０１８１】
既に第１の実施形態でも述べた均等知覚色空間は、上記の演算ビット領域（最大変化量Δｍａｘ）を決定する際にも有用である。
【０１８２】
本実施の形態によれば、原画像の透かし画像形状情報で示される領域において、原画像の特徴を維持しつつ、演算処理によって付加されるノイズ量が画像の階調に従って変化するように、画素値の変更量を設定することが出来る。
【０１８３】
また本第３の実施形態におけるノイズ付加処理は、既に述べたようにノイズ除去処理を行う際に復元されている近傍領域が必要となる。
【０１８４】
従って、透かし画像形状内部に本第３の実施形態で述べたノイズ付加処理を用った後、透かし画像形状外部に第１または第２の実施の形態で述べたノイズ付加処理を用いてノイズを付加するとすれば、透かし画像形状外部のノイズが除去処理されることなしに、透かし画像形状内部のノイズ除去がより困難になる。従って、透かし画像形状内部のノイズ除去には、透かし画像形状内部のノイズ除去のための鍵情報と透かし画像形状外部のノイズ除去のための鍵情報が二重に必要になるという更なる効果も得られる。
【０１８５】
なお、演算ビット領域を決定する際に演算ビット領域決定テーブルＴ＿Ｎを用いたが、数式で表現される演算ビット領域決定関数Ｆを用いることも可能であり、この様な方法も本発明の範疇に含まれる。
【０１８６】
＜第４の実施の形態＞
以上説明した第１乃至第３の実施形態では、画素単位での埋め込みを想定したが、本第４の実施形態では、ＪＰＥＧやＪＰＥＧ２０００等の圧縮符号化画像に対して、可視型電子透かしを埋め込む場合について説明する。
【０１８７】
ＪＰＥＧやＪＰＥＧ２０００等の圧縮符号化方式では、入力色成分は規定されないが、一般にＲ（赤）、Ｇ（緑）、Ｂ（青）の色成分を、Ｙ（輝度）、Ｃｂ（色差）、Ｃｒ（色差）に変換後、離散コサイン変換や離散ウェーブレット変換を実行する場合が多い。
【０１８８】
従って、ＪＰＥＧやＪＰＥＧ２０００圧縮符号化されたカラー画像の輝度成分を表す周波数変換係数を可視型電子透かしの埋め込み成分として用いることで、特別な工夫なしに輝度値に埋め込むことが可能となる。
【０１８９】
ＪＰＥＧ圧縮符号化方式では、ブロック単位の圧縮符号化が行われている。例えば、ＪＰＥＧ圧縮符号化画像には、最小符号化単位（通常８ｘ８画素）が存在し、その単位毎に圧縮符号化の基本的な処理が行われている。従って、ＪＰＥＧ圧縮符号化画像に可視型電子透かしを埋め込む場合、透かし画像形状情報は、画素単位でなく、最小符号化単位とすると、第１乃至第３の実施の形態の手法を応用し易い。
【０１９０】
より詳しく説明すると、画像を８×８画素ブロック単位に周波数成分データに変換するために先ずＤＣＴ変換を行う。そして、この画素ブロックが、画像形状情報の０の領域にあるのか否かを判断し、「１」の領域にあると判断した場合には、ＤＣＴ変換した結果得られる直流成分ＤＣの値を構成するビットに対し、第１乃至第３の実施形態における第１のノイズ付加処理を行ない、「０」の領域にあると判断した場合には第１乃至第３の実施形態における第２のノイズ付加処理を行うことになる。
【０１９１】
図１９の１９０１は、ＪＰＥＧ圧縮符号化方式において、画像を最小符号化単位毎にブロック化して図示した図である。ＪＰＥＧ圧縮符号化画像の場合、最小符号化単位内（１９０１）で離散コサイン変換（ＤＣＴ）が実行される。１９０２は最小符号化単位ごとに得られるＤＣＴ係数を図示したものである。ＤＣＴ係数は直流（平均値）と交流（周波数成分）からなる。
【０１９２】
従って、最小符号化単位ごとに、ＤＣＴ係数のＤＣ成分（直流）に対して、第１の実施の形態で述べた演算ビット領域の演算処理を行えば、最小符号化単位内（１９０１）の平均値を変化させることが出来、ブロック単位の可視型電子透かし埋め込みが実現できる。
【０１９３】
従って透かし画像形状情報は、埋め込み対象最小符号化単位ブロックを指定する情報とすれば、第１乃至第３の実施形態を適用できる。
【０１９４】
なお、第３の実施形態に適用する場合、参照する近傍領域としては、注目画素ブロックの、例えば左隣の画素ブロックのＤＣ成分を参照すると良いであろう。
【０１９５】
以下、各画素ブロックのＤＣＴ変換後の直流成分にノイズを多重化した後は、公知の符号化処理を行ない、符号化することになる。ただし、ＪＰＥＧ圧縮符号化では、量子化処理により或る程度のデータを破棄することになる。通常、この量子化処理では、量子化ステップ（例えば、２とか４等）を決めているので、ＬＳＢもしくは、ＬＳＢから上位に向かう数ビットは除去されることになるので、ノイズ多重化する際には、可視強度Ｓ１、Ｓ２のいずれであっても、この除去されるビットは演算ビット領域から除外することになるであろう。
【０１９６】
一方、ノイズ除去する場合、逆量子化を行うことになるが、これによって得られる直流成分を、第１乃至第３の実施形態における注目画素と見たてて処理すれば良い。第３の実施形態に適用するのであれば、近傍領域は、例えば、復号済みの左隣の画素ブロックにおける直流成分を近傍領域として採用すれば良いであろう。ただし、第３の実施形態に適用する場合、処理単位が８×８画素単位となるので、その画素ブロックが左端にある場合には近傍領域は存在しないので、ノイズ埋め込み時、及び、ノイズ除去時においては、近傍領域を参照せず、自身の画素ブロックに対して固定の演算ビット領域で処理することになるであろう。
【０１９７】
一方、ＪＰＥＧ２０００圧縮符号化画像では、離散ウェーブレット変換（ＤＷＴ）を用いて、画像の形状情報を保持しつつ、画像を低周波から高周波成分に渡って段階的に帯域ごとに分割し、圧縮符号化する。
【０１９８】
図２０はＪＰＥＧ２０００圧縮符号化方式における離散ウェーブレット変換による帯域分割の様子を図的に示した図である。
【０１９９】
離散ウェーブレット変換では、画像に影響が大きい低周波画像成分はＬＬに集まり、ＬＬは原画像の画像特徴をよく保存している。従って、埋め込みに用いる要素を離散ウェーブレット変換（ＤＷＴ）の低周波成分（ＬＬ）とすると、第１乃至第３の実施形態と比較的近い形で可視型電子透かしを埋め込むことが可能である。
【０２００】
図２０の場合、ＬＬは２回ＤＷＴを行った結果の低周波成分のサブブロックであるが、そのサブブロックを構成している全ての画素に対して、第１乃至第３の実施形態と同様の埋め込みを行えばよい。
【０２０１】
なお、第１乃至第３の実施形態で説明した手法をＪＰＥＧ２０００圧縮符号化画像の離散ウェーブレット変換（ＤＷＴ）係数に適用する場合、ＤＷＴ係数は正負の値を取りうることを考慮し、演算ビット領域決定テーブルを設計する必要がある。
【０２０２】
また第１の実施形態をＪＰＥＧ２０００圧縮符号化方式に適応する場合、ノイズ付加処理対象となるＤＷＴ係数をツリー構造内から、可視強度値に応じて、選択してもよい。またサブバンド毎に量子化ステップが設定されている場合、量子化後のＤＷＴ係数を処理対象とすることで可逆のノイズ付加が可能となる。
【０２０３】
また、以下のようにしても良いであろう。
【０２０４】
図２１の２１０１は、図２０から所定の位置のＬＬとツリー構造の関係を満たすＤＷＴ係数を各サブバンドから抜き出したＤＷＴ係数列である。図中、２１０２は、可視強度値とノイズ付加処理対象となるＤＷＴ係数の関係を示す表である。２１０２では、可視強度値の強さに従って、ノイズ付加処理するＤＷＴ係数も変化している。ＬＬに対するノイズ付加処理は、第２の実施の形態のように判定ビット領域を設けて、演算ビット領域を決定した上でノイズ付加処理を施してもよい。その他のサブバンドのＤＷＴ係数は、演算ビット領域を所定の固定値とし、ノイズ付加処理を施してもよい。
【０２０５】
また、ＪＰＥＧ２０００圧縮符号化画像の場合、ＲＯＩ（Ｒｅｇｉｏｎ ｏｆＩｎｔｅｒｅｓｔ：関心領域）のために、画像サイズと同じサイズを持つ１ビットのビットプレーンが用意される（ＪＰＥＧ２０００の基本符号化システムでは、ＲＯＩ領域だけシフトアップして符号化される）。
【０２０６】
従って、関心領域が特に存在せず、透かし画像形状情報を透かし画像として画像閲覧者に示したい場合、透かし画像形状情報をＲＯＩに設定してもよい。
【０２０７】
例えば、著作権情報を示す可視のロゴ情報をＲＯＩに記述しておけば、コンテンツ配信で画像情報を伝送する際に、ロゴ情報から先に視聴者に示すことが出来、利用者にコンテンツの著作権者を明示する効果も実現できる。
【０２０８】
さらに、透かし画像形状情報は既にＲＯＩ情報として、画像と共に符号化されているため、可視型電子透かしを除去する際に必要な鍵情報が僅かですむというメリットがある。
【０２０９】
また、可視型電子透かしの除去に必要な透かし画像形状情報は、画像ファイルのヘッダ等の所定位置に添付することも可能である。このとき、可視型電子透かしが埋め込まれた画像から原画像に戻すときに、画像ファイル以外に必要な鍵情報だけですみ、配送する情報量を削減できる効果もある。
【０２１０】
また、可視型電子透かしの除去に必要な鍵（および透かし画像形状情報）は比較的情報量が少ない為、画像ファイルのヘッダ等の所定位置に添付することも可能である。このとき、特定の利用者のみ可視型電子透かしの除去を可能にするため、上記の鍵（および透かし画像形状情報）を所定の暗号方式（例えば、公開鍵暗号方式）で暗号化し、画像ファイルのヘッダ等の所定位置に添付してもよい。第１乃至第３の実施の形態では、暗号方式としては、非排他的論理和（ＸＯＲ演算）の場合のみ説明したが、所定の処理単位（例えば６４ビット）となるように複数の演算ビット領域を集め、ＤＥＳ等の共通鍵暗号方式や、公開鍵暗号方式を用いることも可能である。
【０２１１】
また、ＤＥＳ等の所定の処理単位毎に処理を行うブロック暗号方式に属する共通鍵暗号方式を用いる場合には、集めたビット数が所定の処理単位に満たない場合には、足りないビット数だけ０や１をパディングし所定単位を満たし、暗号化を行うとよい。このとき元の画素位置に格納できないビットは、ヘッダ等のファイルの所定位置に付加するとよい。
【０２１２】
あるいは、１〜数ビット単位で処理できる（共通鍵暗号方式に属する）ストリーム暗号方式に属する暗号方式を用いてもよい。
【０２１３】
このとき第１〜第３の実施の形態において、初期設定で入力されるのは乱数キーではなく、共通鍵暗号方式では共通鍵、公開鍵暗号方式では埋め込み時には公開鍵、抽出時に秘密鍵となる。
【０２１４】
本実施の形態では、暗号化方式としてＤＥＳを例に挙げたが、他の共通鍵暗号方式、例えば、ＡＥＳ、ＦＥＡＬ、ＩＤＥＡ、ＲＣ２、ＲＣ４、ＲＣ５、ＭＩＳＴＹ、シーザー型暗号、ビジネル暗号・ビューフォート暗号、プレイフェア暗号、ヒル暗号、バーナム暗号等を用いることもできる。
【０２１５】
また本実施の形態では、静止画像についてのみ説明を行ったが、動画像に対しても同様の原理が適用できることは容易に想像できる。例えば、ＭＰＥＧ圧縮符号化方式の場合、中間フレームを埋め込みの対象とすることで、比較的容易に可視型電子透かしを埋め込むことが出来る。ＭｏｔｉｏｎＪＰＥＧ２０００の場合は時間フレーム方向にＪＰＥＧ２０００圧縮符号化画像と同様の手法で可視型電子透かしの埋め込みを繰り返すことで埋め込み可能である。従って、動画像に対して可視型電子透かしを適用することも本発明の範疇であるとする。
【０２１６】
＜ノイズ多重化のサンプルの説明＞
最後に、第１の実施形態で述べたノイズ付加処理を用いて可視型電子透かしを埋め込んだサンプルを添付する。
【０２１７】
ただし、各画像は、本来は多値、すなわち、１画素が多ビットで構成される階調画像であるが、特許出願の明細書に添付する図面は２値画像のみであって多値画像が扱えない。従って、以下に説明する各画像はノイズ多重化した後の多値階調画像そのものではなく、それを更に誤差拡散処理でもって２値化した結果で示していることに留意されたい。
【０２１８】
図２３は横６４０画素×縦４８０画素の８ビットグレースケール画像（原画像）である。図２４は図２３の原画像に対し、図２２で示される透かし画像形状を埋め込んだ画像である。
【０２１９】
図２４では、透かし画像形状内部はＭＳＢ以外のビット位置が演算ビット領域であり、透かし画像形状外部はＭＳＢ， Ｂ６， Ｂ５以外のビット位置が演算ビット領域である。それぞれ鍵から生成されるシリアルビット列との間でＸＯＲ演算を実行し、演算結果で置き換え、可視型電子透かしが埋め込まれた画像（図２４）を生成している。
【０２２０】
図２４の透かし画像形状外部では、上位ビットを多く残していることから、階調性を維持しているが、透かし画像形状内部では、ＭＳＢしか残していないため、非常に階調性が乏しい画像となっている。しかし画像全体が演算処理されているにも関わらず、透かし画像形状内外でノイズ量が異なるため、画像の閲覧者は透かし画像形状情報ははっきりと見ることが出来る。
【０２２１】
上記の実施形態からも明らかなように、本実施形態を実現するその殆どは、ソフトウェアでもって実現できるものである。通常、パーソナルコンピュータ等の汎用情報処理装置にコンピュータプログラムを導入する場合、フロッピー（登録商標）ディスクやＣＤＲＯＭ、或いは、半導体メモリカード等のコンピュータ可読記憶媒体を装置にセットし、インストールプログラムを実行するか、或いは、システムにコピーすることになるので、かかるコンピュータ可読記憶媒体も本願発明の範疇に含まれるのは明らかである。
【０２２２】
また、コンピュータ上で稼働しているＯＳ等が処理の一部又は全部を行う場合、あるいは記憶媒体から読み出されたプログラムコードが、コンピュータに挿入された拡張機能ボードやコンピュータに接続された拡張機能ユニットに備わるメモリに書き込まれた後、そのプログラムコードの指示に基づいて、上記拡張機能ボードや拡張機能ユニットに備わるＣＰＵ等が処理の一部又は全部を行う場合にも、各実施の形態と同等の機能を実現できると共に、同等の効果を得ることができ、本発明の目的を達成することができる。
【０２２３】
以上説明したように、本実施形態によれば、画像及び鍵及び透かし画像形状情報を入力し、画像の構成要素の要素値を選択し、透かし画像形状情報内部であれば、第１のノイズ付加処理を構成要素の要素値に対し実行し、透かし画像形状情報外部であれば、第２のノイズ付加処理を構成要素の要素値に対し実行し、画像に可視の電子透かしを埋め込む。上記の埋め込みにより、可視型電子透かしは、可視型電子透かしが埋め込まれた領域の周囲を解析し、周囲の画素値を元に計算した値で可視型電子透かしの一部を除去することを困難にすると同時に、第１のノイズ付加処理と第２のノイズ付加処理で付加するノイズを視覚的に異ならせることにより、画像の閲覧者に可視型電子透かしを明示することが出来る。
【０２２４】
【発明の効果】
以上説明したように本発明によれば、原画像にノイズを多重化させ、そのノイズを多重化した分布による目視可能な付加情報を可逆的に埋め込むと共に、ノイズ除去に対する高い耐性を有する画像を生成することが可能になる。
【０２２５】
また、他の発明によれば、可視付加情報が多重化された画像から付加情報を除去することでオリジナル画像、もしくはオリジナル画像に近い画像まで復元することが可能になる。
【図面の簡単な説明】
【図１】実施形態における可視型電子透かし埋め込み処理を示すフローチャートである。
【図２】実施形態における可視型電子透かし除去処理を示すフローチャートである。
【図３】第１の実施形態におけるノイズ付加処理のフローチャートである。
【図４】第１の実施の形態における演算ビット領域を表す図である。
【図５】演算処理の具体例を示す図である。
【図６】第１の実施形態におけるノイズ除去処理のフローチャートである。
【図７】第２の実施形態におけるノイズ付加処理のフローチャートである。
【図８】第２の実施形態における判定ビット領域、演算ビット領域、非演算ビット領域の一例を示す図である。
【図９】第２の実施形態における演算ビット領域決定テーブルの一例を示す図である。
【図１０】第２の実施形態における演算ビット領域決定処理を実行する演算ビット領域決定部の内部構成を示すブロック図である。
【図１１】第２の実施形態における演算処理を実行する演算部の内部構成を示すブロック図である。
【図１２】第２の実施形態におけるノイズ除去処理のフローチャートである。
【図１３】第３の実施形態におけるノイズ付加処理のフローチャートである。
【図１４】第３の実施形態におけるノイズ除去処理のフローチャートである。
【図１５】第３の実施形態における原画像とノイズ多重化後の画像の関係を示す図である。
【図１６】第３の実施形態における演算ビット領域決定処理を実行する演算ビット領域決定部の内部構成を示すブロック図である。
【図１７】第３の実施形態における演算ビット領域決定テーブルの一例を示す図である。
【図１８】近傍領域選択解析処理を実行する近傍領域選択解析部の内部構成を示すブロック図である。
【図１９】ＪＰＥＧ圧縮符号化方式における最小符号化単位を図的に示した図である。
【図２０】ＪＰＥＧ２０００圧縮符号化方式における離散ウェーブレット変換による帯域分割の様子を図的に示した図である。
【図２１】所定のＬＬとツリー関係にある離散ウェーブレット変換係数と可視強度値との関係である。
【図２２】本実施形態における透かし画像形状情報の一例を示す図である。
【図２３】第１の実施形態を方法を適用する為の原画像サンプルである。
【図２４】第１の実施形態の方法を適用し、透かし画像形状を埋め込んだ画像サンプルである。
【図２５】実施形態における装置の具体的なブロック構成図である。
Claims (15)
- 多値画像データにノイズを多重化し、ノイズを多重化した分布による目視可能な付加情報を埋め込む画像処理装置であって、
第１の可視強度でもってノイズを多重化する第１のノイズ多重化手段と、
前記第１の可視強度とは異なる第２の可視強度でもってノイズを多重化する第２のノイズ多重化手段と、
所定の画素領域単位に、前記第１、第２の可視強度のいずれで多重化するかを示す情報を、前記付加情報として入力する入力手段と、
前記付加情報に基づき、多値画像データ中の注目画素領域が、前記第１、第２の可視強度のいずれで多重化すべきかを判定する判定手段と、
該判定手段の判定結果に基づいて、注目画素領域に対し、前記第１、第２のノイズ多重化手段のいずれか一方を選択してノイズ多重化を行わせる制御手段と
を備えることを特徴とする画像処理装置。 - 前記第１、第２のノイズ多重化手段それぞれは、
注目画素領域のノイズ多重化する対象となる演算ビット領域を、与えられた可視強度に基づいて決定し、決定された演算ビット領域に、乱数によるノイズのビット情報を可逆的な論理演算を行うことで、ノイズを多重化することを特徴とする請求項１に記載の画像処理装置。 - 前記演算ビット領域は、更に、前記所定画素領域のデータを参照して決定することを特徴とする請求項２に記載の画像処理装置。
- 前記第１、第２のノイズ多重化手段それぞれは、
注目する画素領域のノイズ多重化する対象となる演算ビット領域を、注目画素領域の近傍領域のデータ及び与えられた可視強度に基づいて決定し、決定された演算ビット領域に、乱数によるノイズのビット情報を可逆的な論理演算を行うことで、ノイズを多重化することを特徴とする請求項１に記載の画像処理装置。 - 所定画素領域は単一画素であることを特徴とする請求項１乃至４のいずれか１項に記載の画像処理装置。
- 前記所定画素領域は、所定サイズの画素ブロックであることを特徴とする請求項１乃至４のいずれか１項に記載の画像処理装置。
- 前記所定画素領域は、ＪＰＥＧ圧縮符号化の際の直交変換によって得られた直流成分であることを特徴とする請求項第１項乃至４のいずれか１項に記載の画像処理装置。
- 前記所定画素領域は、ＪＰＥＧ２０００圧縮符号化の際のウェーブレット変換によって得られた低域成分ブロックであることを特徴とする請求項第１項乃至４のいずれか１項に記載の画像処理装置。
- ノイズが可逆的に埋め込まれ、目視可能な付加情報が多重化された多値画像データから、前記付加情報を除去する画像処理装置であって、
第１の可視強度でもって多重化されたノイズを除去する第１のノイズ除去手段と、
前記第１の可視強度とは異なる第２の可視強度でもって多重化されたノイズを除去する第２のノイズ除去手段と、
所定の画素領域単位に、前記第１、第２の可視強度のいずれで多重化されたかを示す情報を、前記付加情報として入力する入力手段と、
前記付加情報に基づき、多値画像データ中の注目画素領域が、前記第１、第２の可視強度のいずれによって多重化されたかを判定する判定手段と、
該判定手段の判定結果に基づいて、注目画素領域に対し、前記第１、第２のノイズ除去手段のいずれか一方を選択して多重化されたノイズを除去する制御手段と
を備えることを特徴とする画像処理装置。 - 多値画像データにノイズを多重化し、ノイズを多重化した分布による目視可能な付加情報を埋め込む画像処理方法であって、
第１の可視強度でもってノイズを多重化する第１のノイズ多重化工程と、
前記第１の可視強度とは異なる第２の可視強度でもってノイズを多重化する第２のノイズ多重化工程と、
所定の画素領域単位に、前記第１、第２の可視強度のいずれで多重化するかを示す情報を、前記付加情報として入力する入力工程と、
前記付加情報に基づき、多値画像データ中の注目画素領域が、前記第１、第２の可視強度のいずれで多重化すべきかを判定する判定工程と、
該判定工程の判定結果に基づいて、注目画素領域に対し、前記第１、第２のノイズ多重化工程のいずれか一方を選択してノイズ多重化を行わせる制御工程と
を備えることを特徴とする画像処理方法。 - ノイズが可逆的に埋め込まれ、目視可能な付加情報が多重化された多値画像データから、前記付加情報を除去する画像処理方法であって、第１の可視強度でもって多重化されたノイズを除去する第１のノイズ除去工程と、
前記第１の可視強度とは異なる第２の可視強度でもって多重化されたノイズを除去する第２のノイズ除去工程と、
所定の画素領域単位に、前記第１、第２の可視強度のいずれで多重化されたかを示す情報を、前記付加情報として入力する入力工程と、
前記付加情報に基づき、多値画像データ中の注目画素領域が、前記第１、第２の可視強度のいずれによって多重化されたかを判定する判定工程と、
該判定工程の判定結果に基づいて、注目画素領域に対し、前記第１、第２のノイズ除去手段のいずれか一方を選択して多重化されたノイズを除去する制御工程と
を備えることを特徴とする画像処理方法。 - 多値画像データにノイズを多重化し、ノイズを多重化した分布による目視可能な付加情報を埋め込む画像処理装置として機能するコンピュータプログラムであって、
第１の可視強度でもってノイズを多重化する第１のノイズ多重化手段と、
前記第１の可視強度とは異なる第２の可視強度でもってノイズを多重化する第２のノイズ多重化手段と、
所定の画素領域単位に、前記第１、第２の可視強度のいずれで多重化するかを示す情報を、前記付加情報として入力する入力手段と、
前記付加情報に基づき、多値画像データ中の注目画素領域が、前記第１、第２の可視強度のいずれで多重化すべきかを判定する判定手段と、
該判定手段の判定結果に基づいて、注目画素領域に対し、前記第１、第２のノイズ多重化手段のいずれか一方を選択してノイズ多重化を行わせる制御手段と
して機能することを特徴とするコンピュータプログラム。 - 請求項１２に記載のコンピュータプログラムを格納することを特徴とするコンピュータ可読記憶媒体。
- ノイズが可逆的に埋め込まれ、目視可能な付加情報が多重化された多値画像データから、前記付加情報を除去する画像処理装置として機能するコンピュータプログラムであって、
第１の可視強度でもって多重化されたノイズを除去する第１のノイズ除去手段と、
前記第１の可視強度とは異なる第２の可視強度でもって多重化されたノイズを除去する第２のノイズ除去手段と、
所定の画素領域単位に、前記第１、第２の可視強度のいずれで多重化されたかを示す情報を、前記付加情報として入力する入力手段と、
前記付加情報に基づき、多値画像データ中の注目画素領域が、前記第１、第２の可視強度のいずれによって多重化されたかを判定する判定手段と、
該判定手段の判定結果に基づいて、注目画素領域に対し、前記第１、第２のノイズ除去手段のいずれか一方を選択して多重化されたノイズを除去する制御手段と
して機能することを特徴とするコンピュータプログラム。 - 請求項１４に記載のコンピュータプログラムを格納することを特徴とするコンピュータ可読記憶媒体。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002191133A JP3919613B2 (ja) | 2002-06-28 | 2002-06-28 | 画像処理装置及び方法、並びにコンピュータプログラム及びコンピュータ可読記憶媒体 |
| US10/600,620 US7200242B2 (en) | 2002-06-28 | 2003-06-23 | Image processing apparatus and method, computer program, and computer-readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002191133A JP3919613B2 (ja) | 2002-06-28 | 2002-06-28 | 画像処理装置及び方法、並びにコンピュータプログラム及びコンピュータ可読記憶媒体 |
Publications (3)
| Publication Number | Publication Date |
|---|---|
| JP2004040237A true JP2004040237A (ja) | 2004-02-05 |
| JP2004040237A5 JP2004040237A5 (ja) | 2005-10-13 |
| JP3919613B2 JP3919613B2 (ja) | 2007-05-30 |
Family
ID=31492009
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2002191133A Expired - Fee Related JP3919613B2 (ja) | 2002-06-28 | 2002-06-28 | 画像処理装置及び方法、並びにコンピュータプログラム及びコンピュータ可読記憶媒体 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US7200242B2 (ja) |
| JP (1) | JP3919613B2 (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009147788A1 (ja) * | 2008-06-04 | 2009-12-10 | パナソニック株式会社 | 暗号化装置及び暗号化システム |
| JP2019220925A (ja) * | 2018-06-22 | 2019-12-26 | パナソニックＩｐマネジメント株式会社 | 撮像装置及び映像出力方法 |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8049933B2 (en) | 2003-09-17 | 2011-11-01 | Canon Kabushiki Kaisha | Copy-forgery-inhibited pattern image generation method and image processing apparatus |
| US20050058476A1 (en) | 2003-09-17 | 2005-03-17 | Canon Kabushiki Kaisha | Copy-forgery-inhibited pattern density parameter determination method, copy-forgery-inhibited pattern image generation method, and image processing apparatus |
| JP4428036B2 (ja) * | 2003-12-02 | 2010-03-10 | ソニー株式会社 | 情報処理装置および方法、プログラム、並びに、情報処理システムおよび方法 |
| US8355018B2 (en) * | 2007-06-15 | 2013-01-15 | Ricoh Co., Ltd. | Independent pixel waveforms for updating electronic paper displays |
| US8416197B2 (en) * | 2007-06-15 | 2013-04-09 | Ricoh Co., Ltd | Pen tracking and low latency display updates on electronic paper displays |
| US8319766B2 (en) * | 2007-06-15 | 2012-11-27 | Ricoh Co., Ltd. | Spatially masked update for electronic paper displays |
| US8913000B2 (en) * | 2007-06-15 | 2014-12-16 | Ricoh Co., Ltd. | Video playback on electronic paper displays |
| US8279232B2 (en) | 2007-06-15 | 2012-10-02 | Ricoh Co., Ltd. | Full framebuffer for electronic paper displays |
| US8203547B2 (en) * | 2007-06-15 | 2012-06-19 | Ricoh Co. Ltd | Video playback on electronic paper displays |
| JP5046047B2 (ja) * | 2008-10-28 | 2012-10-10 | セイコーインスツル株式会社 | 画像処理装置、及び画像処理プログラム |
| BR102014007666B1 (pt) * | 2014-03-28 | 2023-01-10 | Samsung Eletrônica Da Amazônia Ltda | Método para autenticação de transações móveis usando criptografia de vídeo e método para criptografia de vídeo |
| EP4207736A4 (en) | 2020-08-27 | 2024-01-31 | Panasonic Ip Corp America | IMAGE PROCESSING DEVICE, METHOD FOR TRAINING A MACHINE LEARNING MODEL, IDENTIFICATION DEVICE AND IMAGE PROCESSING METHOD |
| CN115439702B (zh) * | 2022-11-08 | 2023-03-24 | 武昌理工学院 | 一种基于频域处理的弱噪声图像分类方法 |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6611607B1 (en) * | 1993-11-18 | 2003-08-26 | Digimarc Corporation | Integrating digital watermarks in multimedia content |
| JP3434088B2 (ja) | 1995-01-17 | 2003-08-04 | 日本電信電話株式会社 | 画像データ変換装置およびその逆変換装置 |
| US5530759A (en) | 1995-02-01 | 1996-06-25 | International Business Machines Corporation | Color correct digital watermarking of images |
| US6233347B1 (en) * | 1998-05-21 | 2001-05-15 | Massachusetts Institute Of Technology | System method, and product for information embedding using an ensemble of non-intersecting embedding generators |
| US6154571A (en) * | 1998-06-24 | 2000-11-28 | Nec Research Institute, Inc. | Robust digital watermarking |
| JP3748495B2 (ja) | 1998-12-11 | 2006-02-22 | キヤノン株式会社 | 画像処理方法、画像処理装置、画像配付システム及び記憶媒体 |
| US6801636B2 (en) * | 2000-01-31 | 2004-10-05 | Canon Kabushiki Kaisha | Image processing apparatus and method, and storage medium |
| JP4218920B2 (ja) * | 2000-04-07 | 2009-02-04 | キヤノン株式会社 | 画像処理装置及び画像処理方法並びに記憶媒体 |
| US6996248B2 (en) * | 2001-06-13 | 2006-02-07 | Qualcomm, Incorporated | Apparatus and method for watermarking a digital image |
-
2002
- 2002-06-28 JP JP2002191133A patent/JP3919613B2/ja not_active Expired - Fee Related
-
2003
- 2003-06-23 US US10/600,620 patent/US7200242B2/en not_active Expired - Fee Related
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009147788A1 (ja) * | 2008-06-04 | 2009-12-10 | パナソニック株式会社 | 暗号化装置及び暗号化システム |
| US8484485B2 (en) | 2008-06-04 | 2013-07-09 | Panasonic Corporation | Encryption device and encryption system |
| JP5346933B2 (ja) * | 2008-06-04 | 2013-11-20 | パナソニック株式会社 | 暗号化装置及び暗号化システム |
| JP2019220925A (ja) * | 2018-06-22 | 2019-12-26 | パナソニックＩｐマネジメント株式会社 | 撮像装置及び映像出力方法 |
| JP7149514B2 (ja) | 2018-06-22 | 2022-10-07 | パナソニックＩｐマネジメント株式会社 | 撮像装置及び映像出力方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| US7200242B2 (en) | 2007-04-03 |
| JP3919613B2 (ja) | 2007-05-30 |
| US20040028256A1 (en) | 2004-02-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3958128B2 (ja) | 画像処理装置及び方法、並びにコンピュータプログラム及びコンピュータ可読記憶媒体 | |
| Kutter et al. | Digital signature of color images using amplitude modulation | |
| Busch et al. | Digital watermarking: From concepts to real-time video applications | |
| US6973195B1 (en) | Method for embedding electronic watermark and method for decoding the same | |
| Yang et al. | A contrast-sensitive reversible visible image watermarking technique | |
| US6697498B2 (en) | Method and computer program product for hiding information in an indexed color image | |
| JP3919613B2 (ja) | 画像処理装置及び方法、並びにコンピュータプログラム及びコンピュータ可読記憶媒体 | |
| JP3960597B2 (ja) | 符号生成装置、画像処理装置、符号生成プログラム、画像処理プログラムおよび記憶媒体 | |
| US20020120849A1 (en) | Parallel processing of digital watermarking operations | |
| Yao et al. | Content-adaptive reversible visible watermarking in encrypted images | |
| JP2004173072A (ja) | 符号生成装置、画像処理装置、符号生成プログラム、画像処理プログラムおよび記憶媒体 | |
| JP3919612B2 (ja) | 画像処理装置及び方法、並びにコンピュータプログラム及びコンピュータ可読記憶媒体 | |
| JP4054619B2 (ja) | 画像処理装置及び方法、並びにコンピュータプログラム及びコンピュータ可読記憶媒体 | |
| Phadikar et al. | Quality access control of compressed color images using data hiding | |
| Ogihara et al. | Data embedding into pictorial images with less distortion using discrete cosine transform | |
| Chen et al. | Robust spatial LSB watermarking of color images against JPEG compression | |
| Negrat et al. | Variable length encoding in multiple frequency domain steganography | |
| Kumar | Steganography using visual cryptography | |
| Verma et al. | Robust color image watermarking scheme using JFIF-YCbCr color space in wavelet domain | |
| Sadh et al. | Dual plane multiple spatial watermarking with self-encryption | |
| Kao et al. | Information Hiding in Lossy Compression Gray Scale Image | |
| Kothari et al. | Video watermarking in spatial domain | |
| JP2004343474A (ja) | 暗号情報埋め込み装置及び暗号情報再生装置 | |
| Muñoz-Ramírez et al. | Steganographic framework for hiding a color image into digital images | |
| Patil et al. | DWT Based Mosaic Image Steganography |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050609 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050609 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20070123 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20070126 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070213 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100223 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110223 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120223 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130223 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140223 Year of fee payment: 7 |
|
| LAPS | Cancellation because of no payment of annual fees |