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);

回复此楼

» 猜你喜欢

同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗已经有6人回复
依托企业入选了国家启明计划青年人才。有无高校可以引进的。已经有6人回复
依托企业入选了国家启明计划青年人才。有无高校可以引进的。已经有10人回复
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）已经有9人回复
有院领导为了换新车，用横向课题经费买了俩车已经有10人回复
遇见不省心的家人很难过已经有24人回复
AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性已经有6人回复
酰胺脱乙酰基已经有13人回复
有时候真觉得大城市人没有县城人甚至个体户幸福已经有10人回复

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

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

智能机器人

Robot (super robot)

我们都爱小木虫

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

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

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

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

返回列表