CN114897657A - 一种数字图像水印嵌入提取方法 - Google Patents

一种数字图像水印嵌入提取方法 Download PDF

Info

Publication number
CN114897657A
CN114897657A CN202210432306.9A CN202210432306A CN114897657A CN 114897657 A CN114897657 A CN 114897657A CN 202210432306 A CN202210432306 A CN 202210432306A CN 114897657 A CN114897657 A CN 114897657A
Authority
CN
China
Prior art keywords
watermark
image
embedding
intermediate frequency
frequency sub
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202210432306.9A
Other languages
English (en)
Other versions
CN114897657B (zh
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.)
Kunming University of Science and Technology
Original Assignee
Kunming University of Science and Technology
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 Kunming University of Science and Technology filed Critical Kunming University of Science and Technology
Priority to CN202210432306.9A priority Critical patent/CN114897657B/zh
Publication of CN114897657A publication Critical patent/CN114897657A/zh
Application granted granted Critical
Publication of CN114897657B publication Critical patent/CN114897657B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T1/00General purpose image data processing
    • G06T1/0021Image watermarking
    • G06T1/005Robust watermarking, e.g. average attack or collusion attack resistant
    • G06T1/0057Compression invariant watermarking
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T1/00General purpose image data processing
    • G06T1/0021Image watermarking
    • G06T1/005Robust watermarking, e.g. average attack or collusion attack resistant

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

本发明涉及一种数字图像水印嵌入提取方法,属于数字图像水印技术领域。本发明首先对原始载体图像进行离散平稳小波变换,提取中频子带LHI、HLI,将提取的中频子带LHI、HLI进行SVD分解。对水印图像进行离散小波变换(DWT),提取中频子带LHW、HLW,将提取的中频子带LHW、HLW进行SVD分解。利用嵌入规则,将水印图像嵌入到原始载体图像中。提取水印图像时,对嵌入水印后的载体图像进行离散平稳小波变换,提取中频子带LHM、HLM,对提取的中频子带LHM、HLM进行SVD分解,利用嵌入规则的反向规则分别提取出水印的两个奇异值矩阵,然后对奇异值矩阵进行反SVD分解得到水印图像两个中频子带,对中频子带进行离散小波反变换得到提取的水印图像。

Description

一种数字图像水印嵌入提取方法
技术领域
本发明涉及一种数字图像水印嵌入提取方法,属于数字图像水印技术领域。
背景技术
随着互联网技术发展的日新月异,数字媒体已经逐步替代传统纸质媒介,逐渐成为人们日常获取信息的重要渠道之一。数字媒体内容的获取方式简单高效,但数字媒体内容大肆传播也带来诸多问题,例如非法复制,伪造,盗版等数字版权保护问题也随之而来,严重侵犯了数字内容出版商的合法权益。为保证数字媒体内容合法传播利用,数字版权保护变得尤其重要。
数字水印技术可以有效地保护数字版权信息。数字水印技术通过将可用于识别版权的信息无形地嵌入到音频、图像、视频中以达到保护数字版权的目的,目前已经得到广泛应用。
发明内容
本发明要解决的技术问题是提供一种数字图像水印嵌入提取方法,解决上述问题。
本发明的技术方案是:一种数字图像水印嵌入提取方法,首先对原始载体图像进行离散平稳小波变换(SWT),提取中频子带LHI、HLI,将提取的中频子带LHI、HLI进行SVD分解。对水印图像进行离散小波变换(DWT),提取中频子带LHW、HLW,将提取的中频子带LHW、HLW进行SVD分解。利用嵌入规则,将水印图像嵌入到原始载体图像中。提取水印图像时,对嵌入水印后的载体图像进行离散平稳小波变换(SWT),提取中频子带LHM、HLM,对提取的中频子带LHM、HLM进行SVD分解,利用嵌入规则的反向规则分别提取出水印的两个奇异值矩阵,然后对奇异值矩阵进行反SVD分解得到水印图像两个中频子带,对中频子带进行离散小波反变换(IDWT)得到提取的水印图像。
具体步骤为:
Step1:水印嵌入过程:
Step1.1:对原始载体图像I进行离散平稳小波变换(SWT),得到LLI、LHI、HLI、HHI四个子带。载体图像I大小为M×M。子带大小为M×M。
Step1.2:提取中频子带LHI,进行SVD分解,得到奇异值矩阵
Figure BDA0003611490030000011
Figure BDA0003611490030000012
Step1.3:提取中频子带HLI,进行SVD分解,得到奇异值矩阵
Figure BDA0003611490030000013
Figure BDA0003611490030000021
Step1.4:水印图像W大小为N×N。原始载体图像与水印图像中M与N的关系:
M=N/2
Step1.5:对水印图像W进行离散小波变换(DWT),得到LLW、LHW、HLW、HHW四个子带。水印图像I大小为N×N。子带大小为
Figure BDA0003611490030000022
Step1.6:提取中频子带LHW,进行SVD分解,得到奇异值矩阵
Figure BDA0003611490030000023
Figure BDA0003611490030000024
Step1.7:提取中频子带HLW,进行SVD分解,得到奇异值矩阵
Figure BDA0003611490030000025
Figure BDA0003611490030000026
Step1.8:对
Figure BDA0003611490030000027
利用下式进行嵌入。
Figure BDA0003611490030000028
Step1.9:对
Figure BDA0003611490030000029
利用下式进行嵌入。
Figure BDA00036114900300000210
Step1.10:对矩阵
Figure BDA00036114900300000211
进行反SVD分解,得到矩阵LHM
Figure BDA00036114900300000212
Step1.11:对矩阵
Figure BDA00036114900300000213
进行反SVD分解,得到矩阵HLM
Figure BDA00036114900300000214
Step1.12:对嵌入水印信息后的中频子带LHM、HLM进行离散平稳小波反变换(ISWT),得到嵌入水印图像后的载体图像I*。
Step2:水印提取过程:
Step2.1:对嵌入水印后的载体图像I*进行离散平稳小波变换(SWT),得到LLI*、LHI*、HLI*、HHI*四个子带。
Step2.2:对中频子带LHI*进行SVD分解,得到奇异值矩阵
Figure BDA00036114900300000215
Figure BDA00036114900300000216
Step2.3:对中频子带HLI*进行SVD分解,得到奇异值矩阵
Figure BDA00036114900300000217
Figure BDA0003611490030000031
Step2.4:利用下式,得到水印图像的一个奇异值矩阵
Figure BDA0003611490030000032
Figure BDA0003611490030000033
Step2.5:利用下式,得到水印图像的另一个奇异值矩阵
Figure BDA0003611490030000034
Figure BDA0003611490030000035
Step2.6:对奇异值矩阵
Figure BDA0003611490030000036
进行反SVD分解,得到矩阵LHW*
Figure BDA0003611490030000037
Step2.7:对奇异值矩阵
Figure BDA0003611490030000038
进行反SVD分解,得到矩阵HLW*
Figure BDA0003611490030000039
Step2.8:对矩阵LHW*、HLW*进行离散平稳小波反变换(ISWT),得到提取的水印图像W*。
本发明的有益效果是:本发明对嵌入水印后的图像进行不同类型的攻击,并对水印图像进行提取,仿真结果表明嵌入水印后的载体图像在遭受到各类攻击后具有较好的不可见性的同时也表现出较好的鲁棒性。
附图说明
图1是本发明的水印嵌入过程图;
图2是本发明水印提取过程图;
图3是本发明三无攻击嵌入水印后的载体图像和提取的水印图像;
图4是本发明四椒盐噪声攻击嵌入水印后的载体图像和提取的水印图像;
图5是本发明JPEG压缩攻击嵌入水印后的载体图像和提取的水印图像;
图6是本发明中值滤波攻击嵌入水印后的载体图像和提取的水印图像。
具体实施方式
下面结合附图和具体实施方式,对本发明作进一步说明。
实施例1:如图1-6所示,一种数字图像水印嵌入提取方法,具体步骤为:
Step1:对原始载体图像进行SWT,提取中频子带LHI、HLI分别进行SVD分解,对水印图像进行DWT,提取中频子带LHW、HLW分别进行SVD分解,利用嵌入规则,将水印图像嵌入到原始载体图像中。
Step2:对嵌入水印后的载体图像进行SWT,提取中频子带LHM、HLM分别进行SVD分解,利用嵌入规则的反向规则分别提取出水印的奇异值矩阵,然后对奇异值矩阵进行反SVD分解得到水印图像两个中频子带,对中频子带进行IDWT得到提取的水印图像。
所述Step1中,利用
Figure BDA0003611490030000041
规则进行嵌入,其中
Figure BDA0003611490030000042
Figure BDA0003611490030000043
为原始载体图像奇异值矩阵,
Figure BDA0003611490030000044
为水印图像奇异值矩阵,
Figure BDA0003611490030000045
为嵌入水印后载体图像的奇异值矩阵,a为嵌入强度。
所述Step2中,利用
Figure BDA0003611490030000046
规则对水印图像进行提取,其中
Figure BDA0003611490030000047
为嵌入水印后载体图像的奇异值矩阵,
Figure BDA0003611490030000048
为提取水印图像的奇异值矩阵,a为嵌入强度。
所述Step1具体为:
对载体图像I进行离散平稳小波变换(SWT),得到LLI、LHI、HLI、HHI四个子带。原始载体图像大小256×256。子带大小256×256。
提取中频子带LHI,进行SVD分解,得到奇异值矩阵
Figure BDA0003611490030000049
Figure BDA00036114900300000410
提取中频子带HLI,进行SVD分解,得到奇异值矩阵
Figure BDA00036114900300000411
Figure BDA00036114900300000412
对水印图像W进行离散小波变换(DWT),得到LLW、LHW、HLW、HHW四个子带。水印图像W大小为512×512。子带大小为256×256。
提取中频子带LHW,进行SVD分解,得到奇异值矩阵
Figure BDA00036114900300000413
Figure BDA00036114900300000414
提取中频子带HLW,进行SVD分解,得到奇异值矩阵
Figure BDA00036114900300000415
Figure BDA00036114900300000416
Figure BDA00036114900300000417
利用下式进行嵌入,a为嵌入强度,取值0.04。
Figure BDA00036114900300000418
Figure BDA0003611490030000051
利用下式进行嵌入,a为嵌入强度,取值0.04。
Figure BDA0003611490030000052
对矩阵
Figure BDA0003611490030000053
进行反SVD分解,得到矩阵LHM
Figure BDA0003611490030000054
对矩阵
Figure BDA0003611490030000055
进行反SVD分解,得到矩阵HLM
Figure BDA0003611490030000056
对嵌入水印信息后的中频子带LHM、HLM进行离散平稳小波反变换(ISWT),得到嵌入水印图像后的载体图像I*。
所述Step1具体为:
对嵌入水印后的载体图像I*进行离散平稳小波变换(SWT),得到LLI*、LHI*、HLI*、HHI*四个子带。嵌入水印后的载体图像I*大小为256×256。
对中频子带
Figure BDA00036114900300000519
进行SVD分解,得到
Figure BDA0003611490030000057
矩阵。
Figure BDA0003611490030000058
对中频子带
Figure BDA00036114900300000520
进行SVD分解,得到
Figure BDA0003611490030000059
矩阵。
Figure BDA00036114900300000510
利用下式,得到水印图像的一个奇异值矩阵
Figure BDA00036114900300000511
a为嵌入强度,取值0.04。
Figure BDA00036114900300000512
利用下式,得到水印图像的另一个奇异值矩阵
Figure BDA00036114900300000513
a为嵌入强度,取值0.04。
Figure BDA00036114900300000514
对奇异值矩阵
Figure BDA00036114900300000515
进行反SVD分解,得到矩阵
Figure BDA00036114900300000521
Figure BDA00036114900300000516
对奇异值矩阵
Figure BDA00036114900300000517
进行反SVD分解,得到矩阵
Figure BDA00036114900300000522
Figure BDA00036114900300000518
对矩阵LHW*、HLW*进行离散平稳小波反变换(ISWT),得到提取的水印图像W*。
为了验证本发明的有效性,本发明使用MATLAB作为实验仿真平台。原始载体图像使用lena灰度图像,图像大小256*256。水印图像为512*512的QR码图像,离散平稳小波(SWT)方案采用Symlets小波基函数,离散小波(DWT)方案采用haar小波基函数。嵌入强度a为0.04。
如图3所示,Host image为原始载体图像lena灰度图,Watermark logo为水印图像,Watermarkedimage为无攻击时嵌入水印后的载体图像,Extractedwatermark为无攻击时提取的水印图像。
本发明使用用峰值信噪比(PSNR)来体现不可感知性。峰值信噪比(PSNR)越高,代表嵌入水印后的载体图像的视觉质量越好。结合相关参考文献,嵌入水印图像后的载体图像峰值信噪比(PSNR)达到30dB或者更高的要求时,人类视觉对嵌入水印后的载体图像不会察觉到明显变化。
本发明使用用相关系数(NC)来体现鲁棒性。相关系数(NC)所表示的含义是提取的水印图像和原水印之间的相似性。
表1是不同攻击对比实验数据。嵌入水印后的载体图像在遭受不同类型的攻击时峰值信噪比(PSNR)以及遭受不同类型的攻击后提取水印与原始水印相关系数(NC)。
Figure BDA0003611490030000061
表1
经过实验,发现本发明提出的算法对等常见的攻击都表现出较好的鲁棒性,同时对于无攻击时也表现出较好的峰值信噪比,具有较好的不可见性。
以上结合附图对本发明的具体实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出各种变化。

Claims (3)

1.一种数字图像水印嵌入提取方法,其特征在于:
Step1:对原始载体图像进行SWT,提取中频子带LHI、HLI分别进行SVD分解,对水印图像进行DWT,提取中频子带LHW、HLW分别进行SVD分解,利用嵌入规则,将水印图像嵌入到原始载体图像中;
Step2:对嵌入水印后的载体图像进行SWT,提取中频子带LHM、HLM分别进行SVD分解,利用嵌入规则的反向规则分别提取出水印的奇异值矩阵,然后对奇异值矩阵进行反SVD分解得到水印图像两个中频子带,对中频子带进行IDWT得到提取的水印图像。
2.根据权利要求1所述的数字图像水印嵌入提取方法,其特征在于:所述Step1中,利用
Figure FDA0003611490020000011
规则进行嵌入,其中
Figure FDA0003611490020000012
为原始载体图像奇异值矩阵,
Figure FDA0003611490020000013
为水印图像奇异值矩阵,
Figure FDA0003611490020000014
为嵌入水印后载体图像的奇异值矩阵,a为嵌入强度。
3.根据权利要求1所述的数字图像水印嵌入提取方法,其特征在于:所述Step2中,利用
Figure FDA0003611490020000015
规则对水印图像进行提取,其中
Figure FDA0003611490020000016
为嵌入水印后载体图像的奇异值矩阵,
Figure FDA0003611490020000017
为提取水印图像的奇异值矩阵,a为嵌入强度。
CN202210432306.9A 2022-04-23 2022-04-23 一种数字图像水印嵌入提取方法 Active CN114897657B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210432306.9A CN114897657B (zh) 2022-04-23 2022-04-23 一种数字图像水印嵌入提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210432306.9A CN114897657B (zh) 2022-04-23 2022-04-23 一种数字图像水印嵌入提取方法

