CN102880997B - Method for embedding watermark image - Google Patents
Method for embedding watermark image Download PDFInfo
- Publication number
- CN102880997B CN102880997B CN201210281560.XA CN201210281560A CN102880997B CN 102880997 B CN102880997 B CN 102880997B CN 201210281560 A CN201210281560 A CN 201210281560A CN 102880997 B CN102880997 B CN 102880997B
- Authority
- CN
- China
- Prior art keywords
- image
- matrix
- image blocks
- watermark
- data
- 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 30
- 239000011159 matrix material Substances 0.000 claims abstract description 80
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 6
- 230000009466 transformation Effects 0.000 claims abstract description 5
- 230000015572 biosynthetic process Effects 0.000 claims description 10
- 230000000739 chaotic effect Effects 0.000 claims description 4
- 230000000903 blocking effect Effects 0.000 abstract 1
- 230000006835 compression Effects 0.000 description 7
- 238000007906 compression Methods 0.000 description 7
- 238000001914 filtration Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 238000010008 shearing Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007850 degeneration Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Landscapes
- Editing Of Facsimile Originals (AREA)
- Image Processing (AREA)
Abstract
The invention relates to a method for embedding a grayscale watermark image. The method is suitable for solving the problem of conflict between the transparency and robustness of a watermarking algorithm. The method comprises the following steps of: blocking an original carrier image, performing discrete cosine transformation (DCT), dividing a coefficient into a low-frequency data block, an intermediate-frequency data block and a high-frequency data block by utilizing a template, performing singular value decomposition (SVD) after a new matrix is formed, calculating the energy ratio of the blocks, and determining an optimal block position into which the watermark image is suitable to be embedded; and adaptively embedding grayscale watermark information into a singular value of the low-frequency block according to the energy ratio.
Description
The application is divisional application, the application number of original application: 201110139206.9, the applying date: 2011-5-26, invention and created name: the embedding of watermarking images and extracting method.
Technical field
The present invention relates to the technical field of image procossing, specifically a kind of embedding grammar of watermarking images.
Background technology
Along with the develop rapidly of Digital object identifier and technique of internet, the safety problem of media product becomes the previous quite important and challenging research topic of order.Digital watermark technology is a kind of potential effective ways realizing copyright protection of digital product; become a study hotspot of information security field at present; also be a kind of important branch of Information hiding research field; digital watermarking is arisen at the historic moment just in this context, and popular direction, the forward position becoming academic circles at present research.
In view of current most of trade mark or copyright marking all represent with gray level image or coloured image, the research therefore how embedding Gray-level Watermarking image has very important practical value, also more meaningful, more realistic demand.Due to the restriction of Information Hiding Capacity, embedding difficulty is high to adopt gray level image to carry out embedding not only as image watermark, and is difficult to realize Blind extracting, while the robustness improving embed watermark, also reduce its not sentience, so the research of this respect is not also a lot.
The transparency and robustness are digital watermark technology two fundamental characteristics.If general embedding is containing the image transparency relatively good (namely PSNR is larger) of watermark, robustness poor (the related coefficient NC of detection is smaller), they are conflicting, the transparency and robustness all with the embedded location of watermark and embedment strength closely bound up.The present invention is exactly that combining image energy spectrometer embeds watermark information strong as far as possible in place in order to solve this technical barrier, improves again and detect related coefficient NC while namely improving Y-PSNR under same embedment strength.
Digital watermarking can be divided into robust watermarking, fragile watermark and semi-fragile watermarking by the effect according to digital watermarking.
(1) fundamental purpose of robust watermarking is the copyright protecting copyright; it requires that the watermark after embedding can stand various conventional signal processing operations; comprise unintentionally or the process of malice, as lossy compression method, filtering, level and smooth, reduce, geometry deformation etc.Robust watermarking is after various process, as long as host image is not corrupted to not spendable degree.All should be able to detect.Therefore the robustness of such watermark requires higher.
(2) fragile watermark is mainly for the protection of the integrality of copyright, the true and false of discriminating digit works.
(3) semi-fragile watermarking requires the useful signal processing operations that can resist to a certain degree, and as JPEG compression etc., than fragile watermark, slightly robust is a little in this watermark.
Summary of the invention
The technical matters that the present invention first will solve is to provide a kind of embedding grammar of Gray-level Watermarking image.
The technical matters of its secondary solution of the present invention solves the transparency in watermarking algorithm to contradict problem with robustness.
For solving the problems of the technologies described above, the invention provides the embedding grammar of watermarking images, comprising the steps:
(1) pre-service of watermarking images: watermarking images is carried out chaos sequence scramble, and the image sequence this watermarking images being deformed into 1 row n row, i.e. n image information;
(2) selection of embedded location: be I (M by matrix by initial carrier image, N) 8 × 8 image blocks (namely length and width are all the rectangular image block of 8 pixels) do not covered mutually of corresponding number are divided into, DCT(discrete cosine transform is carried out to each 8 × 8 image blocks) conversion, low frequency corresponding in each 8 × 8 image blocks and high-frequency data is taken out according to 8 × 8 templates (namely length and width are all the rectangular image template of 8 pixels), and described low frequency and high-frequency data are formed the matrix A of 3 × 3 and the matrix B of 6 × 6 respectively, to matrix A, B carries out svd and asks its matrix norm, solve the energy Ratios g of each 8 × 8 image blocks, energy Ratios
wherein: s
highfor the first singular value of described matrix B, s
lowfor the first singular value of described matrix A, finding out with energy Ratios is corresponding 8 × 8 image blocks of a front n maximal value,
(3) data selection of watermark embedment: watermark information and a described n image information, select correspondingly successively to embed adaptively in the first value of the svd of the described matrix A at the n arranged from big to small by described energy Ratios g 8 × 8 image blocks, embedment strength regulates automatically with the energy Ratios of image block;
(4) DCT inverse transformation is carried out to each 8 × 8 image blocks of embed watermark information, obtain image with watermarked information.
When carrying out the data selection of watermark embedment of above-mentioned steps (3): to the method that the matrix A of n 8 × 8 images carries out svd be one by one: A=USV
t, U ∈ R in formula
m × M, S ∈ R
m × N, V ∈ R
n × N, U and V is orthogonal matrix, and S is diagonal matrix; Take out the first singular value of S diagonal matrix and embed watermark, embedding formula
i=1 Λ n; Embedment strength K, with block energy adjustment more automatic than g, embeds adaptively described in realizing, and its formula is K
i=α+g
i, i=1 Λ m, in formula, α is element factor, and its value condition need ensure that embedding image Y-PSNR (PSNR) with watermarked information is greater than 36db; M is the total block data that original image is divided into described 8 × 8 image blocks; New singular value is substituted into A '=US ' V
t, obtain low-frequency data A ' with watermarked information, and data sent back to the correspondence position in corresponding 8 × 8 image blocks.
Present invention also offers a kind of embedding and extracting method of watermarking images, wherein, the embedding grammar of described watermarking images comprises the steps:
(1) pre-service of watermarking images: watermarking images is carried out chaos sequence scramble, and the image sequence this watermarking images being deformed into 1 row n row, i.e. n image information;
(2) selection of embedded location: be I (M by matrix by initial carrier image, N) 8 × 8 image blocks do not covered mutually of corresponding number are divided into, dct transform is carried out to each 8 × 8 image blocks, low frequency corresponding in each 8 × 8 image blocks and high-frequency data is taken out according to 8 × 8 templates, and described low frequency and high-frequency data are formed the matrix A of 3 × 3 and the matrix B of 6 × 6 respectively, matrix A, B are carried out to svd and ask its matrix norm, solves the energy Ratios g of each 8 × 8 image blocks, energy Ratios
wherein: s
highfor the first singular value of described matrix B, s
lowfor the first singular value of described matrix A; Find out and 8 × 8 image blocks of energy Ratios corresponding to a front n maximal value;
(3) data selection of watermark embedment: watermark information and a described n image information, select correspondingly successively to embed adaptively in the first value of the svd of the described matrix A at n 8 × 8 image blocks by the large minispread of described energy Ratios g, embedment strength regulates automatically with the energy Ratios of image block;
(4) DCT inverse transformation is carried out to each 8 × 8 image blocks of embed watermark information, obtain image with watermarked information;
When carrying out the data selection of watermark embedment of above-mentioned steps (3): to the method that the matrix A of n 8 × 8 images carries out svd be one by one: A=USV
t, U ∈ R in formula
m × M, S ∈ R
m × N, V ∈ R
n × N, U and V is orthogonal matrix, and S is diagonal matrix;
Take out the first singular value of S diagonal matrix and embed watermark, embedding formula
(i=1 Λ n); Embedment strength K, with block energy adjustment more automatic than g, embeds adaptively described in realizing, and its formula is K
i=α+g
i, i=1 Λ m, in formula, α is element factor, and its value condition need ensure that embedding Y-PSNR with watermarked information is greater than 36db; M is the total block data that original image is divided into described 8 × 8 image blocks;
New singular value is substituted into A '=US ' V
t, obtain low-frequency data A ' with watermarked information, and data sent back to the correspondence position in corresponding 8 × 8 image blocks;
The extracting method of described watermarking images, it comprises the steps:
A () is I (M to matrix, N) initial carrier image and matrix are I ' (M, N) 8 × 8 image blocks do not covered mutually being divided into corresponding number containing watermarking images respectively to be checked, and dct transform is carried out to each 8 × 8 image blocks, from each 8 × 8 image blocks, corresponding low frequency and high-frequency data is taken out according to 8 × 8 templates, and the matrix A of 3 × 3 and the matrix B of 6 × 6 are formed respectively to the low frequency in each 8 × 8 image blocks in initial carrier image and high-frequency data, to form respectively containing the low frequency in each 8 × 8 image blocks of watermarking images and high-frequency data the matrix A of matrix 3 × 3 ' with 6 × 6 matrix B ', to described matrix A, A ' carries out svd respectively: A=USV
t, A '=U ' S ' V '
t,
B () solves the watermark information changing value Δ W in each 8 × 8 image blocks
j=S '
j(1,1)-S
j(1,1), j=1 Λ m, m are the total block data that original image is divided into described 8 × 8 image blocks;
C () is in the m that is partitioned in I (M, N) initial carrier image 8 × 8 image blocks from matrix, finding out with energy Ratios is corresponding 8 × 8 image blocks of a front n maximal value, namely finds out individual 8 × 8 image blocks of n;
D the energy Ratios of one of the energy Ratios of one of individual 8 × 8 image blocks of the n described in () and 8 × 8 image blocks of described m is designated as g respectively
n(i) and g
m(j), i=1 Λ n, j=1 Λ m; According to:
If g
n(i) ≡ g
m(j), then image intensity value
draw n image information, K is embedment strength;
E (), by the energy Ratios from big to small sequential of a described n image information according to 8 × 8 image blocks at each image information place, draws 1 row n row image sequence; This image sequence is deformed into the Gray-level Watermarking image of scramble;
F () goes out watermarking images according to the inverse operation of Chaotic Scrambling algorithm from the Gray-level Watermarking Postprocessing technique of above-mentioned scramble;
In described step (a): take out the data A(9 of 1 × 9 from the low frequency and high frequency region coefficient of 8 × 8 template-setups), the data B(36 of 1 × 36) array, this array is deformed into square formation A, the A ' of 3 × 3, again SVD decomposition operation is carried out to square formation A, A ', take out the first value in singular value diagonal matrix.
In described step (c): be I (M to initial carrier image array, N), to be checked is I ' (M containing watermarking images, N), be divided into 8 × 8 image blocks do not covered mutually of corresponding number respectively, and carry out dct transform to each piece, from the low frequency and high frequency region coefficient of 8 × 8 template-setups, take out the data A(9 of 1 × 9), the data B(36 of 1 × 36) array, this array is deformed into the square formation A of 3 × 3,6 × 6 square formation B, again SVD decomposition operation is carried out to square formation A, B, determine energy Ratios
finding out with energy Ratios is corresponding 8 × 8 image blocks of a front n maximal value.
The present invention has positive effect: the embedding grammar of (1) self-adaptation Gray-level Watermarking of the present invention image, utilizes Chaotic Scrambling to remove the correlativity of original watermark and scramble watermark, improves the security of watermark.By making dct transform to initial carrier image I piecemeal and template partitioned matrix SVD decomposes, matrix norm, computing block energy Ratios, looks for the best embedding position of piecemeal, improves the transparency of watermark.In telescopiny, each piecemeal embedment strength regulates automatically with energy Ratios, substantially increases the robustness of watermark.In view of most of trade mark or copyright information are all represent with gray level image or coloured image, utilize the present invention to embed significant Gray-level Watermarking information in initial carrier image I, there is important practical value, more realistic demand.(2) the present invention makes svd in the DCT domain of 8 × 8 image blocks and represents image energy by its matrix norm, utilizes block image energy spectrometer to find method, determines the piecemeal of applicable embed watermark, improve the invisibility of watermark.(3) embedment strength of the present invention's each piecemeal in telescopiny regulates automatically with the energy Ratios g of piecemeal, greatly strengthen the robustness of watermark.(4) piecemeal is divided into basic, normal, high three parts according to 8 × 8 templates by the present invention, by the first value of watermark embedment 3 × 3 Singular Value Decomposition Using, and the robustness of this watermark also improved.
Accompanying drawing explanation
In order to make content of the present invention be more likely to be clearly understood, below basis specific embodiment and by reference to the accompanying drawings, the present invention is further detailed explanation, wherein
Fig. 1 is the watermarking images in embodiment.
Fig. 2 is 8 × 8 templates.
Fig. 3 is the embedding process flow diagram of watermarking images;
Fig. 4 is the extraction process flow diagram of watermarking images;
Image containing watermark is carried out the result schematic diagram that quality factor is the compression of 90,75,50,30,10 by Fig. 5 respectively;
Fig. 6 be make an uproar to adding respectively containing the original image of watermark, the result schematic diagram of Filtering Attacks experiment;
Fig. 7 is to the result schematic diagram sheared respectively containing the original image of watermark, rotate, convergent-divergent attacks experiment.
Embodiment
(embodiment 1)
See Fig. 1-3, the embedded object of the watermarking images of the present embodiment is a matrix is I (M, N) initial carrier image (M=N=512, namely length and width is respectively the image of 512 pixels), its form is unpressed more than 8 high accuracy grey scale bitmaps (bmp or tiff), also can choose the image of JPG form.Gray-level Watermarking image generally chooses trade mark or copyright marking, and the total n of pixel requires less than carrier image (n < M/8*N/8), and in embodiment, selected pixels size is the Gray-level Watermarking of 40 × 49.
Carrier image embedded location is determined and to embed Gray-level Watermarking information Step as follows:
(1) pre-service of watermarking images.For ensureing the security of watermark, watermarking images is carried out some algorithms (as chaos sequence scramble method etc.) scramble, after scramble, image is designated as W, and image 40 × 49 is deformed into 1 row 1960 row image sequence one-dimension array.
(2) embedded location is selected.By initial carrier image I (512,512) 8 × 8 image blocks do not covered mutually of corresponding number are divided into, amount to m=4096 block, dct transform is carried out to each image block, take out low frequency corresponding in each 8 × 8 image blocks and high-frequency data according to 8 × 8 templates, and described low frequency and high-frequency data are formed respectively the B matrix of the A and 6 × 6 of 3 × 3.
One width digital picture can regard the matrix be made up of many real number items as from linear algebra angle.With A ∈ R
m × Nrepresent an image array, its svd is the product of 3 matrixes.
Namely
Or A=USV
t(1)
U ∈ R in formula
m × M, S ∈ R
m × N, V ∈ R
n × N, U and V is orthogonal matrix, and S is diagonal matrix.
S
r=diag(s
1,s
2,Λ,s
r) (2)
On S diagonal line, element meets s
1>=s
2>=Λ>=s
r> s
r+1=Λ=s
n=0, the order of element to be 0, r be matrix A on off-diagonal, formula (1) is exactly the svd of matrix A.
Norm according to matrix norm definition A is
By A=USV
tsubstitution above formula obtains
Because U and V is orthogonal matrix, its inner product is 1, so above formula has again
Represent image energy by image norm, from above-mentioned theory is derived, known image energy is the evolution of singular value quadratic sum, again because that play a decisive role after digital picture svd is s
1(s
1much larger than other singular value), that is other singular value is become 0, only retain s
1the image of rear reconstruct and original image do not have too large degeneration, and this is a feature of digital picture svd.Can be reduced to again accordingly
derive according to above-mentioned theory, when the energy Ratios g of calculating 8 × 8 image block, as long as carry out svd to A, B matrix, solve each image block energy Ratios
(s
highfor the first singular value of described matrix B, s
lowfirst singular value for described matrix A), finding out with energy Ratios is corresponding 8 × 8 image blocks of front 1960 maximal values, can go out which image block and be suitable for embed watermark information by quick position.
(3) data selection of watermark embedment: determine in which image block embed after, watermarking images information is selected to touch in the first value of the svd of the low frequency range (i.e. matrix A) of plate 8 × 8 to embed adaptively, embedment strength regulates than automatic with Block Energy, object is for the image after embedding has comparatively strong robustness, and can resist certain geometric attack.Concrete implementation step is as follows:
1. svd is carried out to matrix A: A=USV
t, take out the first singular value embed watermark of S diagonal matrix, embedding formula
i=1 Λ n; Embedment strength K is with block energy than automatically regulating, and its formula is K
i=α+g
i, i=1 Λ m, in formula, α is element factor, and its value condition need ensure to embed and be greater than 36db(containing the image PSNR of watermark and be less than 36db human eye and almost can perceive image change), generally get 0.1.G is image block energy Ratios, and m is the sum that original image is divided into 8 × 8 pieces, n=1960, m=4096 in the present embodiment.
2. new singular value is substituted into A '=US ' V
t, obtain low-frequency data A ' with watermarked information, and send data back in block correspondence position.
(4) DCT inverse transformation is carried out to each 8 × 8 image blocks of the image of embed watermark, obtain image I ' with watermarked information (512,512).
(embodiment 2)
Recovery and the authenticating step of Gray-level Watermarking image are as follows:
A () is I (512 to initial carrier image array, 512), to be checked is I ' (512 containing watermarking images, 512), be divided into 8 × 8 image blocks do not covered mutually of corresponding number respectively, and dct transform is carried out to each image block, from each 8 × 8 image blocks, corresponding low frequency and high-frequency data is taken out according to 8 × 8 templates, and the matrix A of 3 × 3 and the matrix B of 6 × 6 are formed respectively to the low frequency in each 8 × 8 image blocks in initial carrier image and high-frequency data, to form respectively containing the low frequency in each 8 × 8 image blocks of watermarking images and high-frequency data the matrix A of matrix 3 × 3 ' with 6 × 6 matrix B ', to described matrix A, A ' carries out svd respectively: A=USV
t, A '=U ' S ' V '
t,
B () solves the watermark information changing value Δ W in each 8 × 8 image blocks
j=S '
j(1,1)-S
j(1,1), j=1 Λ m, m are the total block data that original image is divided into described 8 × 8 image blocks;
C initial carrier image array is I (512 by (), 512) 8 × 8 image blocks do not covered mutually of corresponding number are divided into, and carry out dct transform to each piece, from the low frequency and high frequency region coefficient of 8 × 8 template-setups, take out the data A(9 of 1 × 9), the data B(36 of 1 × 36) array, this array is deformed into the square formation A of 3 × 3,6 × 6 square formation B, again SVD decomposition operation is carried out to square formation A, B, determine energy Ratios
find out n maximal value corresponding blocks before energy Ratios.
D the energy Ratios of one of the energy Ratios of one of individual 8 × 8 image blocks of the n described in () and 8 × 8 image blocks of described m is designated as g respectively
n(i) and g
m(j), i=1 Λ n, j=1 Λ m.Because of in 8 × 8 image blocks that information insertion n maximal value before with energy Ratios being is corresponding, the recovery of its watermark is carried out according to following formula:
I=1 Λ n, j=1 Λ m, K in formula
j=α+g
j;
E (), by the energy Ratios from big to small sequential of a described n image information according to 8 × 8 image blocks at each image information place, draws 1 row 1960 row image sequence; This image sequence is deformed into 40 × 49 Gray-level Watermarking images of scramble;
F the inverse operation of () foundation Chaotic Scrambling algorithm goes out the watermarking images of 40 × 49 from 40 × 49 Gray-level Watermarking Postprocessing technique of above-mentioned scramble.
Below simulating, verifying experiment and explanation:
On the block position of all multi-watermark embedding methods is selected (i.e. the process of the watermark transparency), all adopt the mode of computing block variance to find, the present invention have employed Energy Analysis for High in the searching of block position, and experimental result is as table 1:
Intensity | 0.1 | 0.15 | 0.2 | 0.25 |
PSNR1 | 42.4793 | 38.9575 | 36.4587 | 34.5205 |
PSNR2 | 42.5254 | 39.0036 | 36.5048 | 34.5666 |
NC1 | 0.9809 | 0.9847 | 0.9857 | 0.9864 |
NC2 | 0.9922 | 0.9967 | 0.9983 | 0.9995 |
Table 1
(note: the image Y-PSNR that spatial domain variance searching block method embeds containing watermark is PSNR1, and the watermark related coefficient detected is NC2; Frequency domain energy block finding method is designated as PSNR2 and NC2)
Can find out PSNR2>PSNR1, NC2>NC1 each embedment strength from table 1 data, the embedded location that the embedded location that illustration method 2 is selected all wants ratio method 1 to select in invisibility and robustness is more reasonable, superior.
For the checking robustness of algorithm under other common image process is attacked herein, make an uproar to compressing respectively containing watermark original image, adding, filtering, shearing, rotation etc. attack experiment, in experiment, PSNR value be attack the rear calculated value containing watermarking images and original image.
(1) compression attack.Image containing watermark is carried out the compression that quality factor is 90,75,50,30,10, result Fig. 5 respectively.
Algorithm has very strong resistance to compression herein as can be seen from Figure 5, even if the related coefficient detected when quality factor is 10 high compression is all close to 0.9, watermarking images is also clear and legible.
(2) noise, Filtering Attacks.
Algorithm has comparatively strong immunity to gaussian filtering, 3 rank Butterworth low pass (D0=120) attacks herein as can be seen from Figure 6, comparatively responsive to medium filtering.
(3) shearing, rotation, convergent-divergent are attacked.By carry out containing the original image of watermark rotation 45 °, 90 °, amplify 2 times, reduce the attacks such as half, after this kind of attack, image needs to carry out process and returns to and original image same size, and attack algorithm adopts " bilinear interpolation (bilinear) ".Shearing attack is when detecting watermark, and need to do same shearing at original image and detect, experimental result as shown in Figure 7.As long as the one-piece construction of image is not destroyed as can be seen from Figure 7, its watermarking images can detect, and as half-twist, amplification 2 times, reduces half.If picture structure wrecks, the watermark of detection is just more poor, as rotated 45 °, the attack such as shearing.
Above-described embodiment is only for example of the present invention is clearly described, and is not the restriction to embodiments of the present invention.For those of ordinary skill in the field, can also make other changes in different forms on the basis of the above description.Here exhaustive without the need to also giving all embodiments.And these belong to spirit institute's apparent change of extending out of the present invention or change and are still among protection scope of the present invention.
Claims (1)
1. the embedding of watermarking images and an extracting method, is characterized in that:
The embedding grammar of described watermarking images comprises the steps:
(1) pre-service of watermarking images: watermarking images is carried out chaos sequence scramble, and after scramble, image is designated as W, and the image sequence this watermarking images being deformed into 1 row n row, i.e. n image information;
(2) selection of embedded location: by matrix be by initial carrier image
be divided into 8 × 8 image blocks do not covered mutually of corresponding number, dct transform is carried out to each 8 × 8 image blocks, low frequency corresponding in each 8 × 8 image blocks and high-frequency data is taken out according to 8 × 8 templates, and described low frequency and high-frequency data are formed the matrix A of 3 × 3 and the matrix B of 6 × 6 respectively, matrix A, B are carried out to svd and ask its matrix norm, solve the energy Ratios g of each 8 × 8 image blocks, energy Ratios
, wherein:
for the first singular value of described matrix B,
for the first singular value of described matrix A; Find out and 8 × 8 image blocks of energy Ratios corresponding to a front n maximal value;
(3) data selection of watermark embedment: watermark information and a described n image information, select correspondingly successively to embed adaptively in the first value of the svd of the described matrix A at the n arranged from big to small by described energy Ratios g 8 × 8 image blocks, embedment strength regulates automatically with the energy Ratios of image block;
(4) DCT inverse transformation is carried out to each 8 × 8 image blocks of embed watermark information, obtain image with watermarked information;
When carrying out the data selection of watermark embedment of above-mentioned steps (3): to the method that the matrix A of n 8 × 8 images carries out svd be one by one:
, in formula
,
,
, U and V is orthogonal matrix, and S is diagonal matrix;
Take out the first singular value of S diagonal matrix and embed watermark, embedding formula
; Embedment strength K is with block energy adjustment more automatic than g, and embed adaptively described in realizing, its formula is
,
, in formula
for element factor, its value condition need ensure that embedding Y-PSNR with watermarked information is greater than 36db; M is the total block data that original image is divided into described 8 × 8 image blocks;
New singular value is substituted into
, obtain low-frequency data with watermarked information
, and data are sent back to the correspondence position in corresponding 8 × 8 image blocks;
The extracting method of described watermarking images, comprises the steps:
A () to matrix is
initial carrier image and matrix be
be to be checkedly divided into 8 × 8 image blocks respectively containing watermarking images, and dct transform is carried out to each 8 × 8 image blocks, from each 8 × 8 image blocks, take out corresponding low frequency and high-frequency data according to 8 × 8 templates, and the low frequency in each 8 × 8 image blocks in initial carrier image and high-frequency data are formed respectively to the matrix of 3 × 3
with 6 × 6 matrix
, the low frequency in each 8 × 8 image blocks in watermarking images and high-frequency data are formed respectively to the matrix of matrix 3 × 3
with 6 × 6 matrix
, to described matrix
,
carry out svd respectively:
,
;
B () solves the watermark information changing value in each 8 × 8 image blocks
,
, m is the total block data that original image is divided into described 8 × 8 image blocks;
C (), from m 8 × 8 image blocks, finding out with energy Ratios is corresponding 8 × 8 image blocks of a front n maximal value, namely finds out individual 8 × 8 image blocks of n;
D the energy Ratios of one of the energy Ratios of one of individual 8 × 8 image blocks of the n described in () and 8 × 8 image blocks of described m is designated as respectively
with
,
,
; According to:
If
, then image intensity value
, draw n image information, K is embedment strength;
E (), by the energy Ratios size sequential of a described n image information according to 8 × 8 image blocks at each image information place, draws 1 row n row image sequence;
F above-mentioned 1 row n row image sequence is recovered watermarking images according to the inverse operation of Chaotic Scrambling algorithm by ();
In described step (a): take out the data A(9 of 1 × 9 from the low frequency and high frequency region coefficient of 8 × 8 template-setups), the data B(36 of 1 × 36) array, this array is deformed into the square formation of 3 × 3
,
, then to square formation
,
carry out SVD decomposition operation, take out the first value in singular value diagonal matrix.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210281560.XA CN102880997B (en) | 2011-05-26 | 2011-05-26 | Method for embedding watermark image |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210281560.XA CN102880997B (en) | 2011-05-26 | 2011-05-26 | Method for embedding watermark image |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110139206 Division CN102184519B (en) | 2011-05-26 | 2011-05-26 | Method for embedding and extracting watermark images |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102880997A CN102880997A (en) | 2013-01-16 |
CN102880997B true CN102880997B (en) | 2015-02-11 |
Family
ID=47482310
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210281560.XA Active CN102880997B (en) | 2011-05-26 | 2011-05-26 | Method for embedding watermark image |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102880997B (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103295180B (en) * | 2013-05-02 | 2015-12-23 | 陕西师范大学 | The digital watermark method of robustness and transparent optimum balance |
CN105303512B (en) * | 2015-09-18 | 2018-06-29 | 安徽大学 | Video double-watermarking method based on SVD and QR codes |
CN109285106B (en) * | 2018-10-11 | 2021-12-14 | 中南大学 | SLT-SVD-based reversible reduction method for medical images |
CN110163787B (en) * | 2019-04-26 | 2023-02-28 | 江苏水印科技有限公司 | Audio digital robust blind watermark embedding method based on dual-tree complex wavelet transform |
CN111489278B (en) * | 2020-04-07 | 2023-05-16 | 贵州大学 | Text watermark embedding and extracting method based on scrambling diffusion |
CN113190806A (en) * | 2021-04-25 | 2021-07-30 | 北京丁牛科技有限公司 | Data processing method |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1447233A (en) * | 2003-01-23 | 2003-10-08 | 同济大学 | Multi-media data protection method |
CN101493928A (en) * | 2009-02-10 | 2009-07-29 | 国网信息通信有限公司 | Digital watermarking embedding, extracting and quantizing step size coordinating factor optimizing method and device |
CN101847251A (en) * | 2010-06-13 | 2010-09-29 | 电子科技大学 | Digital watermark embedding and extraction method based on chaos sequences |
-
2011
- 2011-05-26 CN CN201210281560.XA patent/CN102880997B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1447233A (en) * | 2003-01-23 | 2003-10-08 | 同济大学 | Multi-media data protection method |
CN101493928A (en) * | 2009-02-10 | 2009-07-29 | 国网信息通信有限公司 | Digital watermarking embedding, extracting and quantizing step size coordinating factor optimizing method and device |
CN101847251A (en) * | 2010-06-13 | 2010-09-29 | 电子科技大学 | Digital watermark embedding and extraction method based on chaos sequences |
Non-Patent Citations (3)
Title |
---|
A geometric distortion resilient image watermarking algorithm based on SVD;ZHOU B,CHEN J.A;《Image and Graphics》;20041231;全文 * |
一种基于离散余弦变换与奇异值分解的数字图像水印算法;刘俊景,蒋华;《微电子学与计算机》;20071231;第24卷(第10期);全文 * |
基于分块奇异值分解的小波域水印算法;龚劬,马素春;《计算机仿真》;20090531;第26卷(第5期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN102880997A (en) | 2013-01-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102880998B (en) | The extracting method of watermarking images | |
CN102880997B (en) | Method for embedding watermark image | |
Keshavarzian et al. | ROI based robust and secure image watermarking using DWT and Arnold map | |
Kumar et al. | Image transformation technique using steganography methods using LWT technique | |
CN106023056A (en) | Zero-watermark embedding and extracting methods and zero-watermark embedding and extracting devices based on DWT and principal component analysis (PCA) compression | |
CN104217388B (en) | A kind of image watermark insertion based on FSSVM, the method and apparatus extracted | |
CN104616244A (en) | Image watermark embedding and extracting method based on back propagation (BP) neural network compressed domain | |
CN104053074B (en) | Video watermarking method based on depth image and Otsu segmentation | |
CN111988492B (en) | Gabor-DCT-based medical image robust watermarking method | |
CN102024249A (en) | Digital image watermarking method based on visual perception characteristics | |
CN107067363B (en) | A kind of non-blind robust Spatial Domain Image Watermark method constant based on pixel order relation | |
CN101719265B (en) | Digital image watermark embedding method and system thereof based on computer-generated hologram | |
CN102930498B (en) | A kind of digital watermark embedding based on Lifting Wavelet and extracting method | |
CN111988490B (en) | Tetrolet-DCT-based medical image robust watermarking method | |
CN1963865A (en) | A safety multifunctional image digital watermark system | |
CN102129657A (en) | Method for embedding multiple watermarks in volume data based on three-dimensional DFT (Delayed-First-Transmission) | |
CN101364300A (en) | Digital watermarking method based on gray theory | |
CN109754358B (en) | Image watermarking method and system based on significance detection and contourlet transformation | |
CN104217390B (en) | A kind of Zero watermarking method, device and watermark extracting method, device | |
CN103049880B (en) | Digital watermarking method based on discrete cosine transformation | |
CN111988491A (en) | Medical image robust watermarking method based on KAZE-DCT | |
CN116012212A (en) | Medical image robust watermarking method based on SUSAN-DCT and Hu moments | |
CN102510492A (en) | Method for embedding multiple watermarks in video based on three-dimensional DWT (Discrete Wavelet Transform) and DFT (Discrete Fourier Transform) | |
Kekre et al. | Robust Watermarking Technique using Hybrid Wavelet Transform Generated from Kekre Transform and Discrete Cosine Transform | |
CN104036450A (en) | Geometric-attack-resistant digital watermark processing method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C53 | Correction of patent of invention or patent application | ||
CB02 | Change of applicant information |
Address after: 213015 Changzhou Province in the Clock Tower District, Jiangsu, Wu Road, No. 1801 Applicant after: Jiangsu University of Technology Address before: 213000 Changzhou Province in the Clock Tower District, Jiangsu, Wu Road, No. 1801 Applicant before: Jiangsu Teachers University of Technology |
|
COR | Change of bibliographic data |
Free format text: CORRECT: APPLICANT; FROM: JIANGSU TECHNOLOGY NORMAL COLLEGE TO: JIANGSU UNIVERSITY OF TECHNOLOGY |
|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |