[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)
blank

相关推荐

  • checkstyle使用_idea checkstyle

    checkstyle使用_idea checkstyleCheckstyle是一款检查java程序代码样式的工具,可以有效的帮助我们检视代码以便更好的遵循代码编写标准,特别适用于小组开发时彼此间的样式规范和统一。Checkstyle提供了高可配置性,以便适

  • 一键生成H5_h5生成器

    一键生成H5_h5生成器H5模板代码一键生成器2021-01技术交流或商业合作:QQ(591033633)微信(15858293899)演示地址见文中目录1、功能介绍22、准备工作53、演示地址及步骤53.1、后台数据接口图形化定义63.2、后台数据接口线上调试63.3、前端页面模板定义73.4、生成页面demo解析11功能介绍笔者作为多年的资深开发,深知调试前端UI是后端开发人员最惧怕的事,并且厌倦了无止境的Ctrl+C和Ctrl+V工作…

  • bitr函数的用法_average函数怎么用

    bitr函数的用法_average函数怎么用BitBlt函数用法  2010-03-0510:25:51|  分类: VC学习|举报|字号 订阅      下载LOFTER我的照片书  |(一)BOOLBitBlt(intx,inty,intnWidth,intnHe

    2022年10月19日
  • MES项目导入-SMT行业解决方案知识

    SMT行业解决方案 随着电子技术的快速发展和电子产品生命周期的不断变短,给PCB板组装制造商提出了越来越高的要求,同时,全球化的市场竞争使企业面临的国际压力越来越大。现代PCB板组装制造企业要想在瞬息万变的市场竞争中站稳脚跟,得以生存发展,必须采用先进的生产模式,快速响应客户订单需求,并提供优质低价的产品。系统集成和信息集成(CIMS)是当前提高企业管理效率的重要技术。

  • 网线之RJ45接口定义及网线线序[通俗易懂]

    网线之RJ45接口定义及网线线序[通俗易懂]RJ45接口定义常见的RJ45接口有两类:用于以太网网卡、路由器以太网接口等的DTE类型,还有用于交换机等的DCE类型。DTE我们可以称做“数据终端设备”,DCE我们可以称做“数据通信设备”。从某种意义来说,DTE设备称为“主动通信设备”,DCE设备称为“被动通信设备”。当两个类型一样的设备使用RJ45接口连接通信时,必须使用交叉线连接。这个可以从如下的RJ45DTE类型引脚定义

  • C++实现堆排序算法[通俗易懂]

    1.实现堆排序算法   用C++实现一个堆排序。2.实现思想① 先将初始文件R[1..n]建成一个大根堆,此堆为初始的无序区 ② 再将关键字最大的记录R[1](即堆顶)和无序区的最后一个记录R[n]交换, 由此得到新的无序区R[1..n-1]和有序区R[n],且满足R[1..n-1].keys≤R[n].key。③ 由于交换后新的根R[1]可能违反堆性质,故应将当前无序…

发表回复

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

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