数字图像处理--算术、几何、谐波、逆谐波均值滤波器Matlab

最新推荐文章于 2022-10-22 08:34:12 发布
最新推荐文章于 2022-10-22 08:34:12 发布
阅读量1.3w 收藏 103
点赞数 14
分类专栏： matlab算法 文章标签： matlab 图像处理
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Dooonald/article/details/78545461
版权

算术均值

close all 
clear all

f=imread('D:/testData/filtering.tif');

[w,h]=size(f);
image= f(:,:);
fsize1=3;
fsize2=5;
fsize3=9;

fssize1=(fsize1-1)/2;
fssize2=(fsize2-1)/2;
fssize3=(fsize3-1)/2;

figure();
subplot(1,2,1);
imshow(image);
xlabel('原图像');

resultImage = image;

for x=1+fssize1:1:w-fssize1
    for y=1+fssize1:1:w-fssize1
        is=f(x-fssize1:1:x+fssize1,y-fssize1:1:y+fssize1);
        resultImage(x,y)=sum(is(:))/numel(is); 
    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('3*3算术均值');

resultImage= f(:,:);
figure();
subplot(1,2,1);
imshow(f);
xlabel('原图像');

for x=1+fssize2:1:w-fssize2
    for y=1+fssize2:1:w-fssize2
        is=f(x-fssize2:1:x+fssize2,y-fssize2:1:y+fssize2);
        resultImage(x,y)=sum(is(:))/numel(is); 

    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('5*5算术均值');

resultImage= f(:,:);
figure();
subplot(1,2,1);
imshow(f);
xlabel('原图像');

for x=1+fssize3:1:w-fssize3
    for y=1+fssize3:1:w-fssize3
        is=f(x-fssize3:1:x+fssize3,y-fssize3:1:y+fssize3);
        resultImage(x,y)=sum(is(:))/numel(is); 

    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('9*9算术均值');

几何均值

close all 
clear all

f=imread('D:/testData/filtering.tif');

[w,h]=size(f);
image= f(:,:);
fsize1=3;
fsize2=5;
fsize3=9;

fssize1=(fsize1-1)/2;
fssize2=(fsize2-1)/2;
fssize3=(fsize3-1)/2;

figure();
subplot(1,2,1);
imshow(image);
xlabel('原图像');

resultImage = image;

for x=1+fssize1:1:w-fssize1
    for y=1+fssize1:1:w-fssize1
        is=f(x-fssize1:1:x+fssize1,y-fssize1:1:y+fssize1);
        resultImage(x,y)=prod(prod(is(:)))^(1/numel(is)); 
    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('3*3几何均值');

resultImage= f(:,:);
figure();
subplot(1,2,1);
imshow(f);
xlabel('原图像');

for x=1+fssize2:1:w-fssize2
    for y=1+fssize2:1:w-fssize2
        %遍历每个点的四周
        is=f(x-fssize2:1:x+fssize2,y-fssize2:1:y+fssize2);
        resultImage(x,y)=prod(prod(is(:)))^(1/numel(is)); 
    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('5*5几何均值');

resultImage= f(:,:);
figure();
subplot(1,2,1);
imshow(f);
xlabel('原图像');

for x=1+fssize3:1:w-fssize3
    for y=1+fssize3:1:w-fssize3
        %遍历每个点的四周
        is=f(x-fssize3:1:x+fssize3,y-fssize3:1:y+fssize3);
        resultImage(x,y)=prod(prod(is(:)))^(1/numel(is)); 
    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('9*9几何均值');

谐波均值

close all 
clear all

f=imread('D:/testData/filtering.tif');

[w,h]=size(f);
image= f(:,:);
fsize1=3;
fsize2=5;
fsize3=9;

fssize1=(fsize1-1)/2;
fssize2=(fsize2-1)/2;
fssize3=(fsize3-1)/2;

figure();
subplot(1,2,1);
imshow(image);
xlabel('原图像');

resultImage = image;

for x=1+fssize1:1:w-fssize1
    for y=1+fssize1:1:w-fssize1
        is=f(x-fssize1:1:x+fssize1,y-fssize1:1:y+fssize1);
        is=1./is;
        resultImage(x,y)=numel(is)/sum(is(:)); 

    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('3*3谐波均值');

resultImage= f(:,:);
figure();
subplot(1,2,1);
imshow(f);
xlabel('原图像');

