【数字图像处理matlab】(HSI变换融合算法）

输入一张高分辨率的全色影像HR，一张低分辨率的多光谱影像MS，采用HSI变换融合算法实现影像融合，其中RGB与HSI影像的相互转换调用自定义函数RGB2HSI()，HSI2RGB()，函数具体代码为。

RGB2HSI()

function Image1=RGB2HSI(image)
%调用代码-------------------------------------------
%imageR=imread('animal.jpg');
%RGB2HSI(image);
%--------------------------------------------------image=im2double(image);
%从三维数组中提取三幅分量图像
R1=image(:,:,1);
G1=image(:,:,2);
B1=image(:,:,3);
I=(R1+G1+B1)/3;        %亮度分量，范围[0,1]
m=min(min(R1,G1),B1);
S=1-3*m./(R1+G1+B1);    %饱和度分量，范围[0,1]
theta=acos(((R1-G1)+(R1-B1))./(2*((R1-G1).^2+((R1-B1).*(G1-B1))).^(1/2))); %弧度
H=theta;   %色度分量，以角度表示，范围是[0,1]（弧度除以2*pi后）
if B1>G1H=2*pi-theta;
end
if S==0H=0;
end
H=H/(2*pi);
Image1=cat(3,H,S,I);
subplot(1,2,1),imshow(image);
title('RGB原图')
subplot(1,2,2),imshow(Image1);
title('转换后的HSI图像')
end

HSI2RGB()

function rgb= HSI2RGB(hsi)
H = hsi(:, :, 1) * 2 * pi;
S = hsi(:, :, 2);
I = hsi(:, :, 3); R = zeros(size(hsi, 1), size(hsi, 2));
G = zeros(size(hsi, 1), size(hsi, 2));
B = zeros(size(hsi, 1), size(hsi, 2)); idx = find( (0 <= H) & (H < 2*pi/3));
B(idx) = I(idx) .* (1 - S(idx));
R(idx) = I(idx) .* (1 + S(idx) .* cos(H(idx)) ./ cos(pi/3 - H(idx)));
G(idx) = 3*I(idx) - (R(idx) + B(idx)); idx = find( (2*pi/3 <= H) & (H < 4*pi/3) );
R(idx) = I(idx) .* (1 - S(idx));
G(idx) = I(idx) .* (1 + S(idx) .* cos(H(idx) - 2*pi/3) ./ cos(pi - H(idx)));
B(idx) = 3*I(idx) - (R(idx) + G(idx)); idx = find( (4*pi/3 <= H) & (H <= 2*pi));
G(idx) = I(idx) .* (1 - S(idx));
B(idx) = I(idx) .* (1 + S(idx) .* cos(H(idx) - 4*pi/3) ./cos(5*pi/3 - H(idx)));
R(idx) = 3*I(idx) - (G(idx) + B(idx)); rgb = cat(3, R, G, B);
rgb = max(min(rgb, 1), 0);
figure,
%subplot(1,2,1),
imshow(rgb);title('转换回的RGB图像')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
H1 = hsi(:, :, 1) * 2 * pi;
S1 = hsi(:, :, 2);
I1 = hsi(:, :, 3);
R1 = zeros(size(hsi, 1), size(hsi, 2));
G1= zeros(size(hsi, 1), size(hsi, 2));
B1 = zeros(size(hsi, 1), size(hsi, 2));
if 0<=H1<2*pi/3;B1=I1.*(1-S1);R1=I1.*(1+(S1.*cos(H1)./cos(pi/3-H1)));G1=3.*I1-(B1+R1);
elseif((2*pi/3<=H1)&&(H1<4*pi/3));R1=I1.*(1-S1);G1=I1.*(1+(S1.*cos(H1-2*pi/3)./cos(pi-H1)));B1=3.*I1-(R1+G1);
else((4*pi/3<=H1)&(H1<2*pi));G1=I1.*(1-S1);B1=I1.*(1+(S1.*cos(H1-4*pi/3)./cos(10*pi/6-H1)));R1=3.*I1-(G1+B1);
end
r=cat(3,R1,G1,B1);
r=max(min(r,1),0);
subplot(1,2,2),imshow(r);
title('转换回的RGB图像') end

基于HSI变换的融合算法：

function F=HSI_melt(HR,MS)
%调用代码------------------------------------------------------
%Image1=imread('HR.jpg');Image2=imread('MS.jpg');
%HSI_melt(Image1,Image2);
%--------------------------------------------------------------
HR=im2double(HR); [row,column,~]=size(HR);
MS=im2double(MS);
MS1=RGB2HSI(MS);     %多光谱影像转换至HSI空间
H=MS1(:,:,1); S=MS1(:,:,2); I=MS1(:,:,3);
MAX2=max(max(I));MIN2=min(min(I));
MAX1=max(max(HR));MIN1=min(min(HR));
%高分辨率影像拉与亮度分量一致,假设变换前[MIN1,MAX2]后的亮度值为[MIN2,MAX2]
A=(MAX2-MIN2)/(MAX1-MIN1);B=(MAX1*MIN2-MIN1*MAX2)/(MAX1-MIN1);
HR1=zeros(row,column);
for i=1:rowfor j=1:columnHR1(i,j)=HR(i,j)*A+B;end
endMS2=cat(3,H,S,HR1);
Image=HSI2RGB(MS2);
figure,imshow(Image);
title('基于HSI变换的影像融合结果');
end

测试：

【数字图像处理matlab】(HSI变换融合算法）相关推荐

数字图像处理MATLAB学习笔记（五）
数字图像处理MATLAB学习笔记(五) Color Image Processing 1 Color Image Representation in MATLAB 这里不多说了,彩色图片在计算机中以R ...
加权平均法融合图像matlab,基于MATLAB的图像融合算法
内容介绍原文档由会员小甜甜发布基于MATLAB的图像融合算法 1.9万字 34页摘要图像融合能够将不同类型传感器获取的同一对象的图像数据进行空间配准.并且采用一定的算法将各图像数据所含的信 ...
《数字图像处理 MATLAB版》学习笔记
学习教材:<数字图像处理 MATLAB版>(第二版) 冈萨雷斯学习过程中的图片代码和及我收集的一些关于数字图像处理的其他学习资料,需要的可以评论留下邮箱(需要购买专栏),加油文章目录 ...
数字图像处理MATLAB学习笔记（一）
数字图像处理MATLAB学习笔记(一) 灰度转换与空间滤波本节主要使用Matlab语言进行灰度转换与空间滤波的使用并对相关数学原理进行总结 1. Intensity Transformer Fun ...
【资源分享】数字图像处理MATLAB版冈萨雷斯+中文高清版+随书源码链接
写在这里的初衷,一是备忘,二是希望得到高人指点,三是希望能遇到志同道合的朋友. 目录 1.数字图像处理MATLAB版冈萨雷斯+中文高清版 2.数字图像处理MATLAB版冈萨雷斯随书源码 1.数字图像处 ...
Win8Metro(C#)数字图像处理--2.31灰度拉伸算法
Win8Metro(C#)数字图像处理--2.31灰度拉伸算法原文:Win8Metro(C#)数字图像处理--2.31灰度拉伸算法 [函数名称] 灰度拉伸函数GrayStretchProces ...
数字图像处理--傅里叶(逆)变换
数字图像处理–傅里叶(逆)变换主要内容 (1)对一副图像进行缩放,显示原始图像和缩放后的图像,分别对其进行傅里叶变换,显示变换后结果: (2)对一副图像进行旋转,显示原始图像和旋转后的图像,分别对其 ...
[数字图像处理Matlab]任选一幅灰度图，自行编写程序，完成直方图均衡化。鼠鼠的数字图像处理实验要求：不能用MATLAB自带的histeq函数实现直方图均衡化
数学理论前提:见数字图像处理(第三版) 李俊山等编著 p50-53 1.利用imhist函数统计像素点个数 Im = imread("C:\Users\鼠鼠\Desktop\数字图像处理m ...
图像处理----入门资料，Matlab r2019最新版，r2008a，《数字图像处理》冈萨雷斯第三版《数字图像处理 Matlab 版》
研究生想研究的方向是图像处理,现在离开学还有四个月,希望能稍微入门吧.同时也希望通过博客的形式记录一下学习的进度,学习的心情.目前学习两周,安装了基本的软件,收集书籍,学习了灰度变换,空间滤波,频率滤 ...

【数字图像处理matlab】(HSI变换融合算法）

【数字图像处理matlab】(HSI变换融合算法）

【数字图像处理matlab】(HSI变换融合算法）相关推荐

最新文章

热门文章