【图像融合】基于matlab加权图像融合【含Matlab源码 722期】

   日期:2021-04-08     浏览:7    评论:0    
核心提示:一、简介基于matlab加权图像融合二、源代码clearg_R=0;g_G=0;g_B=0;h_R=0;h_G=0;h_B=0;fenzi_R=0;fenzi_G=0;fenzi_B=0;fenmu_up_R=0;fenmu_up_G=0;fenmu_up_B=0;fenmu_low_R=0;fenmu_low_G=0;fenmu_low_B=0;tableR=[];tableG=[];tableB=[];up=imread('high.jpg');

一、简介

基于matlab加权图像融合

二、源代码

clear
g_R=0;
g_G=0;
g_B=0;
h_R=0;
h_G=0;
h_B=0;
fenzi_R=0;
fenzi_G=0;
fenzi_B=0;
fenmu_up_R=0;
fenmu_up_G=0;
fenmu_up_B=0;
fenmu_low_R=0;
fenmu_low_G=0;
fenmu_low_B=0;
tableR=[];
tableG=[];
tableB=[];
up=imread('high.jpg');         %读图像
low=imread('low.jpg');

figure(1)
imshow(up);                                     %读RGB数值

[M,N,color]=size(up);

title('加权-RGB表示的高分辨率图像');

figure(2)
imshow(low); 
title('加权-RGB表示的低分辨率图像');
r=double(up(:,:,1));
g=double(up(:,:,2));
b=double(up(:,:,3));
r_low=double(low(:,:,1));
g_low=double(low(:,:,2));
b_low=double(low(:,:,3));
RGB(:,:,1)=0.5*r+0.5*r_low;
RGB(:,:,2)=0.5*g+0.5*g_low;
RGB(:,:,3)=0.5*b+0.5*b_low;
R=RGB(:,:,1);
G=RGB(:,:,2);
B=RGB(:,:,3);
RGB=uint8(round(RGB));   
figure(3)
imshow(RGB)
title('加权-RGB转化后的图像');




              %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
              %                       下面是计算平均梯度G                          %
              %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%     


for ii=1:M-1
    for jj=1:N-1
        g_R=g_R+sqrt((((r(ii+1,jj)-r(ii,jj))^2+(r(ii,jj+1)-r(ii,jj))^2))/2);
        g_G=g_G+sqrt((((g(ii+1,jj)-g(ii,jj))^2+(g(ii,jj+1)-g(ii,jj))^2))/2);
        g_B=g_B+sqrt((((b(ii+1,jj)-b(ii,jj))^2+(b(ii,jj+1)-b(ii,jj))^2))/2);

    end
end
fprintf('\n\n highR的清晰度为:%.4f\n highG的清晰度为:%.4f\n highG的清晰度为:%.4f',...
            g_R/(M-1)/(N-1),g_G/(M-1)/(N-1),g_B/(M-1)/(N-1));              
              
g_R=0;
g_G=0;
g_B=0;
              

for ii=1:M-1
    for jj=1:N-1
        g_R=g_R+sqrt((((r_low(ii+1,jj)-r_low(ii,jj))^2+(r_low(ii,jj+1)-r_low(ii,jj))^2))/2);
        g_G=g_G+sqrt((((g_low(ii+1,jj)-g_low(ii,jj))^2+(g_low(ii,jj+1)-g_low(ii,jj))^2))/2);
        g_B=g_B+sqrt((((b_low(ii+1,jj)-b_low(ii,jj))^2+(b_low(ii,jj+1)-b_low(ii,jj))^2))/2);

    end
end
fprintf('\n\n lowR的清晰度为:%.4f\n lowG的清晰度为:%.4f\n lowG的清晰度为:%.4f',...
         g_R/(M-1)/(N-1),g_G/(M-1)/(N-1),g_B/(M-1)/(N-1));              
              
g_R=0;
g_G=0;
g_B=0;
      

for ii=1:M-1
    for jj=1:N-1
        g_R=g_R+sqrt((((R(ii+1,jj)-R(ii,jj))^2+(R(ii,jj+1)-R(ii,jj))^2))/2);
        g_G=g_G+sqrt((((G(ii+1,jj)-G(ii,jj))^2+(G(ii,jj+1)-G(ii,jj))^2))/2);
        g_B=g_B+sqrt((((B(ii+1,jj)-B(ii,jj))^2+(B(ii,jj+1)-B(ii,jj))^2))/2);

    end

三、运行结果



四、备注

完整代码或者代写添加QQ 1564658423

 
打赏
 本文转载自:网络 
所有权利归属于原作者,如文章来源标示错误或侵犯了您的权利请联系微信13520258486
更多>最近资讯中心
更多>最新资讯中心
更多>相关资讯中心
0相关评论

推荐图文
推荐资讯中心
点击排行
最新信息
新手指南
采购商服务
供应商服务
交易安全
关注我们
手机网站:
新浪微博:
微信关注:

13520258486

周一至周五 9:00-18:00
(其他时间联系在线客服)

24小时在线客服