for x=1+fssize2:1:w-fssize2
    for y=1+fssize2:1:w-fssize2
        %遍历每个点的四周
        is=f(x-fssize2:1:x+fssize2,y-fssize2:1:y+fssize2);
        is=1./is;
        resultImage(x,y)=numel(is)/sum(is(:));  
    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('5*5谐波均值');

resultImage= f(:,:);
figure();
subplot(1,2,1);
imshow(f);
xlabel('原图像');

for x=1+fssize3:1:w-fssize3
    for y=1+fssize3:1:w-fssize3
        %遍历每个点的四周
        is=f(x-fssize3:1:x+fssize3,y-fssize3:1:y+fssize3);
        is=1./is;
        resultImage(x,y)=numel(is)/sum(is(:));   
    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('9*9谐波均值');

逆谐波

close all 
clear all

f=imread('D:/testData/filtering.tif');

[w,h]=size(f);
image= f(:,:);
Q1 = 1.5;
fsize1=3;
fsize2=5;
fsize3=9;

fssize1=(fsize1-1)/2;
fssize2=(fsize2-1)/2;
fssize3=(fsize3-1)/2;

figure();
subplot(1,2,1);
imshow(image);
xlabel('原图像');

resultImage = image;

for x=1+fssize1:1:w-fssize1
    for y=1+fssize1:1:w-fssize1
        is=f(x-fssize1:1:x+fssize1,y-fssize1:1:y+fssize1);
        resultImage(x,y) = sum(is(:).^(Q1+1))/sum(is(:).^(Q1));
    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('3*3逆谐波均值');

resultImage= f(:,:);
figure();
subplot(1,2,1);
imshow(f);
xlabel('原图像');

for x=1+fssize2:1:w-fssize2
    for y=1+fssize2:1:w-fssize2
        %遍历每个点的四周
        is=f(x-fssize2:1:x+fssize2,y-fssize2:1:y+fssize2);
        resultImage(x,y) = sum(is(:).^(Q1+1))/sum(is(:).^(Q1));

    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('5*5逆谐波均值');

resultImage= f(:,:);
figure();
subplot(1,2,1);
imshow(f);
xlabel('原图像');

for x=1+fssize3:1:w-fssize3
    for y=1+fssize3:1:w-fssize3
        %遍历每个点的四周
        is=f(x-fssize3:1:x+fssize3,y-fssize3:1:y+fssize3);
        resultImage(x,y) = sum(is(:).^(Q1+1))/sum(is(:).^(Q1));
    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('9*9逆谐波均值');

中值

close all 
clear all

f=imread('D:/testData/filtering.tif');

[w,h]=size(f);
image= f(:,:);

fsize1=3;
fsize2=5;
fsize3=9;

fssize1=(fsize1-1)/2;
fssize2=(fsize2-1)/2;
fssize3=(fsize3-1)/2;

figure();
subplot(1,2,1);
imshow(image);
xlabel('原图像');

resultImage = image;

for x=1+fssize1:1:w-fssize1
    for y=1+fssize1:1:w-fssize1
        is=f(x-fssize1:1:x+fssize1,y-fssize1:1:y+fssize1);
        temp = sort(is(:));
        resultImage(x,y)= temp((numel(temp)-1)/2);
    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('3*3中值');

resultImage= f(:,:);
figure();
subplot(1,2,1);
imshow(f);
xlabel('原图像');

for x=1+fssize2:1:w-fssize2
    for y=1+fssize2:1:w-fssize2
        %遍历每个点的四周
        is=f(x-fssize2:1:x+fssize2,y-fssize2:1:y+fssize2);
        temp = sort(is(:));
        resultImage(x,y)= temp((numel(temp)-1)/2);
    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('5*5中值');

resultImage= f(:,:);
figure();
subplot(1,2,1);
imshow(f);
xlabel('原图像');

for x=1+fssize3:1:w-fssize3
    for y=1+fssize3:1:w-fssize3
        %遍历每个点的四周
        is=f(x-fssize3:1:x+fssize3,y-fssize3:1:y+fssize3);
        temp = sort(is(:));
        resultImage(x,y)= temp((numel(temp)-1)/2);
    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('9*9中值');

最大值

close all 
clear all

f=imread('D:/testData/filtering.tif');

