【图像融合】基于curvelet变换图像融合matlab代码

【图像融合】基于curvelet变换图像融合matlab代码

1 简介

针对Curvelet分解的不同频率域,分别讨论了低频系数和高频系数的选择原则.在选择低频系数时,定义了局部区域标准方差,采用了"选择"与"平均"相结合的系数选择方案;在选择高频系数时,充分利用Curvelet变换具有方向性的优点,提出了Curvelet域方向对比度的概念,并给出了基于方向对比度的系数选择方案.实验结果表明:本文所给出的融合算法能够很好地保留多幅源图像中的有用信息,得到多个目标聚焦都很清晰的图像.

2 部分代码

clc
clear all
close all% fdct_wrapping_demo_recon.m -- Partial curvelet reconstruction.
X=imread('imgwarp.bmp');
X2=imread('IR1022_3.bmp');
X1=rgb2gray(X);% Forward curvelet transform
disp('Take curvelet transform: fdct_wrapping');
tic; C1 = fdct_wrapping(double(X1),0,2); toc;
tic; C2 &#