CN110570345B - Blind watermarking method for airspace color digital image fused with discrete cosine transform - Google Patents
Blind watermarking method for airspace color digital image fused with discrete cosine transform Download PDFInfo
- Publication number
- CN110570345B CN110570345B CN201910889463.0A CN201910889463A CN110570345B CN 110570345 B CN110570345 B CN 110570345B CN 201910889463 A CN201910889463 A CN 201910889463A CN 110570345 B CN110570345 B CN 110570345B
- Authority
- CN
- China
- Prior art keywords
- watermark
- pixel
- image
- layered
- blue
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000000605 extraction Methods 0.000 claims abstract description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 14
- 230000009466 transformation Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 4
- 238000003384 imaging method Methods 0.000 claims description 3
- 238000013139 quantization Methods 0.000 claims description 3
- 230000006835 compression Effects 0.000 description 8
- 238000007906 compression Methods 0.000 description 8
- 235000002566 Capsicum Nutrition 0.000 description 4
- 239000006002 Pepper Substances 0.000 description 4
- 241000722363 Piper Species 0.000 description 4
- 235000016761 Piper aduncum Nutrition 0.000 description 4
- 235000017804 Piper guineense Nutrition 0.000 description 4
- 235000008184 Piper nigrum Nutrition 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 150000003839 salts Chemical class 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 239000004973 liquid crystal related substance Substances 0.000 description 3
- 238000012795 verification Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T1/00—General purpose image data processing
- G06T1/0021—Image watermarking
- G06T1/005—Robust watermarking, e.g. average attack or collusion attack resistant
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2201/00—General purpose image data processing
- G06T2201/005—Image watermarking
- G06T2201/0065—Extraction of an embedded watermark; Reliable detection
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Editing Of Facsimile Originals (AREA)
- Image Processing (AREA)
Abstract
The invention discloses a blind watermarking method of an airspace color digital image fused with discrete cosine transform by combining the advantages of high running speed of an airspace digital watermarking algorithm and high robustness of a frequency domain digital watermarking algorithm. According to the unique characteristics of the discrete cosine transform direct current coefficients, the invention utilizes the correlation principle of the direct current coefficients of adjacent pixel blocks to finish the embedding and blind extraction of the digital watermark in a space domain without carrying out real discrete cosine transform. The invention embeds the color digital watermark image into the color carrier image, has stronger robustness and higher real-time performance on the premise of better watermark invisibility, solves the problems of low running speed and weak robustness of the digital watermark algorithm of the large-capacity color image, and is suitable for the situation of fast and high-efficiency digital media copyright protection.
Description
Technical Field
The invention belongs to the technical field of network space security, and relates to copyright protection of a color digital image as a digital watermark.
Background
With the rapid development of multimedia technology, efficient storage and transmission of digitized information is possible, but at the same time, a series of serious problems such as piracy, infringement, tampering, etc. are also generated. Therefore, copyright protection is becoming a widespread concern for students at home and abroad. The digital watermark technology is an effective digital product copyright protection and data security maintenance technology, and the successful embedding and extraction of the watermark can effectively protect the digital copyright, thereby solving the problem of copyright protection.
According to different processing modes of carrier images, the digital watermarking algorithm is divided into a spatial watermarking algorithm and a frequency domain watermarking algorithm, wherein the spatial watermarking algorithm has the advantages of simplicity and rapidness in operation, but the robustness is poor, and the frequency domain watermarking algorithm has the advantage of strong robustness, but the forward transformation and the reverse transformation of the corresponding transformation domain enable the operation time of the watermarking algorithm to be longer. Therefore, how to combine the advantages of the two, and on the premise of ensuring good invisibility of the digital watermark, designing a digital watermark algorithm with strong robustness and high instantaneity becomes a current urgent problem to be solved.
Disclosure of Invention
The invention aims to provide a blind watermarking method of an airspace color digital image fused with discrete cosine transform, which comprises two specific watermarking embedding and watermarking extracting processes, wherein the watermarking embedding process is described as follows:
the first step: one width is of the sizeM×MIs a color carrier image of (a)HThree layered carrier images separated into red, green and blueH i And imaging each layered carrierH i Is divided intom×mNon-overlapping pixel blocks of (1), whereini=1, 2, 3 represent red, green, blue three layers, respectively;
and a second step of: one width is of the sizeN×NColor watermark image of (a)WDividing into three layered watermark images of red, green and blue, and simultaneously, carrying out key-based on each layered watermark image in order to improve the security of the watermarkKa i Arnold transformation to obtain three scrambled layered watermark imagesW i The method comprises the steps of carrying out a first treatment on the surface of the Layered watermark imageW i Each decimal pixel value in (2) is converted into 8-bit binary number which is sequentially connected into a length of 8N 2 Watermark bit sequence of (a)SW i Whereini=1, 2, 3 represent red, green, blue three layers, respectively;
and a third step of: sequentially from layered carrier images in longitudinal orderH i Is selected adjacent pixel square blockAAndBwhereini=1, 2, 3 represent red, green, blue three layers, respectively;
fourth step: according to formula (1), pixel blocks are directly calculated in the space domainAAndBdirect current coefficient of (2)dc p ;
Wherein,,p=1, 2 respectively represents pixel blocksAAndB,mis the size of the row (column) of the pixel block,f p (x, y) Is a pixel square blockpFirst, thexLine 1yPixel values of the columns;
fifth step: sequentially from a layered watermark sequenceSW i Selecting watermark bits to be embedded inwThe method comprises the steps of carrying out a first treatment on the surface of the Embedding watermark information according to the correlation of DC coefficients of adjacent pixel blocks by using formulas (2) and (3)wObtaining a pixel blockpModified DC coefficientdc p * ,p=1, 2 respectively represents pixel blocksAAndB;
wherein,,win order for the watermark bits to be embedded,avg=(dc 1 +dc 2 )/2,T i is the firstiThe quantization step size of the layer is determined,i=1, 2, 3 represent red, green, blue three layers, respectively;
sixth step: according to equation (4), using the modified pixel valuef p * (x, y) Replacement of original pixel blockpPixel value of corresponding positionf p (x, y) Obtaining pixel square block containing watermarkA * AndB * whereinp=1, 2;
Seventh step: print pixel blockA * AndB * respectively updated to its layered carrier imageH i Corresponding positions in (a), whereini=1, 2, 3, respectively representing three layers of red, green, blue;
eighth step: repeating the third to seventh steps until all watermark information is embedded, thereby obtaining layered carrier image containing watermarkH i * The method comprises the steps of carrying out a first treatment on the surface of the Finally, three layers of the composite materialLayered carrier image for watermarkingH i * Obtaining colour carrier image containing watermarkH * Whereini=1, 2, 3 represent red, green, blue three layers, respectively;
the watermark extraction process is described as follows:
the first step: image processing of a water-bearing carrier by a dimension-reducing treatmentH * Divided into three layered images of red, green and blueH i * The method comprises the steps of carrying out a first treatment on the surface of the At the same time, each watermark is layered into an imageH i * Is divided into the following sizesm×mNon-overlapping pixel blocks of (1), whereini=1, 2, 3 represent red, green, blue three layers, respectively;
and a second step of: layering images sequentially from a watermark in longitudinal orderH i * Is selected to be adjacent to the watermark pixel blockA * AndB * whereini=1, 2, 3 represent red, green, blue three layers, respectively;
and a third step of: directly computing a watermark pixel block in the airspace according to equation (5)A * AndB * direct current coefficient of (2)dc p * ;
Wherein,,mis the size of the row (column) of the pixel block,f p * (x, y) Is a water-printing pixel squarepFirst, thexLine 1yThe pixel values of the columns are used to determine,p=1, 2 respectively representing a water-containing print pixel blockA * AndB * ;
fourth step: according to a water-bearing pixel blockA * AndB * direct current coefficient of (2)dc p * The size relationship between the pixel blocks is printed from the water by using the formula (6)A * AndB * extracting watermark bits from a watermarkw*Whereinp=1, 2;
Fifth step: repeating the second to fourth steps to obtain the extracted binary watermark sequenceSW i * Watermark binary sequenceSW i * Every 8 bits of binary information are divided into a group and converted into decimal pixel values to form a layered watermark image, whereini=1, 2, 3 represent red, green, blue three layers, respectively;
sixth step: key-based hierarchical watermark imageKa i Is transformed by inverse Arnold to obtain watermark image of each layerW i * The method comprises the steps of carrying out a first treatment on the surface of the Combining watermark extracted images of layersW i * Forming a final extracted watermark imageW * Whereini=1, 2, 3 represent three layers of red, green, and blue, respectively.
The method utilizes the correlation rule of the square direct current coefficients of adjacent images to directly obtain the direct current coefficients of discrete cosine transform in a space domain and complete the embedding and blind extraction of the color digital watermark without carrying out real discrete cosine transform; the method not only has better watermark invisibility, but also has stronger watermark robustness and higher algorithm instantaneity.
Drawings
Fig. 1 (a), 1 (b) are two original color carrier images.
Fig. 2 (a) and 2 (b) are two original color watermark images.
Fig. 3 (a) and 3 (b) are watermark images obtained by embedding the watermarks shown in fig. 2 (a) into the carrier images in sequence in fig. 1 (a) and 1 (b), wherein the structural similarity SSIM values are 0.9343 and 0.9351 in sequence, and the peak signal-to-noise ratio PSNR values are 37.4377dB and 37.2156dB in sequence.
Fig. 4 (a) and 4 (b) are watermarks extracted from fig. 3 (a) and 3 (b) in order, and normalized cross-correlation coefficient NC values thereof are 1.0000 and 1.0000, respectively.
Fig. 5 (a), 5 (b), 5 (c), 5 (d), 5 (e), and 5 (f) are watermarks extracted by sequentially subjecting the watermark image shown in fig. 3 (a) to attacks such as JPEG2000 compression (4:1), JPEG compression (70), salt and pepper noise (0.2%), median filtering (3×3), and scaling (75%), and the normalized cross-correlation coefficient NC values thereof are 0.9684, 0.9573, 0.9950, 0.9142, and 0.9626, respectively.
Fig. 6 (a) and 6 (b) show watermark images obtained by embedding the watermarks shown in fig. 2 (b) into the carrier images in sequence in fig. 1 (a) and 1 (b), wherein the structural similarity SSIM values are 0.9346 and 0.9297 in sequence, and the peak signal-to-noise ratio PSNR values are 37.4091dB and 37.2863dB in sequence.
Fig. 7 (a) and 7 (b) are watermarks extracted from fig. 6 (a) and 6 (b) in order, and normalized cross-correlation coefficient NC values thereof are 1.0000 and 1.0000, respectively.
Fig. 8 (a), 8 (b), 8 (c), 8 (d), 8 (e), and 8 (f) are watermarks extracted by sequentially subjecting the watermark image shown in fig. 6 (a) to attacks such as JPEG2000 compression (4:1), JPEG compression (70), salt and pepper noise (0.2%), median filtering (3×3), and scaling (75%), and the normalized cross-correlation coefficient NC values thereof are 0.9676, 0.9526, 0.9956, 0.9018, 0.9673, respectively.
Detailed Description
The invention aims to provide a blind watermarking method of an airspace color digital image fused with discrete cosine transform, which comprises two specific watermarking embedding and watermarking extracting processes, wherein the watermarking embedding process is described as follows:
the first step: a color carrier image with the size of 512 multiplied by 512 is formedHThree layered carrier images separated into red, green and blueH i And imaging each layered carrierH i Divided into 2 x 2 non-overlapping pixel blocks, whereini=1, 2, 3 represent red, green, blue three layers, respectively;
and a second step of: a color watermark image with the size of 32 multiplied by 32 is displayedWDividing into three layered watermark images of red, green and blue, and simultaneously, carrying out key-based on each layered watermark image in order to improve the security of the watermarkKa i Arnold transformation to obtain three scrambled layered watermark imagesW i The method comprises the steps of carrying out a first treatment on the surface of the Layered watermark imageW i Each decimal pixel value of (a) is converted to an 8-bit binary number (e.g., decimal number 225 may be converted to a binary sequence '11100001'), which in turn is connected to a watermark bit sequenceSW i ,SW i Length of 8 x 32 2 =8192, wherei=1, 2, 3 represent red, green, blue three layers, respectively;
and a third step of: sequentially from layered carrier images in longitudinal orderH i Is selected adjacent pixel square blockAAndBwhereini=1, 2, 3 represent red, green, blue three layers, respectively; here, set upi=1, pixel block selected from red layerAIs thatSelected pixel blockBIs->;
Fourth step: according to formula (1), directly calculating neighboring pixel blocks in the spatial domainAAndBdirect current coefficient of (2)dc p ;
Wherein,,p=1, 2 respectively represents pixel blocksAAndB,mis the size of the row (column) of the pixel block,f p (x, y) Is a pixel square blockpFirst, thexLine 1yPixel values of the columns; at this time, the size of the row (column) of the pixel blockm=2, calculated pixel blockADirect current coefficient of (2)dc 1 Pixel block= 448.5000BDirect current coefficient of (2)dc 2 =448.5000;
Fifth step: sequentially from a layered watermark sequenceSW i Is selected fromTo be embedded with watermark bitswThe method comprises the steps of carrying out a first treatment on the surface of the Embedding watermark information according to the correlation of DC coefficients of adjacent pixel blocks by using formulas (2) and (3)wObtaining a pixel blockpModified DC coefficientdc p * ,p=1, 2 respectively represents pixel blocksAAndB;
wherein,,win order for the watermark bits to be embedded,avg=(dc 1 +dc 2 )/2,T i is the firstiThe quantization step size of the layer is determined,i=1, 2, 3 represent red, green, blue three layers, respectively; at this time, the liquid crystal display device,i=1, from watermark sequenceSW 1 Selected watermark bits to be embeddedw=‘0’,T 1 =24.18,avg448.5000, the modified pixel block is obtained according to the formulas (2), (3)A * Direct current coefficient of (2)dc 1 * Pixel block= 448.5000B * Direct current coefficient of (2)dc 2 * =448.5000;
Sixth step: according to equation (4), using the modified pixel valuef p * (x, y) Replacement of original pixel blockpPixel value of corresponding positionf p (x, y) Obtaining pixel square block containing watermarkA * AndB * ;
wherein the method comprises the steps ofp=1, 2; at this time, the obtained water-containing print pixel blockA * Is thatWater-printing pixel squareB * Is->;
Seventh step: print pixel blockA * AndB * respectively updated to its layered carrier imageH i Corresponding positions in (a), whereini=1, 2, 3, respectively representing three layers of red, green, blue; at this time, the liquid crystal display device,i=1, watered pixel blockA * AndB * is updated to its layered carrier image respectivelyH 1 Corresponding positions in (a);
eighth step: repeating the third to seventh steps until all watermark information is embedded, thereby obtaining layered carrier image containing watermarkH i * The method comprises the steps of carrying out a first treatment on the surface of the Finally, combining three layered watermark-containing layered carrier imagesH i * Obtaining colour carrier image containing watermarkH * Whereini=1, 2, 3 represent red, green, blue three layers, respectively;
the watermark extraction process is described as follows:
the first step: image processing of a water-bearing carrier by a dimension-reducing treatmentH * Divided into three layered images of red, green and blueH i * The method comprises the steps of carrying out a first treatment on the surface of the At the same time, each watermark is layered into an imageH i * Divided into non-overlapping pixel blocks of size 2 x 2, whereini=1, 2, 3 represent red, green, blue three layers, respectively;
and a second step of: layering images sequentially from a watermark in longitudinal orderH i * Is selected to be adjacent to the watermark pixel blockA * AndB * whereini=1, 2, 3 represent red, green, blue three layers, respectively; at this time, set upi=1, a block of water-containing pixels selected from the red layerA * Is thatPixel block containing watermarkB * Is->;
And a third step of: direct computation of watermark pixel blocks in airspace according to equation (5)A * AndB * direct current coefficient of (2)dc p * ;
Wherein,,mis the size of the row (column) of the pixel block,f p * (x, y) Is a water-printing pixel squarepFirst, thexLine 1yThe pixel values of the columns are used to determine,p=1, 2 respectively representing a water-containing print pixel blockA * AndB * the method comprises the steps of carrying out a first treatment on the surface of the At this time, the size of the row (column) of the pixel blockm=2, calculated pixel blockA * Direct current coefficient of (2)dc 1 * Pixel block= 448.5000B * Direct current coefficient of (2)dc 2 * =448.5000;
Fourth step: according to a water-bearing pixel blockA * AndB * direct current coefficient of (2)dc p * The size relationship between the pixel blocks is printed from the water by using the formula (6)A * AndB * extracting watermark bits from a watermarkw * Whereinp=1, 2;
At this time, the liquid crystal display device,dc 1 * -dc 2 * =0, the extracted watermark bits are obtained according to equation (6)w * =‘0’;
Fifth step: repeatedly executing the second to the first steps of the processFour steps, obtaining the extracted binary watermark sequenceSW i * Watermark binary sequenceSW i * Every 8 bits of binary information is divided into a group and converted into decimal pixel values to form a layered watermark image, whereini=1, 2, 3 represent red, green, blue three layers, respectively;
sixth step: key-based hierarchical watermark imageKa i Is transformed by inverse Arnold to obtain watermark image of each layerW i * The method comprises the steps of carrying out a first treatment on the surface of the At the same time, the watermark extraction images of the layers are combinedW i * Forming a final extracted watermark imageW * Whereini=1, 2, 3 represent three layers of red, green, and blue, respectively.
The method has better watermark invisibility, stronger robustness and higher algorithm instantaneity, and is suitable for copyright protection of color digital images as watermarks.
The invention has the effect of verification
To prove the effectiveness of the present invention, two standard 24-bit color images of 512×512 size as shown in fig. 1 (a) and 1 (b) were selected as carrier images, and two 24-bit color images of 32×32 size as shown in fig. 2 (a) and 2 (b) were used as digital watermarks, respectively, for verification.
Fig. 3 (a) and 3 (b) are watermark images obtained by embedding the watermarks shown in fig. 2 (a) into the carrier images in sequence in fig. 1 (a) and 1 (b), wherein the structural similarity SSIM values are 0.9343 and 0.9351 in sequence, and the peak signal-to-noise ratio PSNR values are 37.4377dB and 37.2156dB in sequence; fig. 4 (a) and 4 (b) are watermarks extracted from fig. 3 (a) and 3 (b) in sequence, and normalized cross-correlation coefficient NC values thereof are 1.0000 and 1.0000, respectively; fig. 5 (a), 5 (b), 5 (c), 5 (d), and 5 (e) are watermarks extracted by sequentially subjecting the watermark image shown in fig. 3 (a) to attacks such as JPEG2000 compression (4:1), JPEG compression (70), salt and pepper noise (0.2%), median filtering (3×3), and scaling (75%), and the normalized cross-correlation coefficient NC values thereof are 0.9684, 0.9573, 0.9950, 0.9142, 0.9626, respectively.
Fig. 6 (a) and 6 (b) are watermark images obtained by embedding the watermarks shown in fig. 2 (b) into the carrier images in sequence in fig. 1 (a) and 1 (b), wherein the structural similarity SSIM values are 0.9346 and 0.9297 in sequence, and the peak signal-to-noise ratio PSNR values are 37.4091dB and 37.2863dB in sequence; fig. 7 (a) and 7 (b) are watermarks extracted from fig. 6 (a) and 6 (b) in order, and normalized cross-correlation coefficient NC values thereof are 1.0000 and 1.0000, respectively; fig. 8 (a), 8 (b), 8 (c), 8 (d), 8 (e), and 8 (f) are watermarks extracted after attacks such as JPEG2000 compression (4:1), JPEG compression (70), salt and pepper noise (0.2%), median filtering (3×3), and scaling (4:1) are sequentially performed on the watermark image shown in fig. 6 (a), and normalized cross-correlation coefficient NC values thereof are 0.9676, 0.9526, 0.9956, 0.9018, 0.9673, respectively.
The algorithm runs on the platform 2.60GHZ CPU,4.00GB RAM,Win10,MATLAB (R2017 a) for nearly ten thousand times, the average embedding time of the digital watermark is 0.605655 seconds, the average extraction time is 0.255208 seconds, and the total time is 0.860863 seconds.
In summary, the embedded digital image watermark has higher invisibility, and meets the invisibility requirement of a watermark algorithm; meanwhile, the digital image watermark extracted from various attacked images has better authenticability and higher NC value, which indicates that the method has stronger robustness; in addition, the average running total time of the algorithm is less than 1 second, and the requirement of rapid copyright protection of multimedia big data is met.
Claims (1)
1. A blind watermarking method of a space domain color digital image fused with discrete cosine transform is characterized in that: the method is realized through a specific watermark embedding process and a watermark extraction process, and the watermark embedding process is described as follows:
the first step: one width is of the sizeM×MIs a color carrier image of (a)HThree layered carrier images separated into red, green and blueH i And imaging each layered carrierH i Is divided intom×mNon-overlapping pixel blocks of (1), whereini=1, 2, 3 represent red, green, blue three layers, respectively;
and a second step of: one width is of the sizeN×NColor watermark image of (a)WDividing into three layered watermark images of red, green and blue, and simultaneously, carrying out key-based on each layered watermark image in order to improve the security of the watermarkKa i Arnold transformation to obtain three scrambled layered watermark imagesW i The method comprises the steps of carrying out a first treatment on the surface of the Layered watermark imageW i Each decimal pixel value in (2) is converted into 8-bit binary number which is sequentially connected into a length of 8N 2 Watermark bit sequence of (a)SW i Whereini=1, 2, 3 represent red, green, blue three layers, respectively;
and a third step of: sequentially from layered carrier images in longitudinal orderH i Is selected adjacent pixel square blockAAndBwhereini=1, 2, 3 represent red, green, blue three layers, respectively;
fourth step: according to formula (1), pixel blocks are directly calculated in the space domainAAndBdirect current coefficient of (2)dc p ;
Wherein,,p=1, 2 respectively represents pixel blocksAAndB,mis the size of the row (column) of the pixel block,f p (x, y) Is a pixel square blockpFirst, thexLine 1yPixel values of the columns;
fifth step: sequentially from a layered watermark sequenceSW i Selecting watermark bits to be embedded inwThe method comprises the steps of carrying out a first treatment on the surface of the Embedding watermark information according to the correlation of DC coefficients of adjacent pixel blocks by using formulas (2) and (3)wObtaining a pixel blockpModified DC coefficientdc p * ,p=1, 2 respectively represents pixel blocksAAndB;
wherein,,win order for the watermark bits to be embedded,avg=(dc 1 +dc 2 )/2,T i is the firstiThe quantization step size of the layer is determined,i=1, 2, 3 represent red, green, blue three layers, respectively;
sixth step: according to equation (4), using the modified pixel valuef p * (x, y) Replacement of original pixel blockpPixel value of corresponding positionf p (x, y) Obtaining pixel square block containing watermarkA * AndB * whereinp=1, 2;
Seventh step: print pixel blockA * AndB * respectively updated to its layered carrier imageH i Corresponding positions in (a), whereini=1, 2, 3, respectively representing three layers of red, green, blue;
eighth step: repeating the third to seventh steps until all watermark information is embedded, thereby obtaining layered carrier image containing watermarkH i * The method comprises the steps of carrying out a first treatment on the surface of the Finally, combining three layered watermark-containing layered carrier imagesH i * Obtaining colour carrier image containing watermarkH * Whereini=1, 2, 3 represent red, green, blue three layers, respectively;
the watermark extraction process is described as follows:
the first step: image processing of a water-bearing carrier by a dimension-reducing treatmentH * Divided into three layered images of red, green and blueH i * The method comprises the steps of carrying out a first treatment on the surface of the At the same time, each watermark is layered into an imageH i * Is divided into the following sizesm×mNon-overlapping pixel blocks of (1), whereini=1, 2, 3 represent red, green, blue three layers, respectively;
and a second step of: layering images sequentially from a watermark in longitudinal orderH i * Is selected to be adjacent to the watermark pixel blockA * AndB * whereini=1, 2, 3 represent red, green, blue three layers, respectively;
and a third step of: directly computing a watermark pixel block in the airspace according to equation (5)A * AndB * direct current coefficient of (2)dc p * ;
Wherein,,mis the size of the row (column) of the pixel block,f p * (x, y) Is a water-printing pixel squarepFirst, thexLine 1yThe pixel values of the columns are used to determine,p=1, 2 respectively representing a water-containing print pixel blockA * AndB * ;
fourth step: according to a water-bearing pixel blockA * AndB * direct current coefficient of (2)dc p * The size relationship between the pixel blocks is printed from the water by using the formula (6)A * AndB * extracting watermark bits from a watermarkw*Whereinp=1, 2;
Fifth step: repeating the second to fourth steps to obtain the extracted binary watermark sequenceSW i * Watermark binary sequenceSW i * Every 8 bits of binary information is divided into groups and converted intoDecimal pixel values forming a layered watermark image, whereini=1, 2, 3 represent red, green, blue three layers, respectively;
sixth step: key-based hierarchical watermark imageKa i Is transformed by inverse Arnold to obtain watermark image of each layerW i * The method comprises the steps of carrying out a first treatment on the surface of the Combining watermark extracted images of layersW i * Forming a final extracted watermark imageW * Whereini=1, 2, 3 represent three layers of red, green, and blue, respectively.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910889463.0A CN110570345B (en) | 2019-09-20 | 2019-09-20 | Blind watermarking method for airspace color digital image fused with discrete cosine transform |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910889463.0A CN110570345B (en) | 2019-09-20 | 2019-09-20 | Blind watermarking method for airspace color digital image fused with discrete cosine transform |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110570345A CN110570345A (en) | 2019-12-13 |
CN110570345B true CN110570345B (en) | 2023-07-14 |
Family
ID=68781323
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910889463.0A Active CN110570345B (en) | 2019-09-20 | 2019-09-20 | Blind watermarking method for airspace color digital image fused with discrete cosine transform |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110570345B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111199508B (en) * | 2020-01-03 | 2022-05-17 | 鲁东大学 | Spatial domain color digital image blind watermarking method fusing DCT and DHT |
CN111242828B (en) * | 2020-01-16 | 2023-02-03 | 鲁东大学 | Spatial domain color digital image blind watermarking method fused with discrete Fourier transform |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105761199A (en) * | 2016-03-25 | 2016-07-13 | 鲁东大学 | Robustness color image watermarking method based on combination domain |
CN108648134A (en) * | 2018-05-21 | 2018-10-12 | 鲁东大学 | A kind of spatial domain color digital image blind watermark method of fusion wavelet transform |
CN109035129A (en) * | 2018-08-15 | 2018-12-18 | 鲁东大学 | A kind of color digital image blind watermark method based on two-dimensional discrete sine transform |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2419763A (en) * | 2004-11-01 | 2006-05-03 | Sony Uk Ltd | Image watermarking using weighted code word coefficients |
US10296999B2 (en) * | 2017-02-13 | 2019-05-21 | Macau University Of Science And Technology | Methods and apparatus for color image watermarking |
-
2019
- 2019-09-20 CN CN201910889463.0A patent/CN110570345B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105761199A (en) * | 2016-03-25 | 2016-07-13 | 鲁东大学 | Robustness color image watermarking method based on combination domain |
CN108648134A (en) * | 2018-05-21 | 2018-10-12 | 鲁东大学 | A kind of spatial domain color digital image blind watermark method of fusion wavelet transform |
CN109035129A (en) * | 2018-08-15 | 2018-12-18 | 鲁东大学 | A kind of color digital image blind watermark method based on two-dimensional discrete sine transform |
Non-Patent Citations (1)
Title |
---|
基于四元数的双彩色图像盲水印算法;姚涛;苏庆堂;阙大顺;;计算机工程(第01期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110570345A (en) | 2019-12-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109035129B (en) | Color digital image blind watermarking method based on two-dimensional discrete sine transformation | |
CN108648134B (en) | Spatial domain color digital image blind watermarking method fusing discrete wavelet transform | |
CN110390621B (en) | DCT domain color digital image blind watermarking method based on variable step length | |
CN107895340B (en) | Space-frequency domain color digital image blind watermarking method integrating QR decomposition | |
CN111199508B (en) | Spatial domain color digital image blind watermarking method fusing DCT and DHT | |
CN109102454B (en) | Color QR code digital blind watermarking method integrating fast Fourier transform | |
CN110349073B (en) | Four-system color digital image blind watermarking method based on Schur decomposition | |
CN107274333B (en) | Spatial domain color image blind watermarking method fusing singular value decomposition | |
CN110570345B (en) | Blind watermarking method for airspace color digital image fused with discrete cosine transform | |
CN106204410B (en) | A kind of novel digital watermark method decomposed based on matrix Schur | |
CN112488904B (en) | Matrix singular value decomposition-based quaternion color digital image blind watermarking method | |
CN105023236B (en) | A kind of dual-color image blind watermarking method based on LU Decomposition | |
CN106991636B (en) | Airspace color image blind watermarking method fusing approximate Schur decomposition | |
CN109829846B (en) | Digital image blind watermarking method based on two-dimensional discrete cosine transform | |
CN110415155B (en) | Blind watermarking method for airspace color image fused with haar transformation | |
CN106157233B (en) | A kind of good digital water mark method decomposed based on matrix Schur | |
CN111242828B (en) | Spatial domain color digital image blind watermarking method fused with discrete Fourier transform | |
CN116993567A (en) | Frequency domain blind watermarking method based on Hadamard transform and teaching optimization algorithm | |
CN106169171B (en) | The good digital water mark method decomposed based on Hessenberg | |
CN110415154B (en) | Haer transformation-based quaternary color digital image blind watermarking method | |
CN112488903B (en) | Spatial domain color digital image blind watermarking method fusing multilevel discrete Fourier transform | |
Tomar et al. | A statistical comparison of digital image watermarking techniques | |
CN113191932B (en) | Spatial domain color digital image blind watermarking method fusing discrete Chebyshev transformation | |
CN113191933B (en) | Color digital image blind watermarking method based on approximate Hadamard transform | |
CN117495647A (en) | Color image digital watermarking method based on ULV decomposition |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB03 | Change of inventor or designer information | ||
CB03 | Change of inventor or designer information |
Inventor after: Wang Huanying Inventor after: Yuan Zihan Inventor after: Liu Decheng Inventor after: Zhang Xueting Inventor after: Su Qingtang Inventor before: Su Qingtang Inventor before: Yuan Zihan Inventor before: Liu Decheng Inventor before: Zhang Xueting Inventor before: Wang Huanying |
|
GR01 | Patent grant | ||
GR01 | Patent grant |