[c0,s0] = wavedec2(M1, N, wtype);「建议收藏」

[c0,s0] = wavedec2(M1, N, wtype);「建议收藏」采用小波分解行clear;X1=imread(‘20170901204704.jpg’);imshow(X1);figure;X2=imread(‘20170901204704_副本.jpg’);M1=double(X1)/256;imshow(M1);M2=double(X2)/256;N=4;wtype=’sym4′;[c0,s0]=wavedec2(M1,N,wtype);

大家好,又见面了,我是你们的朋友全栈君。

采用小波分解行

clear ;
X1 = imread(‘20170901204704.jpg’);  imshow(X1); figure;
M1 = double(X1) / 256;   imshow(M1);
N = 4;
wtype = ‘sym4’;
[c0,s0] = wavedec2(M1, N, wtype);

M1导入后为809*803*3,numel(M1) = 1948881;

而length(c0) = 2011320;

两者并非是按照小波变换的定义数获得的。

下例中,作者认为原来程序的低频部分仅为 c0(1:s1(1,1))和c1(1:s1(1,1)是理解错误,这部分是缺少全部的低频信息。他认为低频应该包括c1(1:s1(1,1)*s1(1,2))和c2(1:s1(1,1)*s1(1,2)。

按照小波变换的定义 

[c0,s0] = wavedec2(M1, N, wtype);「建议收藏」

c0(1:s1(1,1)) 包含的是cAn的部分信息。

c1(1:s1(1,1)*s1(1,2))部分则包含的是真正的信息。这个与

 说说wavedec2函数_hugebawu的博客-CSDN博客_wavedec2函数

说明的一致。

使用MATLAB工具wfusimg函数进行图像的融合 – it610.com

x1=imread(‘cs1.jpg’);
x2=imread(‘cs2.jpg’);
x1=double(x1)/256;
x2=double(x2)/256;
xfus=wfusimg(x1,x2,’sym4′,5,’max’,’max’);  %使用wfusimg函数进行图像的融合,其实就是把对低频和高频的处理封装了。
figure;

imshow(xfus);

 和上一个方法相比效果要好得多,这是因为分别对低频个高频进行了处理

wfusimg(x1,x2,’sym4′,5,’max’,’max’)说明

x1,x2为需要融合的图像,采用小波:sym4,分解为5层,近似信号取两幅图中绝对值最大的值,细节信号取两幅图中绝对值最大值,

可选:max,min,mean,img1,img2,rand.

以下是转载的具体方式:

当然,我们也可以不使用小波工具箱中给出的现成函数,而选择自己动手实现基于小波的图像融合算法。这对我们实际了解该算法是如何进行的大有裨益。下面这段示例程序实现了与前面代码段相同的作用,请读者留意注释说明的部分。另需说明的是,下面这段程序仅仅是为了演示算法实现而编写的,因此并没做异常处理的考虑,我们默认待处理的两幅原始图像的尺寸是一样的。

X1 = imread(‘cathe1.bmp’);

X2 = imread(‘cathe2.bmp’);

M1 = double(X1) / 256;

M2 = double(X2) / 256;

N = 4;

wtype = ‘sym4’;

[c0,s0] = wavedec2(M1, N, wtype);

[c1,s1] = wavedec2(M2, N, wtype);

length = size(c1);

Coef_Fusion = zeros(1,length(2));

%低频系数的处理,取平均值

Coef_Fusion(1:s1(1,1)) = (c0(1:s1(1,1))+c1(1:s1(1,1)))/2;

%处理高频系数,取绝对值大者,这里用到了矩阵乘法

MM1 = c0(s1(1,1)+1:length(2));

MM2 = c1(s1(1,1)+1:length(2));

mm = (abs(MM1)) > (abs(MM2));

Y  = (mm.*MM1) + ((~mm).*MM2);

Coef_Fusion(s1(1,1)+1:length(2)) = Y;

%重构

Y = waverec2(Coef_Fusion,s0,wtype);

imshow(Y,[]);

       上面这段代码中用到了小波重构函数waverec2(),函数waverec2()其实就是wavedec2()的相反过程,该函数的常用语法形式如下:

X = waverec2(C,S,’wname’)

X = waverec2(C,S,Lo_R,Hi_R)

       其中,参数C、S、N和’wname’的意义都与二维小波分解函数wavedec2()中定义的一样,这里不再赘言。Lo_R是重构低通滤波器,Hi_R是重构高通滤波器。

自己修改后的代码:(主要是感觉上面的代码没有把低频部分全部计算进去)

clc;
clear all;
p=imread('12.jpg');
q=imread('21.jpg');
p=double(p)/256;
q=double(q)/256;
imshow(p);
figure;
imshow(q);
figure;
[c1,s1]=wavedec2(p,4,'sym4');
[c2,s2]=wavedec2(q,4,'sym4');
length=length(c1);
hecheng=zeros(1,length);
hecheng(1:s1(1,1)*s1(1,2))=(c1(1:s1(1,1)*s1(1,2))+c2(1:s1(1,1)*s1(1,2))/2;
MM1=c1(s1(1,1)*s1(1,2)+1:length);
MM2=c2(s1(1,1)*s1(1,2)+1:length);
mm=(abs(MM1)>abs(MM2));
Y=(mm.*MM1)+(~mm.*MM2);
hecheng(s1(1,1)*s1(1,2)+1:length)=Y;
Y=waverec2(hecheng,s2,'sym4');
imshow(Y);
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/163777.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)


相关推荐

  • django nginx部署_django apache部署

    django nginx部署_django apache部署centos+nginx+uwsgi部署django项目上线

  • 大型电商网站架构案例和技术架构【推荐】

    大型电商网站架构案例和技术架构【推荐】大型网站架构是一个系列文档,欢迎大家关注。本次分享主题:电商网站架构案例。从电商网站的需求,到单机架构,逐步演变为常用的,可供参考的分布式架构的原型。除具备功能需求外,还具备一定的高性能,高可用,可伸缩,可扩展等非功能质量需求(架构目标)。根据实际需要,进行改造,扩展,支持千万PV,是没问题的。本次分享大纲电商案例的原因电商网站需求网站初级架构系统容量估算网站架构分析网站架构优化架构总结电商网站案…

  • 一句话木马怎么连接_js木马源码

    一句话木马怎么连接_js木马源码“EASYNEWS新闻管理系统v1.01正式版”是在企业网站中非常常见的一套整站模版,在该网站系统的留言本组件中就存在着数据过滤不严漏洞,如果网站是默认路径和默认文件名安装的话,入侵者可以利用该漏洞直接上传ASP木马程序控制整个网站服务器。Step1搜索入侵目标使用了“EASYNEWS新闻管理系统v1.01正式版”的网站,在网站页面的底部版权声明处,往往会有关键字符为“WWW.52EAS…

    2022年10月30日
  • python基础(3)列表list「建议收藏」

    python基础(3)列表list「建议收藏」列表列表特点:是一种序列结构,与元组不同,列表具有可变性,可以追加、插入、删除、替换列表中的元素新增元素appendappend添加一个对象,可以是任意类型a=['zhangsa

  • RDIFramework.NET ━ .NET高速信息系统开发框架钜献 V2.9 版本震撼发布

    RDIFramework.NET ━ .NET高速信息系统开发框架钜献 V2.9 版本震撼发布

  • hexdump用法_linux dump命令

    hexdump用法_linux dump命令本文乃fireaxe原创,使用GPL发布,可以自由拷贝,转载。但转载请保持文档的完整性,并注明原作者及原链接。内容可任意使用,但对因使用该内容引起的后果不做任何保证。作者:fireaxe_hq@hotmail.com博客:fireaxe.blog.chinaunix.net转自:http://blog.chinaunix.net/uid-20528014-id-4087756.html开发时经常会…

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号