CN102223561B - Blind watermark embedding and extracting method of stereoscopic video image - Google Patents

Blind watermark embedding and extracting method of stereoscopic video image Download PDF

Info

Publication number
CN102223561B
CN102223561B CN 201110183636 CN201110183636A CN102223561B CN 102223561 B CN102223561 B CN 102223561B CN 201110183636 CN201110183636 CN 201110183636 CN 201110183636 A CN201110183636 A CN 201110183636A CN 102223561 B CN102223561 B CN 102223561B
Authority
CN
China
Prior art keywords
image block
stereoscopic video
video images
coefficient
discrete cosine
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.)
Expired - Fee Related
Application number
CN 201110183636
Other languages
Chinese (zh)
Other versions
CN102223561A (en
Inventor
郁梅
吴爱红
蒋刚毅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Ningbo University
Original Assignee
Ningbo University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Ningbo University filed Critical Ningbo University
Priority to CN 201110183636 priority Critical patent/CN102223561B/en
Publication of CN102223561A publication Critical patent/CN102223561A/en
Application granted granted Critical
Publication of CN102223561B publication Critical patent/CN102223561B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Editing Of Facsimile Originals (AREA)
  • Image Processing (AREA)
  • Compression Or Coding Systems Of Tv Signals (AREA)

Abstract

The invention discloses a blind watermark embedding and extracting method of a stereoscopic video image; the stereoscopic video image is divided into image blocks which are not overlapped mutually, discrete cosine transformation and two-level discrete wavelet transformation are carried out to all image blocks, the relation of numerical values of the discrete cosine transformation direct-current coefficients and two-level discrete wavelet transformation low-frequency coefficients of a left viewpoint view and a right viewpoint view of the stereoscopic video image is analyzed, so as to define the internal relationship of the image blocks and interblock relationship of two image blocks in the stereoscopic video image; a watermark embedding mode is determined though the relationships and a way that whether the watermark information is consistent, wherein the data of the primary stereoscopic video image is not changed when relational embeddedness is carried out, the discrete cosine transformation direct-current coefficient of the image block with small correction magnetite of the discrete cosine transformation direct-current coefficient is corrected when quantization embeddedness is carried out, in this way, the quality of the stereoscopic video image is ensured, when watermark is extracted, no primary stereoscopic video image is needed, so as to realize blind detection.

Description

A kind of blind watermark embedding and extracting method of stereoscopic video images
Technical field
The present invention relates to a kind of digital watermark technology, especially relate to a kind of blind watermark embedding and extracting method of stereoscopic video images.
Background technology
Stereoscopic video images represents is a certain scene or the reference object synchronization information in different angles, and it has increased the depth information of scenery, has strengthened the presence and the sense true to nature of vision.Therefore, the three-dimensional video-frequency system all has broad application prospects in the fields such as advertisement of stereoscopic TV, 3 D video meeting, tele-medicine and special-effect.Along with the appearance of three-dimensional film, the bootlegging of some high-end works and propagation make intellectual property be on the hazard.Therefore, the copyright protection measure towards stereoscopic video images is very necessary.
Digital watermark technology can be well embeds certain secret information in digital products such as digital picture, Voice ﹠ Video, so that the true and reliable property of the copyright of protection digital product, proof digital product, follow the tracks of copy right piracy or the additional information of product is provided.According to different application, digital watermarking can be divided into robust watermarking, fragile watermark and semi-fragile watermarking.Robust watermarking is used for aspects such as copyright protection, and fragile watermark and semi-fragile watermarking are used for the integrity check of video or multi-medium data.Yet existing digital watermark technology is primarily aimed at two dimensional image, becomes the problem that needs are broken through and how digital watermark technology is applied to stereoscopic video images.
Summary of the invention
Technical problem to be solved by this invention provides and is a kind ofly guaranteeing that the stereoscopic video images quality remains unchanged substantially or changes under the little prerequisite, can protect the blind watermark embedding and the extracting method of image copyright effectively.
The present invention solves the problems of the technologies described above the technical scheme that is adopted: a kind of blind watermark embedding method of stereoscopic video images is characterized in that may further comprise the steps:
1. at the watermark built-in end, make S OrgFor being of a size of the original stereoscopic video images of M * N, original stereoscopic video images S OrgComprise left visual point image and right visual point image, remember original stereoscopic video images S OrgLeft visual point image be L Org, remember original stereoscopic video images S OrgRight visual point image be R Org, with original stereoscopic video images S OrgLeft visual point image L OrgWith original stereoscopic video images S OrgRight visual point image R OrgBe divided into respectively (M * N)/(n * n) individual non-overlapping copies is of a size of the image block of n * n, with original stereoscopic video images S OrgLeft visual point image L OrgIn i image block be designated as
Figure BDA0000073150090000021
With original stereoscopic video images S OrgRight visual point image R OrgIn i image block be designated as
Figure BDA0000073150090000022
Original stereoscopic video images S OrgLeft visual point image L OrgIn i image block
Figure BDA0000073150090000023
With original stereoscopic video images S OrgRight visual point image R OrgIn i image block
Figure BDA0000073150090000024
Constitute original stereoscopic video images S OrgIn i image block right, be designated as
Figure BDA0000073150090000025
Wherein, M represents original stereoscopic video images S OrgWide, N represents original stereoscopic video images S OrgHeight, 1≤i≤((M * N)/(n * n));
2. make that W is the binary bitmap that is of a size of m * m, binary bitmap W is carried out the scramble conversion, obtain the binary bitmap behind the scramble, be designated as
Figure BDA0000073150090000026
Wherein, m * m≤((M * N)/(n * n));
3. successively to original stereoscopic video images S OrgIn each image block to adding watermark, its detailed process is as follows:
3.-1, the current original stereoscopic video images S that is handling of definition OrgIn i image block right
Figure BDA0000073150090000027
For current image block right;
3.-2, right to current image block
Figure BDA0000073150090000028
In
Figure BDA0000073150090000029
Carry out discrete cosine transform and secondary wavelet transform respectively, right to current image block
Figure BDA00000731500900000210
In
Figure BDA00000731500900000211
Carry out discrete cosine transform and secondary wavelet transform respectively;
3.-3, make Intra-L represent
Figure BDA00000731500900000212
Internal relations, judge
Figure BDA00000731500900000213
The discrete cosine transform DC coefficient whether greater than
Figure BDA00000731500900000214
The low frequency coefficient of desired location of secondary wavelet transform, if then put
Figure BDA00000731500900000215
Internal relations Intra-L be " 1 ", otherwise, put Internal relations Intra-L be " 0 "; Make Intra-R represent
Figure BDA00000731500900000217
Internal relations, judge
Figure BDA00000731500900000218
The discrete cosine transform DC coefficient whether greater than
Figure BDA00000731500900000219
The low frequency coefficient of desired location of secondary wavelet transform, if then put
Figure BDA00000731500900000220
Internal relations Intra-R be " 1 ", otherwise, put
Figure BDA00000731500900000221
Internal relations Intra-R be " 0 "; Make Intra represent
Figure BDA00000731500900000222
Internal relations, judge
Figure BDA00000731500900000223
Internal relations Intra-L and
Figure BDA00000731500900000224
Internal relations Intra-R whether be " 1 ", if then put
Figure BDA00000731500900000225
Internal relations Intra be " 1 ", otherwise, put
Figure BDA00000731500900000226
Internal relations Intra be " 0 "; Make Inter represent
Figure BDA00000731500900000227
Interblock relation, judge
Figure BDA00000731500900000228
The discrete cosine transform DC coefficient whether greater than
Figure BDA00000731500900000229
The discrete cosine transform DC coefficient, if then put
Figure BDA00000731500900000230
Interblock relations I nter be " 1 ", otherwise, put
Figure BDA00000731500900000231
Interblock relations I nter be " 0 ";
Above-mentioned The secondary wavelet transform desired location low frequency coefficient with
Figure BDA00000731500900000233
The low frequency coefficient of desired location of secondary wavelet transform be the low frequency coefficient of same position of the secondary wavelet transform of image block separately;
3.-4, judge that current image block is right
Figure BDA0000073150090000031
The value of internal relations Intra and the binary bitmap behind the scramble
Figure BDA0000073150090000032
In the pixel value w of i binarized pixel point iWhether identical, if identical, then do not revise
Figure BDA0000073150090000033
With
Figure BDA0000073150090000034
The discrete cosine transform DC coefficient, and the mark current image block is right
Figure BDA0000073150090000035
The watermark embedded mode be 00, execution in step is 3.-7 then; Otherwise execution in step 3.-5;
3.-5, judge that current image block is right
Figure BDA0000073150090000036
The value of interblock relations I nter and the binary bitmap behind the scramble
Figure BDA0000073150090000037
In the pixel value w of i binarized pixel point iWhether identical, if identical, then do not revise
Figure BDA0000073150090000038
With The discrete cosine transform DC coefficient, and the mark current image block is right
Figure BDA00000731500900000310
The watermark embedded mode be 01, execution in step is 3.-7 then; Otherwise execution in step 3.-6;
3.-6, right
Figure BDA00000731500900000311
With
Figure BDA00000731500900000312
The discrete cosine transform DC coefficient quantize respectively according to quantization step S, obtain quantized value Q respectively L, iAnd Q R, i, Q L , i = floor ( DC L , i S ) , Q R , i = floor ( DC R , i S ) , Wherein, DC L, iExpression
Figure BDA00000731500900000315
The discrete cosine transform DC coefficient, DC R, iExpression The discrete cosine transform DC coefficient, floor (x) is for calculating greatest integer function smaller or equal to x; Definition
Figure BDA00000731500900000317
The parity of discrete cosine transform DC coefficient be its quantized value Q L, iParity, the definition
Figure BDA00000731500900000318
The parity of discrete cosine transform DC coefficient be its quantized value Q R, iParity, judge quantized value Q L, iWith quantized value Q R, iWhether be all even number or be all odd number, if then think
Figure BDA00000731500900000319
With The parity unanimity of discrete cosine transform DC coefficient, otherwise, think
Figure BDA00000731500900000321
With
Figure BDA00000731500900000322
The parity of discrete cosine transform DC coefficient inconsistent;
Binary bitmap behind the judgement scramble
Figure BDA00000731500900000323
In the pixel value w of i binarized pixel point iWhether be " 1 ", if it is right then to revise current image block
Figure BDA00000731500900000324
In the discrete cosine transform DC coefficient of any image block, make amended With
Figure BDA00000731500900000326
The parity unanimity of discrete cosine transform DC coefficient, and the mark current image block is right The watermark embedded mode be 10, otherwise it is right to revise current image block
Figure BDA00000731500900000328
In the discrete cosine transform DC coefficient of any image block, make amended
Figure BDA00000731500900000329
With
Figure BDA00000731500900000330
The parity of discrete cosine transform DC coefficient inconsistent, and the mark current image block is right
Figure BDA00000731500900000331
The watermark embedded mode be 10;
3.-7, judge whether i≤m * m sets up, if establishment then adds 1 with the i value, with original stereoscopic video images S OrgIn next image block to right as current image block, and repeated execution of steps 3.-2 to step 3.-7, until handling m * m image block to finishing watermark embed process;
4., the size with binary bitmap, watermark embedded mode and quantization step S be as key, the watermark built-in end extracts end with cipher key delivery to watermark.
Described step 3.-3 in the low frequency coefficient of desired location be the low frequency coefficient of arbitrary position of the secondary wavelet transform of image block.
Described step 3.-6 in that to revise current image block right
Figure BDA0000073150090000041
In the discrete cosine transform DC coefficient of any image block the time, select
Figure BDA0000073150090000042
With
Figure BDA0000073150090000043
Middle discrete cosine transform DC coefficient is revised the discrete cosine transform DC coefficient of that little image block of amplitude and is made amendment.
A kind of blind watermark extracting method that uses the watermarked stereoscopic video images of above-mentioned blind watermark embedding method is characterized in that may further comprise the steps:
1. at the watermark extracting end, make S WatFor being of a size of the stereoscopic video images to be detected of M * N, stereoscopic video images S to be detected WatComprise left visual point image and right visual point image, remember stereoscopic video images S to be detected WatLeft visual point image be L Wat, remember stereoscopic video images S to be detected WatRight visual point image be R Wat, with stereoscopic video images S to be detected WatLeft visual point image L WatWith stereoscopic video images S to be detected WatRight visual point image R WatBe divided into respectively (M * N)/(n * n) individual non-overlapping copies is of a size of the image block of n * n, with stereoscopic video images S to be detected WatLeft visual point image L WatIn j image block be designated as
Figure BDA0000073150090000044
With stereoscopic video images S to be detected WatRight visual point image R WatIn j image block be designated as
Figure BDA0000073150090000045
Stereoscopic video images S to be detected WatLeft visual point image L WatIn j image block
Figure BDA0000073150090000046
With stereoscopic video images S to be detected WatRight visual point image R WatIn j image block
Figure BDA0000073150090000047
Constitute stereoscopic video images S to be detected WatIn j image block right, be designated as
Figure BDA0000073150090000048
Wherein, M represents stereoscopic video images S to be detected WatWide, N represents stereoscopic video images S to be detected WatHeight, 1≤j≤((M * N)/(n * n));
2. successively to stereoscopic video images S to be detected WatIn each image block to extracting watermark, its detailed process is as follows:
2.-1, the current stereoscopic video images S to be detected that is handling of definition WatIn j image block right
Figure BDA0000073150090000049
For current image block right;
2.-2, right to current image block
Figure BDA00000731500900000410
In
Figure BDA00000731500900000411
Carry out discrete cosine transform and secondary wavelet transform respectively, right to current image block
Figure BDA00000731500900000412
In
Figure BDA00000731500900000413
Carry out discrete cosine transform and secondary wavelet transform respectively;
2.-3, make Intra-L ' expression
Figure BDA00000731500900000414
Internal relations, judge
Figure BDA00000731500900000415
The discrete cosine transform DC coefficient whether greater than
Figure BDA00000731500900000416
The low frequency coefficient of desired location of secondary wavelet transform, if then put
Figure BDA00000731500900000417
Internal relations Intra-L ' be " 1 ", otherwise, put Internal relations Intra-L ' be " 0 "; Make Intra-R ' expression
Figure BDA00000731500900000419
Internal relations, judge
Figure BDA00000731500900000420
The discrete cosine transform DC coefficient whether greater than
Figure BDA00000731500900000421
The low frequency coefficient of desired location of secondary wavelet transform, if then put
Figure BDA0000073150090000051
Internal relations Intra-R ' be " 1 ", otherwise, put
Figure BDA0000073150090000052
Internal relations Intra-R ' be " 0 "; Make Intra ' expression
Figure BDA0000073150090000053
Internal relations, judge
Figure BDA0000073150090000054
Internal relations Intra-L ' and
Figure BDA0000073150090000055
Internal relations Intra-R ' whether be " 1 ", if then put
Figure BDA0000073150090000056
Internal relations Intra ' be " 1 ", otherwise, put
Figure BDA0000073150090000057
Internal relations Intra ' be " 0 "; Make Inter ' expression
Figure BDA0000073150090000058
Interblock relation, judge
Figure BDA0000073150090000059
The discrete cosine transform DC coefficient whether greater than
Figure BDA00000731500900000510
The discrete cosine transform DC coefficient, if then put
Figure BDA00000731500900000511
Interblock relations I nter ' be " 1 ", otherwise, put Interblock relations I nter ' be " 0 ";
Above-mentioned The secondary wavelet transform desired location low frequency coefficient with
Figure BDA00000731500900000514
The low frequency coefficient of desired location of secondary wavelet transform be the low frequency coefficient of same position of the secondary wavelet transform of image block separately;
2.-4, right
Figure BDA00000731500900000515
With
Figure BDA00000731500900000516
The discrete cosine transform DC coefficient quantize respectively according to quantization step S, obtain quantized value Q ' respectively L, jAnd Q ' R, j, Q L , j ′ = floor ( DC ′ L , j S ) , Q R , j ′ = floor ( DC ′ R , j S ) , Wherein, DC ' L, jExpression
Figure BDA00000731500900000519
The discrete cosine transform DC coefficient, DC ' R, jExpression
Figure BDA00000731500900000520
The discrete cosine transform DC coefficient, S is a quantization step, floor (x) is for calculating greatest integer function smaller or equal to x;
2.-5, judge that current image block is right
Figure BDA00000731500900000521
The watermark embedded mode whether be 00, if not, then execution in step 2.-6, otherwise, judge again
Figure BDA00000731500900000522
Internal relations Intra-L ' and
Figure BDA00000731500900000523
Internal relations Intra-R ' whether be " 1 ", if be " 1 ", then right from current image block
Figure BDA00000731500900000524
In extract watermark " 1 ", be " 0 " if be " 0 " or one of them, then right from current image block In extract watermark " 0 ", execution in step is 2.-8 then;
2.-6, judge that current image block is right
Figure BDA00000731500900000526
The watermark embedded mode whether be 01, if not, then execution in step 2.-7, otherwise, judge that again current image block is right
Figure BDA00000731500900000527
Interblock relations I nter ' whether be " 1 ", if be " 1 ", then right from current image block
Figure BDA00000731500900000528
In extract watermark " 1 ", if be " 0 ", then right from current image block
Figure BDA00000731500900000529
In extract watermark " 0 ", execution in step is 2.-8 then;
2.-7, judge Q ' L, jAnd Q ' R, jWhether be odd number or be even number, if, then right from current image block
Figure BDA00000731500900000530
In extract watermark " 1 ", otherwise, right from current image block
Figure BDA00000731500900000531
In extract watermark " 0 ", execution in step is 2.-8 then;
2.-8, judge whether j≤m * m sets up, if establishment then adds 1 with the j value, with stereoscopic video images S to be detected WatIn next image block to right as current image block, and repeated execution of steps 2.-2 to step 2.-8, until handling m * m image block to finishing watermark extraction process;
3. the watermarking images that extraction the is obtained random conversion that is inverted, the binary bitmap of finally being extracted.
Described step 2.-3 in the low frequency coefficient of desired location be the low frequency coefficient of arbitrary position of the secondary wavelet transform of image block.
Compared with prior art, the invention has the advantages that by stereoscopic video images being divided into the image block of non-overlapping copies, each image block is implemented discrete cosine transform and secondary wavelet transform, analyze the left side of stereoscopic video images, magnitude relationship between the discrete cosine transform DC coefficient of right visual point image and the low frequency coefficient of secondary wavelet transform, thereby the internal relations of definition image block and the relation of the interblock between two image blocks in stereoscopic video images, by these relations and the watermark information consistent watermark embedded mode that decides whether, do not change the data of original stereoscopic video images when wherein relation embeds, when quantizing to embed the discrete cosine transform DC coefficient of the little image block of discrete cosine transform DC coefficient modification amplitude is made amendment effective like this quality that guarantees stereoscopic video images; Simultaneously when watermark extracting, do not need original stereoscopic video images, realized blind Detecting.What select for use when relation embeds is the relation of the discrete cosine transform DC coefficient of the relation of low frequency coefficient of discrete cosine transform DC coefficient in the piece and secondary wavelet transform and interblock, because the low frequency coefficient of discrete cosine transform DC coefficient and secondary wavelet transform has identical variation tendency, and the pass of the discrete cosine transform DC coefficient of interblock ties up to identical variation relation under the identical attack, therefore the stability of relation is stronger, thereby has effectively improved the robustness of embedding grammar of the present invention.When relation embeds, do not revise the data of original stereo-picture, thereby improved the quality of watermarked stereo-picture effectively.
Description of drawings
Fig. 1 a is the left visual point image in first moment of " puppy " original stereoscopic video images;
Fig. 1 b is the right visual point image in first moment of " puppy " original stereoscopic video images;
Fig. 2 a is original binary bitmap;
Fig. 2 b is the binary bitmap behind the original binary bitmap scramble shown in Fig. 2 a;
Fig. 3 a is first constantly the left visual point image of " puppy " stereoscopic video images after watermarked;
Fig. 3 b is first constantly the right visual point image of " puppy " stereoscopic video images after watermarked;
All watermarks that from the stereoscopic video images of Fig. 3 a and Fig. 3 b formation, extract when Fig. 4 a attacks for not having;
The binary bitmap that from the stereoscopic video images of Fig. 3 a and Fig. 3 b formation, extracts when Fig. 4 b attacks for not having;
Fig. 5 a is for from quality factor being the binary bitmap that extracts the stereoscopic video images after 90 the JPEG compression;
Fig. 5 b is for from quality factor being the binary bitmap that extracts the stereoscopic video images after 70 the JPEG compression;
Fig. 5 c is for from quality factor being the binary bitmap that extracts the stereoscopic video images after 50 the JPEG compression;
Fig. 5 d is for from quality factor being the binary bitmap that extracts the stereoscopic video images after 30 the JPEG compression;
Fig. 5 e is for from quality factor being the binary bitmap that extracts the stereoscopic video images after 20 the JPEG compression;
The binary bitmap of Fig. 6 a for from the stereoscopic video images behind the gaussian filtering (window 3 * 3), extracting;
The binary bitmap of Fig. 6 b for from the stereoscopic video images behind the medium filtering (window 3 * 3), extracting;
The binary bitmap of Fig. 6 c for from the stereoscopic video images behind the medium filtering (window 5 * 5), extracting;
The binary bitmap of Fig. 6 d for from the stereoscopic video images behind the mean filter (window 3 * 3), extracting;
The binary bitmap of Fig. 6 e for from the stereoscopic video images behind the mean filter (window 5 * 5), extracting;
Fig. 7 a is for shearing the binary bitmap that extracts the stereoscopic video images after 1/4 from the top left end;
Fig. 7 b is for shearing the binary bitmap that extracts the stereoscopic video images after 1/8 from the top left end;
Fig. 7 c is for shearing the binary bitmap that extracts the stereoscopic video images after 1/4 from the centre;
Fig. 7 d is for shearing the binary bitmap that extracts the stereoscopic video images after 1/8 from the centre;
Fig. 7 e is for shearing the binary bitmap that extracts the stereoscopic video images after 1/2 from the top left end;
Fig. 8 is for successively to original stereoscopic video images S OrgIn each image block to adding the process schematic diagram of watermark.
Embodiment
Embodiment describes in further detail the present invention below in conjunction with accompanying drawing.
Embodiment one:
The blind watermark embedding method of a kind of stereoscopic video images that present embodiment proposes, it mainly may further comprise the steps:
1. at the watermark built-in end, make S OrgFor being of a size of the original stereoscopic video images of M * N, original stereoscopic video images S OrgComprise left visual point image and right visual point image, remember original stereoscopic video images S OrgLeft visual point image be L Org, remember original stereoscopic video images S OrgRight visual point image be R Org, with original stereoscopic video images S OrgLeft visual point image L OrgWith original stereoscopic video images S OrgRight visual point image R OrgBe divided into respectively (M * N)/(n * n) individual non-overlapping copies is of a size of the image block of n * n, with original stereoscopic video images S OrgLeft visual point image L OrgIn i image block be designated as
Figure BDA0000073150090000071
With original stereoscopic video images S OrgRight visual point image R OrgIn i image block be designated as
Figure BDA0000073150090000072
Original stereoscopic video images S OrgLeft visual point image L OrgIn i image block
Figure BDA0000073150090000073
With original stereoscopic video images S OrgRight visual point image R OrgIn i image block (be left visual point image L OrgWith right visual point image R OrgIn two image blocks of same coordinate position) constitute original stereoscopic video images S OrgIn i image block right, be designated as
Figure BDA0000073150090000075
Wherein, M represents original stereoscopic video images S OrgWide, N represents original stereoscopic video images S OrgHeight, 1≤i≤((M * N)/(n * n)).
In the present embodiment, what the original stereoscopic video images of test was selected for use is the stereoscopic video images in first moment in " puppy " stereoscopic video sequence, resolution sizes is 640 * 480, be M=640, N=480, Fig. 1 a and Fig. 1 b have provided " puppy " original stereoscopic video images first left visual point image and right visual point image constantly respectively.
In the present embodiment, the size n of image block * n can be selected according to the size of stereoscopic video images, as selecting 8 * 8,16 * 16 etc.In actual application, if select too for a short time, then may cause watermarked hiding property bad, and if select too greatly, then can cause the watermark information that embeds very little.In the present embodiment, select n=8.
2. make that W is the binary bitmap that is of a size of m * m, binary bitmap W is carried out the scramble conversion, obtain the binary bitmap behind the scramble, be designated as Wherein, m * m≤((M * N)/(n * n)).In the present embodiment, what the scramble conversion was adopted is the conversion of Arnold scramble, the conversion of Arnold scramble is the one-period conversion, and the data array of conversion must be a square formation, its transformation period is relevant with the size of square formation, for example be 48 the transformation period of 64 * 64 Arnold scramble conversion, and number of transitions L should just can make data reach the purpose of scramble less than 48 so.In the present embodiment, adopt to be of a size of 64 * 64 original binary bitmap shown in Fig. 2 a, the binary bitmap behind its scramble is shown in Fig. 2 b.
3. successively to original stereoscopic video images S OrgIn each image block to adding watermark, its detailed process (Fig. 8 has provided concrete watermark and added process) is as follows:
3.-1, the current original stereoscopic video images S that is handling of definition OrgIn i image block right
Figure BDA0000073150090000082
For current image block right.
3.-2, right to current image block
Figure BDA0000073150090000083
In Carry out respectively discrete cosine transform (Discrete Cosine Transform, DCT) and the secondary wavelet transform (Discrete Wavelet Transform, DWT), right to current image block
Figure BDA0000073150090000085
In
Figure BDA0000073150090000086
Carry out discrete cosine transform and secondary wavelet transform respectively.
3.-3, make Intra-L represent Internal relations, judge
Figure BDA0000073150090000088
The discrete cosine transform DC coefficient whether greater than The low frequency coefficient of desired location of secondary wavelet transform, if then put
Figure BDA00000731500900000810
Internal relations Intra-L be " 1 ", otherwise, put
Figure BDA00000731500900000811
Internal relations Intra-L be " 0 "; Make Intra-R represent Internal relations, judge The discrete cosine transform DC coefficient whether greater than The low frequency coefficient of desired location of secondary wavelet transform, if then put
Figure BDA00000731500900000815
Internal relations Intra-R be " 1 ", otherwise, put
Figure BDA00000731500900000816
Internal relations Intra-R be " 0 "; Make Intra represent
Figure BDA00000731500900000817
Internal relations, judge
Figure BDA00000731500900000818
Internal relations Intra-L and
Figure BDA00000731500900000819
Internal relations Intra-R whether be " 1 ", if then put
Figure BDA00000731500900000820
Internal relations Intra be " 1 ", otherwise, put
Figure BDA00000731500900000821
Internal relations Intra be " 0 "; Make Inter represent
Figure BDA00000731500900000822
Interblock relation, judge
Figure BDA00000731500900000823
The discrete cosine transform DC coefficient whether greater than
Figure BDA00000731500900000824
The discrete cosine transform DC coefficient, if then put
Figure BDA00000731500900000825
Interblock relations I nter be " 1 ", otherwise, put Interblock relations I nter be " 0 ".
Above-mentioned
Figure BDA0000073150090000092
The secondary wavelet transform desired location low frequency coefficient with
Figure BDA0000073150090000093
The low frequency coefficient of desired location of secondary wavelet transform be the low frequency coefficient of same position of the secondary wavelet transform of image block separately.In the present embodiment, the low frequency coefficient of above-mentioned desired location is the low frequency coefficient of arbitrary position of the secondary wavelet transform of image block, as herein The secondary wavelet transform low frequency coefficient and
Figure BDA0000073150090000095
The low frequency coefficient of secondary wavelet transform chooses is the low frequency coefficient of the 0th row the 1st row of the low frequency sub-band of image block separately.
3.-4, judge that current image block is right
Figure BDA0000073150090000096
The value of internal relations Intra and the binary bitmap behind the scramble
Figure BDA0000073150090000097
In the pixel value w of i binarized pixel point iWhether identical, if identical, then do not revise
Figure BDA0000073150090000098
With The discrete cosine transform DC coefficient, and the mark current image block is right
Figure BDA00000731500900000910
The watermark embedded mode be 00, execution in step is 3.-7 then; Otherwise execution in step 3.-5.
At this, the watermark of this step embeds the relation of being called and embeds.
3.-5, judge that current image block is right
Figure BDA00000731500900000911
The value of interblock relations I nter and the binary bitmap behind the scramble
Figure BDA00000731500900000912
In the pixel value w of i binarized pixel point iWhether identical, if identical, then do not revise
Figure BDA00000731500900000913
With
Figure BDA00000731500900000914
The discrete cosine transform DC coefficient, and the mark current image block is right
Figure BDA00000731500900000915
The watermark embedded mode be 01, execution in step is 3.-7 then; Otherwise execution in step 3.-6.
At this, the watermark of this step embeds the relation of being called and embeds.
3.-6, right
Figure BDA00000731500900000916
With The discrete cosine transform DC coefficient quantize respectively according to quantization step S, obtain quantized value Q respectively L, iAnd Q R, i, Q L , i = floor ( DC L , i S ) , Q R , i = floor ( DC R , i S ) , Wherein, DC L, iExpression
Figure BDA00000731500900000920
The discrete cosine transform DC coefficient, DC R, iExpression
Figure BDA00000731500900000921
The discrete cosine transform DC coefficient, S is a quantization step, the value of quantization step S choose according to for watermark information in carrier image hiding property and the compromise of robustness, get S=10 at this, floor (x) is for calculating greatest integer function smaller or equal to x; Definition
Figure BDA00000731500900000922
The parity of discrete cosine transform DC coefficient be its quantized value Q L, iParity, the definition
Figure BDA00000731500900000923
The parity of discrete cosine transform DC coefficient be its quantized value Q R, iParity, judge quantized value Q L, iWith quantized value Q R, iWhether be all even number or be all odd number, if then think
Figure BDA00000731500900000924
With
Figure BDA00000731500900000925
The parity unanimity of discrete cosine transform DC coefficient, otherwise, think
Figure BDA00000731500900000926
With
Figure BDA00000731500900000927
The parity of discrete cosine transform DC coefficient inconsistent.
Binary bitmap behind the judgement scramble In the pixel value w of i binarized pixel point iWhether be " 1 ", if it is right then to revise current image block
Figure BDA00000731500900000929
In the discrete cosine transform DC coefficient of any image block, make amended
Figure BDA00000731500900000930
With
Figure BDA00000731500900000931
The parity unanimity of discrete cosine transform DC coefficient, and the mark current image block is right
Figure BDA0000073150090000101
The watermark embedded mode be 10, otherwise it is right to revise current image block
Figure BDA0000073150090000102
In the discrete cosine transform DC coefficient of any image block, make amended
Figure BDA0000073150090000103
With
Figure BDA0000073150090000104
The parity of discrete cosine transform DC coefficient inconsistent, and the mark current image block is right
Figure BDA0000073150090000105
The watermark embedded mode be 10.
In the present embodiment, right in the modification current image block
Figure BDA0000073150090000106
In the discrete cosine transform DC coefficient of any image block the time, select With
Figure BDA0000073150090000108
Middle discrete cosine transform DC coefficient is revised the discrete cosine transform DC coefficient of that little image block of amplitude and is made amendment.
At this, the watermark of this step embeds to be called and quantizes to embed.
3.-7, judge whether i≤m * m sets up, if establishment then adds 1 with the i value, with original stereoscopic video images S OrgIn next image block to right as current image block, and repeated execution of steps 3.-2 to step 3.-7, until handling m * m image block to finishing watermark embed process.
4., the size with binary bitmap, watermark embedded mode and quantization step S be as key, the watermark built-in end extracts end with cipher key delivery to watermark.
Fig. 3 a and Fig. 3 b have provided left visual point image and the right visual point image after the binary bitmap shown in embedding Fig. 2 b respectively.After watermarked, the left visual point image of " puppy " stereoscopic video images and the Y-PSNR PSNR of right visual point image are respectively 49.34dB and 51.34dB, show that stereoscopic video images and original stereoscopic video images after watermarked are quite similar, the hiding property of watermark information is fine.
Embodiment two:
The blind watermark extracting method of the stereoscopic video images that the blind watermark embedding method that a kind of embodiment of use one of present embodiment proposition provides is watermarked, it mainly may further comprise the steps:
1. at the watermark extracting end, make S WatFor being of a size of the stereoscopic video images to be detected of M * N, stereoscopic video images S to be detected WatComprise left visual point image and right visual point image, remember stereoscopic video images S to be detected WatLeft visual point image be L Wat, remember stereoscopic video images S to be detected WatRight visual point image be R Wat, with stereoscopic video images S to be detected WatLeft visual point image L WatWith stereoscopic video images S to be detected WatRight visual point image R WatBe divided into respectively (M * N)/(n * n) individual non-overlapping copies is of a size of the image block of n * n, with stereoscopic video images S to be detected WatLeft visual point image L WatIn j image block be designated as
Figure BDA0000073150090000109
With stereoscopic video images S to be detected WatRight visual point image R WatIn j image block be designated as
Figure BDA00000731500900001010
Stereoscopic video images S to be detected WatLeft visual point image L WatIn j image block
Figure BDA00000731500900001011
With stereoscopic video images S to be detected WatRight visual point image R WatIn j image block
Figure BDA00000731500900001012
(be left visual point image L WatWith right visual point image R WatIn two image blocks of same coordinate position) constitute stereoscopic video images S to be detected WatIn j image block right, be designated as
Figure BDA00000731500900001013
Wherein, M represents stereoscopic video images S to be detected WatWide, N represents stereoscopic video images S to be detected WatHeight, 1≤j≤((M * N)/(n * n)).
In the present embodiment, stereoscopic video images S to be detected WatSize consistent with the size of the original stereoscopic video images of watermark built-in end, the size of the image block also size with watermark built-in end divided image piece is consistent.
2. successively to stereoscopic video images S to be detected WatIn each image block to extracting watermark, its detailed process is as follows:
2.-1, the current stereoscopic video images S to be detected that is handling of definition WatIn j image block right
Figure BDA0000073150090000111
For current image block right.
2.-2, right to current image block
Figure BDA0000073150090000112
In
Figure BDA0000073150090000113
Carry out discrete cosine transform and secondary wavelet transform respectively, right to current image block
Figure BDA0000073150090000114
In
Figure BDA0000073150090000115
Carry out discrete cosine transform and secondary wavelet transform respectively.
2.-3, make Intra-L ' expression
Figure BDA0000073150090000116
Internal relations, judge
Figure BDA0000073150090000117
The discrete cosine transform DC coefficient whether greater than The low frequency coefficient of desired location of secondary wavelet transform, if then put
Figure BDA0000073150090000119
Internal relations Intra-L ' be " 1 ", otherwise, put
Figure BDA00000731500900001110
Internal relations Intra-L ' be " 0 "; Make Intra-R ' expression
Figure BDA00000731500900001111
Internal relations, judge
Figure BDA00000731500900001112
The discrete cosine transform DC coefficient whether greater than
Figure BDA00000731500900001113
The low frequency coefficient of desired location of secondary wavelet transform, if then put
Figure BDA00000731500900001114
Internal relations Intra-R ' be " 1 ", otherwise, put
Figure BDA00000731500900001115
Internal relations Intra-R ' be " 0 "; Make Intra ' expression
Figure BDA00000731500900001116
Internal relations, judge
Figure BDA00000731500900001117
Internal relations Intra-L ' and
Figure BDA00000731500900001118
Internal relations Intra-R ' whether be " 1 ", if then put
Figure BDA00000731500900001119
Internal relations Intra ' be " 1 ", otherwise, put
Figure BDA00000731500900001120
Internal relations Intra ' be " 0 "; Make Inter ' expression Interblock relation, judge
Figure BDA00000731500900001122
The discrete cosine transform DC coefficient whether greater than
Figure BDA00000731500900001123
The discrete cosine transform DC coefficient, if then put
Figure BDA00000731500900001124
Interblock relations I nter ' be " 1 ", otherwise, put Interblock relations I nter ' be " 0 ".
Above-mentioned
Figure BDA00000731500900001126
The secondary wavelet transform desired location low frequency coefficient with
Figure BDA00000731500900001127
The low frequency coefficient of desired location of secondary wavelet transform be the low frequency coefficient of same position of the secondary wavelet transform of image block separately.In the present embodiment, the low frequency coefficient of above-mentioned desired location is the low frequency coefficient of arbitrary position of the secondary wavelet transform of image block, as herein
Figure BDA00000731500900001128
The secondary wavelet transform low frequency coefficient and
Figure BDA00000731500900001129
The low frequency coefficient of secondary wavelet transform chooses is the low frequency coefficient of the 0th row the 1st row of the low frequency sub-band of image block separately.
2.-4, right
Figure BDA00000731500900001130
With
Figure BDA00000731500900001131
The discrete cosine transform DC coefficient quantize respectively according to quantization step S, obtain quantized value Q ' respectively L, jAnd Q ' R, jAt this, Q L , j ′ = floor ( DC ′ L , j S ) , Q R , j ′ = floor ( DC ′ R , j S ) , Wherein, DC ' L, jExpression
Figure BDA00000731500900001134
The discrete cosine transform DC coefficient, DC ' R, jExpression
Figure BDA00000731500900001135
The discrete cosine transform DC coefficient, S is a quantization step, floor (x) is for calculating greatest integer function smaller or equal to x.
2.-5, judge that current image block is right
Figure BDA0000073150090000121
The watermark embedded mode whether be 00, if not, then execution in step 2.-6, otherwise, judge again Internal relations Intra-L ' and Internal relations Intra-R ' whether be " 1 ", if be " 1 ", then right from current image block
Figure BDA0000073150090000124
In extract watermark " 1 ", be " 0 " if be " 0 " or one of them, then right from current image block
Figure BDA0000073150090000125
In extract watermark " 0 ", execution in step is 2.-8 then.
2.-6, judge that current image block is right The watermark embedded mode whether be 01, if not, then execution in step 2.-7, otherwise, judge that again current image block is right
Figure BDA0000073150090000127
Interblock relations I nter ' whether be " 1 ", if be " 1 ", then right from current image block
Figure BDA0000073150090000128
In extract watermark " 1 ", if be " 0 ", then right from current image block
Figure BDA0000073150090000129
In extract watermark " 0 ", execution in step is 2.-8 then.
2.-7, when current image block right
Figure BDA00000731500900001210
The watermark embedded mode be 10 o'clock, judge Q ' L, jAnd Q ' R, jWhether be odd number or be even number, if, then right from current image block
Figure BDA00000731500900001211
In extract watermark " 1 ", otherwise, right from current image block
Figure BDA00000731500900001212
In extract watermark " 0 ", execution in step is 2.-8 then.
2.-8, judge whether j≤m * m sets up, if establishment then adds 1 with the j value, with stereoscopic video images S to be detected WatIn next image block to right as current image block, and repeated execution of steps 2.-2 to step 2.-8, until handling m * m image block to finishing watermark extraction process.
3. the watermarking images that extraction the is obtained random conversion that is inverted, the binary bitmap of finally being extracted.In the present embodiment, owing to adopt the conversion of L Arnold scramble, therefore can be transferred to the size of the binary bitmap of watermark extracting end according to the watermark built-in end, determine the cycle Z of Arnold scramble conversion, the conversion of (Z-L) inferior Arnold scramble, the binary bitmap of finally being extracted are carried out in all watermarks that extraction is obtained then.
Fig. 4 a and Fig. 4 b have provided does not respectively have all watermarks and the final binary bitmap that obtains of extracting that extracts when attacking from the stereoscopic video images of Fig. 3 a and Fig. 3 b formation.In the present embodiment, adopt normalization similar value HC,
Figure BDA00000731500900001213
Calculate and extract the binary bitmap obtain and the similarity of original binary bitmap, 0≤HC≤1, wherein w (s, t) and
Figure BDA00000731500900001214
Represent respectively original binary bitmap and the binary bitmap coordinate position that extracts for (s, the pixel value of t) locating,
Figure BDA00000731500900001215
The expression xor operator.HC is big more, and the approaching more original binary bitmap of binary bitmap that obtains is extracted in expression, and HC is more little, and expression is extracted the binary bitmap and the original binary bitmap that obtain and differed bigger.Can estimate the robustness of watermarking algorithm by the size of HC value.Do not having under the situation of attacking, the normalization similar value HC of the binary bitmap that present embodiment extracts from the stereoscopic video images of Fig. 3 a and Fig. 3 b formation is 1, and is promptly identical with original binary bitmap.
Further, the stereoscopic video images that is embedded with watermark that obtains is carried out distortion handle, comprise JPEG compression, the gaussian filtering of different windows size, medium filtering, mean filter under the different quality factors, and the shearing of position in various degree.Fig. 5 a, Fig. 5 b, Fig. 5 c, Fig. 5 d and Fig. 5 e have provided the binary bitmap that extracts the stereoscopic video images after quality factor q is respectively 90,70,50,30,20 JPEG compression respectively; Fig. 6 a, Fig. 6 b, Fig. 6 c, Fig. 6 d and Fig. 6 e have provided the binary bitmap that extracts respectively from the stereoscopic video images behind gaussian filtering (window 3 * 3), medium filtering (window 3 * 3), medium filtering (window 5 * 5), mean filter (window 3 * 3), the mean filter (window 5 * 5); Fig. 7 a, Fig. 7 b, Fig. 7 c, Fig. 7 d and Fig. 7 e have provided the binary bitmap that extracts respectively from the stereoscopic video images after top left end shearing 1/4, top left end shearing 1/8, middle shearing 1/4, middle shearing 1/8, the top left end shearing 1/2.Table 1 has provided above-mentioned various distortion and has handled the normalization similar value HC that the binary bitmap that obtains is extracted in the back.From Fig. 5 a to Fig. 7 e and table 1 as seen, under distortion dissimilar, in various degree was attacked, the watermark that the inventive method embeds had identification preferably, illustrates that embedding of the present invention and extracting method have robustness preferably.
Table 1 is dissimilar, the normalization similar value HC of watermark detection is attacked down in distortion in various degree
JPEG compresses Q=90 JPEG compresses Q=70 JPEG compresses Q=50 JPEG compresses Q=30 JPEG compresses Q=20
HC=0.9258 HC=0.9063 HC=0.8601 HC=0.7324 HC=0.7048
The strange window of gaussian filtering Medium filtering 3 * 3 Medium filtering 5 * 5 Mean filter 3 * 3 Mean filter 5 * 5
HC=0.9932 HC=0.8965 HC=0.8091 HC=0.9216 HC=0.7986
The top left end shears 1/4 The top left end shears 1/8 Middle shearing 1/4 Middle shearing 1/8 The top left end shears 1/2
HC=0.9198 HC=0.9819 HC=0.8369 HC=0.9216 HC=0.8438

Claims (3)

1. the blind watermark embedding method of a stereoscopic video images is characterized in that may further comprise the steps:
1. at the watermark built-in end, make S OrgFor being of a size of the original stereoscopic video images of M * N, original stereoscopic video images S OrgComprise left visual point image and right visual point image, remember original stereoscopic video images S OrgLeft visual point image be L Org, remember original stereoscopic video images S OrgRight visual point image be R Org, with original stereoscopic video images S OrgLeft visual point image L OrgWith original stereoscopic video images S OrgRight visual point image R OrgBe divided into respectively (M * N)/(n * n) individual non-overlapping copies is of a size of the image block of n * n, with original stereoscopic video images S OrgLeft visual point image L OrgIn i image block be designated as
Figure FDA00002958474800011
With original stereoscopic video images S OrgRight visual point image R OrgIn i image block be designated as
Figure FDA00002958474800012
Original stereoscopic video images S OrgLeft visual point image L OrgIn i image block With original stereoscopic video images S OrgRight visual point image R OrgIn i image block
Figure FDA00002958474800014
Constitute original stereoscopic video images S OrgIn i image block right, be designated as
Figure FDA00002958474800015
Wherein, M represents original stereoscopic video images S OrgWide, N represents original stereoscopic video images S OrgHeight, 1≤i≤((M * N)/(n * n));
2. make that W is the binary bitmap that is of a size of m * m, binary bitmap W is carried out the scramble conversion, obtain the binary bitmap behind the scramble, be designated as
Figure FDA00002958474800016
Wherein, m * m≤((M * N)/(n * n));
3. successively to original stereoscopic video images S OrgIn each image block to adding watermark, its detailed process is as follows:
3.-1, the current original stereoscopic video images S that is handling of definition OrgIn i image block right
Figure FDA00002958474800017
For current image block right;
3.-2, right to current image block
Figure FDA00002958474800018
In
Figure FDA00002958474800019
Carry out discrete cosine transform and secondary wavelet transform respectively, right to current image block
Figure FDA000029584748000110
In Carry out discrete cosine transform and secondary wavelet transform respectively;
3.-3, make Intra-L represent
Figure FDA000029584748000112
Internal relations, judge
Figure FDA000029584748000113
The discrete cosine transform DC coefficient whether greater than
Figure FDA000029584748000114
The low frequency coefficient of desired location of secondary wavelet transform, if then put
Figure FDA000029584748000115
Internal relations Intra-L be " 1 ", otherwise, put Internal relations Intra-L be " 0 "; Make Intra-R represent
Figure FDA000029584748000117
Internal relations, judge
Figure FDA000029584748000118
The discrete cosine transform DC coefficient whether greater than
Figure FDA000029584748000119
The low frequency coefficient of desired location of secondary wavelet transform, if then put
Figure FDA000029584748000120
Internal relations Intra-R be " 1 ", otherwise, put
Figure FDA000029584748000121
Internal relations Intra-R be " 0 "; Make Intra represent Internal relations, judge Internal relations Intra-L and
Figure FDA00002958474800023
Internal relations Intra-R whether be " 1 ", if then put
Figure FDA00002958474800024
Internal relations Intra be " 1 ", otherwise, put
Figure FDA00002958474800025
Internal relations Intra be " 0 "; Make Inter represent
Figure FDA00002958474800026
Interblock relation, judge The discrete cosine transform DC coefficient whether greater than
Figure FDA00002958474800028
The discrete cosine transform DC coefficient, if then put
Figure FDA00002958474800029
Interblock relations I nter be " 1 ", otherwise, put
Figure FDA000029584748000210
Interblock relations I nter be " 0 ";
Above-mentioned
Figure FDA000029584748000211
The secondary wavelet transform desired location low frequency coefficient with
Figure FDA000029584748000212
The low frequency coefficient of desired location of secondary wavelet transform be the low frequency coefficient of same position of the secondary wavelet transform of image block separately;
Described step 3.-3 in the low frequency coefficient of desired location be the low frequency coefficient of arbitrary position of the secondary wavelet transform of image block;
3.-4, judge that current image block is right The value of internal relations Intra and the binary bitmap behind the scramble
Figure FDA000029584748000214
In the pixel value w of i binarized pixel point iWhether identical, if identical, then do not revise
Figure FDA000029584748000215
With
Figure FDA000029584748000216
The discrete cosine transform DC coefficient, and the mark current image block is right
Figure FDA000029584748000217
The watermark embedded mode be 00, execution in step is 3.-7 then; Otherwise execution in step 3.-5;
3.-5, judge that current image block is right
Figure FDA000029584748000218
The value of interblock relations I nter and the binary bitmap behind the scramble
Figure FDA000029584748000219
In the pixel value w of i binarized pixel point iWhether identical, if identical, then do not revise With
Figure FDA000029584748000221
The discrete cosine transform DC coefficient, and the mark current image block is right
Figure FDA000029584748000222
The watermark embedded mode be 01, execution in step is 3.-7 then; Otherwise execution in step 3.-6;
3.-6, right
Figure FDA000029584748000223
With
Figure FDA000029584748000224
The discrete cosine transform DC coefficient quantize respectively according to quantization step S, obtain quantized value Q respectively L, iAnd Q R, i, Q L , i = floor ( DC L , i S ) , Q R , i = floor ( DC R , i S ) , Wherein, DC L, iExpression
Figure FDA000029584748000227
The discrete cosine transform DC coefficient, DC R, iExpression The discrete cosine transform DC coefficient, floor (x) is for calculating greatest integer function smaller or equal to x; Definition
Figure FDA000029584748000229
The parity of discrete cosine transform DC coefficient be its quantized value Q L, iParity, the definition
Figure FDA000029584748000230
The parity of discrete cosine transform DC coefficient be its quantized value Q R, iParity, judge quantized value Q L, iWith quantized value Q R, iWhether be all even number or be all odd number, if then think
Figure FDA000029584748000231
With
Figure FDA000029584748000232
The parity unanimity of discrete cosine transform DC coefficient, otherwise, think
Figure FDA000029584748000233
With
Figure FDA000029584748000234
The parity of discrete cosine transform DC coefficient inconsistent;
Binary bitmap behind the judgement scramble
Figure FDA000029584748000235
In the pixel value w of i binarized pixel point iWhether be " 1 ", if it is right then to revise current image block
Figure FDA000029584748000236
In the discrete cosine transform DC coefficient of any image block, make amended
Figure FDA00002958474800031
With
Figure FDA00002958474800032
The parity unanimity of discrete cosine transform DC coefficient, and the mark current image block is right
Figure FDA00002958474800033
The watermark embedded mode be 10, otherwise it is right to revise current image block
Figure FDA00002958474800034
In the discrete cosine transform DC coefficient of any image block, make amended
Figure FDA00002958474800035
With
Figure FDA00002958474800036
The parity of discrete cosine transform DC coefficient inconsistent, and the mark current image block is right
Figure FDA000029584748000313
The watermark embedded mode be 10;
3.-7, judge whether i≤m * m sets up, if establishment then adds 1 with the i value, with original stereoscopic video images S OrgIn next image block to right as current image block, and repeated execution of steps 3.-2 to step 3.-7, until handling m * m image block to finishing watermark embed process;
4., the size with binary bitmap, watermark embedded mode and quantization step S be as key, the watermark built-in end extracts end with cipher key delivery to watermark.
2. the blind watermark embedding method of a kind of stereoscopic video images according to claim 1 is characterized in that during described step 3.-6 in that to revise current image block right
Figure FDA00002958474800037
In the discrete cosine transform DC coefficient of any image block the time, select With
Figure FDA000029584748000314
Middle discrete cosine transform DC coefficient is revised the discrete cosine transform DC coefficient of that little image block of amplitude and is made amendment.
3. blind watermark extracting method that uses the watermarked stereoscopic video images of the described blind watermark embedding method of claim 1 is characterized in that may further comprise the steps:
1. at the watermark extracting end, make S WatFor being of a size of the stereoscopic video images to be detected of M * N, stereoscopic video images S to be detected WatComprise left visual point image and right visual point image, remember stereoscopic video images S to be detected WatLeft visual point image be L Wat, remember stereoscopic video images S to be detected WatRight visual point image be R Wat, with stereoscopic video images S to be detected WatLeft visual point image L WatWith stereoscopic video images S to be detected WatRight visual point image R WatBe divided into respectively (M * N)/(n * n) individual non-overlapping copies is of a size of the image block of n * n, with stereoscopic video images S to be detected WatLeft visual point image L WatIn j image block be designated as With stereoscopic video images S to be detected WatRight visual point image R WatIn j image block be designated as
Figure FDA00002958474800039
Stereoscopic video images S to be detected WatLeft visual point image L WatIn j image block
Figure FDA000029584748000310
With stereoscopic video images S to be detected WatRight visual point image R WatIn j image block Constitute stereoscopic video images S to be detected WatIn j image block right, be designated as
Figure FDA000029584748000312
Wherein, M represents stereoscopic video images S to be detected WatWide, N represents stereoscopic video images S to be detected WatHeight, 1≤j≤((M * N)/(n * n));
2. successively to stereoscopic video images S to be detected WatIn each image block to extracting watermark, its detailed process is as follows:
2.-1, the current stereoscopic video images S to be detected that is handling of definition WatIn j image block right
Figure FDA00002958474800041
For current image block right;
2.-2, right to current image block In
Figure FDA00002958474800043
Carry out discrete cosine transform and secondary wavelet transform respectively, right to current image block
Figure FDA00002958474800044
In
Figure FDA00002958474800045
Carry out discrete cosine transform and secondary wavelet transform respectively;
2.-3, make Intra-L' represent
Figure FDA00002958474800046
Internal relations, judge
Figure FDA00002958474800047
The discrete cosine transform DC coefficient whether greater than
Figure FDA00002958474800048
The low frequency coefficient of desired location of secondary wavelet transform, if then put
Figure FDA00002958474800049
Internal relations Intra-L' be " 1 ", otherwise, put
Figure FDA000029584748000410
Internal relations Intra-L' be " 0 "; Make Intra-R' represent
Figure FDA000029584748000411
Internal relations, judge
Figure FDA000029584748000412
The discrete cosine transform DC coefficient whether greater than
Figure FDA000029584748000413
The low frequency coefficient of desired location of secondary wavelet transform, if then put
Figure FDA000029584748000414
Internal relations Intra-R' be " 1 ", otherwise, put
Figure FDA000029584748000415
Internal relations Intra-R' be " 0 "; Make Intra' represent
Figure FDA000029584748000416
Internal relations, judge
Figure FDA000029584748000417
Internal relations Intra-L' and
Figure FDA000029584748000418
Internal relations Intra-R' whether be " 1 ", if then put
Figure FDA000029584748000419
Internal relations Intra' be " 1 ", otherwise, put
Figure FDA000029584748000420
Internal relations Intra' be " 0 "; Make Inter' represent Interblock relation, judge
Figure FDA000029584748000422
The discrete cosine transform DC coefficient whether greater than The discrete cosine transform DC coefficient, if then put
Figure FDA000029584748000424
Interblock relations I nter' be " 1 ", otherwise, put
Figure FDA000029584748000425
Interblock relations I nter' be " 0 ";
Above-mentioned
Figure FDA000029584748000426
The secondary wavelet transform desired location low frequency coefficient with
Figure FDA000029584748000427
The low frequency coefficient of desired location of secondary wavelet transform be the low frequency coefficient of same position of the secondary wavelet transform of image block separately;
Described step 2.-3 in the low frequency coefficient of desired location be the low frequency coefficient of arbitrary position of the secondary wavelet transform of image block;
2.-4, right
Figure FDA000029584748000428
With
Figure FDA000029584748000429
The discrete cosine transform DC coefficient quantize respectively according to quantization step S, obtain quantized value Q ' respectively L, jAnd Q ' R, j, Q L , j ′ = floor ( DC ′ L , j S ) , Q R , j ′ = floor ( DC ′ R , j S ) , Wherein, DC' L, jExpression
Figure FDA000029584748000432
The discrete cosine transform DC coefficient, DC' R, jExpression
Figure FDA000029584748000433
The discrete cosine transform DC coefficient, S is a quantization step, floor (x) is for calculating greatest integer function smaller or equal to x;
2.-5, judge that current image block is right
Figure FDA000029584748000434
The watermark embedded mode whether be 00, if not, then execution in step 2.-6, otherwise, judge again
Figure FDA000029584748000435
Internal relations Intra-L' and
Figure FDA000029584748000436
Internal relations Intra-R' whether be " 1 ", if be " 1 ", then right from current image block
Figure FDA000029584748000437
In extract watermark " 1 ", be " 0 " if be " 0 " or one of them, then right from current image block
Figure FDA000029584748000438
In extract watermark " 0 ", execution in step is 2.-8 then;
2.-6, judge that current image block is right
Figure FDA000029584748000439
The watermark embedded mode whether be 01, if not, then execution in step 2.-7, otherwise, judge that again current image block is right
Figure FDA00002958474800051
Interblock relations I nter' whether be " 1 ", if be " 1 ", then right from current image block
Figure FDA00002958474800052
In extract watermark " 1 ", if be " 0 ", then right from current image block
Figure FDA00002958474800053
In extract watermark " 0 ", execution in step is 2.-8 then;
2.-7, judge Q ' L, jAnd Q ' R, jWhether be odd number or be even number, if, then right from current image block
Figure FDA00002958474800054
In extract watermark " 1 ", otherwise, right from current image block In extract watermark " 0 ", execution in step is 2.-8 then;
2.-8, judge whether j≤m * m sets up, if establishment then adds 1 with the j value, with stereoscopic video images S to be detected WatIn next image block to right as current image block, and repeated execution of steps 2.-2 to step 2.-8, until handling m * m image block to finishing watermark extraction process;
3. the watermarking images that extraction the is obtained random conversion that is inverted, the binary bitmap of finally being extracted.
CN 201110183636 2011-07-01 2011-07-01 Blind watermark embedding and extracting method of stereoscopic video image Expired - Fee Related CN102223561B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110183636 CN102223561B (en) 2011-07-01 2011-07-01 Blind watermark embedding and extracting method of stereoscopic video image

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110183636 CN102223561B (en) 2011-07-01 2011-07-01 Blind watermark embedding and extracting method of stereoscopic video image

Publications (2)

Publication Number Publication Date
CN102223561A CN102223561A (en) 2011-10-19
CN102223561B true CN102223561B (en) 2013-07-31

Family

ID=44779964

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110183636 Expired - Fee Related CN102223561B (en) 2011-07-01 2011-07-01 Blind watermark embedding and extracting method of stereoscopic video image

Country Status (1)

Country Link
CN (1) CN102223561B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103366337A (en) * 2013-08-06 2013-10-23 山东大学 Digital watermarking method capable of being applied to 2D-3D (Two Dimensional-Three Dimensional) conversion
CN103533343B (en) * 2013-09-29 2015-05-20 宁波大学 Stereo image quality objective evaluation method based on digital watermarking
CN103763578A (en) * 2014-01-10 2014-04-30 北京酷云互动科技有限公司 Method and device for pushing program associated information
CN106507216B (en) * 2016-09-23 2019-08-02 曙光信息产业(北京)有限公司 The method that three-dimensional video-frequency is embedded in the method, apparatus of watermark and extracts watermark
CN106803230A (en) * 2017-01-16 2017-06-06 浙江工业大学 A kind of image encryption method based on digital watermarking
CN108648130B (en) * 2018-04-04 2022-05-27 宁波大学 Totally-blind digital watermarking method with copyright protection and tampering positioning functions
CN111260533B (en) * 2020-01-17 2023-04-25 山东师范大学 Image watermarking method and system for fusing texture rule features in image blocks and between blocks

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101699508B (en) * 2009-09-03 2012-01-11 中兴通讯股份有限公司 Image digital watermark embedding and extracting method and system

Also Published As

Publication number Publication date
CN102223561A (en) 2011-10-19

Similar Documents

Publication Publication Date Title
CN102223561B (en) Blind watermark embedding and extracting method of stereoscopic video image
Wang et al. Blind image data hiding based on self reference
Keshavarzian et al. ROI based robust and secure image watermarking using DWT and Arnold map
Jayalakshmi et al. Digital watermarking in contourlet domain
Divecha et al. Implementation and performance analysis of DCT-DWT-SVD based watermarking algorithms for color images
Ernawan et al. An improved imperceptibility and robustness of 4x4 dct-svd image watermarking with a modified entropy
CN103974144A (en) Video digital watermarking method based on characteristic scale variation invariant points and microscene detection
Khalifa et al. A robust non-blind algorithm for watermarking color images using multi-resolution wavelet decomposition
CN101583033A (en) Method for protecting H.264 video data by using robust watermarks
Thajeel et al. A NEW COLOR IMAGE WATERMARKING TECHNIQUE USING MULTIPLE DECOMPOSITIONS.
Ernawan et al. Block-based Tchebichef image watermarking scheme using psychovisual threshold
CN1276387C (en) Synchronous self-adaptable watermark method based on image continuity
CN109829846B (en) Digital image blind watermarking method based on two-dimensional discrete cosine transform
Yadav et al. Transform domain robust image-adaptive watermarking: prevalent techniques and their evaluation
Hana et al. A robust multiple watermarking scheme based on the DWT
Yu et al. A new stereo image watermarking method for 3D media
Shanthakumari et al. Video Steganography using LSB matching revisited algorithm
Al-Gindy et al. A graphical user interface watermarking technique for the copyright protection of colour images using colour watermarks
Azizi et al. Hybrid image watermarking using local complexity variations
Moniruzzaman et al. Robust RGB color image watermarking scheme based on DWT-SVD and chaotic system
Gao et al. A video dual watermarking algorithm against geometric attack based on integer wavelet and SIFT
Keshavarzian A new ROI and block based watermarking scheme using DWT
Kadu et al. Secure data hiding using robust firefly algorithm
Leng et al. High payload data hiding based on just noticeable distortion profile and LSB substitution
CN104143173A (en) Image self-adaption blind watermarking algorithm based on DWT-DCT

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130731

Termination date: 20210701