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

最新推荐文章于 2023-06-14 17:34:39 发布
最新推荐文章于 2023-06-14 17:34:39 发布
阅读量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 581
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均值滤波器比较_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对一幅图像使用各类滤波器进行处理
12-11
总的来说，MATLAB提供了强大的图像处理功能，使得对图像应用算术、几何、谐波和逆谐波滤波变得简单易行。理解这些滤波器的工作原理和如何在MATLAB中实现它们，对于图像处理和分析的专业人士来说是非常有价值的。通过...
小波图像分割matlab代码-dip-project:数字图像处理项目
05-24
小波图像分割matlab代码数字图像处理项目 2017年秋季课程-数字图像处理 一些数字图像主题的简单matlab实现 话题 直方图均衡 空间增强方法复制填充，比例尺，均值，拉普拉斯算子，Sobel 频域滤波理想，巴特沃思，高斯...
空间滤波-反谐波平均滤波器
Henry的博客
10-22 1327
反谐波平均滤波器适用于降低或者消除椒盐噪声，它在处理的时候，要满足一些条件Q = 0 的时候，反谐波滤波器简化为算术平均滤波器，也就是邻域均值滤波器Q > 0 的时候，处理胡椒噪声Q < 0 的时候，处理盐粒噪声Q = -1 的时候，简化为谐波平均滤波器2. 代码实现算法实现的部分为：实现最困难的部分在pad和如果手动卷积，具体的过程可以看代码讲解部分，这里不做多余介绍在处理被kernel 包围区域的时候，我们不是让被处理的区域与滤波器加权平均，这里是非线性的操作。
[数字图像处理]图像去噪初步(1)--均值滤波器
weixin_34244102的博客
07-25 536
1.图像去噪的前言     上一篇博文中，我对噪声的类型进行了介绍，也使用的Matlab对各种噪声进行了实现。旧话重提，一幅图像，甚至是一个信号的老化，能够使用下面模型来表示。 能够使用下面算式来表示 这里，因为退化函数的作用，使得原图像产生退化(比方，运动模糊)，然后在加上一个加性噪声项。       本博文，主要对去除加性噪声的线性滤波器的性能进行了比較...
[数字图像处理] 逆谐波均值去噪
吴政亿的博客
05-17 8004
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实现图像空间域滤波复原
最新发布
anxinbuxinan的博客
06-14 745
自适应空间域滤波器：自适应空间域滤波器是指根据局部像素的统计特征，自适应地选择滤波器模板的大小和权值，并计算该像素的值。例如，自适应中值滤波器就是一种常见的自适应空间域滤波器，它能够根据局部像素的方差和中值等信息，自适应地选择滤波器模板的大小和性质，从而更好地处理图像中的噪声。顺序统计滤波器：顺序统计滤波器是指对某个像素周围的相邻像素按照一定顺序进行排序，并选择其中某个位置的像素值作为该像素的值。该滤波器的原理是用某个像素周围的相邻像素的平均值替代该像素的值，从而实现去噪的效果。
基于matlab的图像处理——算术均值、中值、几何滤波器
L_superhero的博客
05-24 3968
问题： 如下图所示，白条的宽度和高度分别是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逆谐波均值滤波
临风暖阳的博客
04-13 1429
% 逆谐波均值滤波 clc,clear,close all % 清理命令区、清理工作区、关闭显示图形 warning off % 消除警告 feature jit off % 加速代码运行 im = imread(‘brain.bmp’); % 原图像 im = imnoise(im,‘gaussian’,0,1e-3); % 原图像 + 白噪声 im1 = conharmmean_filter(im,3,3,3/2); % 应用逆谐波均值滤波 figure.
Matlab滤波器设计示例
热门推荐
chenxy_bwave的专栏
11-20 1万+
Matlab提供了非常丰富的滤波器设计工具，其中最基本的是designfilt()函数。designfilt()函数是通用的滤波器设计工具，可以用于设计各种类型的数字滤波器，以及模拟滤波器的数字近似。Matlab还提供了其它很多专用类型的滤波器设计函数，还有基于图形化界面进行滤波器设计的工具filterDesigner(旧版本为fdatool)。本文综合但简要地介绍各种matlab滤波器设计工具，顺便代码示例介绍用于滤波器频谱观测用的工具、以及如何使用所生成的滤波器对象。
【OpenCV 例程200篇】96. 谐波平均滤波器
youcans的博客
02-11 3561
谐波平均滤波器既能处理盐粒噪声（白色噪点），又能处理类似于高斯噪声的其他噪声，但不能处理胡椒噪声（黑色噪点）
图像恢复及滤波处理
火雨(Nick)
12-26 1万+
图像恢复及滤波处理 1.基本概念：       a.图像恢复是通过计算机处理，对质量下降的图像加以重建或恢复的处理过程。因摄像机与物体相对运动、系统误差、畸变、噪声等因素的影响，使图像往往不是真实景物的完善映像。在图像恢复中，需建立造成图像质量下降的退化模型，然后运用相反过程来恢复原来图像，并运用一定准则来判定是否得到图像的最佳恢复。在遥感图像处理中，为消除遥感图像的失真、畸变，恢复目标的反射
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值