位置: IT常识 - 正文

数字图像处理之matlab实验(三):空间滤波器(数字图像处理-应用篇)

编辑:rootadmin
数字图像处理之matlab实验(三):空间滤波器 1、空间滤波原理

推荐整理分享数字图像处理之matlab实验(三):空间滤波器(数字图像处理-应用篇),希望有所帮助,仅作参考,欢迎阅读内容。

文章相关热门搜索词:数字图像处理之父,数字图像处理之父,数字图像处理之MATLAB图像基本处理,数字图像处理之怀旧特效,数字图像处理之MATLAB图像基本处理,数字图像处理之MATLAB图像基本处理,数字图像处理之怀旧特效,数字图像处理之父,内容如对您有帮助,希望把文章链接给更多的朋友!

空间滤波,就是在原图像上,用一个固定尺寸的模板去做卷积运算,得到的新图像就是滤波结果。滤波,就是过滤某种信号的意思。过滤哪种信号取决于模板设计,如果是锐化模板,处理后就保留高频信号,如果是平滑模板,处理后就保留低频信号。

(1)模板运算

图像处理中模板能够看作是n*n(n通常是奇数)的窗体。模板连续地运动于整个图像中,对模板窗体范围内的像素做相应处理。

模板运算主要分为:

①模板卷积。模板卷积是把模板内像素的灰度值和模板中相应的灰度值相乘,求平均值赋给当前模板窗体的中心像素。作为它的灰度值;

②模板排序。模板排序是把模版内像素的灰度值排序,取某个顺序统计量作为模板中心像素灰度值。

Matlab中做模版卷积十分高效,取出模版内子矩阵和模版权重点乘求平均就可以,已图示为例,3X3的模板在图像上滑动,原图像f(x,y) 经过模板处理后变成了g(x,y)。

(2)边界处理

处理边界有非常多种做法:

①重复图像边缘上的行和列。

②卷绕输入图像(假设第一列紧接着最后一列)。

③在输入图像外部填充常数(例如零)。

④去掉不能计算的行列。仅对可计算的像素计算卷积。

(3)空间域滤波

把模板运算运用于图像的空间域增强的技术称为空间域滤波。依据滤波频率空间域滤波分为平滑滤波(减弱和去除高频分量)和锐化滤波(减弱和去除低频分量),依据滤波计算特点又分为线性滤波和非线性滤波。

因此空间域滤波可分为:

数字图像处理之matlab实验(三):空间滤波器(数字图像处理-应用篇)

分类     线性       非线性

平滑   线性平滑   非线性平滑

锐化   线性锐化   非线性锐化

2、平滑滤波器(1)添加噪声

噪声主要分类为两类,高斯噪声和椒盐噪声。

高斯噪声在每个像素上都会出现,赋值服从高斯分布。

椒盐噪声出现位置随机,所以可以控制椒盐噪声的密度,椒盐噪声的幅度确定,椒噪声偏暗,盐噪声偏亮。

Image = mat2gray( imread('original_pattern.jpg') ,[0 255]);noiseIsp=imnoise(Image,'salt & pepper',0.1); %添加椒盐噪声,密度为0.1 imshow(noiseIsp,[0 1]); title('椒盐噪声图像');noiseIg=imnoise(Image,'gaussian'); %添加高斯噪声,默认均值为0,方差为0.01figure;imshow(noiseIg,[0 1]); title('高斯噪声图像'); (2)平滑滤波器

 平滑滤波器可以去除图像的噪声,使图像变得模糊。包括:中值滤波、均值滤波、高斯滤波。

高斯滤波、均值滤波去除高斯噪声。

(3)均值滤波

Image=imread('Letters-a.jpg');noiseI=imnoise(Image,'gaussian'); %添加高斯噪声subplot(221),imshow(Image),title('原图');subplot(222),imshow(noiseI),title('高斯噪声图像');result1=filter2(fspecial('average',3),noiseI); %3×3均值滤波result2=filter2(fspecial('average',7),noiseI); % 7×7均值滤波subplot(223),imshow(uint8(result1)),title('3×3均值滤波');subplot(224),imshow(uint8(result2)),title('7×7均值滤波');(4)中值滤波

 

Image=rgb2gray(imread('lotus.bmp'));noiseI=imnoise(Image,'salt & pepper',0.1);result=medfilt2(noiseI); %3×3中值滤波subplot(121),imshow(noiseI),title('椒盐噪声图像');subplot(122),imshow(uint8(result)),title('3×3中值滤波');(5)自编程实现高斯滤波

