• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    迪恩网络公众号

一个matlab数字图像处理程序的解释

原作者: [db:作者] 来自: [db:来源] 收藏 邀请
clc;                                                   %clc是清除command window里的内容
clear all;                                             %clear是清除workspace里的变量
close all;                                             %close all来关闭所有已经打开的图像窗口
image_ori = imread(\'skeleton_orig.bmp\');               %读取图像数据
image_ori = rgb2gray(image_ori);                       % 将rgb模式转换成灰度图
figure;                                                %画图
subplot(2,4,1);                                        %把区域分成2行4列,并把图像显示到第一个
imshow(image_ori);                                     %显示图像
title(\'a.original image\');                             %显示标题
image_ext = [zeros(1,500);image_ori;zeros(1,500)]\';    %分别在上下两边加一行,并转置
image_ext = [zeros(1,802);image_ext;zeros(1,802);]\';   %分别在左右两边加一列,并转置
image_double = im2double(image_ext);                   %将imdata_add转换成双精度

% lapa_modelace Operating 
lapa_mode = [-1,-1,-1;-1,8,-1;-1,-1,-1];                %拉普拉斯变换用模板[-1,-1,-1;-1,8,-1;-1,-1,-1]; 
lapa_mode = lapa_mode(:);                               %转换成列矩阵
for x = 2:1:801
    for y = 2:1:501
    A =image_double([x-1:x+1],[y-1:y+1]);               %从image_double中以(x,y)为中心取出一个3*3的矩阵 
    image_lapa(x-1,y-1) = lapa_mode\'*A(:);              %将取出的矩阵与模板进行拉普拉斯变换取得的值放在image_lapa的(x-1,y-1)位置
    end
end
subplot(2,4,2);
imshow(image_lapa);
title(\'b.lapalace image\');
subplot(2,4,3);

%Sharpening Operating
image_ab_sha = image_lapa+im2double(image_ori);           
imshow(image_ab_sha);
title(\'c.a+b sharpening image\');

%Grads Operating
sobel_mode_1= [-1,-2,-1;0,0,0;1,2,1];                     %定义梯度处理的一个掩膜
sobel_mode_1= sobel_mode_1(:);
sobel_mode_2 = [-1,0,1;-2,0,2;-1,0,1];                    %定义梯度处理的另一个掩膜
sobel_mode_2 = sobel_mode_2(:);
for x = 2:1:801
    for y = 2:1:501
    A =image_double([x-1:x+1],[y-1:y+1]); 
    mid_1(x-1,y-1) = sobel_mode_1\'*A(:); 
    mid_2(x-1,y-1) = sobel_mode_2\'*A(:);
    end
end
image_grad = abs(mid_1)+abs(mid_2);                        %对gx、gy绝对值化,然后求和
subplot(2,4,4);
imshow(image_grad);
title(\'d.grads image\');

%Smoothing Operating
Mxy_Ext = [zeros(2,500);image_grad;zeros(2,500)];           %分别在上下两边加两行
Mxy_Ext = [zeros(2,804);Mxy_Ext\';zeros(2,804)]\';            %先将Mxy_Ext转置,然后分别在左右两边加一列,最后将所得结果再转置
Linear_smooth = ones(5,5);                                  %创建一个5*5矩阵,初始值全是1
Linear_smooth = Linear_smooth(:);                           %将矩阵Linear_smooth转换成列矩阵
for x = 3:1:802                                             %大循环,x从3--802,步进是1
    for y = 3:1:502                                         %小循环,y从3--502,步进是1
        A = Mxy_Ext([x-2:x+2],[y-2:y+2]);                   %以Emx_Ext(x,y)为中心,截取5*5矩阵复制到A中
        image_smooth(x-2,y-2) = Linear_smooth\'*A(:);         %将获取的矩阵值各个值求和然后放在新矩阵image_smooth的(x-2,y-2)位置
    end
end
image_smooth = image_smooth/25;                             %求平均值
subplot(2,4,5);                                             %把区域分成1行4列,并把图像显示到第一个
imshow(image_smooth);                                       %显示图像
title(\'e.smoothing image\');                                 %显示标题

%Masking Operating 
image_mask = image_ab_sha.*image_smooth;                    %点乘
subplot(2,4,6);
imshow(image_mask);
title(\'f.Masking image\');

%Sharpening Operating
image_af_sha = im2double(image_ori) + image_mask;
subplot(2,4,7);
imshow(image_af_sha);
title(\'g.a+f sharpening image\');

%Exponential  Transform
image_fin = image_af_sha.^0.5;                              %幂律变换 
subplot(2,4,8);
imshow(image_fin);
title(\'h.final result\');

 


鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
上一篇:
[转]常用Delphi开发资料网址发布时间:2022-07-18
下一篇:
比较R语言、perl语言,matlab中for循环和while循环的使用发布时间:2022-07-18
热门推荐
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap