“程序设计与算法训练”课程设计“二值图像数字水印技术的实践”

数据结果课程设计（题目抄袭自某牛逼985高校）。

项目地址（含报告等）

课程编号: 课程性质:必修

程序设计与算法训练课程设计报告

院 系： 计算机与信息系

班 级：

姓 名：

学 号：

指导教师：

选题名称： 二值图像数字水印技术实践

2019 年 1 月 20 日 至 2019 年 5 月 5 日

目录

程序设计与算法训练课程设计报告

一、 实验概述

1.1 课程设计题目

1.2 课程设计目的

1.3 系统主要内容与功能

1.3.1 设计内容

1.3.2 设计功能

1.3.3 系统类图

1.3.4 属性和方法定义

1.3.5 实验环境与工具

二、 实验原理

2.1 图像水印技术简述

2.2 图像处理技术简述

2.3 Qt 及 Qt

2.4 相关算法简介

2.5 技术流程

三、 实验结果

四、 算法分析

五、 总结

附录(源代码与其他)

源代码

其他

实验代码

实验结果

一、 实验概述 1.1 课程设计题目

题目要求对给定的一种简单的二值图像的数字水印算法编程实现。

1.2 课程设计目的

对数字水印技术建立一定的认识,能建立位矩阵、位向量等 ADT,并能用这些 ADT 完成给定二值图像数字水印的嵌入和抽取。

1.3 系统主要内容与功能 1.3.1 设计内容

具体包括以下内容:

(1)图像的读取与保存,及相应的矩阵和向量的运算;

(2)二值图像水印算法的实现;

(3)软件的界面和接口设计,信号发送和槽位的设计;

(4)软件鲁棒性分析、算法鲁棒性分析和相关总结。

1.3.2 设计功能

(1)设计合理的数据结构,编程实现算法;

(2)给定测试图片,按照指定的 bmp 格式,保存于外存中。

1.3.3 系统类图

( &array);

( &array);

( &);

(uchar* , uchar* dst, const int width, const int , const int );

uchar* (uchar* , const int width, const int );

( src, key);

(const int );

( &dir);

uchar* (const char *, int& , int& );

( &str);

uchar* (uchar* , uchar* , const int size);

uchar* (uchar* , const int width, const int , int x_off, int y_off);

uchar* (uchar* , uchar* edge, const int size, code);

bool (const char* , uchar* , const , const width);

BMIH;

BMFH;

int ;

int ;

int ;

int ;

* ;

( * = 0);

~();

void ();

void ged(const &arg1);

void ();

void ();

Ui:: *ui;

image;

key;

uchar* dst;

表 1.1 系统类图

1.3.4 属性和方法定义

类名

成员类型

类型

成员名

描述

属性

BMIH

BMP图像信息头

BMFH

BMP图像文件头

int

BMP图像宽

int

BMP图像高

int

图像类型，每像素位数

int

图像数据每行字节数

*

BMP图像调色板信息

方法

秘钥序列转换为二进制信息

秘钥序列转换为字符串

二进制信息转秘钥序列

使用水印秘钥编码原图

uchar*

获得从提取到的边缘图像，用作添加水印位置的参考

使用秘钥序列加密二进制

生成秘钥序列

uchar*

读取BMP图片

字符串转换为秘钥序列

uchar*

获得和相减得到的结果

uchar*

获得原图像右移x_off个单位，下移y_off个单位后得到的图像

uchar*

由边缘图像和原图获得水印编码后的图像

bool

保存编码水印信息的图像

属性

Ui::*

ui

维护图形界面

image

维护读入的BMP图像

key

秘钥序列

uchar*

dst'

维护编码水印的目标图像

方法

无

( * = 0)

构造函数，初始化界面

无

~()

析构函数，释放内存（Ui）

槽位

void

()

浏览、选择图像

void

ged(const &arg1)

修改水印信息

void

()

水印信息编码BMP图像

void

()

解码水印BMP图像，获取水印信息

表 1.2 属性和方法定义

1.3.5 实验环境与工具

(1)操作系统: 16.04 LTS;

(2)开发工具:Qt、Qt ;

(3)实现语言:C++。

二、 实验原理 2.1 图像水印技术简述

随着互联网和信息技术的快速发展,近年来数字内容的未授权获取,传输,操纵和分发的问题变得越来越严重。信息安全研究引

起了人们的广泛关注。除了一般采用的数字加密算法之外,近年来用于信息安全的影像视觉算法包括光学图像加密、认证和水印算

法被广泛地研究和应用起来。影像视觉信息安全算法通常拥有并行高速处理和多维能力的优势。信息隐藏技术,即图像水印技术,也

称隐写技术是一种隐蔽性地改变载波信号以嵌入隐藏消息,即水印信号的技术。可以对各种各样类别的信号执行信息隐藏,其中包括

但不限于:音频信号、图像信号和视频信号。信息隐藏技术(图像水印技术)允许将特定的信息加入到需要保护的媒体信息中,加入的

信息一般为具有特定意义的内容,如版权所有者信息、发行标志、特定代码等。而图像水印技术也因而成为了数字图像处理专业当

下或未来重要的研究领域,在知识产权的保护等方面有着广泛的应用前景。为了确保大规模在线分发的多媒体内容的版权和知识产

权,我们需要通过有效的保护来控制分发和传播,控制来自盗版用户或未经授权普通用户的恶意操纵和恶意拷贝传播。

为了提高效率,水印需要良好的隐蔽性,并且拥有嵌入高容量和有效载荷,能够在确保有效载荷的安全传输的同时,对最常见的图像处理(恶意或非恶意)进行鲁棒性处理。此外水印技术还有以下的特点:

(1)水印后的主图像不应存在显著地信息损坏和分辨率降低;

(2)水印应具有在主图像中良好的隐蔽性,即不可见性,一般用图像的峰值信噪比 (PSNR) 与结构相似形（SSIM）来描述;

(3)水印应具有良好的鲁棒性,不易从主图像中被损坏。水印应可以承受不同类型的信息损坏打击,例如 JPEG 压缩、裁剪、旋

转、缩放、噪声、滤波运算和模糊运算。应该特别指出,该特点仅适用于一部分鲁棒性极好的水印算法中;

(4)未经授权的用户/盗版用户应很难非法访问到该水印信息。

2.2 图像处理技术简述

在图像水印技术实现时,需要大量运用到图像处理知识和技术。数字图像处理( Image )是通过计算机对图像

进行去除噪声、增强、复原、分割、提取特征等处理的方法和技术。数字图像处理( Image )又称为计算机图像

处理,它是指将图像信号转换成数字信号并利用计算机对其进行处理的过程。数字图像处理的产生和迅速发展主要受三个因素的影

响:一是计算机的发展;二是数学的发展(特别是离散数学理论的创立和完善);三是广泛的农牧业、林业、环境、军事、工业和医学等

方面的应用需求的增长。

在图像水印技术中,本着不损坏原有图像质量的原则,该算法需要对图像快速的读写能力、对分块图像稳定和鲁棒的运算能力,

以及优秀的边缘提取能力。

2.3 Qt 及 Qt

Qt 是一个跨平台的 C++图形用户界面应用程序框架。它为应用程序开发者提供建立艺术级图形用户界面所需的所有功能。

它是完全面向对象的,很容易扩展,并且允许真正的组件编程。

作为一个优秀的 C++框架,由于其优秀的信号与槽机制,Qt 被广泛地应用于软件编写中。本次课程设计将采用 Qt 和 Qt

作为程序界面的编写工具。

2.4 相关算法简介

数字图像水印算法是将一段信息附加在图像上的算法,其中被附加信息的图像被称为主机图像(host image),简称主图像。根据

水印算法的鲁棒性,水印算法可以分为鲁棒/稳健的水印算法( )和脆弱的水印算法(),其中绝

大多数的数字水印算法都属于前者。如上文所提到,鲁棒的水印在主图像受到攻击和扭曲时,仍能保持完整。由于鲁棒的水印很难从

主机图像中移除,因此常常用于保护版权;另一方面,脆弱的水印一般仅用于验证,即主机图像的完整性检查。完整的脆弱水印表示主

机图像处于其原始形式,并没有收到编辑、损坏或更改。在本次课程设计中,为了简洁起见,我们采用脆弱的水印算法。

经典的数字图像水印算法构建的系统包括双随机相位编码(DRPE)系统,离轴全息系统,相移全息系统,优化的仅相位掩模结构,

联合变换相关器(JTC),重影成像系统和 系统,等等。它们在各自的领域都具有相当不错的效果。在本次课程设计中,

为了兼顾效率和效果,我们采用一种简单的二值图像数字水印算法,其基本思路如下:由于水印算法需要改变主机图像的像素值,从而

一定程度上改变主机图像的信息。而被改变像素值,因而保存着水印信息的像素点被称为隐藏点,它决定了隐藏了怎样的信息。而简

单的处理方法容易导致图像质量的下降,例如,在全白的图像块中插入一个黑色的噪声点。另一方面,隐藏点也应避免选取在图像中

的细线区域、直线边的中间像素、孤立像素等图像信息熵较大的点。由此,二值图像的数字水印嵌入算法的关键是隐藏点的选取,以

下是隐藏点的选取规则:

二值图像数据应隐藏域图像中黑/白色区域的边界上,但上述诸如细线区域等仍不适于隐藏数据的边界。因此隐藏点应具有以

下的特点:它是边界像素,并且不同时是左边界和右边界/上边界和下边界,以确保避免将信息隐藏在细线区域等。这需要一个优秀的

边界提取算法。

在隐藏点被确定后,水印信息应转换为二进制序列,并保存在隐藏点中。同理,因此所有可被转换为二进制序列的诸如文字信

号、音频信号、图像信号和视频信号,并可进行反转换的信号都可以作为该算法的水印信息。

在本次课程设计中,我们简单地以字符串信息和二进制序列信息为例,来保存相关的水印。值得一提的是,为了保证水印信息更

好地隐蔽,水印信息被加入前,需要用一段密钥进行加密,并在提取时,使用该密钥进行解密。

2.5 技术流程

本次课程设计采用的主要框架包括 C++与 Qt 等,其中 C++作为程序实现语言,Qt 作为软件实现框架,使用 C++中的指针工具作

为图像处理工具,BMP 图像的读写采用文件头信息解析模式。

本次编程是在 16.04 LTS 上进行的,在该环境下,常用的图像解码工具包括 Qt 的 模块以及 库。其

中使用 C++进行 bmp 文件解码过程较为复杂,需要鲁棒的文件读写设计以及接口设计,同时,对于不同类型的文件需要不同的解码方

式,不广泛适用于.png,.img,.jpg,.giff 等等图像格式的使用,但减少了第三方库的利用。而 框架的引入对于文件数据对象的

空间开辟和释放过于频繁,降低了程序的效率,并需要引入标准模板库框架。作为对比, 库拥有极高的鲁棒性、延展性和高

效性,但需要引入第三方库,软件打包较为复杂。在“尽量引入较少的第三方库”的原则下,我们采用了纯C++语言,根据 BMP 图像的特

点设计合理的数据结构进行编码、存储图像,以提高本软件的实用性。

本次课程设计的技术流程如下:首先设计并实现边缘信息提取的算法()函数,以确定隐藏点,再设计数字隐藏的算法

()函数。在软件流程中,首先点击“浏览”按钮,选取二值图像,再输入需要加入的水印,最后点击“编码”按钮,则生成并显示

了附有水印信息的图像。该图像被保存在.pro 的 Qt 项目文件同名文件夹下,被命名为 .bmp。此外,对于已编码的

.bmp,点击“解码”按钮,会以对话框形式弹出被解码出的水印信息。

三、 实验结果

本次课程设计生成的软件结果如下:

软件打开后的界面如下,为避免用户每次打开软件后都需要输入水印,过于麻烦,本软件内置了水印二进制序列,如图3.1所示。

图 3.1 软件开启界面

输入二值图像后,会被展示在软件中,如图3.2所示。

图 3.2 显示所输入二值图像

点击“编码”按钮后,生成编码图像,如图3.3所示。

图 3.3 生成编码图像

点击“解码”按钮,弹出水印信息,如图3.4所示。

图 3.4 显示水印信息

本程序也支持字符串的水印编码,将水印附带的 选取为“”,然后编码并解码后的结果,如图3.5所示。

图 3.5 字符串水印编码图像

四、 算法分析

本程序经过第三方测试发现并无 bug 存在,鲁棒性优秀,可以完整的实现简单的二进制序列的二值数字图像水印算法功能。

在第二章节中提到,优秀的水印算法拥有上述的四个特点:a.不改变原图信息;b.隐蔽性;c.鲁棒性;d.不可访问性。本报告将从这四个

角度评判该算法的性能。

1. 由第三章节的内容易见,该算法具有从肉眼中几乎无法分辨与原图像的区别,同时如表 3 所示,通过计算原图与被水印编码图

的峰值信噪比与结构相似比,两者均处于合理的置信区间,因而完美地符合了不改变原图信息和隐蔽性的原则;

2. 该算法经过轻微地噪声扰动后就无法保存水印信息,如图 6 所示,因此鲁棒性极低,然而这符合脆弱的水印算法的特点,可以用

于验证图像的完整性;

3. 从 .bmp 和原图之间的差分影像可以轻易地获得加密后的二进制序列。但是由于有密钥的存在,无密钥的非授权用

户和盗版用户无法访问到原始水印的信息,具有不可访问性。关于不可访问性,在附录中有所展示。

对比

属性

原图与水印图

原图与加噪声水印图

brain

hippo

brain

hippo

峰值信噪比(PSNR)

39.6375

27.5873

44.7908

13.7145

7.6078

13.621

结构相似性(SSIM)

0.9999

0.9963

0.9998

0.4152

0.3125

0.3331

表 4.1峰值信噪比与结构相似性实验数据

图 4.1原图、水印图、加噪声水印图对比

综上所述,本算法在效率极高的同时,完美地满足了优秀的脆弱水印算法的四个特点,同时保证了效果和效率。

五、 总结

在本次课程设计中,按照题目要求与提示,针对给定的二值图像数字水印算法完成了编程实现、算法分析与相关内容的简单介

绍。先将水印信息转换为二进制序列,然后,算法利用密钥对序列进行转换,从而降低了第三方获取水印的可能性。通过查阅有关资

料,对数字水印建立了一定的认识,并构建了位矩阵、位相量等 ADT。综述,本次课程设计完成了既定的各个要求。

附录(源代码与其他) 源代码

头文件 .h:

#ifndef ADDDITIONAL_UTILITY_H
#define ADDDITIONAL_UTILITY_H#include
#include
#include
#include
#include
#include
#include
#pragma pack(2)// In the Linux environmentusing byteArray = QVector;//******************************************top******************************************************
//In the Linux environment
typedef struct BITMAPFILEHEADER
{u_int16_t bfType;u_int32_t bfSize;u_int16_t bfReserved1;u_int16_t bfReserved2;u_int32_t bfOffBits;
}BITMAPFILEHEADER;
typedef struct BITMAPINFOHEADER
{u_int32_t biSize;u_int32_t biWidth;u_int32_t biHeight;u_int16_t biPlanes;u_int16_t biBitCount;u_int32_t biCompression;u_int32_t biSizeImage;u_int32_t biXPelsPerMeter;u_int32_t biYPelsPerMeter;u_int32_t biClrUsed;u_int32_t biClrImportant;struct BITMAPINFOHEADER &operator=( struct BITMAPINFOHEADER &BMIH ){biSize = BMIH.biSize;biWidth = BMIH.biWidth;biHeight = BMIH.biHeight;biPlanes = BMIH.biPlanes;biClrUsed = BMIH.biClrUsed;biSizeImage = BMIH.biSizeImage;biCompression = BMIH.biCompression;biClrImportant = BMIH.biClrImportant;biXPelsPerMeter = BMIH.biXPelsPerMeter;biYPelsPerMeter = BMIH.biYPelsPerMeter;//qDebug() << "214235423" << endl;return *this;}}BITMAPINFODEADER;typedef unsigned char BYTE;
typedef unsigned short WORD;
typedef unsigned int DWORD;
typedef long LONG;typedef struct tagRGBQUAD
{BYTE rgbBlue;BYTE rgbGreen;BYTE rgbRed;BYTE rgbReserved;
}RGBQUAD;typedef struct tagIMAGEDATA
{BYTE blue;BYTE green;//cancel the annotation by DFZBYTE red;//cancel the annotation by DFZ
}IMAGEDATA;//********************************************bottom*******************************************#endif // ADDDITIONAL_UTILITY_H

头文件 .h:

#ifndef BMPUTIL_H
#define BMPUTIL_H#include "addditional_utility.h"class watermark
{
public:QString array2byte(byteArray &array);QString array2str(byteArray &array);byteArray byte2Array(QString &number);byteArray decodeImg(uchar* buffer, uchar* dst, const int width, const int height, const int length);uchar* edgeExtract(uchar* buffer, const int width, const int height);byteArray encode(byteArray src, byteArray key);byteArray generateKey(const int length);byteArray img2Array(QString &dir);uchar* readBmp(const char *bmpName, int& bmpWidth, int& bmpHeight);byteArray str2Array(QString &str);uchar* substract(uchar* buffer1, uchar* buffer2, const int size);uchar* translation(uchar* buffer, const int width, const int height, int x_off, int y_off);uchar* watermarkImg(uchar* buffer, uchar* edge, const int size, byteArray code);bool savebmp(const char* filename, uchar* buffer, const u_int32_t height, const u_int32_t width);
private:BITMAPINFODEADER BMIH;BITMAPFILEHEADER BMFH;int biWidth;int biHeight;int biBitCount;int lineByte;RGBQUAD* pColorTable;
protected:};#endif // BMPUTIL_H

头文件 .h:

#ifndef MAINWINDOW_H
#define MAINWINDOW_H#include 
#include namespace Ui
{
class MainWindow;
}using byteArray = QVector;class MainWindow : public QMainWindow
{Q_OBJECT
public:explicit MainWindow(QWidget *parent = 0);~MainWindow();private slots:void on_pushButtonBrowse_clicked();void on_lineEdit_textChanged(const QString &arg1);void on_pushButtonEncode_clicked();void on_pushButtonDecode_clicked();private:Ui::MainWindow *ui;QPixmap image;byteArray key;uchar* dst;
};#endif // MAINWINDOW_H

源文件 .cpp:

#include "addditional_utility.h"

源文件 .cpp:

#include "bmputil.h"uchar* watermark::readBmp(const char *bmpName, int& bmpWidth, int& bmpHeight)
{FILE *fp = fopen(bmpName, "rb");if(fp == Q_NULLPTR){QMessageBox::warning(Q_NULLPTR, "Error", "Error in Open File!");return Q_NULLPTR;}// skip the fileheaderfseek(fp, sizeof(BITMAPFILEHEADER), SEEK_CUR);// read the infoheader//BITMAPINFOHEADER* head = new BITMAPINFOHEADER;watermark WT;//infohead = &(WT.BMIH);BITMAPINFODEADER infohead = WT.BMIH;fread(&infohead,sizeof(BITMAPINFOHEADER),1,fp);//fread(infoHead, sizeof(BITMAPINFOHEADER), 1, fp);bmpWidth = infohead.biWidth;bmpHeight = infohead.biHeight;biBitCount = infohead.biBitCount;lineByte = (bmpWidth*biBitCount/8+3)/4*4;//qDebug() << infohead.biSize << endl;//qDebug() << infohead.biWidth << endl;//qDebug() << infohead.biHeight << endl;//qDebug() << infohead.biBitCount << endl;if (biBitCount == 8){pColorTable = new RGBQUAD[256];fread(pColorTable, sizeof(RGBQUAD), 256, fp);uchar* pBmpBuf = new uchar[ bmpWidth * bmpHeight ];fread(pBmpBuf, sizeof(uchar), bmpWidth * bmpHeight, fp);fclose(fp);uchar* buffer = new uchar[bmpWidth * bmpHeight];for(int i = 0; i < bmpHeight; i++){for(int j = 0; j 0.5){res.append(true);}else{res.append(false);}}return res;
}//获得原buffer图像右移x_off个单位，下移y_off个单位后得到的图像
uchar* watermark::translation(uchar* buffer, const int width, const int height, int x_off, int y_off)
{uchar* res = new uchar[width*height];for(int i = 0; i buffer2[i] ? buffer1[i] - buffer2[i] : 0;}return res;
}// 获得从buffer提取到的边缘图像，用作添加水印位置的参考
uchar* watermark::edgeExtract(uchar* buffer, const int width, const int height)
{uchar* BL = substract(translation(buffer, width, height, 1, 0), buffer, width*height);uchar* BR = substract(translation(buffer, width, height, -1, 0), buffer, width*height);uchar* BT = substract(translation(buffer, width, height, 0, 1), buffer, width*height);uchar* BB = substract(translation(buffer, width, height, 0, -1), buffer, width*height);uchar* B1 = new uchar[height*width];for(int i = 0; i= height){continue;}for(int b = -1; b<2; b++){if(b +j <0 || b+j>=width){continue;}sum1 += buffer[(a + i) *width + (b + j)]/255;sum2 += B1[(a + i) *width + (b + j)];}}if(sum1 == sum2){res[i*width + j] = 0;}}}}return res;
}// 由边缘图像和原图获得水印编码后的图像
uchar* watermark::watermarkImg(uchar* buffer, uchar* edge, const int size, byteArray code)
{uchar* res = new uchar[size];for(int i = 0; i

#include "mainwindow.h"
#include "ui_mainwindow.h"#include 
#include 
#include 
#include #include "bmputil.h"
//#include "bmputil.cpp"MainWindow::MainWindow(QWidget *parent) :QMainWindow(parent),ui(new Ui::MainWindow)
{ui->setupUi(this);image = QPixmap();QStringList bands = QStringList() << "QString" << "byteArray";ui->comboBoxWaterMark->setModel(new QStringListModel(bands));ui->comboBoxWaterMark->setCurrentIndex(1);ui->lineEditWaterMark->setText("01010101010101010101010101010101");
}MainWindow::~MainWindow()
{delete ui;
}void MainWindow::on_pushButtonBrowse_clicked()
{auto file = QFileDialog::getOpenFileName(this, tr("Append selected images"));ui->lineEdit->setText(file);
}void MainWindow::on_lineEdit_textChanged(const QString &arg1)
{Q_UNUSED(arg1);watermark WT;if(QFileInfo(ui->lineEdit->displayText()).exists()){int height, width;uchar* buffer = WT.readBmp(ui->lineEdit->displayText().toStdString().data(), width, height);if(buffer){QPixmap img = QPixmap::fromImage(QImage(buffer, width, height, QImage::Format_Grayscale8));QGraphicsScene *scene = new QGraphicsScene;scene->addPixmap(img);ui->graphicsViewPrevious->setScene(scene);ui->graphicsViewPrevious->show();ui->graphicsViewPrevious->fitInView(img.rect(), Qt::KeepAspectRatio);}}
}void MainWindow::on_pushButtonEncode_clicked()
{byteArray code;watermark WT;if(ui->comboBoxWaterMark->currentIndex() == 1){QRegExp regx("[0-1]+$");QValidator *validator = new QRegExpValidator(regx, this );ui->lineEditWaterMark->setValidator( validator );QString str1 = ui->lineEditWaterMark->text();code = WT.byte2Array(str1);delete validator;}else{QRegExp regx(".+\n");QValidator *validator = new QRegExpValidator(regx, this );ui->lineEditWaterMark->setValidator( validator );QString str2 = ui->lineEditWaterMark->text();//code = str2Array(ui->lineEditWaterMark->text());code = WT.str2Array(str2);delete validator;}key = WT.generateKey(code.length());int width, height;uchar* buffer = WT.readBmp(ui->lineEdit->displayText().toStdString().data(), width, height);uchar* edge = WT.edgeExtract(buffer, width, height);dst = WT.watermarkImg(buffer, edge, width*height, WT.encode(code, key));//imwrite("encode.bmp", dst*255);image = QPixmap::fromImage(QImage(dst, width, height, QImage::Format_Grayscale8));QGraphicsScene *scene = new QGraphicsScene;scene->addPixmap(image);ui->graphicsViewAfter->setScene(scene);ui->graphicsViewAfter->show();ui->graphicsViewAfter->fitInView(image.rect(), Qt::KeepAspectRatio);//WT.saveBmp(dst);//qDebug() << "In void MainWindow::on_pushButtonEncode_clicked(), height = " << height << endl;//qDebug() << "void MainWindow::on_pushButtonEncode_clicked(), width = " << width << endl;WT.savebmp("encode.bmp", dst, height, width);
}void MainWindow::on_pushButtonDecode_clicked()
{int width, height;watermark WT;uchar* buffer = WT.readBmp(ui->lineEdit->displayText().toStdString().data(), width, height);byteArray code = WT.encode(WT.decodeImg(buffer, dst, width, height, key.length()), key);if(ui->comboBoxWaterMark->currentIndex() == 1){QMessageBox::warning(this, "Decode", "The watermark is " + WT.array2byte(code)+"!");}if(ui->comboBoxWaterMark->currentIndex() == 0){QMessageBox::warning(this, "Decode", "The watermark is " + WT.array2str(code)+"!");}
}

#include "mainwindow.h"
#include int main(int argc, char *argv[])
{QApplication a(argc, argv);MainWindow w;w.setWindowTitle("8 bits Binary BMP picture watermarking program");w.show();return a.exec();
}

IDI_ICON1 ICON DISCARDABLE t14.ico

#-------------------------------------------------
#
# Project created by QtCreator 2019-02-03T14:45:42
#
#-------------------------------------------------QT       += core guigreaterThan(QT_MAJOR_VERSION, 4): QT += widgetsTARGET = ImageEncodeAdvance
TEMPLATE = app# The following define makes your compiler emit warnings if you use
# any feature of Qt which has been marked as deprecated (the exact warnings
# depend on your compiler). Please consult the documentation of the
# deprecated API in order to know how to port your code away from it.
DEFINES += QT_DEPRECATED_WARNINGS# You can also make your code fail to compile if you use deprecated APIs.
# In order to do so, uncomment the following line.
# You can also select to disable deprecated APIs only up to a certain version of Qt.
#DEFINES += QT_DISABLE_DEPRECATED_BEFORE=0x060000    # disables all the APIs deprecated before Qt 6.0.0#INCLUDEPATH += $$PWD/opencv/include
INCLUDEPATH += /usr/local/include \/usr/local/include/opencv4 \/usr/local/include/opencv4/opencv2#LIBS += $$PWD/opencv/lib/opencv_world320.lib
LIBS += /usr/local/lib/libopencv_calib3d.so.4.0.1 \/usr/local/lib/libopencv_core.so.4.0.1 \/usr/local/lib/libopencv_features2d.so.4.0.1 \/usr/local/lib/libopencv_flann.so.4.0.1 \/usr/local/lib/libopencv_highgui.so.4.0.1 \/usr/local/lib/libopencv_imgcodecs.so.4.0.1 \/usr/local/lib/libopencv_imgproc.so.4.0.1 \/usr/local/lib/libopencv_ml.so.4.0.1 \/usr/local/lib/libopencv_objdetect.so.4.0.1 \/usr/local/lib/libopencv_photo.so.4.0.1 \/usr/local/lib/libopencv_stitching.so.4.0.1 \/usr/local/lib/libopencv_videoio.so.4.0.1 \/usr/local/lib/libopencv_video.so.4.0.1SOURCES += \main.cpp \mainwindow.cpp \bmputil.cpp \addditional_utility.cppHEADERS += \mainwindow.h \bmputil.h \addditional_utility.hFORMS += \mainwindow.uiTARGET = "The Binary BMP Picture watermarking Program"RC_FILE = my_app.rc

#include   
#include   
#include   
#include 
#include 
#include 
#include 
#include #include   
#pragma pack(2)
using namespace std;  typedef struct BITMAPFILEHEADER  
{   u_int16_t bfType;   u_int32_t bfSize;   u_int16_t bfReserved1;   u_int16_t bfReserved2;   u_int32_t bfOffBits;   
}BITMAPFILEHEADER;   
typedef struct BITMAPINFOHEADER  
{   u_int32_t biSize;   u_int32_t biWidth;   u_int32_t biHeight;   u_int16_t biPlanes;   u_int16_t biBitCount;   u_int32_t biCompression;   u_int32_t biSizeImage;   u_int32_t biXPelsPerMeter;   u_int32_t biYPelsPerMeter;   u_int32_t biClrUsed;   u_int32_t biClrImportant;   
}BITMAPINFODEADER;  class OperateBMP
{
public:int readBmp();int saveBmp();void work();void compareBMP();~OperateBMP();private:BITMAPFILEHEADER BMFH;BITMAPINFODEADER BMIH;int biWidth;            //图像宽  int biHeight;           //图像高  int biBitCount;         //图像类型，每像素位数  unsigned char *pBmpBuf; //存储图像数据  int lineByte;           //图像数据每行字节数string originFileName_string;string newFileName_string;const char *originFileName_char;const char *newFileName_char;string pictureName_without_type;//string resultName_string;char *resultName_char;int lengthOfResultName_char;string compareResult;
};OperateBMP::~OperateBMP()
{/*delete pBmpBuf;delete originFileName_char;delete newFileName_char;delete resultName_char;*/
}int OperateBMP::readBmp()  
{  FILE *fp;  cout << "Please input the name of the origin BMP file:" << endl;cin >> originFileName_string;originFileName_char = new char[ originFileName_string.length() ];originFileName_char = originFileName_string.c_str();if( (fp = fopen(originFileName_char,"rb")) == NULL)  //以二进制的方式打开文件  {  cout<<"The file "<

clc;
close all;path1 = '/home/u25th_engineer/lena256.bmp';
path2 = '/home/u25th_engineer/lena256_encode.bmp';
path3 = '/home/u25th_engineer/lena256_encode_noise.bmp';X = imread(path1);Y = imread(path2);%KKK = filter2( fspecial('sobel'), Y );
%X1 = mat2gray(Y)
%Z=X1;
Z = imnoise(Y, 'salt & pepper');%添加椒盐噪声，也可以改成其他噪声
%Z = imnoise(Y,'gaussian',0.1,0.004);
imwrite( Z, path3 );
%A = fspecial('average',3); %生成系统预定义的3X3滤波器  figure;
subplot(1, 3, 1); imshow(X); title('a.Lena256原图像');
subplot(1, 3, 2); imshow(Y); title('b.加水印图像');
subplot(1, 3, 3); imshow(Z); title('c.水印图像加噪声');X = double(X); 
Y = double(Y); D = Y - X;
MSE = sum(D(:).*D(:)) / numel(Y);%均方根误差MSEd = 0;
e = 0;
file_name = path1;       
cover_object = double(imread(file_name));
Mc = size(cover_object, 1);     %原图像行数
Nc = size(cover_object, 2);     %原图像列数
file_name = path2;
watermarked_image = double(imread(file_name));
%计算信噪比
for i = 1 : Mcfor j = 1 : Nca(i, j) = cover_object(i, j) ^ 2;b(i, j) = cover_object(i, j) - watermarked_image(i, j);d = d + a(i, j);e = e + b(i, j) ^ 2;end
end
PSNR = 10 * log10(d / e);MAE = mean(mean(abs(D)));%平均绝对误差w = fspecial('gaussian', 11, 1.5); %window 加窗  
K(1) = 0.01; 
K(2) = 0.03; 
L = 255; 
Y = double(Y); 
X = double(X); 
C1 = (K(1)*L) ^ 2;
C2 = (K(2)*L) ^ 2; 
w = w/sum(sum(w)); 
ua = filter2(w, Y, 'valid');%对窗口内并没有进行平均处理，而是与高斯卷积，  
ub = filter2(w, X, 'valid'); % 类似加权平均  
ua_sq = ua.*ua;
ub_sq = ub.*ub;
ua_ub = ua.*ub; 
siga_sq = filter2(w, Y.*Y, 'valid') - ua_sq; 
sigb_sq = filter2(w, X.*X, 'valid') - ub_sq; 
sigab = filter2(w, Y.*X, 'valid') - ua_ub;
ssim_map = ((2*ua_ub + C1).*(2*sigab + C2))./((ua_sq + ub_sq + C1).*(siga_sq + sigb_sq + C2)); 
MSSIM = mean2(ssim_map); display(MSE);%均方根误差MSE
display(PSNR);%峰值信噪比
display(MAE);%平均绝对误差
display(MSSIM);%结构相似性SSIM%display(d);
%display(e);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
X = double(X); 
Y = double(Z); D = Y - X;
MSE1 = sum(D(:).*D(:)) / numel(Y);%均方根误差MSEd = 0;
e = 0;
file_name = path1;       
cover_object = double(imread(file_name));
Mc = size(cover_object, 1);     %原图像行数
Nc = size(cover_object, 2);     %原图像列数
file_name = path3;
watermarked_image = double(imread(file_name));
%计算信噪比
for i = 1 : Mcfor j = 1 : Nca(i, j) = cover_object(i, j) ^ 2;b(i, j) = cover_object(i, j) - watermarked_image(i, j);d = d + a(i, j);e = e + b(i, j) ^ 2;end
end
PSNR1 = 10 * log10(d / e);MAE1 = mean(mean(abs(D)));%平均绝对误差w = fspecial('gaussian', 11, 1.5); %window 加窗  
K(1) = 0.01; 
K(2) = 0.03; 
L = 255; 
Y = double(Y); 
X = double(X); 
C1 = (K(1)*L) ^ 2;
C2 = (K(2)*L) ^ 2; 
w = w/sum(sum(w)); 
ua = filter2(w, Y, 'valid');%对窗口内并没有进行平均处理，而是与高斯卷积，  
ub = filter2(w, X, 'valid'); % 类似加权平均  
ua_sq = ua.*ua;
ub_sq = ub.*ub;
ua_ub = ua.*ub; 
siga_sq = filter2(w, Y.*Y, 'valid') - ua_sq; 
sigb_sq = filter2(w, X.*X, 'valid') - ub_sq; 
sigab = filter2(w, Y.*X, 'valid') - ua_ub;
ssim_map = ((2*ua_ub + C1).*(2*sigab + C2))./((ua_sq + ub_sq + C1).*(siga_sq + sigb_sq + C2)); 
MSSIM1 = mean2(ssim_map); display(MSE1);%均方根误差MSE
display(PSNR1);%峰值信噪比
display(MAE1);%平均绝对误差
display(MSSIM1);%结构相似性SSIM%display(d);
%display(e);