Publications (2)

Publication Number Publication Date
CN114897657A true CN114897657A (zh) 2022-08-12
CN114897657B CN114897657B (zh) 2024-10-01

Family

ID=82717183

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210432306.9A Active CN114897657B (zh) 2022-04-23 2022-04-23 一种数字图像水印嵌入提取方法

Country Status (1)

Country Link
CN (1) CN114897657B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116542840A (zh) * 2023-07-06 2023-08-04 云南日报报业集团 一种基于新闻文本图像的鲁棒水印方法
CN116563081A (zh) * 2023-07-06 2023-08-08 云南日报报业集团 一种空间域盲图像水印方法
CN116596737A (zh) * 2023-07-06 2023-08-15 云南日报报业集团 一种图像水印嵌入与提取方法
CN116645259A (zh) * 2023-07-18 2023-08-25 云南日报报业集团 一种基于rdwt-hd-svd的水印嵌入提取方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103955880A (zh) * 2014-04-11 2014-07-30 杭州电子科技大学 基于Zernike矩的DWT-SVD鲁棒盲水印方法
CN104680470A (zh) * 2013-11-26 2015-06-03 上海宝钢包装股份有限公司 一种基于dwt-svd的强鲁棒性全息水印嵌入和提取方法
CN110264390A (zh) * 2019-06-24 2019-09-20 上海海事大学 一种基于双水印嵌入的数字水印方法
CN112132731A (zh) * 2020-09-10 2020-12-25 郑州轻工业大学 采用预设psnr的dwt-svd域自适应鲁棒水印算法
CN112801846A (zh) * 2021-02-09 2021-05-14 腾讯科技(深圳)有限公司 一种水印嵌入、提取方法、装置、计算机设备和存储介质

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104680470A (zh) * 2013-11-26 2015-06-03 上海宝钢包装股份有限公司 一种基于dwt-svd的强鲁棒性全息水印嵌入和提取方法
CN103955880A (zh) * 2014-04-11 2014-07-30 杭州电子科技大学 基于Zernike矩的DWT-SVD鲁棒盲水印方法
CN110264390A (zh) * 2019-06-24 2019-09-20 上海海事大学 一种基于双水印嵌入的数字水印方法
CN112132731A (zh) * 2020-09-10 2020-12-25 郑州轻工业大学 采用预设psnr的dwt-svd域自适应鲁棒水印算法
CN112801846A (zh) * 2021-02-09 2021-05-14 腾讯科技(深圳)有限公司 一种水印嵌入、提取方法、装置、计算机设备和存储介质

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
ADNAN MUSTAFA CHEEMA: ""A Novel Optimized Semi-Blind Scheme for Color Image Watermarking"", 《IEEE ACCESS》, vol. 8, 15 September 2020 (2020-09-15), pages 169525 - 169547, XP011810550, DOI: 10.1109/ACCESS.2020.3024181 *
王莹: ""Contourlet域的盲数字水印算法研究"", 《中国优秀硕士学位论文全文数据库 信息科技辑》, no. 2020, 15 March 2020 (2020-03-15), pages 138 - 185 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116542840A (zh) * 2023-07-06 2023-08-04 云南日报报业集团 一种基于新闻文本图像的鲁棒水印方法
CN116563081A (zh) * 2023-07-06 2023-08-08 云南日报报业集团 一种空间域盲图像水印方法
CN116596737A (zh) * 2023-07-06 2023-08-15 云南日报报业集团 一种图像水印嵌入与提取方法
CN116542840B (zh) * 2023-07-06 2023-09-22 云南日报报业集团 一种基于新闻文本图像的鲁棒水印方法
CN116563081B (zh) * 2023-07-06 2023-09-26 云南日报报业集团 一种空间域盲图像水印方法
CN116645259A (zh) * 2023-07-18 2023-08-25 云南日报报业集团 一种基于rdwt-hd-svd的水印嵌入提取方法
CN116645259B (zh) * 2023-07-18 2023-09-26 云南日报报业集团 一种基于rdwt-hd-svd的水印嵌入提取方法

Also Published As

Publication number Publication date
CN114897657B (zh) 2024-10-01

Similar Documents

Publication Publication Date Title
CN114897657B (zh) 一种数字图像水印嵌入提取方法
Divecha et al. Implementation and performance analysis of DCT-DWT-SVD based watermarking algorithms for color images
Al-Haj Combined DWT-DCT digital image watermarking
CN105512999B (zh) 一种双变换的彩色图像全息水印方法
Malonia et al. Digital image watermarking using discrete wavelet transform and arithmetic progression technique
Patil et al. DWT based invisible watermarking technique for digital images
Salama et al. A robust digital image watermarking technique based on wavelet transform
Kumar et al. A Robust Algorithm for Digital Image copyright protection
CN116503230A (zh) 一种基于双通道的鲁棒图像水印算法
Maity et al. An image watermarking scheme using HVS characteristics and spread transform
Shahraeini et al. A robust digital image watermarking approach against jpeg compression attack based on hybrid fractal-wavelet
Singh et al. Spatial and frequency domain for grey level digital images
Lin et al. A robust watermark scheme for copyright protection
CN111340675B (zh) 一种基于稀疏表示的彩色图案水印嵌入和提取方法
CN116645259B (zh) 一种基于rdwt-hd-svd的水印嵌入提取方法
Al-Asmari et al. A DWT ordering scheme for hiding data in images using pixel value difference
Al-Haj A hybrid digital image watermarking algorithm
Zhu et al. A novel blind watermarking scheme in Contourlet domain based on singular value decomposition
Prajapati Transform based digital image watermarking techniques for image authentication
Fazli et al. Advanced wavelet based steganography for colored images
Bairagee et al. Review of digital image watermarking techniques
Kumari et al. Robust image watermarking algorithm DWT-SVD-AES for telemedicine application
Channapragada et al. Digital watermarking algorithm based on complete complementary code
Xiong An improved DCT based color image watermarking scheme
CN116596737A (zh) 一种图像水印嵌入与提取方法

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
GR01 Patent grant
GR01 Patent grant