【图像去噪】基于matlab GUI HSI彩色图像去噪【含Matlab源码 1786期】

举报
海神之光 发表于 2022/05/28 22:28:23 2022/05/28
【摘要】 一、HSI彩色图像去噪简介 针对彩色图像中噪声难以去除的问题，根据HSI空间独特的色彩分离特点，对受高噪声污染的彩色图像的噪声去除进行了研究。首先将彩色图像投影到色彩特征空间HSI中，将色彩信息与亮度特...

一、HSI彩色图像去噪简介

针对彩色图像中噪声难以去除的问题，根据HSI空间独特的色彩分离特点，对受高噪声污染的彩色图像的噪声去除进行了研究。首先将彩色图像投影到色彩特征空间HSI中，将色彩信息与亮度特征信息进行分离操作，然后对该空间中的色彩分量H和S提出极坐标下距离阈值去噪方法进行处理，在保持色彩不失真的情况下去除噪声。

HSI空间中，色调H与饱和度S均表示颜色信息，H分量通常以角度度量，S通常表示离I轴的远近程度，人类感知的色彩是这两个分量的合成。文献[2]中建立的色度模型保持了图像的色彩特性。文献[3]中基于这种模型，提出色度中值滤波器去除色彩中的噪声，虽然保真色彩，但滤掉的噪声点很少。本文在色度模型的基础上，在保持颜色不失真的情况下，构造了平均最小距离阈值，尽可能滤掉更多的噪声点，提出了极坐标中距离阈值去噪（distance threshold denoising in polar coordinates,DST），对H和S去除噪声，具体处理过程如下：

步骤1分块预处理。将彩色图像进行HSI分解后，对饱和度分量矩阵S和色调分量矩阵H进行相同的分块处理，块大小为r×c。将两个分量矩阵中任意相同位置的块记为块矩阵s、块矩阵h，并寻找块中像素点数值的最小值，分别记为sm in和hm in。

步骤2极坐标系下块中颜色量的定义。设块中任一点(m,n)处表示颜色的量为Cmn，且满足m=1,2,⋯,r;  n=1,2,⋯,c，若smn、hmn分别表示点(m,n)处的饱和度值和色调值，则对于块矩阵s和h，有smn=s(m,n),hmn=h(m,n)，那么在极坐标系下可定义：
在这里插入图片描述
其中，i为虚数单位。这说明用块矩阵s和h联合定义的块中每一个表示颜色的量都是个复数，点(i,j)处的颜色量Cij也是个复数。

步骤3距离阈值的设置。块中的点(i,j)与其他任一点(m,n)的距离dmnij可由下式计算。
在这里插入图片描述
其中，在块矩阵s和h中，sij、hij是点(i,j)处的饱和度值和色调值，即sij=s(i,j),hij=h(i,j);smn、hmn是点(m,n)处的饱和度值和色调值，即smn=s(m,n),hmn=h(m,n)。

当m=1,2,⋯,r;n=1,2,⋯,c时，对于固定点(i,j),所有的dmnij可以构成像素点(i,j)与块中其他任意像素点的距离矩阵Dij∈Rr×c，设Dij中非0的最小值为dij，则块中每一个点对应一个非0的最小距离值。所有最小距离值的平均值可通过下式计算得到。
在这里插入图片描述
步骤4噪声的消除。若在块中点(i,j)处有dij>d，则说明同等条件下该像素点与周围其他像素点的差距较大，相似性太小，在此处产生了突变。根据噪声点的突变型，可判断该点是噪声点，需要去除掉。设块中去除噪声后的饱和度和色调块矩阵分别为ends和endh，具体处理过程如下所示。
在这里插入图片描述
其中，i=1,2,⋯,r;j=1,2,⋯,c。对ends和endh进行分块逆变换，得到最终的去除噪声后的饱和度矩阵end S和色调矩阵end H。

二、部分源代码

function varargout = myHSI(varargin)
% MYHSI MATLAB code for myHSI.fig
%      MYHSI, by itself, creates a new MYHSI or raises the existing
%      singleton*.
%
%      H = MYHSI returns the handle to a new MYHSI or the handle to
%      the existing singleton*.
%
%      MYHSI('CALLBACK',hObject,eventData,handles,...) calls the local
%      function named CALLBACK in MYHSI.M with the given input arguments.
%
%      MYHSI('Property','Value',...) creates a new MYHSI or raises the
%      existing singleton*.  Starting from the left, property value pairs are
%      applied to the GUI before myHSI_OpeningFcn gets called.  An
%      unrecognized property name or invalid value makes property application
%      stop.  All inputs are passed to myHSI_OpeningFcn via varargin.
%
%      *See GUI Options on GUIDE's Tools menu.  Choose "GUI allows only one
%      instance to run (singleton)".
%
% See also: GUIDE, GUIDATA, GUIHANDLES

% Edit the above text to modify the response to help myHSI

% Last Modified by GUIDE v2.5 04-Mar-2022 13:01:39

% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name',       mfilename, ...
                   'gui_Singleton',  gui_Singleton, ...
                   'gui_OpeningFcn', @myHSI_OpeningFcn, ...
                   'gui_OutputFcn',  @myHSI_OutputFcn, ...
                   'gui_LayoutFcn',  [] , ...
                   'gui_Callback',   []);
if nargin && ischar(varargin{1})
    gui_State.gui_Callback = str2func(varargin{1});
end

