实验6-阶段复习

已于 2023-04-17 09:38:53 修改
阅读量189 收藏 1
点赞数
文章标签： 计算机视觉 opencv 图像处理
于 2023-03-27 10:36:38 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_46311585/article/details/129788253
版权
该文详细描述了如何使用MATLAB进行图像处理，包括位平面分离来展示图像的各个位平面，RGB图像的颜色分量分离与重组，以及使用LSB(LeastSignificantBit)顺序和随机隐藏技术隐藏与提取信息的过程。还涉及到了4bit直接替换法的信息隐藏方法，并提供了相应的MATLAB代码示例。
摘要由CSDN通过智能技术生成
  • 实验题目（此栏实验前由老师填写）

  • 位平面分离

编码实现：


X=double(rgb2gray(imread('D:\1751\pic\airplane.bmp')));%将原始图像转换为灰度图像，如果本身就是灰度图可以去掉rgb2gray函数
[h w] = size(X);
subplot(3,3,1);
imshow(X,[]);
title('原始图像');
for k=1:8
 for i=1:h
     for j=1:w
      tmp(i,j) = bitget(X(i,j),k);%bitget函数首先将X(i,j)处灰度值分解为二进制串，然后取第k位
      
     end
 end
 subplot(3,3,k+1);
 imshow(tmp,[]);
 ind = num2str(k);
 imwrite(tmp,['D:\1751\pic\',ind,'.bmp'],'bmp')
 imti = ['第',ind,'个位平面'];
 title(imti);
end
    1. 读入一个灰度图。

    1. 按位分离出8个位平面，并保存到磁盘上。

    1. 显示各个位平面。

***4、分别读取前述的8个位平面图，合成为一个灰度图。

 A = imread('D:\1751\pic\houses.bmp');
Array=cell(1,8);
% 显示8个位平面图像
for i=8:-1:1Array{i} = bitshift(bitget(A,i),i-1);
end %%%复原
%%定义全0矩阵
[M,N]=size(Array{1});
B=zeros(M,N,'uint8');for i=1:8B=bitset(B,i,bitget(Array{i},i));
end
figure;
subplot(1,2,1);
imshow(A);title('原始图像');
subplot(1,2,2)
imshow(B);title('复原后图像')

  • RGB图像颜色分量分离

编写代码实现：

1、读入一个RGB图像。

2、按RGB图像颜色分量分离。分别保存到磁盘中。

3、显示各个颜色分量。

 P=imread('D:\1751\pic\flowers.bmp'); %载入图片
figure;
%提取图像三通道信息
R1=P;G1=P;B1=P;
R1(:,:,2)=0;R1(:,:,3)=0;G1(:,:,1)=0;G1(:,:,3)=0;B1(:,:,1)=0;B1(:,:,2)=0;
%保存图像
imwrite(R1,'D:\1751\pic\R1.png','png');
imwrite(G1,'D:\1751\pic\G1.png','png');
imwrite(B1,'D:\1751\pic\B1.png','png');
% 显示图像
subplot(1,3,1);imshow(R1,[]);title('R通道');
subplot(1,3,2);imshow(G1,[]);title('G通道');
subplot(1,3,3);imshow(B1,[]);title('B通道');

4、分量合成：

   编辑代码实现：

（1）读取上述的三个颜色分量对应的图片

（2）合成为RGB图像并显示出来。

R=imread('D:\1751\pic\R1.png');
G=imread('D:\1751\pic\G1.png');
B=imread('D:\1751\pic\B1.png');

RGB1(:,:,1)=R(:,:,1);
RGB1(:,:,2)=G(:,:,2);
RGB1(:,:,3)=B(:,:,3);


imwrite(RGB1,'D:\1751\pic\RGB1.png');
imshow(RGB1);

三、 LSB顺序隐藏技术

1、创建文本文件f1.txt,保存你的学号及姓名全拼。

2、在所给的图像中任选一个（灰度图或RGB图像），编写代码，使用LSB顺序隐藏技术，把secure01.txt隐藏到图像中去，并保存隐藏后的图像。

3、显示隐藏前后的图像。

4、显示隐藏前后图像的差别图像。

%文件名:lsbhide.m
%程序员:李巍
%编写时间:2004.2.29
%函数功能:本函数将完成在LSB上的顺序信息隐秘
[ste_cover,len_total]=lsbhide('D:\1751\\pic\lenna.bmp','D:\1751\pic\secure01.txt','D:\1751\pic\lena02.bmp');
%参数说明:
%input是信息隐秘载体图像
%file是秘密消息文件
%output是信息隐秘后生成图像
%ste_cover是信息隐秘后图像矩阵 
%len_total是秘密消息的长度，即容量

function  [ste_cover,len_total]=lsbhide(input,file,output)
%读入图像矩阵
cover=imread(input);
ste_cover=cover;
ste_cover=double(ste_cover);
%将文本文件转换为二进制序列
f_id=fopen(file,'r');
[msg,len_total]=fread(f_id,'ubit1'); 

%判断嵌入消息量是否过大
[m,n]=size(ste_cover);
if len_total>m*n
    error('嵌入消息量过大，请更换图像');
end
%p作为消息嵌入位数计数器   
p=1; 
for f2=1:n
    for f1=1:m
        ste_cover(f1,f2)=ste_cover(f1,f2)-mod(ste_cover(f1,f2),2)+msg(p,1);
        if p==len_total
            break;
         end
         p=p+1;
     end
     if p==len_total
         break;
     end
 end
ste_cover=uint8(ste_cover);

imwrite(ste_cover,output);
%显示实验结果 
subplot(1,2,1);imshow(cover); title('  原始图像  ');
subplot(1,2,2);imshow(output);title('隐藏信息的图像');


end

5、提取隐藏的信息，保存为f2.txt.显示文件内容。

%文件名:lsbget.m
%程序员:李巍
%编写时间:2004.2.29
%函数功能:本函数将完成提取隐秘于LSB上的秘密消息
%输入格式举例:
result=lsbget('D:\1751\pic\lena02.bmp',156,'D:\1751\pic\f2.txt');
%参数说明:
%output是信息隐秘后的图像
%len_total是秘密消息的长度
%goalfile是提取出的秘密消息文件
%result是提取的消息 
function result=lsbget(output,len_total,goalfile)
ste_cover=imread(output);
ste_cover=double(ste_cover);
%判断嵌入消息量是否过大
[m,n]=size(ste_cover);
frr=fopen(goalfile,'a');
%p作为消息嵌入位数计数器,将消息序列写回文本文件
p=1;
for f2=1:n
    for f1=1:m
         if bitand(ste_cover(f1,f2),1)==1
             fwrite(frr,1,'ubit1');
             result(p,1)=1;
         else  
             fwrite(frr,0,'ubit1');
             result(p,1)=0;
         end
         if p==len_total
             break;
         end    
         p=p+1;
    end
    if p==len_total
        break;
    end
end
fclose(frr);
end

***6、LSB随机隐藏的实现：对上述问题，采用随机隐藏的算法再实现一次，并比较隐藏前后图像的差别图像。

%文件名:randlsbhide.m
%程序员:李巍
%编写时间:2004.3.2
%函数功能:本函数将完成随机选择LSB的信息隐秘
%输入格式举例:
[ste_cover,len_total]=randlsbhide('D:\1751\pic\lenna.bmp','D:\1751\pic\secure01.txt','D:\1751\pic\lena03.bmp',213);
%参数说明:
%input是信息隐蔽载体图像
%file是秘密消息文件
%output是信息隐秘后生成图像
%key是随机间隔函数的密钥 
%ste_cover是信息隐秘后图像矩阵 
%len_total是秘密消息的长度，即容量
function  [ste_cover,len_total]=randlsbhide(input,file,output,key)
%读入图像矩阵
cover=imread(input);
ste_cover=cover;
ste_cover=double(ste_cover);
%将文本文件转换为二进制序列
f_id=fopen(file,'r');
[msg,len_total]=fread(f_id,'ubit1'); 
%判断嵌入消息量是否过大
[m,n]=size(ste_cover);
if len_total>m*n
    error('嵌入消息量过大，请更换图像');
end
%p作为消息嵌入位数计数器
p=1;
%调用随机间隔函数选取像素点 
[row,col]=randinterval(ste_cover,len_total,key);
%在LSB隐秘消息 
for i=1:len_total
   ste_cover(row(i),col(i))=ste_cover(row(i),col(i))-mod(ste_cover(row(i),col(i)),2)+msg(p,1);
 if p==len_total
   break;
 end
 p=p+1;
end
ste_cover=uint8(ste_cover);
imwrite(ste_cover,output);
%显示实验结果 
subplot(1,2,1);imshow(cover); title('  原始图像  ');
subplot(1,2,2);imshow(output);title('隐藏信息的图像');
end
%文件名:randlsbget.m
%程序员:李巍
%编写时间:2004.2.29
%函数功能:本函数将完成提取隐秘于LSB上的秘密消息
%输入格式举例:
result=randlsbget('D:\1751\pic\lena03.bmp',156,'D:\1751\pic\secure03.txt',1988)
 
%参数说明:
%output是信息隐秘后的图像
%len_total是秘密消息的长度
%goalfile是提取出的秘密消息文件
%key是随机间隔函数的密钥 
%result是提取的消息 
function result=randlsbget(output,len_total,goalfile,key)
ste_cover=imread(output);
ste_cover=double(ste_cover);
%判断嵌入消息量是否过大
[m,n]=size(ste_cover);
frr=fopen(goalfile,'a');
%p作为消息嵌入位数计数器,将消息序列写回文本文件
p=1;
%调用随机间隔函数选取像素点 
[row,col]=randinterval(ste_cover,len_total,key);
for i=1:len_total
     if bitand(ste_cover(row(i),col(i)),1)==1
         fwrite(frr,1,'ubit1');
         result(p,1)=1;
     else  
         fwrite(frr,0,'ubit1');
         result(p,1)=0;
     end
     if p==len_total
         break;
     end    
     p=p+1;
end
fclose(frr);
end

四、4bit直接替换法

编码实现：

1、读取一个灰度图（cover）（如果是RGB图像，则取其中一个颜色分量），每个像素的低4位置0.

2、读取一个待隐藏的灰度图(secure)，每个像素低4位清0，再将高4位移至低4位

3、将cover的每个像素的低4位用secure的像素低4位替换。

4、保存替换结果，并显示出来（ste-cover）。

cover = imread('D:\1751\pic\houses.bmp') ; 
ste_cover = cover; 
ste_cover=double(ste_cover); 


f_id = fopen('D:\1751\pic\sec.txt', 'r') ; 
[ msg, len_total] = fread( f_id, 'ubit1');
[ m, n] = size( ste_cover) ; 
if len_total > m*n
error( '嵌入的信息量太大') ;
 end

p = 1; 
for f2 = 1:n
 for f1 = 1:m
%ste_cover( f1, f2) = ste_cover( f1, f2) -mod( ste_cover( f1, f2) , 2) + msg( p, 1) ; 
ste_cover( f1, f2) = bitand(ste_cover( f1, f2), 254) ;
 ste_cover( f1, f2) = ste_cover( f1, f2) + msg( p, 1) ;
if p == len_total
 break; 
end 
p = p + 1; 
end 
if p == len_total
 break; 
end
end
 ste_cover = uint8( ste_cover) ; %生成隐密后的图像;
 imwrite( ste_cover, 'D:\1751\pic\houses02.bmp') ; %显示隐藏结果; 
subplot( 1, 2, 1) ; imshow( cover) ; title( '原始图像') ;
 subplot( 1, 2, 2) ; imshow( 'houses02.bmp') ; title( '隐藏信息后的图像');

5、编写代码实现提取功能

（1）读取ste-cover.

（2）定义图像变量，把ste-cover中提取的低4位，保存至该图像变量并移动至高4位。显示该图像变量。

%文件名:lsbget.m

%函数功能:本函数将完成提取隐秘于LSB上的秘密消息
%输入格式举例:result=lsbgetnew('d.bmp',1736,'secret.txt');
%参数说明:
%output是信息隐秘后的图像
%len_total是秘密消息的长度
%goalfile是提取出的秘密消息文件
%result是提取的消息 
%function result=lsbgetnew(output,len_total,goalfile)
ste_cover=imread('D:\1751\pic\houses02.bmp');
ste_cover=double(ste_cover);
%判断嵌入消息量是否过大
[m,n]=size(ste_cover);
frr=fopen('D:\1751\pic\sec02.txt','w');
%p作为消息嵌入位数计数器,将消息序列写回文本文件
p=1;
for f2=1:n
    for f1=1:m
         if bitand(ste_cover(f1,f2),1)==1
             fwrite(frr,1,'ubit1');
             result(p,1)=1;
         else  
             fwrite(frr,0,'ubit1');
             result(p,1)=0;
         end
         if p==len_total
             break;
         end    
         p=p+1;
    end
    if p==len_total
        break;
    end
end
fclose(frr);

五、exer6.bmp、exer7.bmp图像中隐藏有一段文字，请尝试编写提取程序提取隐藏的信息。

%文件名:lsbget.m
%程序员:李巍
%编写时间:2004.2.29
%函数功能:本函数将完成提取隐秘于LSB上的秘密消息
%输入格式举例:
result=lsbgetrgb('D:\1754\pic\exer6.bmp',944,'D:\1754\secret.txt',2);
%参数说明:
%output是信息隐秘后的图像
%len_total是秘密消息的长度
%goalfile是提取出的秘密消息文件
%result是提取的消息 
function result=lsbgetrgb(output,len_total,goalfile,level)
ste_cover=imread(output);
ste_cover=ste_cover(:,:,level);
ste_cover=double(ste_cover);
%判断嵌入消息量是否过大
[m,n]=size(ste_cover);
frr=fopen(goalfile,'a');
%p作为消息嵌入位数计数器,将消息序列写回文本文件
p=1;
for f2=1:n
    for f1=1:m
         if bitand(ste_cover(f1,f2),1)==1
             fwrite(frr,1,'ubit1');
             result(p,1)=1;
         else  
             fwrite(frr,0,'ubit1');
             result(p,1)=0;
         end
         if p==len_total
             break;
         end    
         p=p+1;
    end
    if p==len_total
        break;
    end
end
fclose(frr);
end

 六、蓝盾平台隐写题练习

 登录蓝盾平台：10.10.84.82

（1） 完成“基于LSB的数字水印嵌入与提取”实验

   

我是小豆子
关注
关于MATLAB中使用的UBIT1问题
Carry On
10-21 2900
最近在看基于LSB方法的信息隐藏的时候
信息隐藏与数字水印实验4-LSB信息隐藏的卡方分析
taylorswift的博客
07-07 3722
一、 实验目的 了解什么是隐写分析（steganalysis），隐写分析与信息隐藏和数字水印的关 系。掌握基于图像的 LSB 隐写的分析方法，设计并实现一种基于图像的 LSB 卡方隐写分析方法。二、 实验环境 （1）Windows10 或 11 系统 （2）Python 3.7 或以上 （3）图像文件 man.bmp三、 实验要求 （1）在灰度图 man.bmp 利用 LSB 方法随机嵌入 0 和 1 的比特，嵌入率为 100% ， 后 生 成 mangraystego.bmp ， 在 一 行 两 列 中
实验8-基于DCT的数字水印隐藏及攻击提取技术
m0_46311585的博客
04-17 929
3、读取“huawei.png（或者du.png,或者kgy.jpg）” ，变换为灰度图，尺寸大小设置为64x64（即前一图像的尺寸/8的大小）。1、参考所附代码，读取一个RGB图‘baboon.png’，取得该图的R分量。4、参考所给代码，设置嵌入因子为0.1，对mark.bmp进行DCT变换。6、嵌入后的图像进行DCT逆变换，并保存为withmark.bmp。二、参考所给代码，从withmark.bmp中提取嵌入的水印图像。4、提取隐藏的信息，保存为secure02.txt.显示文件内容。
信息隐藏将txt文件合并到jpg文件中_在图像中隐藏秘密消息Steghide Kali Linux隐写术...
weixin_39708557的博客
12-04 389
在计算机科学中，将信息隐藏在图像，文档，程序，有效载荷，消息，音乐，HTML页面，可移动媒体等文件内部的技术被称为隐写术，其做法是在其他非秘密文本中隐藏消息或信息。或数据。Steghide支持以下图像格式：JPEG，BMP，WAV，AU文件。本教程使用Kali linux进行演示。安装Steghide要安装Steghide，请在Linux中打开一个终端，然后使用以下命令。apt-get insta...
信息隐藏（2）---直接4bit替换
maff的专栏
11-23 4178
所谓直接4bit 替换法, 就是直接用秘密图像像素值的高4bit 去替换载体图像像素值的低4bit % 函数功能: 直接将秘密图像的高4bit 隐藏在RGB 载体图像的R, G, B 层中所选的那一层的低4bit, % 并将秘密图像提取出来, 最后显示。要求载体图像的大小大于等于秘密图像的大小, 且秘密图像是二值或灰度图像 % 输入格式: % data = imagehide( ‘c: \ l
华为-HCNP-ospf-综合实验-大图-知识点回顾-复习-路由交换-by LeeSir.rar
02-09
【华为 HCNP OSPF 综合实验】是华为认证网络专家（Huawei Certified Network Professional，简称HCNP）课程中的...在复习阶段，不断实践和回顾这些知识点至关重要，因为只有将理论与实践相结合，才能真正掌握网络技术。
吉林大学微机实验代码实验报告及上机复习2021年.zip
12-02
1. 复习实验过程：理解实验的每一个阶段，加深对实验原理的理解。 2. 学习问题解决方法：看到他人在实验中遇到的困难及应对策略，提高自己的问题解决能力。 3. 提升分析能力：通过对实验数据的分析，学习如何从...
生态学-期末复习资料整理.pdf
12-30
它们具有综合作用、主导因子作用、不可替代性和可补偿性、阶段性、直接和间接作用以及限制性作用。例如，Liebig最小因子定律指出，生物生长受限于最稀缺的资源，而Shelford耐性定律则说明生物对环境因子有耐受范围，...
教育心理学--复习资料全.docx
10-10
量的研究则涉及实验法（如实验室实验法和自然实验法）、问卷调查法和测验法，旨在量化数据并进行统计分析。 进入中学生阶段，他们的心理发展具有连续性、顺序性、阶段性、不均衡性和差异性等特点。在认知方面，...
2017届高三化学二轮复习课件：第1部分-专题十三　化学实验基础-考点2.ppt
10-18
6. **实验报告撰写**：学习如何撰写完整的实验报告，包括实验目的、原理、步骤、数据记录、结果分析和讨论等部分。这是将实验操作转化为理论知识的关键环节。 7. **特定考点2**：虽然具体考点未明确，但通常可能...
Matlab相关函数
hxxjxw的博客
02-28 1016
imread() 从图形文件读取图像 img=imread('su.bmp') matlab读入图像的数据是uint8，而matlab中数值一般采用double型（64位）存储和运算。所以要先将图像转为double格式的才能运算 img=double(img) imwrite() 将图像写入图形文件 %将一个 50×50 的灰度值数组写入myGray.pn...
实验12-信息隐藏技术综合复习
最新发布
m0_46311585的博客
05-08 273
1、编写代码exer03_1.m,实现功能：采用4bit替换法“huawei.bmp”隐藏到“kgy.bmp”中（R、G、B层中任选一个），隐藏后的图像保存为“ kgy01.bmp”,显示该图像（标题为你的学号）。二、编写代码exer02.m, 实现功能：读取“pepper.bmp”, 将其R、G、B颜色分量的最低2位清0，再组合成新彩色图像，保存为“pepper01.bmp”，显示该图像（标题为你的学号）***九、继续复习基于LSB的顺序隐藏技术（图像至图像的嵌入，二值化水印图像）（综合复习9的第9题）
基于DCT的信息隐藏
毕业作品网站
10-13 1240
在图片嵌入程序中，首先读取水印图片sss.bmp，通过rgb2gray函数使其变成灰度图，然后通过im2bw函数是水印图像变成二值图，同时记录下水印图像的长rm和宽cm，以同样的方法读取载体图片cover_image.bmp。id=13912 DCT变换把水印图片分成8X8的矩阵之后，通过after(x+i,y+j)=before(x+i,y+j)+alpha*k(1)，在结果矩阵的中频段嵌入水印，然后通过blkproc函数做DCT逆变换后输出保存在markresule.bmp图片中。
LSB图像信息隐藏实验（附源代码）
热门推荐
cheeseandcake的博客
11-01 4万+
LSB图像信息隐藏 1.掌握对图像的基本操作。 2.能够用 LSB 算法对图像进行信息隐藏 3.能够用 LSB 提取算法提取隐藏进图像的信息 4.能够反映 jpeg 压缩率与误码率之间的关系 三、【实验过程】 实验算法 1：信息隐藏 实验算法2：信息提取 实验算法3：图片压缩率与误码率关系
LSB信息隐藏算法
HUANGliang_的博客
12-31 5446
将一张图片藏入另一张图片的算法 3.使用LSB信息隐藏算法，分别嵌入1bit/pixel、2bit/pixel、3bit/pixel，并显示伪装图像。2.制作水印图片，可以嵌入1个位平面的、2个位平面的、3个位平面的。1.将载体图片分成8个位平面，并显示。4.提取嵌入的秘密信息。...
在Linux上实现将文件隐藏到图像中
09-01 367
Steganography(隐写术)是一个宽泛的话题。这个教程只涉及如何将文件掩藏到图像文件中，这也是Steganography的核心。不要认为安全专家不能破解这个方法。 这里描述的方法是非常简单的，因此即使是一个中级的安全专家也可以在几小时内将其破...
MIDI信息隐藏（MATLAB，信息隐藏与数字水印实验）
Hardworking666的博客
12-19 3914
文章目录一、MIDI简介二、MIDI隐藏代码三、取出MIDI隐藏信息 一、MIDI简介 MIDI(Musical Instrument Digital Interface，乐器数字接口)，电子乐器制造商们建立的通信标准，辅助音乐创作、确定电脑音乐程序，用音符的数字控制信号来记录音乐。 一个标准MIDI文件基本上由两部分组成：头块和音轨块。头块用来描述整个MIDI文件基本信息。音轨块则包含一系列由MIDI消息构成的MIDI数据流。原则上，可为某种声音、某种乐谱或某种乐器等分配一个音轨块。 MIDI文件中前..
ZBIT、ZBIT2、UBIT是什么？
emowuyi
09-22 464
下面是我最近遇到的问题，然后再TI官方论坛中找到的回答。 Q1:what is thisZBIT,ZBIT2,UBIT? code:in osal.c void osal_start_system( void ){#if !defined (ZBIT) for(;;) // Forever Loop#endif { uint8 idx = 0; ... A1:ZBITis a c...
项目验收关键：华为HCIE-R&S复习与实验手册
"华为HCIE-R&S学习备考必备复习实验手册" 本文主要探讨了项目管理和项目验收的关键点，尤其强调了在项目收尾阶段的重要工作和准备验收材料的必要性。项目管理是确保项目成功的核心环节，对于HCIE-R&S（路由与交换）...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值