利用Matlab实现基于小波变换的遥感图像融合
使用Matlab进行多传感器数据融合与信息融合的方法
使用Matlab进行多传感器数据融合与信息融合的方法引言：随着科技的飞速发展，多传感器系统已经成为很多领域中不可或缺的一部分。
多传感器系统可以通过融合来自不同传感器的信息，提高系统的准确性和可靠性。
在这方面，Matlab作为一种强大的计算工具，提供了丰富的函数和算法，为多传感器数据融合与信息融合提供了便捷的解决方案。
一、多传感器数据融合的概念及应用领域1.1 多传感器数据融合的概念多传感器数据融合即利用不同传感器采集到的信息，通过适当的算法和方法，将这些信息结合起来，形成更准确、更全面的信息作为输出。
通过融合来自不同传感器的数据，可以弥补各个传感器本身的局限性，提高系统的感知和决策能力。
1.2 多传感器数据融合的应用领域多传感器数据融合技术在许多领域都有广泛的应用。
例如，在智能交通系统中，多传感器数据融合可以实现交通流量监测、车辆行为分析和交通事故预测等功能。
在军事领域，多传感器数据融合可以用于目标跟踪、情报获取和作战决策等方面。
此外，多传感器数据融合还被广泛应用于环境监测、安防监控、无人机控制等领域。
二、多传感器数据融合的方法和算法2.1 传感器数据预处理在进行数据融合之前，首先需要对传感器采集到的数据进行预处理。
常见的预处理方法包括噪声滤波、数据校正和数据对齐等。
噪声滤波可以有效地去除传感器数据中的噪声，提高数据的质量。
数据校正可以将不同传感器的数据转化为统一的参考坐标系，便于数据融合。
数据对齐可以解决不同传感器采样率不一致的问题，将数据对齐到统一的时间轴上。
2.2 传感器数据融合方法传感器数据融合方法根据融合的目标和需求可以分为两种：低层数据融合和高层数据融合。
低层数据融合主要是将多个传感器的原始数据融合在一起，形成一个多源数据集，如加权平均、最大值和最小值融合等。
高层数据融合则是在低层融合的基础上，进一步分析和处理融合后的数据，提取更高级别的信息，如卡尔曼滤波、粒子滤波和神经网络等。
2.3 传感器数据融合算法在Matlab中，有丰富的函数和算法可用于多传感器数据融合。
brovey变换融合算法matlab
Brovey变换融合算法MATLAB实现一、引言Brovey变换是一种常用于遥感图像融合的方法，它能够有效地融合多光谱和全色图像，生成既有全色图像的空间细节，又保留多光谱图像的光谱信息的融合图像。
以下是在MATLAB中实现Brovey变换融合算法的步骤。
二、算法原理Brovey变换的基本思想是对每一个像素进行归一化处理，然后再将归一化后的图像进行乘法运算。
这样可以增强图像的对比度，提高图像的空间分辨率，同时保持原有的光谱特性。
具体算法步骤如下：1. 对多光谱图像（MSI）和全色图像（PAN）进行归一化处理。
2. 对归一化后的MSI和PAN进行乘法运算。
3. 对得到的乘积图像进行反归一化处理，得到最终的融合图像。
三、MATLAB实现以下是基于上述算法原理的MATLAB实现代码：```Matlabfunction fused_img = Brovey_Transform(MSI, PAN)% 输入：MSI - 多光谱图像，PAN - 全色图像% 输出：fused_img - 融合后的图像% 将MSI和PAN转换为double类型，方便后续计算MSI = double(MSI);PAN = double(PAN);% 获取MSI和PAN的大小[M, N, ~] = size(MSI);[P, Q] = size(PAN);% 确保MSI和PAN大小一致，如果不一致，则需要进行重采样或其他处理if M ~= P || N ~= Qerror('MSI和PAN的大小不一致，无法进行Brovey变换');end% 对MSI和PAN进行归一化处理MSI_norm = MSI ./ repmat(sum(MSI, 3), [1, 1, size(MSI, 3)]);PAN_norm = PAN ./ sum(PAN(:));% 对归一化后的MSI和PAN进行乘法运算fused_img = zeros(size(MSI));for i = 1:size(MSI, 3)fused_img(:, :, i) = MSI_norm(:, :, i) .* PAN_norm;end% 对得到的乘积图像进行反归一化处理，得到最终的融合图像min_val = min(min(min(fused_img)));max_val = max(max(max(fused_img)));fused_img = (fused_img - min_val) / (max_val - min_val); % 归一化到[0,1]区间fused_img = uint8(fused_img * 255); % 转换为uint8类型，方便显示和存储end```请注意，这只是一个基础的Brovey变换实现，实际应用中可能需要进行一些优化和调整，例如对输入图像进行预处理（如去噪、增强等），对输出图像进行后处理（如颜色校正、对比度拉伸等）。
小波变换-图像融合matlab代码
%对图像进行放大算子的运算
PIC3 = conv2(conv2(Y, 2*fw, 'valid'),2*fw','valid');
%第i1级图像重构；
PIC1 = PIC3 + L{i};
%选取图像范围
PIC1 = PIC1(1:k1(i),1:k2(i));
% %end;
% %end;
% %%
% %for k=1:256,
% % p1(k)=p1(k)/(d);
% % p2(k)=p2(k)/(d);
% %end;
%
% %for i=1:256
% % for j=1:256
% % p3(i,j) =p3(i,j)/(d);
%图像隔行隔列插值扩展恢复到原尺寸图像
[c d] = size(Y4);
Y6 = zeros(2*c, 2*d);
Y6(1:2:2*c,1:2:2*d) = Y4;
Y7 = zeros(2*c, 2*d);
%PIC2 = X2;
X1 = PIC1;
X2 = PIC2;
%定义滤波窗口；
fw = 1/16.*[1 4 6 4 1];
z =3;
L = cell(1,z);
L1 = cell(1,z);
for i = 1:z % N1
G3 = conv2(conv2(Y3, fw, 'valid'),fw', 'valid');
%将图像进行隔行隔列减半
[a b] = size(G2);
Y4 = G2(1:2:a, 1:2:b);
Matlab中的遥感图像处理与分析方法
Matlab中的遥感图像处理与分析方法遥感图像处理与分析是遥感技术的重要组成部分，它包括对获取的遥感图像进行预处理、增强、分类、信息提取等一系列操作。
Matlab作为一种强大的科学计算软件，提供了丰富的工具箱和函数，广泛应用于遥感图像处理与分析领域。
本文将介绍Matlab中一些常用的遥感图像处理与分析方法，并探讨其在实际应用中的价值。
一、遥感图像的读取与显示在进行遥感图像处理与分析之前，首先需要将遥感图像读取到Matlab中。
Matlab提供了多种读取图像的函数，如imread、multibandread等。
通过这些函数，可以将遥感图像以矩阵的形式存储在Matlab的变量中，方便后续的处理。
读取遥感图像后，我们可以使用imshow函数在Matlab中显示图像。
通过调整imshow函数的参数，可以实现对图像的缩放、亮度、对比度等的调整。
此外，Matlab还提供了imtool函数，可以在一个窗口中同时显示多幅图像，方便进行比较和分析。
二、遥感图像的预处理遥感图像的预处理是遥感图像处理与分析的重要步骤之一。
预处理的目的是消除图像中的噪声、增强图像的对比度、调整图像的亮度等，为后续的处理提供更好的数据基础。
在Matlab中，可以使用多种函数实现遥感图像的预处理。
例如，imadjust函数可以调整图像的亮度和对比度，imnoise函数可以在图像中添加噪声，medfilt2函数可以进行中值滤波，去除图像中的椒盐噪声等。
此外，Matlab还提供了一些专门用于遥感图像处理的工具箱，如Image Processing Toolbox、Computer Vision Toolbox等，这些工具箱提供了丰富的函数和工具，便于进行图像的预处理操作。
三、遥感图像的增强与融合遥感图像的增强与融合是遥感图像处理与分析的重要任务之一。
增强可以使图像中的细节更加清晰，对于提取图像中的信息非常有帮助。
融合可以将来自不同传感器或不同时刻的遥感图像融合在一起，得到更全面的信息。
MATLAB中的遥感图像处理方法解析
MATLAB中的遥感图像处理方法解析遥感图像处理是一项重要的技术，广泛应用于农业、环境保护、城市规划等领域。
MATLAB作为一种强大的科学计算软件，提供了许多有效的图像处理工具和算法，使得遥感图像的处理更加简便高效。
本文将通过几个实例，介绍MATLAB 中常用的遥感图像处理方法。
一、图像预处理遥感图像通常存在一些噪声和失真。
为了提高图像质量和后续分析的精确性，需要对图像进行预处理。
MATLAB提供了各种滤波器和降噪算法，如中值滤波、高斯滤波和小波变换。
这些方法可以降低图像中的噪声，并使细节更加清晰。
二、图像增强图像增强是提升图像视觉效果的重要方法。
在遥感图像处理中，一般采用直方图均衡化和对比度拉伸等方法。
直方图均衡化可以使图像的亮度分布更均匀，增强图像的视觉效果。
对比度拉伸则通过扩展图像的动态范围，使得图像中的细节更加丰富。
三、影像分割影像分割是将图像分割成不同的区域或目标的过程。
MATLAB提供了多种分割算法，如基于阈值的分割、基于区域的分割和基于边缘的分割。
这些方法可以帮助我们从遥感图像中提取出感兴趣的目标，为后续的分析提供有效的数据。
四、特征提取特征提取是从遥感图像中提取出有意义的特征信息的过程。
常用的特征包括纹理特征、形状特征和光谱特征等。
MATLAB提供了一系列用于特征提取的函数和工具箱，如灰度共生矩阵、哈尔小波变换和主成分分析等。
这些方法可以帮助我们从遥感图像中提取出有价值的特征，用于后续的分类和识别任务。
五、图像分类图像分类是将图像分成不同的类别或类别的过程。
在遥感图像处理中，一般采用监督学习和无监督学习的方法。
监督学习需要样本标注数据，可以通过支持向量机和随机森林等算法进行分类。
无监督学习则不需要标注数据，常用的方法有k均值聚类和自组织映射网络等。
MATLAB提供了这些算法的实现和函数，方便我们进行遥感图像的分类和识别。
六、图像融合图像融合是将多个传感器或多个波段的图像进行融合，得到更全面、更丰富的信息的过程。
在Matlab中进行图像融合与图像叠加的方法与技巧
在Matlab中进行图像融合与图像叠加的方法与技巧引言：随着数字图像处理和计算机视觉领域的发展，图像融合和图像叠加变得越来越重要。
图像融合是指将多幅图像合成为一幅具有更清晰、更丰富信息的图像，而图像叠加则是在保留所叠加图像的原始信息的同时，使图像更加丰富和易于理解。
Matlab作为一种强大的科学计算工具，提供了丰富的图像处理函数和工具箱，可以很方便地进行图像融合与图像叠加。
一、图像融合的方法与技巧1. 融合算法图像融合的基本方法有加权平均法、空间域融合法、频域融合法、小波融合法等。
加权平均法是最简单的方法，通过计算图像像素的平均值来融合。
空间域融合法是通过对直接融合的图像进行空间域操作来提取融合结果。
频域融合法则是通过将图像转换到频域，然后进行频域操作来实现融合。
小波融合法是基于小波变换的方法，利用小波分析的多尺度分解能力对图像进行分析和融合。
根据具体需求和图像的特点，选择合适的融合算法是非常重要的。
2. 图像预处理在进行图像融合之前，通常需要进行图像预处理，以提高融合结果的质量。
常用的图像预处理方法包括灰度拉伸、直方图均衡化、滤波等。
灰度拉伸是通过对图像的像素值进行线性变换，将图像像素值的范围拉伸到合适的范围内，从而增加图像的对比度。
直方图均衡化则是将图像的像素值在灰度直方图上均匀分布，以增强图像的细节。
滤波是通过对图像进行滤波操作，如低通滤波、高通滤波等，以去除图像中的噪声和不需要的细节。
3. 图像融合的策略图像融合的策略可以根据具体需求来选择。
常见的策略包括全局融合和局部融合。
全局融合是将所有图像的信息进行融合，得到整体的融合结果。
而局部融合则是将不同图像的不同区域进行融合，以保留更多的细节和纹理。
根据具体应用和需求，选择合适的融合策略可以使融合结果更加符合实际需求。
4. 参数设置与调整在进行图像融合过程中，不同的算法和方法有各自的参数，根据不同的图像和具体应用，需要适时地进行参数的设置和调整。
图像处理matlab及图像融合图像镶嵌图像拼接
图像处理matlab及图像融合图像镶嵌图像拼接在实际的对图像处理过程中，由于我们读出的图像是unit8型，⽽在MATLAB的矩阵运算中要求所有的运算变量为double型（双精度型）。
因此读出的图像数据不能直接进⾏相加求平均，因此必须使⽤⼀个函数将图像数据转换成双精度型数据。
MATLAB中提供了这样的函数：im2double函数，其语法格式为：I2 = im2double(I1)其中I1是输⼊的图像数据，它可能是unit8或unit16型数据，通过函数的变化输出I2为⼀个double型数据，这样两图像数据就可以⽅便的进⾏相加等代数运算.要把double的图像（范围是0到1）再次转化为256灰度值的，可以这样Igrey= uint8(I2*255)图像类型转换函数：dither() 通过颜⾊抖动，把真彩图像转换成索引图像或灰度图象转换成⼆值图像gray2ind() 将灰度图像（或⼆值图像）转换成索引图像grayslice() 通过设定的阈值将灰度图象转换成索引图像im2bw() 通过设定亮度阈值将灰度、真彩、索引图象转换成⼆值图像ind2gray() 将索引图象转换成灰度图象ind2rgb() 将索引图象转换成真彩⾊图像mat2gray() 将⼀个数据矩阵转换成⼀幅灰度图象rgb2gray() 将真彩转换成灰度图象rgb2ind() 将真彩转换成索引图象图像类型与类型间的转换1。
索引图像：包括⼀个数据矩阵X和⼀个⾊图阵MAP。
矩阵元素值指向MAP中的特定颜⾊向量。
2。
灰度图像：数据矩阵I，I中的数据代表了颜⾊灰度值。
矩阵中的元素可以是double类型、8位或16位⽆符号的整数类型。
3。
RGB图像：即真彩图像。
矩阵中每个元素为⼀个数组，数组的元素定义了像素的红、绿、蓝颜⾊值。
RGB数组可以是double类型、8位或16位⽆符号的整数类型。
4。
⼆值图像：⼀个数据阵列，每个象素只能取0或1。
矩阵的基本运算⾏列式求值：det(A)矩阵加减：+、-矩阵相乘：*矩阵左除：A/B %相当于inv(A)*B矩阵右除：A\B %相当于A*inv(B)矩阵的幂：^矩阵转置：'矩阵求共轭(实部相同，虚部相反)：conj(X)矩阵求逆：inv(X)级数的求和与收敛symsum(fun,var,a,b):其中fun是通项表达式，var为求和变量，a为求和起点，b为求和终点例如：I为1/[n*(2n+1)]从1到正⽆穷的和，求Isyms n;f1=1/(n*(2*n+1));I=symsum(f1,n,1,inf)计算结果为：I =2-2*log(2)空间曲⾯mesh()函数语法：mesh(Z):mesh(X,Y,Z,C):其中C是⽤来定义相应点颜⾊等属性的数组例：求x^2+y^2=z的空间曲⾯x=-4:4;y=x;[X,Y]=meshgrid(x,y);%⽣成x,y坐标Z=X.^2+Y.^2;mesh(X,Y,Z)曲⾯图[x,y]=meshgrid(xa,ya) 当xa,ya分别为m维和n维⾏向量，得到x和y均为n⾏m列矩阵。
Matlab中的图像融合和多模态图像分析技术
Matlab中的图像融合和多模态图像分析技术图像处理是一项非常重要的技术，在许多领域都有广泛的应用，如医学影像分析、计算机视觉、遥感图像处理等。
在图像处理中，图像融合和多模态图像分析技术是两个非常重要的方面。
本文将介绍在Matlab中实现图像融合和多模态图像分析的方法和技术。
一、图像融合技术图像融合是指将多个不同模态或不同源的图像融合为一个具有更丰富信息的图像。
在图像融合技术中，常用的方法有像素级融合和特征级融合。
1.1 像素级融合像素级融合是指将多幅图像的像素按照一定的规则进行融合。
在Matlab中，可以使用imfuse函数来实现像素级融合。
该函数可以通过设置不同的融合模式来实现不同的效果，如加权平均、最大值、最小值等。
通过调整各个模态的权重，可以获得不同的融合效果。
1.2 特征级融合特征级融合是指将多幅图像的特征进行融合。
在Matlab中，可以使用特征提取和特征匹配的方法来实现特征级融合。
首先，使用不同的特征提取方法，如SIFT、SURF等，提取多幅图像的特征点。
然后，使用特征匹配的方法，如RANSAC算法，将多幅图像的特征点进行匹配和融合。
最后，根据匹配结果，可以生成一幅具有更丰富信息的图像。
二、多模态图像分析技术多模态图像分析是指对多模态图像进行分析和处理，以获得更全面和准确的信息。
在Matlab中，可以使用多种方法和技术来实现多模态图像分析。
2.1 图像配准图像配准是多模态图像分析的基础，它是将多幅图像进行准确的空间或特征对齐。
在Matlab中，可以使用imregister函数来实现图像配准。
该函数可以通过设置不同的配准方法和参数，如相位相关、归一化互相关等，来实现不同的配准效果。
2.2 图像分割图像分割是将图像中的目标或区域进行划分和提取的过程。
在多模态图像分析中，图像分割可以用来提取不同模态之间的特征。
在Matlab中，可以使用多种图像分割算法，如阈值分割、区域生长、边缘检测等，来实现图像分割。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
ＸＵ Ｓ ＨＥＮＧＸＩＡＮＧ ＨＵ ＣＨＡＯ ＸＵ ＹＵＮＱＩＮＧ
摘要：为 了 更 好 地 进 行 不 同 分 辨 率 图 像 的 融 合 ， 提 出 了 一 种 在 Ｍａｔｌａｂ 平 台 上 实 现 基 于 小 波 变 换 的 多 光 谱 图 像 与 高 分 辨 率 图 像融合方法。实验结果分析表明， 得到的融合图像与原多光谱图像相比， 细节信息更为突出， 整体信息更为丰富， 基本达到 了提高融合增强的目的。 关键词：小波变换； 图像融合； 遥感 中图分类号：ＴＰ７５１ 文献标识码：Ｂ
ｆｅｃｔｓ． Ｉｔ ｎｏｔ ｏｎｌｙ ｐｒｅｓｅｒｖｅｓ ｓｐｅｃｔｒａｌ ｉｎｆｏｒｍａｔｉｏｎ ｏｆ ｔｈｅ ｏｒｉｇｉｎａｌ ｍｕｌｔｉ－ ｓｐｅｃｔｒａｌ ｉｍａｇｅ ｗｅｌｌ， ｂｕｔ ａｌｓｏ ｅｎｈａｎｃｅｓ ｓｐａｔｉａｌ ｄｅｔａｉｌ ｉｎｆｏｒｍａｔｉｏｎ．
技 Ｋｅｙ Ｗｏｒ ｄｓ：ｗａｖｅｌｅｔ ｔｒ ａｎｓｆｏｒ ｍ， ｉｍａｇｅ ｆｕｓｉｏｎ， ｒ ｅｍｏｔｅ ｓｅｎｓｉｎｇ
表 １ 图像处理及小波变换函数
Ｍａｔｌａｂ 软件提供的工具箱 （如图像处理 工 具 、小波 工 具 箱 等）都是由该领域内的学术水平很高的专家编写的， 用户无需编 写 自 己 学 科 范 围 内 的 基 础 程 序 ， 而 直 接 进 行 高 、精 、尖 的 研 究 。 它提供了一个非常方便快捷的算法研究平台， 让用户把精力集 中 在算 法 而 不是 编 程 上， 极 大 地 方便 了 问 题 的 研 究 。Ｍａｔｌａｂ 工 具箱提供了实现图像处理和小波变换的一些函数， 下面仅介绍
令， 其分析处理结果便以数值或图形方式显示出来。
本文在分析遥感图像融合方法和小波变换的基础上， 提出
了一种基于小波变换的遥感图像融合方 法 ， 并 借助 Ｍａｔｌａｂ 图 像
处理和小波分析工具箱进行了仿真， 同时依据一些图像融合的
评价标准对融合结果图像进行了评价。
１ 图像处理和小波变换工具箱实现
遥感图像融合
术
近年来， 小波变换在遥感图像融合中的应用已有诸多报
导， 但 都 是 采用 ＦＯＲＴＲＡＮ 语 言 和 Ｃ 等 高 级 语 言 编 程 实 现 的 。
创 小波变换的数学公式很多， 对数学知识要求很严， 要自己编程
实现 并 不 容易 。Ｍａｔｌａｂ 的 问 世则 简 化 了这 部 分 工 作 ， 它 由 主 程
新 序及功能丰富的工具箱组成。其中小波工具箱提供了大量功能
完 善 的 小 波 工 具（包 括 各 种 小 波 分 析 函 数 等）， 是 进 行 图 像 处 理
的理想工具。在 Ｍａｔｌａｂ 环境下， 对图像的 分 析和 处 理 可采 用 人
机 交 互 的 方 式 ， 用 户 只 需 按 Ｍａｔｌａｂ 的 格 式 要 求 给 出 相 应 的 命
创 征 信 息 ， 根 据 不 同 融 合 目 的 ， 分 别 针 对 图 像 的 高 频 域 和 低 频 域
从上述系数融合方法中选取合适的融合准则进行处理， 以达到
更好的融合效果， 从整体上提高融合质量。
４ 结束语
文中将小波变换应用于遥感图像融合处理， 取得了良好的
效果， 不仅保留了多光谱图像的光谱信息， 而且有效地提高了
融 合 图 像 的 空 间 信 息 量 ， 为 遥 感 图 像 提 高 清 晰 度 、可 靠 性 和 解
译度提供了有力的支持。另外 ， 在 Ｍａｔｌａｂ 平台 上 非 常方 便 快 捷
Ａｂｓｔｒ ａｃｔ：Ｉｎ ｔｈｉｓ ｐａｐｅｒ， ａ ｍｅｔｈｏｄ ｏｆ ｉｍａｇｅ ｆｕｓｉｏｎ ｍｅｔｈｏｄ ｗａｓ ｐｒｏｐｏｓｅｄ ｏｎ Ｍａｔｌａｂ ｐｌａｔｆｏｒｍ ｂａｓｅｄ ｏｎ ｗａｖｅｌｅｔ ｔｒａｎｓｆｏｒｍ ｔｏ ｆｕｓｅ ａ ｌｏｗ－
ｒｅｓｏｌｕｔｉｏｎ ｍｕｌｔｉ － ｓｐｅｃｔｒａｌ ｉｍａｇｅ ａｎｄ ｈｉｇｈ － ｒｅｓｏｌｕｔｉｏｎ ｐａｎｃｈｒｏｍａｔｉｃ ｉｍａｇｅ． Ｅｘｐｅｒｉｍｅｎｔａｌ ｒｅｓｕｌｔｓ ｓｈｏｗ ｔｈａｔ ｔｈｉｓ ｓｃｈｅｍｅ ｈａｓ ｂｅｔｔｅｒ ｅｆ－
合图像＇ ）；
２ 实验结果及评价
由于没有标准图像可以比较， 对融合结果的评价一直是困 扰人们的一个问题。考虑到图像大小的一致性， 本文一方面选 择 ＳＰＯＴ 全 色波 段 图 像和 以 之 配准 过 的 ＴＭ 多波 段 图 像作 为 结 果判定的参考图像； 另一方面， 从数理统计角度和视觉量测方 面进行分析和比较， 以评价小波变换方法的融合性能。图 １ 中 （ａ）和（ｂ）为融合前遥感图像， （ｃ）为小波变换融合结果图像。
（２）用 ＴＭ 的低频图像替代 ＳＰＯＴ 全色图像的低频图像； （３）用替 换 后的 ＴＭ 低 频 图像 与 ＳＰＯＴ 全 色图 像 的 细节 图 像 进行小波逆变换得到融合结果图像。 １．２ Ｍａｔｌａｂ 实现 为了验证小波变换图像融合方法的正确性和有效性， 作者 选取低分辨率的 ＴＭ 多光 谱 图 像 ｄｍｔｍ．ｉｍｇ 和 高 分 辨率 的 ＳＰＯＴ 全 色 图 像 ｓｐｏｔｓ．ｉｍｇ 在 Ｍａｔｌａｂ 平 台 上 进 行 融 合 实 验 。 两 图 像 均 已经过几何校正， 并已精确配准， 选取的实验区大小 为 ４００×４００ 像素。Ｍａｔｌａｂ 中实现小波变换融合的主要程序如下： ％读取并显示原始多光谱和全色图像数据 ｉｍｇＰａｎ ＝ ｉｍｒｅａｄ（＇ ｓｐｏｔｓ．ｂｍｐ＇ ）； ｉｍｇＭｕｌ ＝ ｉｍｒｅａｄ（＇ ｄｍｔｍ．ｂｍｐ＇ ）； ｓｕｂｐｌｏｔ（１，３，１）， ｉｍｓｈｏｗ（ｉｍｇＭｕｌ）， ｘｌａｂｅｌ （＇ （ａ）低 分 辨 率 多 光
表 ２ 几种融合结果图像统计参数对比
从表 ２ 的 实验 数 据 可以 看 出 ， 融合 图 像 的标 准 差 、熵 、清 晰 度、空间频率等四项指标均比原始多光谱图像有了较大的提 高， 这表明融合图像的信息量比原始图像有明显增加， 而且图 像 的 细 节 反 差 、纹 理 和 清 晰 度 得 到 较 大 的 提 高 ， 融 合 图 像 质 量 较 好 。上 述 指 标 的 统 计 参 数 的 变 化 趋 势 与 融 合 方 法 的 理 论 分 析 结果一致， 与人的视觉效果的变化一致， 充分说明了小波变换 融合方法是可行的。
地完成了相关融合算法及效果评价， 大大提高了工作效率， 方
便了问题的研究。
技
然而遥感图像融合， 因涉及到 同 ， 很难 有 一 个具 体 模 型的 情 况 ， 就
要 用 基 于 弹 性 模 型 的 方 法 。基 于 图 像 本 身 所 存 在 的 相 关 性 和 特
图 １ 源图像及小波变换融合结果图像 从主观目视判读可以看出， 融合结果图像的空间分辨率 比原多光谱图像有较大程度的提高， 有效地增强了图像的空 间 细 节 信 息 量 ， 融 合 图 像 中 道 路 、山 脉 的 纹 理 清 楚 ， 边 缘 轮 廓 也清晰可见。 为了对小波变换法从量上进行评价， 这里采用了保留光谱 信 息 和 空 间 细 节 信 息 两 个 方 面 的 评 价 指 标 ：图 像 灰 度 均 值 、标 准 差 、信 息 熵 、交 叉 熵 、清 晰 度 、扭 曲 程 度 、偏 差 指 数 、空 间 频 率 和 峰值信噪比。借助 Ｍａｔｌａｂ 的图像处理工具 箱 中 的相 关 函 数， 上 述指标均可容易得到。为简便起见， 取融合前后图像的 ＲＧＢ 三 个光谱分量统计得到的评价指标数据的平均值来比较， 具体统 计数据如表 ２ 所示。
图像处理 文章编号：１００８－ ０５７０（２００８）０１－ ３－ ０３０２－ ０２
中 文 核 心 期 刊 《 微 计 算 机 信 息 》（ 管 控 一 体 化 ）２００８ 年 第 ２４ 卷 第 １－３ 期
利用 Ｍａ ｔｌａ ｂ 实现基于小波变换的遥感图像融合
Ｒｅ ｍ ｏ ｔｅ Ｓ ｅ ｎ ｓ ｉｎ ｇ Ｉｍ ａ ｇ ｅ Ｆｕ ｓ ｉｏ ｎ Ｂａ ｓ ｅ ｄ ｏ ｎ Ｗａ ｖｅ ｌｅ ｔ Ｔｒａ ｎ ｓ ｆｏ ｒｍ Ａｃｈ ｉｅ ｖｅ ｄ ｂ ｙ Ｍａ ｔｌａ ｂ
（１）对配 准 后 ＴＭ 和 ＳＰＯＴ 全 色 图像 分 别 进行 ｎ 次 小 波 变 换 （ｎ 通 常取 ２ 或 ３）， 得 到各 自 相 应分 辨 率 的低 频 图 像和 高 频 细 节 纹 理 图像 ， 本 文 中选 用 的 小波 基 是 长度 为 １３ 的 Ｄａｕｂｅｃｈｉｅｓ 小 波基， 小波变换级数 ｎ＝３；
徐胜祥： 硕士 讲师 基 金 项 目 ： 孝 感 学 院 自 然 科 学 立 项 项 目 资 助 （ｚ２００７０１０）
程序中涉及的几个主要的函数， 见表 １ 所示。 １．１ 基于小波变换的图像融合方法 小波分析是近年来图像分析的热点， 有关小波应用于融合
的研究已有不少报道， 其中比较经典的基于小波变换的融合方 法是以 Ｍａｌｌａｔ 算法为核心的离散正交小波变换法。对于二维的 图像信号来说， 经过一次离散正交小波变换后， 图像被分解为 ４ 幅， 其中左 上 角 一幅 是 原 图像 的 平 滑逼 近（低 频）， 左 下 角 为垂 直 细 节 ， 右 上 角 为 水 平 细 节 ， 右 下 角 为 原 图 像 的 细 节 分 量（高 频）。 然后， 可以继续对左上角的低频分量进行同样的第 ２ 次、第 ３ 次……小波变换。利用其分频特征和重构特性， 可以用于以非 线性的对数映射方式融合不同类型的图像数据， 使融合后的图 像既保留原高分辨率遥感图像的结构信息， 又融合多光谱图像 丰富的光谱信息， 提高遥感图像的解译能力、分类精度。采用离 散正 交 小波 变 换 进行 ＴＭ 和 ＳＰＯＴ 两 图像 的 数 据融 合 ， 具 体 步 骤如下：
－ ３０２ － ３６０元 ／ 年 邮局订阅号： ８２－９４６
《现场总线技术应用 ２００ 例》
您的论文得到两院院士关注
图像处理
谱图像＇ ）； ｓｕｂｐｌｏｔ（１，３，２）， ｉｍｓｈｏｗ（ｉｍｇＰａｎ）， ｘｌａｂｅｌ （＇ （ｂ）高 分 辨 率 全 色