| 查看: 1683 | 回复: 0 | ||
[求助]
hilbert图像置乱算法
|
|
以下是我编的hilbert图像置乱算法，不能运行出置乱后的图像，求高手指点 clear all; I=imread('lena.bmp'); figure,imshow(I),title('原图'); I1=double(I); img=I1; [N,N]=size(img); temp1=zeros(1,N); temp2=zeros(1,N); for X=1:1:N for Y=1:1:N temp1(end+1)=mod(X,2); temp1=(temp1-temp1(end))/2; temp1=temp1(end:-1:1); temp2(end+1)=mod(Y,2); temp2=(temp2-temp2(end))/2; temp2=temp2(end:-1:1); Xj=temp1; Yj=temp2; for k=1:100 Xjk=temp1(k); Yjk=temp2(k); XXjk=Xjk; YYjk=Yjk; Ak=sum((1-XXjk)*(1-YYjk)); Bk=sum(XXjk*YYjk); if mod(Ak,2)>0 XXj(k+1)=Yj(k+1); YYj(k+1)=Xj(k+1); else XXj(k+1)=Xj(k+1); YYj(k+1)=Yj(k+1); end if mod(Bk,2)>0 XXj(k+1)=1-Xj(k+1); YYj(k+1)=1-Yj(k+1); else XXj(k+1)=XXj(k+1); YYj(k+1)=YYj(k+1); end tjk=XXjk+YYjk+2*XXjk*(1-YYjk); j=sum(tjk/4^k); I1(X,Y)=j; I2=I1(X,Y); end end end figure imshow(I2); title('编码图像'); |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
fileCode有新解读？
已经有16人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
基金中了
已经有14人回复
-
关于filecode
已经有7人回复
-
娱乐
已经有6人回复
找到一些相关的精华帖子，希望有用哦~
求教有没有什么算法判断一个网络中节点的拓扑位置是无差异的？
已经有4人回复
- Ni-Cu-P 三元合金的TEM图像是晶粒与非晶的混合网格结构，怎么解释？
已经有26人回复
- 固体废物处理与处置怎么讲更好
已经有10人回复
- mtDNA图像检测不理想，求原因
已经有18人回复
- 【求助】zeiss显微镜CCD （型号为MRC）相应的图像采集卡
已经有5人回复
- 想在广东省内寻求可以调剂的博导，人工智能、模式识别、图像处理、算法研究、机器学习
已经有6人回复
- 关于离职后住房公积金的算法
已经有10人回复
- 【交流】为什么拉伸试验, 断裂位置在平行段之外不好?
已经有10人回复
- origin8画图导出eps图像中文乱码
已经有3人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
100