[w,h]=size(f);
image= f(:,:);

fsize1=3;
fsize2=5;
fsize3=9;

fssize1=(fsize1-1)/2;
fssize2=(fsize2-1)/2;
fssize3=(fsize3-1)/2;

figure();
subplot(1,2,1);
imshow(image);
xlabel('原图像');

resultImage = image;

for x=1+fssize1:1:w-fssize1
    for y=1+fssize1:1:w-fssize1
        is=f(x-fssize1:1:x+fssize1,y-fssize1:1:y+fssize1);
        temp = is(:);
        resultImage(x,y)= max(temp);
    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('3*3最大值');

resultImage= f(:,:);
figure();
subplot(1,2,1);
imshow(f);
xlabel('原图像');

for x=1+fssize2:1:w-fssize2
    for y=1+fssize2:1:w-fssize2
        %遍历每个点的四周
        is=f(x-fssize2:1:x+fssize2,y-fssize2:1:y+fssize2);
        temp = is(:);
        resultImage(x,y)= max(temp);
    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('5*5最大值');

resultImage= f(:,:);
figure();
subplot(1,2,1);
imshow(f);
xlabel('原图像');

for x=1+fssize3:1:w-fssize3
    for y=1+fssize3:1:w-fssize3
        %遍历每个点的四周
        is=f(x-fssize3:1:x+fssize3,y-fssize3:1:y+fssize3);
        temp = is(:);
        resultImage(x,y)= max(temp);
    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('9*9最大值');

最小值

close all 
clear all

f=imread('D:/testData/filtering.tif');

[w,h]=size(f);
image= f(:,:);

fsize1=3;
fsize2=5;
fsize3=9;

fssize1=(fsize1-1)/2;
fssize2=(fsize2-1)/2;
fssize3=(fsize3-1)/2;

figure();
subplot(1,2,1);
imshow(image);
xlabel('原图像');

resultImage = image;

for x=1+fssize1:1:w-fssize1
    for y=1+fssize1:1:w-fssize1
        is=f(x-fssize1:1:x+fssize1,y-fssize1:1:y+fssize1);
        temp = is(:);
        resultImage(x,y)= min(temp);
    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('3*3最小值');

resultImage= f(:,:);
figure();
subplot(1,2,1);
imshow(f);
xlabel('原图像');

for x=1+fssize2:1:w-fssize2
    for y=1+fssize2:1:w-fssize2
        %遍历每个点的四周
        is=f(x-fssize2:1:x+fssize2,y-fssize2:1:y+fssize2);
        temp = is(:);
        resultImage(x,y)= min(temp);
    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('5*5最小值');

resultImage= f(:,:);
figure();
subplot(1,2,1);
imshow(f);
xlabel('原图像');

for x=1+fssize3:1:w-fssize3
    for y=1+fssize3:1:w-fssize3
        %遍历每个点的四周
        is=f(x-fssize3:1:x+fssize3,y-fssize3:1:y+fssize3);
        temp = is(:);
        resultImage(x,y)= min(temp);
    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('9*9最小值');

中点

close all 
clear all

f=imread('D:/testData/filtering.tif');

[w,h]=size(f);
image= f(:,:);

fsize1=3;
fsize2=5;
fsize3=9;

fssize1=(fsize1-1)/2;
fssize2=(fsize2-1)/2;
fssize3=(fsize3-1)/2;

figure();
subplot(1,2,1);
imshow(image);
xlabel('原图像');

resultImage = image;

for x=1+fssize1:1:w-fssize1
    for y=1+fssize1:1:w-fssize1
        is=f(x-fssize1:1:x+fssize1,y-fssize1:1:y+fssize1);
        temp = is(:);
        resultImage(x,y)= (max(temp) + min(temp))/2;
    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('3*3中点');

resultImage= f(:,:);
figure();
subplot(1,2,1);
imshow(f);
xlabel('原图像');

for x=1+fssize2:1:w-fssize2
    for y=1+fssize2:1:w-fssize2
        %遍历每个点的四周
        is=f(x-fssize2:1:x+fssize2,y-fssize2:1:y+fssize2);
        temp = is(:);
        resultImage(x,y)= (max(temp) + min(temp))/2;
    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('5*5中点');

resultImage= f(:,:);
figure();
subplot(1,2,1);
imshow(f);
xlabel('原图像');

