24小时热门版块排行榜

返回列表

西贝墨将军

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 5.3
帖子: 10
在线: 2.9小时
虫号: 3754560
注册: 2015-03-21
专业: 集成电路设计与测试

[求助] 图像处理中的关于可视密码的问题——秘密图像的分存与恢复

这是一个（3,4）门限的拉格朗日差值方法的秘密图像分存与恢复的程序，怎么能改成（2,4）门限的呢？
clear all:
close all;
M=imread('0.jpg');
ss=rgb2gray(M);
figure;
imshow(ss);
%生成子图像
[m n]=size(ss);
for i=1:m*n
if ss(i)>250
ss(i)=250;
end
end
s=double(ss)+1;
x=[1 2 3 4];
g1=zeros(m,n);
g2=zeros(m,n);
g3=zeros(m,n);
g4=zeros(m,n);
yy1=zeros(m,n);
yy2=zeros(m,n);
yy3=zeros(m,n);
yy4=zeros(m,n);
y1=zeros(m,n);
y2=zeros(m,n);
y3=zeros(m,n);
y4=zeros(m,n);
for j=1:m*n
a1=mod(2*j,251);
a2=mod(3*j,251);
f=[a1 a2 s(j)];
g1(j)=polyval(f,x(1));
yy1(j)=mod(g1(j),251);
g2(j)=polyval(f,x(2));
yy2(j)=mod(g2(j),251);
g3(j)=polyval(f,x(3));
yy3(j)=mod(g3(j),251);
g4(j)=polyval(f,x(4));
yy4(j)=mod(g4(j),251);
end
y1=uint8(yy1-1);
y2=uint8(yy2-1);
y3=uint8(yy3-1);
y4=uint8(yy4-1);
figure,imshow(y1);
figure,imshow(y2);
figure,imshow(y3);
figure,imshow(y4);
%恢复
l1=x(2)*x(3)/[(x(1)-x(2))*(x(1)-x(3))];
l2=x(1)*x(3)/[(x(2)-x(1))*(x(2)-x(3))];
l3=x(1)*x(2)/[(x(3)-x(1))*(x(3)-x(2))];
rr1=zeros(m,n);
r=zeros(m,n);
for j=1:m*n
rr1(j)=mod(yy1(j)*l1+yy2(j)*l2+yy3(j)*l3,251);
end
r=uint8(rr1-1);
figure,imshow(r);

回复此楼

» 猜你喜欢

计算机、0854电子信息（085401-058412）调剂已经有4人回复
基金申报已经有3人回复
国自然申请面上模板最新2026版出了吗？已经有9人回复
溴的反应液脱色已经有6人回复
纳米粒子粒径的测量已经有7人回复
常年博士招收(双一流，工科) 已经有4人回复
推荐一本书已经有10人回复
参与限项已经有5人回复
有没有人能给点建议已经有5人回复
假如你的研究生提出不合理要求已经有12人回复

1楼 2015-04-07 22:12:48

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

点击这里搜索更多相关资源

科研从小木虫开始，人人为我，我为人人

相关版块跳转我要订阅楼主西贝墨将军的主题更新

返回列表