Image=imread('Letters-a.jpg');sigma1=0.6; sigma2=10; r=3; % 高斯模板的参数NoiseI= imnoise(Image,'gaussian'); %加噪gausFilter1=fspecial('gaussian',[2*r+1 2*r+1],sigma1); gausFilter2=fspecial('gaussian',[2*r+1 2*r+1],sigma2); result1=imfilter(NoiseI,gausFilter1,'conv');result2=imfilter(NoiseI,gausFilter2,'conv');subplot(231);imshow(Image);title('原图');subplot(232);imshow(NoiseI);title('高斯噪声图像');subplot(233);imshow(result1);title('sigma1 =0.6高斯滤波');subplot(234);imshow(result2);title('sigma2 =10高斯滤波');%imwrite(uint8(NoiseI),'gr.bmp');%imwrite(uint8(result1),'gr1.bmp');%imwrite(uint8(result2),'gr2.bmp');%编写高斯滤波函数实现[height,width]=size(NoiseI); for x=-r:r for y=-r:r H(x+r+1,y+r+1)=1/(2*pi*sigma1^2).*exp((-x.^2-y.^2)/(2*sigma1^2)); endend H=H/sum(H(:)); %归一化高斯模板H result3=zeros(height,width); %滤波后图像 midimg=zeros(height+2*r,width+2*r); %中间图像 midimg(r+1:height+r,r+1:width+r)=NoiseI; for ai=r+1:height+r for aj=r+1:width+r temp_row=ai-r; temp_col=aj-r; temp=0; for bi=1:2*r+1 for bj=1:2*r+1 temp= temp+(midimg(temp_row+bi-1,temp_col+bj-1)*H(bi,bj)); end end result3(temp_row,temp_col)=temp; end end subplot(235);imshow(uint8(result3));title('myself高斯滤波'); 3、锐化滤波器

 (1)梯度算子

Image=im2double(rgb2gray(imread('lotus.jpg')));subplot(131),imshow(Image),title('原图像');[h,w]=size(Image);edgeImage=zeros(h,w);for x=1:w-1 for y=1:h-1 edgeImage(y,x)=abs(Image(y,x+1)-Image(y,x))+abs(Image(y+1,x)-Image(y,x)); endendsubplot(132),imshow(edgeImage),title('梯度图像');sharpImage=Image+edgeImage;subplot(133),imshow(sharpImage),title('锐化图像');(2)Robert算子 

Image=im2double(rgb2gray(imread('lotus.jpg')));subplot(221),imshow(Image),title('原图像');BW= edge(Image,'roberts');subplot(222),imshow(BW),title('边缘检测');H1=[1 0; 0 -1];H2=[0 1;-1 0];R1=imfilter(Image,H1);R2=imfilter(Image,H2);edgeImage=abs(R1)+abs(R2);subplot(223),imshow(edgeImage),title('Robert梯度图像');sharpImage=Image+edgeImage;subplot(224),imshow(sharpImage),title('Robert锐化图像');(3)Sobel算子 

Image=im2double(rgb2gray(imread('lotus.jpg')));subplot(221),imshow(Image),title('原图像');BW= edge(Image,'sobel');subplot(222),imshow(BW),title('边缘检测');H1=[-1 -2 -1;0 0 0;1 2 1];H2=[-1 0 1;-2 0 2;-1 0 1];R1=imfilter(Image,H1);R2=imfilter(Image,H2);edgeImage=abs(R1)+abs(R2);subplot(223),imshow(edgeImage),title('Sobel梯度图像');sharpImage=Image+edgeImage;subplot(224),imshow(sharpImage),title('Sobel锐化图像');(4)多个模板边缘检测 

 

clear,clc,close all;Image=im2double(rgb2gray(imread('lotus.jpg')));H1=[-1 -1 -1;0 0 0;1 1 1];H2=[0 -1 -1;1 0 -1; 1 1 0];H3=[1 0 -1;1 0 -1;1 0 -1];H4=[1 1 0;1 0 -1;0 -1 -1];H5=[1 1 1;0 0 0;-1 -1 -1];H6=[0 1 1;-1 0 1;-1 -1 0];H7=[-1 0 1;-1 0 1;-1 0 1];H8=[-1 -1 0;-1 0 1;0 1 1];R1=imfilter(Image,H1);R2=imfilter(Image,H2);R3=imfilter(Image,H3);R4=imfilter(Image,H4);R5=imfilter(Image,H5);R6=imfilter(Image,H6);R7=imfilter(Image,H7);R8=imfilter(Image,H8);edgeImage1=abs(R1)+abs(R7);sharpImage1=edgeImage1+Image;f1=max(max(R1,R2),max(R3,R4));f2=max(max(R5,R6),max(R7,R8));edgeImage2=max(f1,f2);sharpImage2=edgeImage2+Image;subplot(221),imshow(edgeImage1),title('两个模板边缘检测');subplot(222),imshow(edgeImage2),title('八个模板边缘检测');subplot(223),imshow(sharpImage1),title('两个模板边缘锐化');subplot(224),imshow(sharpImage2),title('八个模板边缘锐化');(5)Laplacian算子

Image=im2double(rgb2gray(imread('lotus.jpg')));subplot(131),imshow(Image),title('原图像');H=fspecial('laplacian',0);R=imfilter(Image,H);edgeImage=abs(R);subplot(132),imshow(edgeImage),title('Laplacian梯度图像');H1=[0 -1 0;-1 5 -1;0 -1 0];sharpImage=imfilter(Image,H1);subplot(133),imshow(sharpImage),title('Laplacian锐化图像');
本文链接地址:https://www.jiuchutong.com/zhishi/295182.html 转载请保留说明!

上一篇:【2】YOLOv8原理解析:重新定义实时目标检测的速度和精度(yolov3原理及代码解析)

下一篇:uni-app 从0 到 1 制作一个项目,收藏等于学会(uni-app实战教程)

  • 苹果13用的什么屏幕(苹果13用的什么充电接口)

    苹果13用的什么屏幕(苹果13用的什么充电接口)

  • 三星buds2怎么开启游戏模式(三星buds2怎么开机)

    三星buds2怎么开启游戏模式(三星buds2怎么开机)

  • iphone6更新ios12卡吗(苹果6手机升级到了12系统后卡了怎么办)

    iphone6更新ios12卡吗(苹果6手机升级到了12系统后卡了怎么办)

  • 华为nova5的文字扫描功能在哪(华为nova5文字识别怎么使用)

    华为nova5的文字扫描功能在哪(华为nova5文字识别怎么使用)

  • soulmate关系可以解除吗(soulmate关系有啥用)

    soulmate关系可以解除吗(soulmate关系有啥用)

  • 无线鼠标怎么看充满电(无线鼠标怎么看电量)

    无线鼠标怎么看充满电(无线鼠标怎么看电量)

  • oppoa59怎么分屏(oppoa59手机怎样分屏)

    oppoa59怎么分屏(oppoa59手机怎样分屏)

  • excel背景图片怎么设置(excel背景图片怎么设置平铺)

    excel背景图片怎么设置(excel背景图片怎么设置平铺)

  • 修改微信号条件(修改微信号满足条件)

    修改微信号条件(修改微信号满足条件)

  • 华为matepad可以装Windows吗(华为matepad可以做ppt吗)

    华为matepad可以装Windows吗(华为matepad可以做ppt吗)

  • 微信语音对方一直无应答(微信语音对方一直忙线中是不是设置了)

    微信语音对方一直无应答(微信语音对方一直忙线中是不是设置了)

  • 笔记本可以开热点给手机用吗(笔记本可以开热点吗win7)

    笔记本可以开热点给手机用吗(笔记本可以开热点吗win7)

  • qq匹配聊天失败是什么原因(qq匹配聊天显示匹配失败)

    qq匹配聊天失败是什么原因(qq匹配聊天显示匹配失败)

  • 拍立得为什么拍出来是黑色的(拍立得为什么拍出来偏的)

    拍立得为什么拍出来是黑色的(拍立得为什么拍出来偏的)

  • creative cloud卸载不掉(creative cloud卸载后果)

    creative cloud卸载不掉(creative cloud卸载后果)

  • 抖音显示超级乐迷什么意思(抖音显示超级乐队怎么弄)

    抖音显示超级乐迷什么意思(抖音显示超级乐队怎么弄)

  • 天猫国际和天猫有什么区别(天猫国际和天猫国际直营店有区别吗)

    天猫国际和天猫有什么区别(天猫国际和天猫国际直营店有区别吗)

  • 掌上生活大学生认证有风险吗(掌上生活大学生额度)

    掌上生活大学生认证有风险吗(掌上生活大学生额度)

  • id全称是什么(id是什么)

    id全称是什么(id是什么)

  • 光驱是多媒体信息吗(光驱是多媒体信号吗)

    光驱是多媒体信息吗(光驱是多媒体信号吗)

  • z370和z390主板区别(z370和z390主板区别知乎)

    z370和z390主板区别(z370和z390主板区别知乎)

  • 中国移动微信公众号怎么关注(中国移动微信公众号)

    中国移动微信公众号怎么关注(中国移动微信公众号)

  • 怎么查找我创建的微信群(怎么查找我创建的小程序)

    怎么查找我创建的微信群(怎么查找我创建的小程序)

  • 快手直播黄钻是快币吗(快手直播黄钻是谁给的)

    快手直播黄钻是快币吗(快手直播黄钻是谁给的)

  • 华硕笔记本连不了网怎么办(华硕笔记本连不上5g无线网络)

    华硕笔记本连不了网怎么办(华硕笔记本连不上5g无线网络)

  • 小米手环怎么调亮度(小米手环怎么调时间日期)

    小米手环怎么调亮度(小米手环怎么调时间日期)

  • 筛选状态下如何粘贴(筛选状态下如何快速下拉填充)

    筛选状态下如何粘贴(筛选状态下如何快速下拉填充)

  • 收取违约金是否交税
  • 5月份企业所得税
  • 收到个人保险费会计分录
  • 企业所得税季度申报时间
  • 合并报表营业收入怎么算出来的
  • 企业预缴的增值税收入算不算企业的销售收入
  • 小规模纳税人取的专票能抵扣吗
  • 小额贷款公司的钱从哪里来
  • 公交车收购
  • 研究开发费用加计扣除多少
  • 销售人员的提成属于什么工资
  • 计提工资数大于实际支付数怎么办?
  • 工会账户销户划拨流程
  • 特殊销售方式的计税依据
  • 财产租赁所得个人所得税计算方法2020
  • 退休人员是否还继续享受社保
  • 小规模纳税人咨询费税率
  • 增值税报表附表三如何填写
  • 不适用研究开发费用税前加计扣除政策的有
  • 进项发票已入账抵扣后开负数红冲怎么处理
  • 待抵扣进项税额转出会计分录
  • 来料加工出口退税申报操作
  • 费用转入本年利润
  • 每个季度零申报需要做什么
  • 附加税的税率是多少2023年
  • 往来账乱如何处理
  • 个人综合所得适用什么税率?
  • thinkpad x230笔记本
  • window10系统有几种
  • 合同资产减值怎么填
  • php中字符串函数
  • 股权转让协议的注意事项
  • 固定资产清理是三栏式还是多栏式
  • 塞里雅兰瀑布图片高清
  • 固定资产非正常中断没超过三个月
  • 悉尼奥林匹克公园
  • php构造函数重载
  • php sha1加密 解密
  • 销售费用凭证
  • 人工智能大模型体验报告3.0
  • deepsort 跟踪
  • 融资租赁的会计科目
  • 出口退税率和进口税率的区别
  • 装修公司开的劳务发票
  • js怎么存储数据
  • mysql备份原理
  • 销售成本是付现成本吗
  • 生育津贴现在可以打个人账户么
  • 小企业会计准则和企业会计准则的区别
  • 发票上的货物要怎么填?
  • 收到税务局退回税款会计分录
  • 增值税进项税额加计抵减政策
  • 归纳投入外币资产怎么写
  • 不得免征和抵扣税额抵减额分录
  • 员工提前预支工资合法吗
  • 预缴土地增值税的会计处理
  • 进项发票超过期限不认证怎么办
  • 事业单位私车公用如何处罚
  • 会计分录要怎么写
  • 营业外支出增加说明了什么
  • 苹果序列号查询
  • 微软出win9了吗
  • linux系统汉化教程
  • win8系统如何查看电脑mac物理地址
  • awk中执行命令
  • win8系统桌面图标
  • eclipse4.9.0安装windowbuilder
  • javascript typeof id==='string'?document.getElementById(id):id解释 原创
  • unity smooth
  • scrollview不滚动
  • sockaddr_in和sockaddr
  • Unity3d MonoDevelop格式化代码
  • docker编写dockerfile
  • ubuntu如何下载
  • jquery 树结构
  • 苏州封闭式高中
  • 达州房管局备案查询系统
  • 注销税务登记存货如何处理
  • 小规模纳税人的认定标准是什么
  • 明星征税多少
  • 免责声明:网站部分图片文字素材来源于网络,如有侵权,请及时告知,我们会第一时间删除,谢谢! 邮箱:opceo@qq.com

    鄂ICP备2023003026号

    网站地图: 企业信息 工商信息 财税知识 网络常识 编程技术

    友情链接: 武汉网站建设