kuwahara_filter滤波器----MATLAB

%% kuwahara_filter滤波器
clc,clear,close all  % 清理命令区、清理工作区、关闭显示图形
warning off       % 消除警告
feature jit off      % 加速代码运行
[filename ,pathname]=...uigetfile({'*.bmp';'*.tif';'*.jpg';},'选择图片'); %选择图片路径
str=[pathname filename]; % 合成路径+文件名
im = imread(str);        % 读图
% 转化为灰度图像
if size(im,3)==1im = im;
elseim = rgb2gray(im);
end
im = imnoise(im,'gaussian',0,1e-3); % 原图像 + 白噪声figure,
subplot(121),imshow(im);title('原始图像')
colormap(jet)  % 颜色
shading interp % 消隐
im1 = kuwahara(im,5);
subplot(122),imshow(im1);title('kuwahara滤波图像')
colormap(jet)  % 颜色
shading interp % 消隐
function [Y,Xpad] = kuwahara(X,WINSZ)
% kuwahara_filter滤波器
% 对图像边界轮廓有较强的保护作用
% kuwahara nonlinear edge-preserving filtering
% 函数输入：
%          X：二维图像矩阵
%          WINSZ： window size
% 函数输出：
%          Y：滤波图像
%          Xpad：点扩展矩阵
% The Kuwahara filter 4块，(最直观的如下所示，5x5 pixels).
%
%    ( a  a  ab   b  b)
%    ( a  a  ab   b  b)
%    (ac ac abcd bd bd)
%    ( c  c  cd   d  d)
%    ( c  c  cd   d  d)if nargin < 2WINSZ = 5;  % 默认值
endif ~isa(X,'double')X = im2double(X);  % double类型
endPad = floor(WINSZ/2);  % 向下取整Xpad = padarray(X,[Pad, Pad],'replicate');  % 点扩展
% A = [ 1 2;
%       3 4 ];
% B = padarray(A,[3 2],'replicate','post')
% B =
%      1     2     2     2
%      3     4     4     4
%      3     4     4     4
%      3     4     4     4
%      3     4     4     4
[padRows,padCols] = size(Xpad);                    % 求维数
Y = zeros(size(X));nRowIters = length((Pad+1):(padRows-Pad));
count = 1;
for i = (Pad+1):(padRows-Pad)for j = (Pad+1):(padCols-Pad)% 分成每一个小窗，Q1~Q4w = Xpad((i-Pad):(i+Pad),(j-Pad):(j+Pad));wnw = w(1:(Pad+1),1:(Pad+1));wne = w(1:(Pad+1),(Pad+1):WINSZ);wsw = w((Pad+1):WINSZ,1:(Pad+1));wse = w((Pad+1):WINSZ,(Pad+1):WINSZ);% 计算方差s = var([wnw(:),wne(:),wsw(:),wse(:)]);m = mean([wnw(:),wne(:),wsw(:),wse(:)]);[y,k] = min(s);  % 计算均值Y(i,j) = m(k);endend

kuwahara_filter滤波器----MATLAB相关推荐

常用滤波器Matlab程序设计
常用滤波器Matlab程序设计 (低通滤波器.高通滤波器.带通滤波器.带阻滤波器) 以下四个滤波器都是切比雪夫I型数字滤波器 1.低通滤波器低通滤波(Low-pass filter) 是一种过滤 ...
并联型APF/有源电力滤波器/Matlab/Simulink仿真 *dq/FBD谐波/无功检测
并联型APF/有源电力滤波器/Matlab/Simulink仿真 *dq/FBD谐波/无功检测 *两相旋转坐标系(dq).两相静止坐标系(αβ)下的PI控制 *SVPWM调制方式 (含仿真介绍文档,不 ...
并联型APF/有源电力滤波器/Matlab/Simulink仿真 dq/FBD谐波/无功检测
并联型APF/有源电力滤波器/Matlab/Simulink仿真 *dq/FBD谐波/无功检测 *两相旋转坐标系(dq).两相静止坐标系(αβ)下的PI控制 *SVPWM调制方式 (含仿真介绍文档) ...
有源电力滤波器matlab仿真，并联型apf仿真fft分析谐波电流检测ipiq法
有源电力滤波器matlab仿真, 并联型apf仿真fft分析谐波电流检测ipiq法跟踪电流控制(传统滞环控制空间电压矢量滞环控制) 总谐波畸变率降至3%以下 ID:695064569089802 ...
fir抽取滤波器matlab,Xilinx FPGA之抽取FIR滤波器的实现
过去半年有几位客户请我帮助他们设计和实现数字下变频器所用的降采样(即"抽取")滤波器,这种滤波器在软件无线电与数据采集类应用中都很常见.这项工作即便对于经验丰富的设计师也不是一件小 ...
语音怎么进入滤波器matlab,基于Matlab的语音信号滤波器的设计与实现
(4)双线性变换法设计Butterworth滤波器程序如下: fs=22050; x1=wavread('wo.wav'); t=0:1/22050:(size(x1)-1)/22050; Au=0 ...
fir滤波器matlab实现_关于FIRamp;IIR系统的算法说明以及结果验证（1）
首先,做一个简短的开场白. 本贴主要是一个关于信号处理方面的学习笔记.主要目的有三, 作为研究笔记留存分享我个人的理解与专业人士进行意见交换另外由于楼主不是教课员,因此无法保证算法分析以及个人理 ...
matlab中多边形滤波器,几种常见空间滤波器MATLAB实现
本文链接:https://blog.csdn.net/LYduring/article/details/80443573 一.目的实现算术均值滤波器.几何均值滤波器.中值滤波器.修正的阿尔法均值滤波 ...
双边滤波器—— Matlab实现
例:先用双边滤波器(BF)对原图像进行滤波得到低频部分,原图和低频作差后得到高频分量,高频分量和低频分量分别增强后再进行合成. 双边滤波的特点是保边去噪,相较于高斯滤波,在平滑图像的同时,增加了对图像 ...

kuwahara_filter滤波器----MATLAB

kuwahara_filter滤波器----MATLAB相关推荐

最新文章

热门文章