找回密码
 注册
关于网站域名变更的通知
查看: 630|回复: 1
打印 上一主题 下一主题

基于matlab自适应形态学之图像去噪

[复制链接]
  • TA的每日心情

    2019-11-20 15:22
  • 签到天数: 2 天

    [LV.1]初来乍到

    跳转到指定楼层
    1#
    发表于 2021-4-29 17:53 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式

    EDA365欢迎您登录!

    您需要 登录 才可以下载或查看,没有帐号?注册

    x
    9 ]- k0 m1 h; J( Y2 S. O3 V/ a- A
    一、简介
    # |/ k' Y# {% u; t5 a7 }. o数字图像在进行数学形态滤波去噪时,根据噪声特点可以尝试采用维数由小到大的结构元素来进行处理,进而达到滤除不同噪声的目的。采用数学形态学的多结构元素,可以更多地保持数字图像的几何特征.  s9 H8 g# ]% F" A- ]9 }
    ( c2 Z; j* n/ c2 @+ T
    ! a% c5 A" m) N' w& {
    二、源代码
    " _. q  f& s- x/ P) Q5 H7 C9 j- c7 n$ B
    • %% Figure 1
    • A = double(imread('retina.png'));
    • A = A(:,:,2);
    • x1=361;y1=493;
    • x2=289;y2=207;
    • Z1 = GAN(A,18,y1,x1,'CLIP');
    • Z2 = GAN(A,18,y2,x2,'CLIP');
    • Z = Z1+Z2;
    • figure
    • subplot(121);
    • hold on;
    • imshow(A,[]);
    • plot(x1,y1,'w*');title('points (x_1,y_1) and (x_2,y_2)')
    • plot(x2,y2,'w*');
    • subplot(122);
    • imshow(max(A,255*Z),[]);title('GANs of (x_1,y_1) and (x_2,y_2)')
    • %% Figure 2
    • A = double(imread('lena256.jpg'));
    • se = strel('disk',2);
    • A_dilation = imdilate(A,se);
    • A_erosion = imerode(A,se);
    • mtol = 20;
    • A_GANdilation = GANdilation(A,A,mtol,'CLIP');
    • A_GANerosion = GANerosion(A,A,mtol,'CLIP');
    • figure
    • subplot(231);imshow(A,[]);title('original');
    • subplot(232);imshow(A_dilation,[]);title('classical dilation');
    • subplot(233);imshow(A_erosion,[]);title('classical erosion');
    • subplot(235);imshow(A_GANdilation,[]);title('GAN dilation');
    • subplot(236);imshow(A_GANerosion,[]);title('GAN erosion');
    • %% Figure 3
    • A = double(imread('vaisseaux3.png'));
    • C = floor(imContrasteLIP(A));
    • B_CLIP = toggleContrastCLIP(A,5);
    • B_LIP = GANtoggleContrastLIP(A,255-C,30);
    • figure
    • subplot(221);imshow(A,[]);title('original');
    • subplot(222);imshow(C,[]);title('LIP contrast');
    • subplot(223);imshow(B_CLIP,[]);title('classical toggle contrast');
    • subplot(224);imshow(B_LIP,[]);title('GAN toggle contrast');
    • %% Figure 4
    • A = double(imread('lena256.jpg'));
    • Anoisy = double(uint8(255*imnoise(A/255,'salt & pepper',0.02)));
    • fun = @(x) median(x(:));
    • B = nlfilter(Anoisy,[3 3],fun);
    • C = GANmedian(Anoisy,Anoisy,20,'CLIP');
    • D = ZGANmedian(Anoisy,Anoisy,20,'CLIP',4);
    • figure
    • subplot(231);imshow(A,[]);title('noise free'); title('original');
    • subplot(232);imshow(Anoisy,[]);title('noisy');title('noisy');
    • subplot(234);imshow(B,[]);title('classical median');title('classical median');
    • subplot(235);imshow(C,[]);title('GAN median');title('GAN median');
    • subplot(236);imshow(D,[]);title('Z-GAN median');title('Z-GAN median');

    • * v8 N# Z9 t& i0 v! t, z, q
       
    4 O7 y9 U/ m2 W% v  a/ I1 Z9 l: Z0 J" x  R* B' \
    0 l* ]5 Y: }, o0 p8 R1 g) v
    三、运行结果6 X( `- y- r3 I: s5 }8 p
    1 U" T0 T- H9 j2 @* q# Y1 o

    9 z9 h& l% b, B0 x& r# O9 ?, o
    5 B1 M5 `& x& I) ]
  • TA的每日心情
    慵懒
    2020-8-28 15:16
  • 签到天数: 1 天

    [LV.1]初来乍到

    2#
    发表于 2021-4-30 09:49 | 只看该作者
    采用数学形态学的多结构元素
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

    推荐内容上一条 /1 下一条

    EDA365公众号

    关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

    GMT+8, 2025-10-29 23:14 , Processed in 0.125000 second(s), 26 queries , Gzip On.

    深圳市墨知创新科技有限公司

    地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

    快速回复 返回顶部 返回列表