function rate=RSana(input) %输入为512*512的图片--RSana('lenna.bmp')
W = imread(input);
W = double(W);
W=W(:)'; % 将像素矩阵转化为一维数组
[m,n] = size(W);
k = n/4;
asd = W;
R = zeros(2, 4); %用来存储两组RS的四个值
M = [1;0;0;1]; %用来实现随即翻转
bloc = zeros(4, 1);
for j = 1 : k %做非负和非正翻转,再计算相关性
bloc=asd((j-1)*4+1:j*4); %每四个像素放在一起
summ(1) = Abf(bloc);
summ(2) = Abf(posturn(bloc,M));
summ(3) = Abf(negturn(bloc,M));
if summ(2) > summ(1) %Rm
R(1,1) = R(1,1) + 1;
end
if summ(2) < summ(1) %Sm
R(1,2) = R(1,2) + 1;
end
if summ(3) > summ(1) %R_m
R(1,3) = R(1,3) + 1;
end
if summ(3) < summ(1) %S_m
R(1,4) = R(1,4) + 1;
end
end
asd = posturn(asd, ones(n, 1)); %先对载体全做正翻转,接着做非负和非正翻转,再计算相关性
for j = 1 : k
bloc=asd((j-1)*4+1:j*4);
summ(1) = Abf(bloc);
summ(2) = Abf(posturn(bloc,M));
summ(3) = Abf(negturn(bloc,M));
if summ(2) > summ(1) %Rm
R(2,1) = R(2,1) + 1;
end
if summ(2) < summ(1) %Sm
R(2,2) = R(2,2) + 1;
end
if summ(3) > summ(1) %R_m
R(2,3) = R(2,3) + 1;
end
if summ(3) < summ(1) %S_m
R(2,4) = R(2,4) + 1;
end
end
R = R/k;
y = (R(1,3)-R(1,4))-(R(1,1)-R(1,2));
dpz = R(1,1) - R(1,2); dpo = R(2,1) - R(2,2);
dnz = R(1,3) - R(1,4); dno = R(2,3) - R(2,4);
C = [2 * (dpo + dpz), (dnz - dno - dpo - 3 * dpz), (dpz - dnz)];
r = roots(C);
p = r./(r - 0.5);
rate=p(2);
RS.rar_RS LSB_RS 隐藏分析_rs分析_rs隐写分析_隐写分析
版权申诉
197 浏览量
2022-09-22
20:01:50
上传
评论 2
收藏 1KB RAR 举报
JonSco
- 粉丝: 75
- 资源: 1万+
最新资源
- bj计算机思维citel
- 课程设计作业-SpringBoot企业产品仓库管理系统(文档 + 源码 + 数据库文件)
- python相关知识点说明
- 基于 STM32 实现交通灯控制程序源码+详细文档+全部资料(优秀项目).zip
- linux简介及基础教程及特点阐述.txt
- 数据结构课程设计课作业资料
- java简介及基础教程及特点阐述.txt
- 一个简单的数据分析脚本示例,它使用Python和pandas库来读取CSV文件、执行一些基本的数据清洗和转换,并进行简单的可视化
- AD7124+PT100 基于STM32f103的代码 检测多路温度(100~-200)并输出
- CHROME webdriver Version: 125.0.6422.78 win64
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
评论0