for x=1+fssize3:1:w-fssize3
    for y=1+fssize3:1:w-fssize3
        %遍历每个点的四周
        is=f(x-fssize3:1:x+fssize3,y-fssize3:1:y+fssize3);
        temp = is(:);
        resultImage(x,y)= (max(temp) + min(temp))/2;
    end
end


subplot(1,2,2);
imshow(resultImage);
xlabel('9*9中点');
bllddee
关注
专栏目录
MATLAB图像均值滤波
03-29
对图像添加一个方差为0.05的椒盐噪声（自己可以添加高斯噪声），然后对图像进行3*3算数均值滤波，几何均值滤波，Q=1.5大小的逆谐波滤波处理，可以观察处理后的结果
MATLAB谐波均值滤波
临风暖阳的博客
04-13 2249
% 谐波均值滤波 clc,clear,close all % 清理命令区、清理工作区、关闭显示图形 warning off % 消除警告 feature jit off % 加速代码运行 im = imread(‘brain.bmp’); % 原图像 im = imnoise(im,‘gaussian’,0,1e-3); % 原图像 + 白噪声 im1 = harmonymean_filter(im,3,3); % 应用谐波均值滤波 figure(‘colo.
谐波均值滤波matlab,谐波均值滤波器.PPT
weixin_35486456的博客
03-19 1094
谐波均值滤波器数字图像处理Digital Image Processing 第6章 图像复原(Image Restoration) 第6章 图像复原(Image Restoration) 6.1 Fundamentals of Image Restoration and Degradation Model 6.1 Fundamentals of Image Restoration and De...
css的伪类选择器来实现一些图形
weixin_47450807的博客
09-20 580
css的伪类选择器来实现一些图形 爱心 图片为： <!DOCTYPE html> <html lang="cn"> <head> <meta charset="UTF-8"> <meta name="viewport" content="width=device-width, initial-scale=1.0"> <title>Document</title> <style>
matlab调和均值滤波_MATLAB--数字图像处理 图像噪声与滤波处理
weixin_39541693的博客
01-04 970
一、实验名称图像的噪声与滤波处理二、实验目的1.熟悉MATLAB软件的使用。2.掌握图像的噪声与滤波处理。三、实验内容1.对一张图片添加不同强度的高斯噪声、椒盐噪声2.对图像进行不同模板的均值滤波、高斯加权滤波、中值滤波，对比结果并分析不同滤波方式的优劣四、实验仪器与设备Win10 64位电脑MATLAB2017a五、实验原理高斯噪声 所谓高斯噪声是指它的概率密度函数服从高斯分布（即正态分布...
matlab调和均值滤波_数字图像处理基础 — 高斯滤波
weixin_31066503的博客
01-04 218
高斯滤波，本文主要讲其如何通过C语言实现。不太擅长写理论性质的文章，这里仅仅阐述自己怎么实现以及简单的优化过程。通常我们对获取的图像进行进一步处理时，往往需要先进行一次降噪，而通常我们选择的是高斯滤波也称高斯模糊，其因为较好的降噪效果而被广泛使用，具体的原因大家可以做实验以及看相关的论文去理解。如果你经常使用MATLAB或者Python，可能只需要调用一下库函数即可实现。以matlab为例，其高斯...
MATLAB几种均值滤波器的滤波效果比较_MATLAB均值滤波器比较_matlab_均值滤波器_
10-03
本篇将详细介绍几种均值滤波器——算术均值滤波器、几何均值滤波器、谐波均值滤波器和逆谐波均值滤波器，并通过MATLAB进行滤波效果的比较。 首先，均值滤波器是最基础的线性滤波器之一，其工作原理是用图像窗口内...
给我一个算术均值滤波器 几何均值滤波器 谐波均值滤波器 逆谐波均值滤波器这4个滤波器同时处理同一张图片的MATLAB代码
最新发布
05-14
以下是MATLAB代码示例： ```matlab % 读取图片 img = imread('test.jpg');...其中，算术均值滤波器和几何均值滤波器使用了 `imfilter` 函数，谐波均值滤波器和逆谐波均值滤波器使用了 `nlfilter` 函数。
1710121401-赵奕翔-数字图像处理-实验项目06实验报告1
08-08
此外，还尝试了几何均值滤波器，以及逆谐波滤波器，其中Q值分别设置为-1.5和1.5，观察它们对噪声图像的复原效果。 【基于顺序统计滤波器的图像复原】 顺序统计滤波器是一种非线性滤波器，它考虑了像素邻域内的像素...
算术滤波器：处理图像处理中的算术滤波器-matlab开发
05-30
可用于图像处理
用matlab处理 算术平均法 滤波
12-06
b1 是采样数据点 1024个，采样频率5000Hz 算术平均法 滤波 是m文件，没装matlab也可以用写字板打开
用matlab对一幅图像使用各类滤波器进行处理
12-11
总的来说，MATLAB提供了强大的图像处理功能，使得对图像应用算术、几何、谐波和逆谐波滤波变得简单易行。理解这些滤波器的工作原理和如何在MATLAB中实现它们，对于图像处理和分析的专业人士来说是非常有价值的。通过...
小波图像分割matlab代码-dip-project:数字图像处理项目
05-24
小波图像分割matlab代码数字图像处理项目 2017年秋季课程-数字图像处理 一些数字图像主题的简单matlab实现 话题 直方图均衡 空间增强方法复制填充，比例尺，均值，拉普拉斯算子，Sobel 频域滤波理想，巴特沃思，高斯...
空间滤波-反谐波平均滤波器
Henry的博客
10-22 1322
反谐波平均滤波器适用于降低或者消除椒盐噪声，它在处理的时候，要满足一些条件Q = 0 的时候，反谐波滤波器简化为算术平均滤波器，也就是邻域均值滤波器Q > 0 的时候，处理胡椒噪声Q < 0 的时候，处理盐粒噪声Q = -1 的时候，简化为谐波平均滤波器2. 代码实现算法实现的部分为：实现最困难的部分在pad和如果手动卷积，具体的过程可以看代码讲解部分，这里不做多余介绍在处理被kernel 包围区域的时候，我们不是让被处理的区域与滤波器加权平均，这里是非线性的操作。
[数字图像处理]图像去噪初步(1)--均值滤波器
weixin_34244102的博客
07-25 535
1.图像去噪的前言     上一篇博文中，我对噪声的类型进行了介绍，也使用的Matlab对各种噪声进行了实现。旧话重提，一幅图像，甚至是一个信号的老化，能够使用下面模型来表示。 能够使用下面算式来表示 这里，因为退化函数的作用，使得原图像产生退化(比方，运动模糊)，然后在加上一个加性噪声项。       本博文，主要对去除加性噪声的线性滤波器的性能进行了比較...
[数字图像处理] 逆谐波均值去噪
吴政亿的博客
05-17 8001
151220129 计科 吴政亿 在ch09上，老师对于逆谐波均值提出一个问题: f^(x,y)=∑(s,t)∈Sxyg(s,t)Q+1∑(s,t)∈Sxyg(s,t)Qf^(x,y)=∑(s,t)∈Sxyg(s,t)Q+1∑(s,t)∈Sxyg(s,t)Q\hat{f}(x,y)=\frac{\sum_{(s,t)\in S_{xy}}{g(s,t)^{Q+1}}}{\sum_{(s,t...
基于matlab的图像处理——算术均值、中值、几何滤波器
L_superhero的博客
05-24 3959
问题： 如下图所示，白条的宽度和高度分别是7像素和210像素，白条之间的间隔是17像素，采用以下滤波器对测试图像进行滤波，并输出结果 3x3算术均值滤波器 5x5算术均值滤波器 7x7算术均值滤波器 同理使用中值滤波器 同理使用几何滤波器 3x3算术均值滤波器 I=imread("511.png"); subplot(2,1,1); imshow(I) title("原图") I1= imfilter(I,fspecial('average',3)); subplot(2,1,2); imsh
MATLAB--空域滤波学习笔记
小席的学习日志
05-30 2724
空域滤波 空域滤波是指对图像中每个像素为中心的邻域进行预定义操作，然后用得到的结果代替邻域中心像素点的像素值。滤波器的中心访问输入图像中的每个像素，就生成了滤波后的图像。 空域滤波器可分为线性空域滤波器和非线性空域滤波器。如果在图像像素上执行的是线性操作，则该滤波器称为线性空间滤波器，否则，滤波器就称为非线性空间滤波器。 1、线性空域滤波 我们关注奇数尺寸m×n（假设m=2a+1且n=2b+1）的...
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值