if nargout
    [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
    gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT


% --- Executes just before myHSI is made visible.
function myHSI_OpeningFcn(hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject    handle to figure
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
% varargin   command line arguments to myHSI (see VARARGIN)

% Choose default command line output for myHSI
handles.output = hObject;
axes(handles.axes1);
set(gca, 'vis', 'off') 
axes(handles.axes2);
set(gca, 'vis', 'off') 
axes(handles.axes3);
set(gca, 'vis', 'off') 
axes(handles.axes4);
set(gca, 'vis', 'off') 
% Update handles structure
guidata(hObject, handles);

% UIWAIT makes myHSI wait for user response (see UIRESUME)
% uiwait(handles.figure1);


% --- Outputs from this function are returned to the command line.
function varargout = myHSI_OutputFcn(hObject, eventdata, handles) 
% varargout  cell array for returning output args (see VARARGOUT);
% hObject    handle to figure
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Get default command line output from handles structure
varargout{1} = handles.output;


% --- 打开图像.
function pushbutton1_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton1 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
filter = { ...
'*.bmp;*.jpg;*.gif;*.emf', '所有图像文件 (*.bmp; *.jpg; *.gif; *emf)'; ...
'*.bmp', '位图文件 (*.bmp)'; ...
'*.jpg', 'JPEG文件 (*.jpg)'; ...
'*.gif', 'GIF文件 (*.gif)'; ...
'*.emf', '图元文件 (*.emf)'; ...
'*.*', '所有文件 (*.*)' ...
};

% 选择文件
[filename, pathname] = uigetfile( filter, '打开...');
if isequal(filename,0) | isequal(pathname,0), return, end

% 读入图像并显示
[X, map] = imread([pathname filename]);
% colormap(map)
% image(X)
axes(handles.axes1); cla reset;
imshow(X);title('原图');
% 设置坐标系为适合显示图像的方式
% axis image
set(gca, 'vis', 'off') 
handles.I=X;
guidata(hObject, handles);


  
 
    
  
  • 1
    • 
  
  • 2
    • 
  
  • 3
    • 
  
  • 4
    • 
  
  • 5
    • 
  
  • 6
    • 
  
  • 7
    • 
  
  • 8
    • 
  
  • 9
    • 
  
  • 10
    • 
  
  • 11
    • 
  
  • 12
    • 
  
  • 13
    • 
  
  • 14
    • 
  
  • 15
    • 
  
  • 16
    • 
  
  • 17
    • 
  
  • 18
    • 
  
  • 19
    • 
  
  • 20
    • 
  
  • 21
    • 
  
  • 22
    • 
  
  • 23
    • 
  
  • 24
    • 
  
  • 25
    • 
  
  • 26
    • 
  
  • 27
    • 
  
  • 28
    • 
  
  • 29
    • 
  
  • 30
    • 
  
  • 31
    • 
  
  • 32
    • 
  
  • 33
    • 
  
  • 34
    • 
  
  • 35
    • 
  
  • 36
    • 
  
  • 37
    • 
  
  • 38
    • 
  
  • 39
    • 
  
  • 40
    • 
  
  • 41
    • 
  
  • 42
    • 
  
  • 43
    • 
  
  • 44
    • 
  
  • 45
    • 
  
  • 46
    • 
  
  • 47
    • 
  
  • 48
    • 
  
  • 49
    • 
  
  • 50
    • 
  
  • 51
    • 
  
  • 52
    • 
  
  • 53
    • 
  
  • 54
    • 
  
  • 55
    • 
  
  • 56
    • 
  
  • 57
    • 
  
  • 58
    • 
  
  • 59
    • 
  
  • 60
    • 
  
  • 61
    • 
  
  • 62
    • 
  
  • 63
    • 
  
  • 64
    • 
  
  • 65
    • 
  
  • 66
    • 
  
  • 67
    • 
  
  • 68
    • 
  
  • 69
    • 
  
  • 70
    • 
  
  • 71
    • 
  
  • 72
    • 
  
  • 73
    • 
  
  • 74
    • 
  
  • 75
    • 
  
  • 76
    • 
  
  • 77
    • 
  
  • 78
    • 
  
  • 79
    • 
  
  • 80
    • 
  
  • 81
    • 
  
  • 82
    • 
  
  • 83
    • 
  
  • 84
    • 
  
  • 85
    • 
  
  • 86
    • 
  
  • 87
    • 
  
  • 88
    • 
  
  • 89
    • 
  
  • 90
    • 
  
  • 91
    • 
  
  • 92
    • 
  
  • 93
    • 
  
  • 94
    • 
  
  • 95
    • 
  
  • 96
    • 
  
  • 97
    • 
  
  • 98
    • 
  
  • 99
    • 
  
  • 100
    • 
  
  • 101
    • 
  
  • 102
    • 
  
  • 103
    • 
  
  • 104
    • 
  
  • 105
    • 
  
  • 106
    • 
  
  • 107
    • 
  
  • 108
    • 
  
  • 109
    • 
  
  • 110
    • 
  
  • 111
    • 
  
  • 112
    •

三、运行结果

在这里插入图片描述

四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1] 蔡利梅.MATLAB图像处理——理论、算法与实例分析[M].清华大学出版社，2020.
[2]杨丹,赵海滨,龙哲.MATLAB图像处理实例详解[M].清华大学出版社，2013.
[3]周品.MATLAB图像处理与图形用户界面设计[M].清华大学出版社，2013.
[4]刘成龙.精通MATLAB图像处理[M].清华大学出版社，2015.
[5]杨培,高雷阜,訾玲玲.HSI空间上高噪声彩色图像去噪方法研究[J].计算机科学与探索. 2020,14(11)

文章来源: qq912100926.blog.csdn.net，作者：海神之光，版权归原作者所有，如需转载，请联系作者。

原文链接：qq912100926.blog.csdn.net/article/details/123534112

【版权声明】本文为华为云社区用户转载文章，如果您发现本社区中有涉嫌抄袭的内容，欢迎发送邮件进行举报，并提供相关证据，一经查实，本社区将立刻删除涉嫌侵权内容，举报邮箱： cloudbbs@huaweicloud.com
UI MATLAB
  • 点赞
  • 收藏
  • 关注作者

评论（0）

0/1000
抱歉，系统识别当前为高风险访问，暂不支持该操作

全部回复

上滑加载中

设置昵称

在此一键设置昵称，即可参与社区互动！

*长度不超过10个汉字或20个英文字符，设置后3个月内不可修改。

*长度不超过10个汉字或20个英文字符，设置后3个月内不可修改。

确认 取消

加入云驻计划，成为创作者

  • 华为云周边好礼
  • 免费体验产品
  • 特殊身份标识
  • 线下官方门票
  • 内部专家零距离
  • 与10000+优质创作者共同成长
立即加入