JP2006014171A - Electronic watermark inserting method, electronic watermark inserting apparatus, electronic watermark inserting program, electronic watermark detection method, electronic watermark detecting apparatus, electronic watermark detection program, data transmission method, data transmission system, and data transmission program - Google Patents
Electronic watermark inserting method, electronic watermark inserting apparatus, electronic watermark inserting program, electronic watermark detection method, electronic watermark detecting apparatus, electronic watermark detection program, data transmission method, data transmission system, and data transmission program Download PDFInfo
- Publication number
- JP2006014171A JP2006014171A JP2004191408A JP2004191408A JP2006014171A JP 2006014171 A JP2006014171 A JP 2006014171A JP 2004191408 A JP2004191408 A JP 2004191408A JP 2004191408 A JP2004191408 A JP 2004191408A JP 2006014171 A JP2006014171 A JP 2006014171A
- Authority
- JP
- Japan
- Prior art keywords
- circle
- image
- positions
- frequency component
- predetermined number
- 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
Landscapes
- Image Processing (AREA)
- Editing Of Facsimile Originals (AREA)
- Facsimile Image Signal Circuits (AREA)
- Image Analysis (AREA)
Abstract
Description
本発明は、画像データに電子透かしを挿入するための電子透かし挿入方法、その装置及びそのプログラム、電子透かしが挿入された画像データから電子透かしを検出するための電子透かし検出方法、その装置及びそのプログラム並びに電子透かしを用いてデータを伝送するためのデータ伝送方法、データ伝送システム及びデータ伝送プログラムに関する。 The present invention relates to a digital watermark insertion method for inserting a digital watermark into image data, its apparatus and its program, a digital watermark detection method for detecting a digital watermark from image data into which the digital watermark has been inserted, its apparatus and its The present invention relates to a program, a data transmission method for transmitting data using digital watermarking, a data transmission system, and a data transmission program.
近年、画像データに電子透かしを挿入する技術が開発され、著作権等を有する者等を特定する情報、著作物等の利用を許諾する場合の利用方法及び条件に関する情報等を電子透かしの形態で画像データに挿入することが実用化されてきている。 In recent years, a technique for inserting a digital watermark into image data has been developed, and information for identifying a person who has a copyright, etc., information on usage methods and conditions for permitting the use of a copyrighted work, etc. in the form of a digital watermark. Insertion into image data has been put into practical use.
通常は電子透かしは電子データの形態で挿入され、検出されるが、印刷にも対応した電子透かしも登場してきている。これは、電子透かしが挿入された画像データが印刷された場合であっても、印刷された画像データから電子透かしを検出できるようにしたものである。 Normally, a digital watermark is inserted and detected in the form of electronic data, but a digital watermark corresponding to printing has also appeared. This makes it possible to detect a digital watermark from the printed image data even when the image data into which the digital watermark is inserted is printed.
ところが、従来からの印刷にも対応した電子透かしは、スキャナを用いて印刷物を読み取ることを前提としている。このような前提のもとでは、読み取られた画像の縮尺は予め知られており、画像が回転していることはない。また、スキャナは6400dpi程度の高解像度を有するので高周波成分に埋め込まれている電子透かしを検出することができる。更に、カメラで読み取る場合と異なり手振れが生ずる場合もない。更に、スキャナにより一定の光源により画像が照射されるので、信号対雑音比、コントラストもよい。 However, conventional digital watermarking that also supports printing is premised on reading a printed matter using a scanner. Under such a premise, the scale of the read image is known in advance, and the image is never rotated. Further, since the scanner has a high resolution of about 6400 dpi, it can detect a digital watermark embedded in a high frequency component. Further, unlike the case of reading with a camera, there is no case of camera shake. Further, since the image is emitted from the scanner with a constant light source, the signal-to-noise ratio and the contrast are good.
一方、近年カメラを内蔵し、電子メール機能及びWWW(World Wide Web)ブラウザ機能を有する携帯電話が普及してきている。電子透かしが挿入された画像をこのような携帯電話のカメラで読み取り、携帯電話からそのような画像を電子メールの添付ファイル又はホームページのフォーム情報として送れば、WWWサーバ側は電子透かしを検出し、それに基づいてデータベースから情報を検索し、携帯電話はWWWサーバ、メールサーバ等からその様な情報の提供等の様々なサービスを受けることができる。 On the other hand, cellular phones with built-in cameras and electronic mail functions and WWW (World Wide Web) browser functions have recently become widespread. If an image with a digital watermark inserted is read by the camera of such a mobile phone, and the image is sent from the mobile phone as an attached file of an e-mail or form information on a homepage, the WWW server side detects the digital watermark, Based on this, information is retrieved from the database, and the mobile phone can receive various services such as provision of such information from a WWW server, a mail server, or the like.
従って、携帯電話の内蔵カメラ等で読み取った画像から電子透かしを検出することを可能とする電子透かし挿入方法、その装置及びそのプログラム並びに電子透かし検出方法、その装置及びそのプログラムの登場が望まれていた。 Therefore, the advent of a digital watermark insertion method, apparatus and program thereof, and digital watermark detection method, apparatus and program thereof that enable detection of a digital watermark from an image read by a built-in camera or the like of a mobile phone is desired. It was.
しかしながら、携帯電話の内蔵カメラ等で読み取った画像は、画像とカメラの間の距離等が変動要因となり、縮尺が不定である。また、携帯電話の内蔵カメラ等で読み取った画像は、画像とカメラの間の相対回転角等が変動要因となり、少なからず回転している場合が多い。更に、大多数の携帯電話の内蔵カメラで読み取った画像は、解像度が低い。例えば、3cm×4cmの画像を6400dpiのスキャナで読み取った場合の画素数は、約7600万であるのに対し、現在普及している大多数の携帯電話の内蔵カメラの画素数は30万程度であることより明らかである。更に、携帯電話の内蔵カメラ等で読み取った画像は、手振れが生じている場合がある。更に、携帯電話の内蔵カメラ等で読み取った画像は、一般の外光のもとで読み取られるので、信号対雑音比及びコントラストが十分でない場合が多い。 However, the scale of an image read by a built-in camera of a mobile phone is indefinite due to factors such as the distance between the image and the camera. In addition, an image read by a built-in camera of a mobile phone is often rotated not a little because of a variation factor such as a relative rotation angle between the image and the camera. Furthermore, the resolution of the images read by the built-in cameras of most mobile phones is low. For example, the number of pixels when a 3 cm × 4 cm image is read by a 6400 dpi scanner is about 76 million, whereas the number of built-in cameras of the majority of mobile phones currently popular is about 300,000. It is clearer than there is. Further, an image read by a built-in camera of a mobile phone may cause camera shake. Furthermore, since an image read with a built-in camera of a mobile phone is read under ordinary external light, the signal-to-noise ratio and contrast are often insufficient.
そこで、縮尺が不定であり、回転していて、解像度が低く、手振れが生じていて、信号対雑音比が悪く、コントラストも悪い画像であっても電子透かしの形態で埋め込まれたデータを安定して検出することを可能とするロバストな電子透かし挿入方法、その装置及びそのプログラム並びに電子透かし検出方法、その装置およびそのプログラムを提供することを目的とした発明が発明され、特許文献1で公開されている。
しかしながら、特許文献1に記載されている発明のように、拡大縮小率を検出するための非ゼロ周波数成分を周波数領域に設けると、これを逆フーリエ変換して原画像に加算した場合に、原画像の画質がある程度落ちてしまう。また、回転角度を検出するための点及びデータを表す非ゼロ周波数成分も周波数領域に設けられ、これを逆フーリエ変換したものも原画像に加算するので、これによっても原画像の画質がある程度落ちてしまう。従って、周波数領域に設ける非ゼロ周波数成分の数を、むやみに増やさずに、できるだけ減らすことが、画質を低下させないという観点からみると、望ましい。
However, when a non-zero frequency component for detecting the enlargement / reduction ratio is provided in the frequency domain as in the invention described in
一方、同じ数の非ゼロ周波数成分を周波数領域に設けるのであれば、できる限り多くの数の非ゼロ周波数成分をデータを表すために利用することが、埋め込むべきデータのとりうる値の数を増やすという観点からみると、望ましい。そうすれば、より多くの情報量を有するデータを画像に埋め込むことが可能となる。 On the other hand, if the same number of non-zero frequency components are provided in the frequency domain, using as many non-zero frequency components as possible to represent the data increases the number of possible values of the data to be embedded. From this point of view, it is desirable. Then, data having a larger amount of information can be embedded in the image.
また、特許文献1に記載されている発明によると、データを表す非ゼロ周波数成分が挿入される周波数領域における円の半径が予め決まっているので、撮影した画像の拡大縮小率が判ると、データを表す非ゼロ周波数成分が挿入されている円の半径位置が予測できてしまうため、データの秘匿性が弱く、従って、例えば、悪意を持った第三者に悪用されるおそれがあった。
In addition, according to the invention described in
そこで、本発明は、周波数領域に設ける点のうちデータを表すために利用する非ゼロ周波数成分を多くすることを可能とする電子透かし挿入方法、その装置及びそのプログラム、電子透かし検出方法、その装置及びそのプログラム並びにデータ伝送方法、そのシステム及びそのプログラムを提供することを目的とする。 Therefore, the present invention provides a digital watermark insertion method, apparatus and program thereof, digital watermark detection method and apparatus capable of increasing non-zero frequency components used for representing data among points provided in the frequency domain. And a program thereof, a data transmission method, a system thereof, and a program thereof.
また、本発明は、データの秘匿性を高めることを可能とする電子透かし挿入方法、その装置及びそのプログラム、電子透かし検出方法、その装置及びそのプログラム並びにデータ伝送方法、そのシステム及びそのプログラムを提供することを目的とする。 The present invention also provides a digital watermark insertion method, apparatus and program thereof, digital watermark detection method, apparatus and program thereof, data transmission method, system and program thereof that can enhance data confidentiality. The purpose is to do.
本発明の第1の観点によれば、周波数領域において、半径が不定な第1の円上に並ぶ第1の所定数の候補位置のうちの挿入すべきデータにより相互に組み合わされる位置であって前記第1の所定数より少ない第2の所定数のものに非ゼロ周波数成分を生成するステップと、前記周波数領域において、前記第1の円とは異なる第2の円上の位置に画像の回転角度を検出するための非ゼロ周波数成分を生成するステップと、生成された前記非ゼロ周波数成分を合わせて前記周波数領域から空間領域に変換するステップと、原画像の全部又は一部のブロックに前記変換により得た画像を加算又は減算するステップと、を備えることを特徴とする電子透かし挿入方法が提供される。 According to the first aspect of the present invention, in the frequency domain, the first predetermined number of candidate positions arranged on a first circle having an indefinite radius are positions that are combined with each other by data to be inserted. Generating a non-zero frequency component for a second predetermined number less than the first predetermined number, and rotating the image to a position on a second circle different from the first circle in the frequency domain Generating a non-zero frequency component for detecting an angle; transforming the generated non-zero frequency component together from the frequency domain to a spatial domain; and adding all or a part of a block of an original image to the block And a step of adding or subtracting an image obtained by the conversion.
上記の電子透かし挿入方法において、前記第2の円上の位置とは、周波数領域の原点を中心として180度以外の角度を回転させた後の全ての該位置が回転させる前の全ての該位置と重なるようなことが、該角度を180度以外のいかなる値としても、ないようなものであってもよい。 In the above-described digital watermark insertion method, the positions on the second circle are all the positions before all the positions after being rotated by an angle other than 180 degrees around the origin of the frequency domain. The angle may not be any value other than 180 degrees.
上記の電子透かし挿入方法において、前記第1の円及び前記第2の円とは異なる第3の円上に並ぶ第3の所定数の候補位置のうちの挿入すべき前記データにより相互に組合わされる位置であって前記第3の所定数より少ない第4の所定数のものに非ゼロ周波数成分を生成するステップを更に備え、前記変換のステップでは、前記第3の円上に生成された非ゼロ周波数成分も前記周波数領域から前記空間領域に変換してもよい。 In the above-described electronic watermark insertion method, the data to be inserted are combined with each other among the third predetermined number of candidate positions arranged on a third circle different from the first circle and the second circle. And generating a non-zero frequency component at a fourth predetermined number less than the third predetermined number, wherein the converting step includes generating a non-zero frequency component on the third circle. Zero frequency components may also be transformed from the frequency domain to the spatial domain.
本発明の第2の観点によれば、画像より1以上の画像ブロックを得るステップと、前記1以上の画像ブロックのデータを空間領域から周波数領域に変換するステップと、前記1以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得るステップと、前記周波数領域において、前記各周波数成分の振幅を基に、前記画像の回転角度を検出するための非ゼロ周波数成分を検出することにより前記画像の回転角度を検出するステップと、前記周波数領域において、円周上の第1の所定数の位置であって検出された前記回転角度を回転基準としたもののうちの前記第1の所定数よりも少ない第2の所定数の何れかの位置に非ゼロ周波数成分があるか否かを異なる半径の各円毎に検出するステップと、前記第2の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せを基に、前記画像に挿入されているデータを検出するステップと、を備えることを特徴とする電子透かし検出方法が提供される。 According to a second aspect of the present invention, obtaining one or more image blocks from an image, converting data of the one or more image blocks from a spatial domain to a frequency domain, Obtaining an amplitude of each frequency component based on each frequency component; and detecting a non-zero frequency component for detecting a rotation angle of the image based on the amplitude of each frequency component in the frequency domain. And detecting the rotation angle of the image by the first predetermined number of rotations at the first predetermined number of positions on the circumference and using the detected rotation angles as a rotation reference in the frequency domain. Detecting whether there is a non-zero frequency component at any of a second predetermined number of positions less than the number for each circle of different radii, and at any of the second predetermined number of positions Said has been detected based on the combination of the position of the non-zero frequency component, digital watermark detection method characterized by comprising the steps of: detecting the data inserted in the image is provided.
上記の電子透かし検出方法において、前記画像の回転角度を検出するための非ゼロ周波数成分は、円上の位置に分布していてもよい。 In the digital watermark detection method, the non-zero frequency component for detecting the rotation angle of the image may be distributed at positions on a circle.
上記の電子透かし検出方法において、前記円上の位置とは、周波数領域の原点を中心として180度以外の角度を回転させた後の全ての該位置が回転させる前の全ての該位置と重なるようなことが、該角度を180度以外のいかなる値としても、ないようなものであってもよい。 In the above-described digital watermark detection method, the positions on the circle are such that all the positions after rotating an angle other than 180 degrees centering on the origin of the frequency domain overlap all the positions before the rotation. However, the angle may not be any value other than 180 degrees.
上記の電子透かし検出方法において、データに係る前記非ゼロ周波数成分があるか否かを異なる半径毎に検出するステップは、半径が短い円から半径が長い円に向かって行われ、データに係る前記非ゼロ周波数成分があるか否かを異なる半径毎に検出するステップにより、データに係る前記非ゼロ周波数成分が検出された後に、前記周波数領域において、円周上の第3の位置であって検出された前記回転画像を回転基準としたもののうちの前記第3の所定数よりも少ない第4の所定数の何れかの位置に非ゼロ周波数成分があるか否かを、前記データに係る非ゼロ周波数成分が検出された円よりも半径が長い各円毎に検出するステップを更に備え、前記画像に挿入されているデータを検出するステップでは、前記第2の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せのみならず、前記第4の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せも基にして、前記画像に挿入されているデータを検出してもよい。 In the above-described digital watermark detection method, the step of detecting whether there is the non-zero frequency component related to data for each different radius is performed from a circle with a short radius toward a circle with a long radius, and the data related to the data By detecting whether or not there is a non-zero frequency component for each different radius, after the non-zero frequency component related to data is detected, it is detected at a third position on the circumference in the frequency domain. Whether or not there is a non-zero frequency component at any position of the fourth predetermined number smaller than the third predetermined number of the rotation images that are used as a rotation reference is determined based on the non-zero value related to the data The method further comprises a step of detecting each circle having a radius longer than that of the circle in which the frequency component is detected, and the step of detecting data inserted in the image has a position at any one of the second predetermined number. Not only the combination of the positions of the non-zero frequency components detected to be detected, but also the combination of the positions of the non-zero frequency components detected to be in any one of the fourth predetermined number of positions. The data inserted in the image may be detected.
上記の電子透かし検出方法において、データに係る前記非ゼロ周波数成分があるか否かを異なる半径毎に検出するステップは、半径が長い円から半径が短い円に向かって行われ、データに係る前記非ゼロ周波数成分があるか否かを異なる半径毎に検出するステップにより、データに係る前記非ゼロ周波数成分が検出された後に、前記周波数領域において、円周上の第3の位置であって検出された前記回転画像を回転基準としたもののうちの前記第3の所定数よりも少ない第4の所定数の何れかの位置に非ゼロ周波数成分があるか否かを、前記データに係る非ゼロ周波数成分が検出された円よりも半径が短い各円毎に検出するステップを更に備え、前記画像に挿入されているデータを検出するステップでは、前記第2の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せのみならず、前記第4の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せも基にして、前記画像に挿入されているデータを検出してもよい。 In the above digital watermark detection method, the step of detecting whether or not there is the non-zero frequency component related to the data for each different radius is performed from a circle with a long radius toward a circle with a short radius, and the data related to the data By detecting whether or not there is a non-zero frequency component for each different radius, after the non-zero frequency component related to data is detected, it is detected at a third position on the circumference in the frequency domain. Whether or not there is a non-zero frequency component at any position of the fourth predetermined number smaller than the third predetermined number of the rotation images that are used as a rotation reference is determined based on the non-zero value related to the data The method further comprises a step of detecting each circle having a radius shorter than that of the circle in which the frequency component is detected, and the step of detecting the data inserted in the image has a position at any one of the second predetermined number. Not only the combination of the positions of the non-zero frequency components detected to be detected, but also the combination of the positions of the non-zero frequency components detected to be in any one of the fourth predetermined number of positions. The data inserted in the image may be detected.
上記の電子透かし検出方法において、前記画像を空間フィルタに通すステップを更に備え、前記画像より1以上の画像ブロックを得るステップは、前記空間フィルタを通った画像より1以上の画像ブロックを得るステップであってもよい。 In the above-described digital watermark detection method, the method further includes the step of passing the image through a spatial filter, and the step of obtaining one or more image blocks from the image is a step of obtaining one or more image blocks from the image passed through the spatial filter. There may be.
本発明の第3の観点によれば、周波数領域において、半径が不定な第1の円上に並ぶ第1の所定数の候補位置のうちの挿入すべきデータにより相互に組み合わされる位置であって前記第1の所定数より少ない第2の所定数のものに非ゼロ周波数成分を生成するステップと、前記周波数領域において、前記第1の円とは異なる第2の円上の位置に画像の回転角度を検出するための非ゼロ周波数成分を生成するステップと、生成された前記非ゼロ周波数成分を合わせて前記周波数領域から空間領域に変換するステップと、原画像の全部又は一部のブロックに前記変換により得た画像を加算又は減算するステップと、前記加算又は減算を経た画像より1以上の画像ブロックを得るステップと、前記1以上の画像ブロックのデータを空間領域から周波数領域に変換するステップと、前記1以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得るステップと、前記周波数領域において、前記各周波数成分の振幅を基に、前記画像の回転角度を検出するための非ゼロ周波数成分を検出することにより前記画像の回転角度を検出するステップと、前記周波数領域において、円周上の第1の所定数の位置であって検出された前記回転角度を回転基準としたもののうちの前記第1の所定数よりも少ない第2の所定数の何れかの位置に非ゼロ周波数成分があるか否かを異なる半径の各円毎に検出するステップと、前記第2の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せを基に、前記画像に挿入されているデータを検出するステップと、を備えることを特徴とするデータ伝送方法が提供される。 According to the third aspect of the present invention, in the frequency domain, the first predetermined number of candidate positions arranged on a first circle having an indefinite radius are positions that are combined with each other by data to be inserted. Generating a non-zero frequency component for a second predetermined number less than the first predetermined number, and rotating the image to a position on a second circle different from the first circle in the frequency domain Generating a non-zero frequency component for detecting an angle; transforming the generated non-zero frequency component together from the frequency domain to a spatial domain; and adding all or a part of a block of an original image to the block A step of adding or subtracting the image obtained by the conversion, a step of obtaining one or more image blocks from the image that has undergone the addition or subtraction, and the frequency of the data of the one or more image blocks from a spatial domain Converting to a region; obtaining an amplitude of each frequency component based on each frequency component of the one or more image blocks; and rotating the image based on the amplitude of each frequency component in the frequency region. Detecting a rotation angle of the image by detecting a non-zero frequency component for detecting an angle; and the rotation detected at a first predetermined number of positions on a circumference in the frequency domain. Detecting for each circle of different radii whether or not there is a non-zero frequency component at any of a second predetermined number less than the first predetermined number of angles based on rotation. Detecting data inserted in the image based on a combination of positions of the non-zero frequency components detected to be at any one of the second predetermined number of positions. Data transmission method according to claim is provided.
本発明によれば、拡大縮小率を検出するための非ゼロ周波数成分を配置する必要が無いので、非ゼロ周波数成分の数を所定数とするという前提のもとでは、データを表すために非ゼロ周波数成分を有効に利用することがでる。 According to the present invention, since it is not necessary to arrange non-zero frequency components for detecting the enlargement / reduction ratio, the non-zero frequency components are represented in order to represent data under the assumption that the number of non-zero frequency components is a predetermined number. The zero frequency component can be used effectively.
また、本発明によれば、拡大縮小率を検出するための非ゼロ周波数成分を配置する必要が無いので、データを表すための非ゼロ周波数成分の数及び回転角度を検出するための非ゼロ周波数成分の数を所定数とするという前提のもとでは、画質を向上させることができる。 In addition, according to the present invention, since it is not necessary to arrange non-zero frequency components for detecting the enlargement / reduction ratio, the number of non-zero frequency components for representing data and the non-zero frequency for detecting the rotation angle are used. Under the assumption that the number of components is a predetermined number, the image quality can be improved.
更に、本発明によれば、データを表す非ゼロ周波数成分が挿入される周波数領域における円の半径が不定であるので、データの秘匿性を高めることができる。 Furthermore, according to the present invention, since the radius of the circle in the frequency region where the non-zero frequency component representing the data is inserted is indefinite, the confidentiality of the data can be improved.
以下、図面を参照して本発明を実施するための最良の形態について詳細に説明する。 The best mode for carrying out the present invention will be described below in detail with reference to the drawings.
図1は、本発明の実施形態で利用する周波数領域におけるパターンを示す。図1において、黒丸は非ゼロ周波数成分を示し、白抜きの丸は、ゼロ周波数成分を示す。ここで、非ゼロ周波数成分とは、振幅がゼロでない周波数成分のことであり、ゼロ周波数成分とは、振幅がゼロである周波数成分のことである。 FIG. 1 shows a pattern in the frequency domain used in the embodiment of the present invention. In FIG. 1, black circles indicate non-zero frequency components, and white circles indicate zero frequency components. Here, the non-zero frequency component is a frequency component whose amplitude is not zero, and the zero frequency component is a frequency component whose amplitude is zero.
図1に示すように、最外円は、20等分され、20個の位置に非ゼロ周波数成分又はゼロ周波数成分が配置される。20個の位置の何れもがX軸又はY軸と重複しないように、角度のオフセットを設けて最外円は20等分されている。最内円も、20等分され、20個の位置のうちの10個の位置に非ゼロ周波数成分が配置される。 As shown in FIG. 1, the outermost circle is divided into 20 equal parts, and non-zero frequency components or zero frequency components are arranged at 20 positions. The outermost circle is divided into 20 equal parts by providing an angle offset so that none of the 20 positions overlaps the X axis or the Y axis. The innermost circle is also divided into 20 equal parts, and non-zero frequency components are arranged at 10 positions out of the 20 positions.
図1に示すパターンにおいて、非ゼロ周波数成分又はゼロ周波数成分が配置される位置が、X軸又はY軸と重複しないように選択されているのは、電子透かしを検出する際に、画像をフーリエ変換して得られた周波数領域において、X軸上及びY軸上にノイズに起因する非ゼロ周波数成分が多数発生し、このような非ゼロ周波数成分が電子透かし挿入装置により挿入した非ゼロ周波数成分と区別が付かなくなる可能性が高くなるからである。 In the pattern shown in FIG. 1, the position where the non-zero frequency component or the zero frequency component is arranged is selected so as not to overlap with the X axis or the Y axis. In the frequency domain obtained by conversion, a large number of non-zero frequency components due to noise are generated on the X-axis and the Y-axis, and such non-zero frequency components are inserted by the digital watermark insertion device. This is because there is a high possibility that they cannot be distinguished from each other.
図1に示すパターンは、図2及び3に示すパターンに分解することができる。図2に示すパターンは、データを表すパターンであり、図3に示すパターンは、回転角度を検出するためのパターンである。 The pattern shown in FIG. 1 can be decomposed into the patterns shown in FIGS. The pattern shown in FIG. 2 is a pattern representing data, and the pattern shown in FIG. 3 is a pattern for detecting a rotation angle.
図2を参照すると、データを表すパターンは、最外円を20等分して得られる20個の位置に形成される。20個の位置のうちの10個の位置に非ゼロ周波数成分が配置され、残りの10個の位置にゼロ周波数成分が配置される(非ゼロ周波数成分が配置されない。)。周波数成分の振幅は原点を中心として点対称に分布するので、20個の位置から10個の位置を選択することは、10個の位置から5個の位置を選択することと等価である。従って、データを表すパターンは10C5=252通りのデータを表すことができる。 Referring to FIG. 2, patterns representing data are formed at 20 positions obtained by dividing the outermost circle into 20 equal parts. Non-zero frequency components are arranged at 10 of the 20 positions, and zero frequency components are arranged at the remaining 10 positions (non-zero frequency components are not arranged). Since the amplitude of the frequency component is distributed point-symmetrically around the origin, selecting 10 positions from 20 positions is equivalent to selecting 5 positions from 10 positions. Therefore, the pattern representing data can represent 10 C 5 = 252 kinds of data.
図3を参照すると、回転角度を検出するためのパターンは、最内円を20等分して得られる20個の位置のうちの10個の位置に非ゼロ周波数成分を配置することにより形成される。 Referring to FIG. 3, the pattern for detecting the rotation angle is formed by arranging non-zero frequency components at 10 positions out of 20 positions obtained by dividing the innermost circle into 20 equal parts. The
画像の回転角度を検出するためには、円上に10個の非ゼロ周波数成分が存在することが検出されるまで、円の半径をゼロから徐々に増やし、検出された時の10個の非ゼロ周波数成分の周波数領域における極座標で見た位置の角度を検出する。 In order to detect the rotation angle of the image, the radius of the circle is gradually increased from zero until it is detected that there are 10 non-zero frequency components on the circle, and the 10 non-zero when detected. The angle of the position seen in polar coordinates in the frequency domain of the zero frequency component is detected.
画像の回転角度が検出されたならば、非ゼロ周波数成分又はゼロ周波数成分が存在するはずである20個の位置の周波数領域の極座標で見た角度を知ることができるので、円上のそれらの角度の位置にウィンドウを設け、そのウィンドウの位置に10個の非ゼロ周波数成分があることが検出されるまで、円の半径を最内円の半径から徐々に増やし、検出された時の10個の非ゼロ周波数成分の周波数領域における極座標で見た位置を検出する。そして、それらの位置の組合せ及び回転角度よりデータを検出する。 If the rotation angle of the image is detected, we can know the angle seen in polar coordinates in the frequency domain of the 20 positions where there should be non-zero frequency components or zero frequency components, so those on the circle A window is provided at the angle position, and the radius of the circle is gradually increased from the radius of the innermost circle until it is detected that there are 10 non-zero frequency components at the position of the window. The position seen in the polar coordinate in the frequency domain of the non-zero frequency component of is detected. Data is detected from the combination of the positions and the rotation angle.
回転角度を検出するためのパターンは、図3に示されるようなものには限られない。例えば、最内円を2分割して得た2つの位置に非ゼロ周波数成分を挿入しても回転角度を検出することができる。また、図4に示すように最内円に分布する非ゼロ周波数成分を挿入しても回転角度を検出することができる。ここで、回転角度が0度である場合と回転角度が180度である場合には、周波数成分の振幅の分布は同一であるので、0度の回転角度と180度の回転角度を区別する必要はない。従って、回転角度を検出するためのパターンは、一般には、周波数領域の原点を中心として180度以外の角度を回転させた後の全てのパターン中の位置が回転させる前のパターン中の全ての位置と重なるようなことが、角度を180度以外のいかなる値としても、ないようなものである。 The pattern for detecting the rotation angle is not limited to that shown in FIG. For example, even if a non-zero frequency component is inserted into two positions obtained by dividing the innermost circle into two, the rotation angle can be detected. Further, as shown in FIG. 4, the rotation angle can be detected even when non-zero frequency components distributed in the innermost circle are inserted. Here, when the rotation angle is 0 degree and when the rotation angle is 180 degrees, the amplitude distribution of the frequency components is the same, so it is necessary to distinguish between the 0 degree rotation angle and the 180 degree rotation angle. There is no. Therefore, the pattern for detecting the rotation angle is generally all positions in the pattern before rotating the positions in all patterns after rotating an angle other than 180 degrees around the origin of the frequency domain. The angle does not overlap with any value other than 180 degrees.
また、図4に示すように、図1に示す最外円よりも半径が長い円を設け、その円上の20個の位置のうちの10個の位置に非ゼロ周波数成分を挿入し、残りの10個の位置にゼロ周波数成分を挿入するように(非ゼロ周波数成分を挿入しないように)してもよい。この場合、その円上の非ゼロ周波数成分の位置の組合せの数は、10C5=252であるので、図1に示す最外円上の非ゼロ周波数成分の位置の組合せと掛け合わせて、252×252=63504通りのデータを表すことができる。 Further, as shown in FIG. 4, a circle having a radius longer than the outermost circle shown in FIG. 1 is provided, and non-zero frequency components are inserted into 10 positions out of 20 positions on the circle, and the rest The zero frequency components may be inserted at the ten positions (not to insert non-zero frequency components). In this case, since the number of combinations of the positions of the non-zero frequency components on the circle is 10 C 5 = 252, it is multiplied by the combination of the positions of the non-zero frequency components on the outermost circle shown in FIG. 252 × 252 = 63504 data can be represented.
更に、図1及び4の例では、全ての円を20分割したが、各円毎に分割数が異なるようにしてもよい。また、分割位置の周波数領域の極座標で見た角度を円毎に変えても良い。但し、電子透かし挿入側と電子透かし検出側との間で、分割数及び分割位置の円相互間での相対角度を揃えておく必要がある。 Furthermore, in the example of FIGS. 1 and 4, all the circles are divided into 20 parts, but the number of divisions may be different for each circle. Moreover, you may change the angle seen in the polar coordinate of the frequency domain of a division position for every circle. However, it is necessary to make the relative angle between the number of divisions and the circles of the division positions uniform between the digital watermark insertion side and the digital watermark detection side.
図5は、本発明の実施形態による電子透かし挿入装置の構成を示すブロック図である。 FIG. 5 is a block diagram showing the configuration of the digital watermark insertion apparatus according to the embodiment of the present invention.
図5を参照すると、本発明の実施形態による電子透かし挿入装置は、回転角度検出用パターン生成部101、データ保持部102、データパターン生成部103、合成部104、逆離散フーリエ変換部105、二値化部106及び加算部107を備える。なお、これらの部分はハードウェアによって実現することもできるが、コンピュータをこれらの部分として機能させるためのプログラムをコンピュータが読み込んで実行することによっても実現することができる。
Referring to FIG. 5, the digital watermark insertion apparatus according to the embodiment of the present invention includes a rotation angle detection
次に、電子透かし挿入装置により行われる電子透かし挿入方法について、図5及び6を参照して説明する。 Next, a digital watermark insertion method performed by the digital watermark insertion apparatus will be described with reference to FIGS.
回転角度検出用パターン生成部101は、図3に示すような回転角度検出用パターンを生成する(ステップS151)。データ保持部102は、挿入するべきデータを保持する。データパターン生成部103は、データ保持部102に保持されているデータの値に応じて、図2に示すようなデータパターンを生成する(ステップS152)。合成部104は、回転角度検出用パターン生成部101が生成した回転角度検出用パターン及びデータパターン生成部103が生成したデータパターンを合成して、図1に示すようなパターンを生成する(ステップS153)。逆フーリエ変換部105は、合成部104が生成した図1に示すようなパターンを逆フーリエ変換して、空間領域におけるパターン(空間領域パターン)を生成する(ステップS154)。二値化部106は、逆フーリエ変換部105が生成した空間領域パターンを二値化する(ステップS155)。加算部107は、二値化部106により二値化された空間領域パターンを原画像の黄成分の全部又は一部のブロックに加算する(ステップS156)。
The rotation angle detection
なお、二値化部106を通さない空間領域パターンを原画像の黄成分に加算してもよい。また、加算部107を減算部に置き換え、二値化された空間領域パターン又は二値化されていない空間領域パターンを原画像の黄成分から減算してもよい。
Note that a spatial region pattern that does not pass through the
図7は、本発明の実施形態による電子透かし検出装置の構成を示すブロック図である。 FIG. 7 is a block diagram showing the configuration of the digital watermark detection apparatus according to the embodiment of the present invention.
図7を参照すると、本発明の実施形態による電子透かし検出装置は、空間フィルタ201、ブロック化部202、離散フーリエ変換部203、振幅計算部204、振幅積算部205、回転角度検出部206及びデータ検出部207を備える。なお、これらの部分はハードウェアによって実現することもできるが、コンピュータをこれらの部分として機能させるためのプログラムをコンピュータが読み込んで実行することによっても実現することができる。
Referring to FIG. 7, a digital watermark detection apparatus according to an embodiment of the present invention includes a
次に、図7に示す電子透かし検出装置により行われる電子透かし検出方法について、図7及び8等を参照して説明する。 Next, a digital watermark detection method performed by the digital watermark detection apparatus shown in FIG. 7 will be described with reference to FIGS.
空間フィルタ201は、図9に示すような重みを有することから明らかなように、帯域通過特性を有す。電子透かしが挿入された画像の青成分(電子透かし挿入画像)に、空間フィルタ201が適用される(ステップS251)。これにより、図1に示すパターン以外の非ゼロ周波数成分は抑圧される。ブロック化部202は、空間フィルタが適用された電子透かし挿入画像を1以上のブロックにブロック化する(ステップS252)。離散フーリエ変換部203は、ブロック化部202が出力する各ブロックの電子透かし挿入画像を離散フーリエ変換して、各ブロックの各周波数成分を得る(ステップS253)。振幅計算部204は、各ブロックの各周波数成分の振幅を計算する(ステップS254)。振幅積算部205は、ブロック間で各周波数毎に振幅を積算する(ステップS255)。回転角度検出部206は、振幅積算部205が出力する各周波数成分の振幅に基づいて、図3に示すような最内円の非ゼロ周波数成分を検出し、これに基づいて画像の回転角度を検出する(ステップS256)。この検出の詳細については後述する。データ検出部208は、振幅積算部205が出力する振幅及び回転角度検出部206が出力する回転角度に基づいてデータを検出する(ステップS257)。この検出の詳細についても後述する。
As apparent from the fact that the
図10は、本発明の他の形態による電子透かし検出装置の構成を示すブロック図である。 FIG. 10 is a block diagram showing a configuration of a digital watermark detection apparatus according to another embodiment of the present invention.
図10を参照すると、本発明の実施形態による電子透かし検出装置は、空間フィルタ201、ブロック化部202、離散フーリエ変換部203、周波数成分積算部304、振幅計算部305、回転角度検出部206及びデータ検出部207を備える。なお、これらの部分はハードウェアによって実現することもできるが、コンピュータをこれらの部分として機能させるためのプログラムをコンピュータが読み込んで実行することによっても実現することができる。
Referring to FIG. 10, a digital watermark detection apparatus according to an embodiment of the present invention includes a
次に、図10に示す電子透かし検出装置により行われる電子透かし検出方法について、図10及び11等を参照して説明する。 Next, a digital watermark detection method performed by the digital watermark detection apparatus shown in FIG. 10 will be described with reference to FIGS.
空間フィルタ201、ブロック化部202、離散フーリエ変換部203の動作は、図7に示す電子透かし検出装置のものと同様であるので、説明を省略する。周波数成分積算部304は、離散フーリエ変換部203が出力する各ブロックの周波数成分を複素数のまま積算する(ステップS354)。振幅計算部305は、周波数成分積算部304が出力する複素数の積算値の振幅を計算する(ステップS355)。回転角度検出部206は、振幅計算部305が出力する各周波数成分の振幅に基づいて、図3に示すような最内円の非ゼロ周波数成分を検出し、これに基づいて画像の回転角度を検出する(ステップS256)。この検出の詳細については後述する。データ検出部208は、振幅積算部205が出力する振幅及び回転角度検出部206が出力する回転角度に基づいてデータを検出する(ステップS257)。この検出の詳細についても後述する。
The operations of the
次に、回転角度を検出するステップ(ステップS256)の詳細について、図12を参照して説明する。 Next, details of the step of detecting the rotation angle (step S256) will be described with reference to FIG.
まず、パラメータである半径を周波数領域における最小周波数に対応するゼロに初期化する(ステップS401)。次に、現在の半径の円周上にN1個(図1の例では、10個)の非ゼロ周波数成分が存在するか否かを判断する(ステップS402)。存在しなければ(ステップS402でNO)、半径を周波数領域における離散周波数単位間隔に対応するΔRだけ長くして(ステップS403)、再度、ステップS402を行う。 First, the parameter radius is initialized to zero corresponding to the minimum frequency in the frequency domain (step S401). Next, it is determined whether N1 (10 in the example of FIG. 1) non-zero frequency components exist on the circumference of the current radius (step S402). If it does not exist (NO in step S402), the radius is increased by ΔR corresponding to the discrete frequency unit interval in the frequency domain (step S403), and step S402 is performed again.
ステップS402においては、以下のようなことをする。すなわち、360度を2×N1に分割して得た角度毎に2×N1個のウィンドウを設け、全周波数成分をそのウィンドウを通した後にしきい値を超える振幅の非ゼロ周波数成分がN1個検出できるか否かをウィンドウの角度をずらしながらみていく。検出できた場合には、そのN1個の非ゼロ周波数成分が検出しようとしている非ゼロ周波数成分である。 In step S402, the following is performed. That is, 2 × N1 windows are provided for each angle obtained by dividing 360 degrees into 2 × N1, and N1 nonzero frequency components having amplitudes exceeding the threshold value after passing all frequency components through the window. Check if it can be detected while shifting the window angle. If it can be detected, the N1 non-zero frequency components are non-zero frequency components to be detected.
存在すれば(ステップS402でYES)、現在の半径の円が最内円であることを意味するので、検出されたN1個の非ゼロ周波数成分の分布を基に回転角度を検出する(ステップS404)。 If it exists (YES in step S402), it means that the circle with the current radius is the innermost circle, and therefore the rotation angle is detected based on the distribution of the detected N1 non-zero frequency components (step S404). ).
図12には示していないが、半径がゼロになってもN1個の非ゼロ周波数成分が検出できない場合には、ウィンドウの半径方向の幅を広げてステップS402及びステップS403を繰り返す。それでもどの半径でもN1個の非ゼロ周波数成分が検出できない場合には、ウィンドウの半径方向の幅を戻した上で上記しきい値を下げて、ステップS402及びS403を繰り返す。それでもどの半径でもN1個の非ゼロ周波数成分が検出できない場合には、再度ウィンドウの半径方向の幅を変えてステップS402及びS403を繰り返す。以下同様なことを繰り返す。 Although not shown in FIG. 12, if N1 non-zero frequency components cannot be detected even when the radius becomes zero, steps S402 and S403 are repeated with the width in the radial direction of the window widened. If N1 non-zero frequency components cannot be detected at any radius, the threshold value is lowered after returning the radial width of the window, and steps S402 and S403 are repeated. If N1 non-zero frequency components still cannot be detected at any radius, steps S402 and S403 are repeated again by changing the radial width of the window. The same is repeated thereafter.
次に、データを検出するステップ(ステップS257)の詳細について、図13を参照して説明する。 Next, details of the step of detecting data (step S257) will be described with reference to FIG.
まず、パラメータである半径を回転角度検出用の円である最内円の半径に初期化する(ステップS421)。次に、現在の半径の円周上にN1個(図1の例では、10個)の非ゼロ周波数成分が存在するか否かを判断する(ステップS422)。存在しなければ(ステップS422でNO)、半径を周波数領域における離散周波数単位間隔に対応するΔRだけ長くして(ステップS423)、再度、ステップS422を行う。 First, a radius as a parameter is initialized to a radius of an innermost circle that is a circle for detecting a rotation angle (step S421). Next, it is determined whether or not N1 (10 in the example of FIG. 1) non-zero frequency components exist on the circumference of the current radius (step S422). If it does not exist (NO in step S422), the radius is increased by ΔR corresponding to the discrete frequency unit interval in the frequency domain (step S423), and step S422 is performed again.
ステップS422における検出方法はステップS402における検出方法と同様なものである。また、検出できない場合には上記と同様な再試行を行う。 The detection method in step S422 is the same as the detection method in step S402. If it cannot be detected, a retry similar to the above is performed.
存在すれば(ステップS422でYES)、現在の半径の円が最外円であることを意味するので、検出されたN1個の非ゼロ周波数成分の分布及びステップS256で検出された回転角度を基にデータを検出する(ステップS424)。 If it exists (YES in step S422), it means that the circle with the current radius is the outermost circle, and therefore based on the distribution of the detected N1 non-zero frequency components and the rotation angle detected in step S256. The data is detected (step S424).
なお、図4に示すように、図1に示す最外円の外側に更にデータ円がある場合には、ステップS421乃至S423と同様なことをそのデータ円についても行い、ステップS424では両円の非ゼロ周波数成分の分布及びステップS256で検出された回転角度を基にデータを検出する。 As shown in FIG. 4, when there is a data circle outside the outermost circle shown in FIG. 1, the same processing as in steps S421 to S423 is also performed for the data circle. Data is detected based on the distribution of non-zero frequency components and the rotation angle detected in step S256.
なお、図9に示す空間フィルタの通過帯域は、ナイキスト周波数の半分の周波数を中心とするので、カメラにより撮影した場合に、約80%乃至約130%の範囲で拡大又は縮小されていても、空間フィルタを通した画像においては、電子透かしが強調され、原画像が抑圧される効果を得ることができる。カメラにより撮影した画像が前記範囲よりを超えて拡大又は縮小されている場合は、図4に示す空間フィルタを通した画像においては、電子透かしが強調され、原画像が抑圧される効果を得ることができず、従って、電子透かしを検出できなくなる場合もありうるが、そのような場合には、空間フィルタの通過帯域を変化させて、再度、電子透かしの検出を再試行するようにすればよい。 Note that the passband of the spatial filter shown in FIG. 9 is centered on a half of the Nyquist frequency, so that when captured by a camera, even if it is enlarged or reduced in a range of about 80% to about 130%, In the image that has passed through the spatial filter, the digital watermark is emphasized and the effect of suppressing the original image can be obtained. When the image taken by the camera is enlarged or reduced beyond the above range, the digital watermark is emphasized and the original image is suppressed in the image that has passed through the spatial filter shown in FIG. Therefore, the digital watermark may not be detected, but in such a case, the passband of the spatial filter may be changed and the detection of the digital watermark may be retried again. .
本発明は、印刷物に印刷される画像に電子透かしを挿入し、その画像をカメラ等で撮影し、撮影された画像から電子透かしを検出することにより、電子透かしに対応するデータを取得するために利用することができる。 In order to obtain data corresponding to a digital watermark by inserting a digital watermark into an image printed on a printed material, photographing the image with a camera or the like, and detecting the digital watermark from the captured image Can be used.
101 回転角度検出用パターン生成部
102 データ保持部
103 データパターン生成部
104 合成部
105 逆離散フーリエ変換部
106 二値化部
107 加算部(減算部)
201 空間フィルタ
202 ブロック化部
203 離散フーリエ変換部
204 振幅計算部
205 振幅積算部
206 回転角度検出部
207 データ検出部
304 周波数成分積算部
305 振幅計算部
DESCRIPTION OF
201
Claims (23)
前記周波数領域において、前記第1の円とは異なる第2の円上の位置に画像の回転角度を検出するための非ゼロ周波数成分を生成するステップと、
生成された前記非ゼロ周波数成分を合わせて前記周波数領域から空間領域に変換するステップと、
原画像の全部又は一部のブロックに前記変換により得た画像を加算又は減算するステップと、
を備えることを特徴とする電子透かし挿入方法。 In the frequency domain, the first predetermined number of candidate positions arranged on a first circle with an indefinite radius are positions that are combined with each other by data to be inserted, and the second number smaller than the first predetermined number Generating a non-zero frequency component for a predetermined number;
Generating a non-zero frequency component for detecting a rotation angle of the image at a position on a second circle different from the first circle in the frequency domain;
Combining the generated non-zero frequency components from the frequency domain to the spatial domain;
Adding or subtracting the image obtained by the conversion to all or a part of the blocks of the original image;
An electronic watermark insertion method comprising:
前記第2の円上の位置とは、周波数領域の原点を中心として180度以外の角度を回転させた後の全ての該位置が回転させる前の全ての該位置と重なるようなことが、該角度を180度以外のいかなる値としても、ないようなものであることを特徴とする電子透かし挿入方法。 The digital watermark insertion method according to claim 1,
The positions on the second circle are such that all the positions after rotating at an angle other than 180 degrees around the origin of the frequency domain overlap with all the positions before rotating. An electronic watermark insertion method characterized in that the angle is not set to any value other than 180 degrees.
前記第1の円及び前記第2の円とは異なる第3の円上に並ぶ第3の所定数の候補位置のうちの挿入すべき前記データにより相互に組合わされる位置であって前記第3の所定数より少ない第4の所定数のものに非ゼロ周波数成分を生成するステップを更に備え、
前記変換のステップでは、前記第3の円上に生成された非ゼロ周波数成分も前記周波数領域から前記空間領域に変換することを特徴とする電子透かし挿入方法。 The digital watermark insertion method according to claim 1,
Of the third predetermined number of candidate positions arranged on a third circle different from the first circle and the second circle, the positions are combined with each other by the data to be inserted, and the third circle Generating a non-zero frequency component in a fourth predetermined number less than the predetermined number of
In the converting step, a non-zero frequency component generated on the third circle is also converted from the frequency domain to the spatial domain.
前記1以上の画像ブロックのデータを空間領域から周波数領域に変換するステップと、
前記1以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得るステップと、
前記周波数領域において、前記各周波数成分の振幅を基に、前記画像の回転角度を検出するための非ゼロ周波数成分を検出することにより前記画像の回転角度を検出するステップと、
前記周波数領域において、円周上の第1の所定数の位置であって検出された前記回転角度を回転基準としたもののうちの前記第1の所定数よりも少ない第2の所定数の何れかの位置に非ゼロ周波数成分があるか否かを異なる半径の各円毎に検出するステップと、
前記第2の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せを基に、前記画像に挿入されているデータを検出するステップと、
を備えることを特徴とする電子透かし検出方法。 Obtaining one or more image blocks from the image;
Converting the data of the one or more image blocks from a spatial domain to a frequency domain;
Obtaining the amplitude of each frequency component based on each frequency component of the one or more image blocks;
Detecting the rotation angle of the image in the frequency domain by detecting a non-zero frequency component for detecting the rotation angle of the image based on the amplitude of each frequency component;
Any one of a second predetermined number smaller than the first predetermined number among the first predetermined number of positions on the circumference in the frequency domain and the rotation angle detected as a rotation reference. Detecting for each circle of different radii whether there is a non-zero frequency component at the position of
Detecting data inserted in the image based on a combination of positions of the non-zero frequency components detected to be in any of the second predetermined number of positions;
An electronic watermark detection method comprising:
前記画像の回転角度を検出するための非ゼロ周波数成分は、円上の位置に分布することを特徴とする電子透かし検出方法。 The digital watermark detection method according to claim 4,
The digital watermark detection method according to claim 1, wherein the non-zero frequency components for detecting the rotation angle of the image are distributed at positions on a circle.
前記円上の位置とは、周波数領域の原点を中心として180度以外の角度を回転させた後の全ての該位置が回転させる前の全ての該位置と重なるようなことが、該角度を180度以外のいかなる値としても、ないようなものであることを特徴とする電子透かし検出方法。 The digital watermark detection method according to claim 5, wherein
The position on the circle means that all the positions after rotating an angle other than 180 degrees around the origin of the frequency domain overlap with all the positions before rotating. A digital watermark detection method characterized in that it does not have any value other than degrees.
データに係る前記非ゼロ周波数成分があるか否かを異なる半径毎に検出するステップは、半径が短い円から半径が長い円に向かって行われ、
データに係る前記非ゼロ周波数成分があるか否かを異なる半径毎に検出するステップにより、データに係る前記非ゼロ周波数成分が検出された後に、前記周波数領域において、円周上の第3の位置であって検出された前記回転画像を回転基準としたもののうちの前記第3の所定数よりも少ない第4の所定数の何れかの位置に非ゼロ周波数成分があるか否かを、前記データに係る非ゼロ周波数成分が検出された円よりも半径が長い各円毎に検出するステップを更に備え、
前記画像に挿入されているデータを検出するステップでは、前記第2の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せのみならず、前記第4の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せも基にして、前記画像に挿入されているデータを検出することを特徴とする電子透かし検出方法。 The digital watermark detection method according to claim 4,
The step of detecting whether there is the non-zero frequency component related to the data for each different radius is performed from a circle with a short radius toward a circle with a long radius,
A third position on the circumference in the frequency domain after the non-zero frequency component related to the data is detected by detecting whether the non-zero frequency component related to the data exists for each different radius. Whether or not there is a non-zero frequency component at any position of the fourth predetermined number smaller than the third predetermined number of the detected rotation images as a rotation reference. Detecting for each circle having a longer radius than the circle in which the non-zero frequency component is detected,
In the step of detecting data inserted in the image, not only the combination of the positions of the non-zero frequency components detected to be in any of the second predetermined number of positions, but also the fourth predetermined An electronic watermark detection method, comprising: detecting data inserted in the image based on a combination of positions of the non-zero frequency components detected to be in any position of the number.
データに係る前記非ゼロ周波数成分があるか否かを異なる半径毎に検出するステップは、半径が長い円から半径が短い円に向かって行われ、
データに係る前記非ゼロ周波数成分があるか否かを異なる半径毎に検出するステップにより、データに係る前記非ゼロ周波数成分が検出された後に、前記周波数領域において、円周上の第3の位置であって検出された前記回転画像を回転基準としたもののうちの前記第3の所定数よりも少ない第4の所定数の何れかの位置に非ゼロ周波数成分があるか否かを、前記データに係る非ゼロ周波数成分が検出された円よりも半径が短い各円毎に検出するステップを更に備え、
前記画像に挿入されているデータを検出するステップでは、前記第2の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せのみならず、前記第4の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せも基にして、前記画像に挿入されているデータを検出することを特徴とする電子透かし検出方法。 The digital watermark detection method according to claim 4,
The step of detecting whether there is the non-zero frequency component related to the data for each different radius is performed from a circle with a long radius toward a circle with a short radius,
A third position on the circumference in the frequency domain after the non-zero frequency component related to the data is detected by detecting whether the non-zero frequency component related to the data exists for each different radius. Whether or not there is a non-zero frequency component at any position of the fourth predetermined number smaller than the third predetermined number of the detected rotation images as a rotation reference. Detecting for each circle having a shorter radius than the circle in which the non-zero frequency component is detected,
In the step of detecting data inserted in the image, not only the combination of the positions of the non-zero frequency components detected to be in any of the second predetermined number of positions, but also the fourth predetermined An electronic watermark detection method, comprising: detecting data inserted in the image based on a combination of positions of the non-zero frequency components detected to be in any position of the number.
前記画像を空間フィルタに通すステップを更に備え、
前記画像より1以上の画像ブロックを得るステップは、前記空間フィルタを通った画像より1以上の画像ブロックを得るステップであることを特徴とする電子透かし検出方法。 The digital watermark detection method according to claim 4,
Further comprising passing the image through a spatial filter;
The step of obtaining one or more image blocks from the image is a step of obtaining one or more image blocks from an image that has passed through the spatial filter.
前記周波数領域において、前記第1の円とは異なる第2の円上の位置に画像の回転角度を検出するための非ゼロ周波数成分を生成するステップと、
生成された前記非ゼロ周波数成分を合わせて前記周波数領域から空間領域に変換するステップと、
原画像の全部又は一部のブロックに前記変換により得た画像を加算又は減算するステップと、
前記加算又は減算を経た画像より1以上の画像ブロックを得るステップと、
前記1以上の画像ブロックのデータを空間領域から周波数領域に変換するステップと、
前記1以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得るステップと、
前記周波数領域において、前記各周波数成分の振幅を基に、前記画像の回転角度を検出するための非ゼロ周波数成分を検出することにより前記画像の回転角度を検出するステップと、
前記周波数領域において、円周上の第1の所定数の位置であって検出された前記回転角度を回転基準としたもののうちの前記第1の所定数よりも少ない第2の所定数の何れかの位置に非ゼロ周波数成分があるか否かを異なる半径の各円毎に検出するステップと、
前記第2の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せを基に、前記画像に挿入されているデータを検出するステップと、
を備えることを特徴とするデータ伝送方法。 In the frequency domain, the first predetermined number of candidate positions arranged on a first circle with an indefinite radius are positions that are combined with each other by data to be inserted, and the second number smaller than the first predetermined number Generating a non-zero frequency component for a predetermined number;
Generating a non-zero frequency component for detecting a rotation angle of the image at a position on a second circle different from the first circle in the frequency domain;
Combining the generated non-zero frequency components from the frequency domain to the spatial domain;
Adding or subtracting the image obtained by the conversion to all or a part of the blocks of the original image;
Obtaining one or more image blocks from the image that has undergone the addition or subtraction;
Converting the data of the one or more image blocks from a spatial domain to a frequency domain;
Obtaining the amplitude of each frequency component based on each frequency component of the one or more image blocks;
Detecting the rotation angle of the image in the frequency domain by detecting a non-zero frequency component for detecting the rotation angle of the image based on the amplitude of each frequency component;
Any one of a second predetermined number smaller than the first predetermined number among the first predetermined number of positions on the circumference in the frequency domain and using the detected rotation angle as a rotation reference. Detecting for each circle of different radii whether there is a non-zero frequency component at the position of
Detecting data inserted in the image based on a combination of positions of the non-zero frequency components detected to be in any of the second predetermined number of positions;
A data transmission method comprising:
前記周波数領域において、前記第1の円とは異なる第2の円上の位置に画像の回転角度を検出するための非ゼロ周波数成分を生成する手段と、
生成された前記非ゼロ周波数成分を合わせて前記周波数領域から空間領域に変換する手段と、
原画像の全部又は一部のブロックに前記変換により得た画像を加算又は減算する手段と、
を備えることを特徴とする電子透かし挿入装置。 In the frequency domain, the first predetermined number of candidate positions arranged on a first circle having an indefinite radius are positions that are combined with each other by data to be inserted, and are less than the first predetermined number. Means for generating a non-zero frequency component in a predetermined number;
Means for generating a non-zero frequency component for detecting a rotation angle of the image at a position on a second circle different from the first circle in the frequency domain;
Means for combining the generated non-zero frequency components into a spatial domain from the frequency domain;
Means for adding or subtracting the image obtained by the conversion to all or a part of the blocks of the original image;
An electronic watermark insertion apparatus comprising:
前記第2の円上の位置とは、周波数領域の原点を中心として180度以外の角度を回転させた後の全ての該位置が回転させる前の全ての該位置と重なるようなことが、該角度を180度以外のいかなる値としても、ないようなものであることを特徴とする電子透かし挿入装置。 The digital watermark insertion apparatus according to claim 11, wherein
The positions on the second circle are such that all the positions after rotating at an angle other than 180 degrees around the origin of the frequency domain overlap with all the positions before rotating. An electronic watermark insertion apparatus characterized in that the angle is not set to any value other than 180 degrees.
前記第1の円及び前記第2の円とは異なる第3の円上に並ぶ第3の所定数の候補位置のうちの挿入すべき前記データにより相互に組合わされる位置であって前記第3の所定数より少ない第4の所定数のものに非ゼロ周波数成分を生成する手段を更に備え、
前記変換の手段では、前記第3の円上に生成された非ゼロ周波数成分も前記周波数領域から前記空間領域に変換することを特徴とする電子透かし挿入装置。 The digital watermark insertion apparatus according to claim 11, wherein
Of the third predetermined number of candidate positions arranged on a third circle different from the first circle and the second circle, the positions are combined with each other by the data to be inserted, and the third circle Means for generating a non-zero frequency component in a fourth predetermined number less than the predetermined number of
The digital watermark insertion apparatus characterized in that the converting means converts a non-zero frequency component generated on the third circle from the frequency domain to the spatial domain.
前記1以上の画像ブロックのデータを空間領域から周波数領域に変換する手段と、
前記1以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得る手段と、
前記周波数領域において、前記各周波数成分の振幅を基に、前記画像の回転角度を検出するための非ゼロ周波数成分を検出することにより前記画像の回転角度を検出する手段と、
前記周波数領域において、円周上の第1の所定数の位置であって検出された前記回転角度を回転基準としたもののうちの前記第1の所定数よりも少ない第2の所定数の何れかの位置に非ゼロ周波数成分があるか否かを異なる半径の各円毎に検出する手段と、
前記第2の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せを基に、前記画像に挿入されているデータを検出する手段と、
を備えることを特徴とする電子透かし検出装置。 Means for obtaining one or more image blocks from the image;
Means for transforming the data of the one or more image blocks from a spatial domain to a frequency domain;
Means for obtaining the amplitude of each frequency component based on each frequency component of the one or more image blocks;
Means for detecting the rotation angle of the image in the frequency domain by detecting a non-zero frequency component for detecting the rotation angle of the image based on the amplitude of each frequency component;
Any one of a second predetermined number smaller than the first predetermined number among the first predetermined number of positions on the circumference in the frequency domain and the rotation angle detected as a rotation reference. Means for detecting for each circle of different radii whether there is a non-zero frequency component at the position of
Means for detecting data inserted in the image based on a combination of positions of the non-zero frequency components detected to be in any of the second predetermined number of positions;
An electronic watermark detection apparatus comprising:
前記画像の回転角度を検出するための非ゼロ周波数成分は、円上の位置に分布することを特徴とする電子透かし検出装置。 The digital watermark detection apparatus according to claim 14,
The digital watermark detection apparatus, wherein the non-zero frequency component for detecting the rotation angle of the image is distributed at positions on a circle.
前記円上の位置とは、周波数領域の原点を中心として180度以外の角度を回転させた後の全ての該位置が回転させる前の全ての該位置と重なるようなことが、該角度を180度以外のいかなる値としても、ないようなものであることを特徴とする電子透かし検出装置。 The digital watermark detection apparatus according to claim 15, wherein
The position on the circle means that all the positions after rotating an angle other than 180 degrees around the origin of the frequency domain overlap with all the positions before rotating. An electronic watermark detection apparatus characterized by having no value other than degrees.
データに係る前記非ゼロ周波数成分があるか否かを異なる半径毎に検出する手段は、半径が短い円から半径が長い円に向かって行われ、
データに係る前記非ゼロ周波数成分があるか否かを異なる半径毎に検出する手段により、データに係る前記非ゼロ周波数成分が検出された後に、前記周波数領域において、円周上の第3の位置であって検出された前記回転画像を回転基準としたもののうちの前記第3の所定数よりも少ない第4の所定数の何れかの位置に非ゼロ周波数成分があるか否かを、前記データに係る非ゼロ周波数成分が検出された円よりも半径が長い各円毎に検出する手段を更に備え、
前記画像に挿入されているデータを検出する手段では、前記第2の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せのみならず、前記第4の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せも基にして、前記画像に挿入されているデータを検出することを特徴とする電子透かし検出装置。 The digital watermark detection apparatus according to claim 14,
The means for detecting whether there is the non-zero frequency component related to the data for each different radius is performed from a circle with a short radius toward a circle with a long radius,
A third position on the circumference in the frequency domain after the non-zero frequency component related to the data is detected by means for detecting for each different radius whether or not the non-zero frequency component related to the data exists. Whether or not there is a non-zero frequency component at any position of the fourth predetermined number smaller than the third predetermined number of the detected rotation images as a rotation reference. Further comprising means for detecting each circle having a longer radius than the circle in which the non-zero frequency component is detected,
In the means for detecting data inserted in the image, not only the combination of the positions of the non-zero frequency components detected to be in any of the second predetermined number of positions, but also the fourth predetermined An electronic watermark detection apparatus for detecting data inserted in the image based on a combination of positions of the non-zero frequency components detected to be in any position of the number.
データに係る前記非ゼロ周波数成分があるか否かを異なる半径毎に検出する手段は、半径が長い円から半径が短い円に向かって行われ、
データに係る前記非ゼロ周波数成分があるか否かを異なる半径毎に検出する手段により、データに係る前記非ゼロ周波数成分が検出された後に、前記周波数領域において、円周上の第3の位置であって検出された前記回転画像を回転基準としたもののうちの前記第3の所定数よりも少ない第4の所定数の何れかの位置に非ゼロ周波数成分があるか否かを、前記データに係る非ゼロ周波数成分が検出された円よりも半径が短い各円毎に検出する手段を更に備え、
前記画像に挿入されているデータを検出する手段では、前記第2の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せのみならず、前記第4の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せも基にして、前記画像に挿入されているデータを検出することを特徴とする電子透かし検出装置。 The digital watermark detection apparatus according to claim 14,
The means for detecting whether there is the non-zero frequency component related to the data for each different radius is performed from a circle with a long radius toward a circle with a short radius,
A third position on the circumference in the frequency domain after the non-zero frequency component related to the data is detected by means for detecting for each different radius whether or not the non-zero frequency component related to the data exists. Whether or not there is a non-zero frequency component at any position of the fourth predetermined number smaller than the third predetermined number of the detected rotation images as a rotation reference. Further comprising means for detecting each circle having a shorter radius than the circle in which the non-zero frequency component is detected,
In the means for detecting data inserted in the image, not only the combination of the positions of the non-zero frequency components detected to be in any of the second predetermined number of positions, but also the fourth predetermined An electronic watermark detection apparatus for detecting data inserted in the image based on a combination of positions of the non-zero frequency components detected to be in any position of the number.
前記画像を空間フィルタに通す手段を更に備え、
前記画像より1以上の画像ブロックを得る手段は、前記空間フィルタを通った画像より1以上の画像ブロックを得る手段であることを特徴とする電子透かし検出装置 The digital watermark detection apparatus according to claim 4,
Further comprising means for passing the image through a spatial filter;
The digital watermark detection apparatus characterized in that the means for obtaining one or more image blocks from the image is means for obtaining one or more image blocks from the image that has passed through the spatial filter.
前記周波数領域において、前記第1の円とは異なる第2の円上の位置に画像の回転角度を検出するための非ゼロ周波数成分を生成する手段と、
生成された前記非ゼロ周波数成分を合わせて前記周波数領域から空間領域に変換する手段と、
原画像の全部又は一部のブロックに前記変換により得た画像を加算又は減算する手段と、
前記加算又は減算を経た画像より1以上の画像ブロックを得る手段と、
前記1以上の画像ブロックのデータを空間領域から周波数領域に変換する手段と、
前記1以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得る手段と、
前記周波数領域において、前記各周波数成分の振幅を基に、前記画像の回転角度を検出するための非ゼロ周波数成分を検出することにより前記画像の回転角度を検出する手段と、
前記周波数領域において、円周上の第1の所定数の位置であって検出された前記回転角度を回転基準としたもののうちの前記第1の所定数よりも少ない第2の所定数の何れかの位置に非ゼロ周波数成分があるか否かを異なる半径の各円毎に検出する手段と、
前記第2の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せを基に、前記画像に挿入されているデータを検出する手段と、
を備えることを特徴とするデータ伝送装置。 In the frequency domain, the first predetermined number of candidate positions arranged on a first circle having an indefinite radius are positions that are combined with each other by data to be inserted, and the second predetermined number is less than the first predetermined number. Means for generating a non-zero frequency component in a predetermined number;
Means for generating a non-zero frequency component for detecting a rotation angle of the image at a position on a second circle different from the first circle in the frequency domain;
Means for combining the generated non-zero frequency components into a spatial domain from the frequency domain;
Means for adding or subtracting the image obtained by the conversion to all or a part of the blocks of the original image;
Means for obtaining one or more image blocks from the image that has undergone the addition or subtraction;
Means for transforming the data of the one or more image blocks from a spatial domain to a frequency domain;
Means for obtaining the amplitude of each frequency component based on each frequency component of the one or more image blocks;
Means for detecting the rotation angle of the image in the frequency domain by detecting a non-zero frequency component for detecting the rotation angle of the image based on the amplitude of each frequency component;
Any one of a second predetermined number smaller than the first predetermined number among the first predetermined number of positions on the circumference in the frequency domain and using the detected rotation angle as a rotation reference. Means for detecting for each circle of different radii whether there is a non-zero frequency component at the position of
Means for detecting data inserted into the image based on a combination of positions of the non-zero frequency components detected to be in any of the second predetermined number of positions;
A data transmission device comprising:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2004191408A JP4172587B2 (en) | 2004-06-29 | 2004-06-29 | Digital watermark insertion method, digital watermark insertion device, digital watermark insertion program, digital watermark detection method, digital watermark detection device, digital watermark detection program, data transmission method, data transmission system, and data transmission program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2004191408A JP4172587B2 (en) | 2004-06-29 | 2004-06-29 | Digital watermark insertion method, digital watermark insertion device, digital watermark insertion program, digital watermark detection method, digital watermark detection device, digital watermark detection program, data transmission method, data transmission system, and data transmission program |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2006014171A true JP2006014171A (en) | 2006-01-12 |
JP4172587B2 JP4172587B2 (en) | 2008-10-29 |
Family
ID=35780818
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2004191408A Expired - Fee Related JP4172587B2 (en) | 2004-06-29 | 2004-06-29 | Digital watermark insertion method, digital watermark insertion device, digital watermark insertion program, digital watermark detection method, digital watermark detection device, digital watermark detection program, data transmission method, data transmission system, and data transmission program |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4172587B2 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011150519A (en) * | 2010-01-21 | 2011-08-04 | Kyodo Printing Co Ltd | Rotation angle detector |
-
2004
- 2004-06-29 JP JP2004191408A patent/JP4172587B2/en not_active Expired - Fee Related
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011150519A (en) * | 2010-01-21 | 2011-08-04 | Kyodo Printing Co Ltd | Rotation angle detector |
Also Published As
Publication number | Publication date |
---|---|
JP4172587B2 (en) | 2008-10-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kutter et al. | Fair benchmark for image watermarking systems | |
US8890950B2 (en) | Methods and arrangements for processing image data | |
Poljicak et al. | Discrete Fourier transform–based watermarking method with an optimal implementation radius | |
US8694049B2 (en) | Fast signal detection and distributed computing in portable computing devices | |
US7970166B2 (en) | Steganographic encoding methods and apparatus | |
US20040001608A1 (en) | Image processor and image processing method | |
Kutter et al. | Fair evaluation methods for image watermarking systems | |
JP4503679B2 (en) | Document protection methods | |
US9697578B2 (en) | Assessment of camera phone distortion for digital watermarking | |
JP2002247344A (en) | Method and apparatus for processing image | |
Huang et al. | Attack LSB matching steganography by counting alteration rate of the number of neighbourhood gray levels | |
JP2003032473A (en) | Embedding and detection of electronic watermark | |
JP4199540B2 (en) | Halftone watermarking and related applications | |
Li et al. | Wavelet-based fragile watermarking scheme for image authentication | |
JP4204263B2 (en) | Digital watermark insertion method and apparatus, and digital watermark detection method and apparatus | |
JP4172587B2 (en) | Digital watermark insertion method, digital watermark insertion device, digital watermark insertion program, digital watermark detection method, digital watermark detection device, digital watermark detection program, data transmission method, data transmission system, and data transmission program | |
JP4172588B2 (en) | Digital watermark insertion method, digital watermark insertion apparatus, digital watermark insertion program, digital watermark detection method, digital watermark detection apparatus, and digital watermark detection program | |
JP4270504B2 (en) | Digital watermark insertion method, digital watermark insertion apparatus, digital watermark insertion program, digital watermark detection method, digital watermark detection apparatus, and digital watermark detection program | |
JP6920719B2 (en) | Digital watermarking device and method | |
JP2005123974A (en) | System and method for watermark detection system and method for providing information | |
JP4541213B2 (en) | Digital watermark insertion method, digital watermark detection method, digital watermark insertion device, and digital watermark detection device | |
JP2005039603A (en) | Method and system for electronic watermark detection, information providing system, and information providing method | |
JP2005117406A (en) | System and method for detecting electronic watermark, and system and method for providing information | |
D’Angelo et al. | Watermark-based authentication | |
JP2005295458A (en) | Information detecting apparatus, information embedding apparatus, image processing system, information detecting method, and information embedding method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20070605 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080513 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080709 |
|
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: 20080806 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20080806 |
|
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: 20110822 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110822 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120822 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120822 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130822 Year of fee payment: 5 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |