图像处理基于matlab GUI算子+滤波器+数米粒含Matlab源码 007期

Posted 2021-05-23 唐门佛怒唐莲

一、简介

基于matlab GUI算子+滤波器+数米粒

二、源代码

function varargout = mygui(varargin)


% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
global im;
gui_State = struct('gui_Name',       mfilename, ...
                   'gui_Singleton',  gui_Singleton, ...
                   'gui_OpeningFcn', @mygui_OpeningFcn, ...
                   'gui_OutputFcn',  @mygui_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 mygui is made visible.
function mygui_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 mygui (see VARARGIN)

% Choose default command line output for mygui
handles.output = hObject;

% Update handles structure
guidata(hObject, handles);

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


% --- Outputs from this function are returned to the command line.
function varargout = mygui_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;


% --- Executes on button press in open_pushbutton.
function open_pushbutton_Callback(hObject, eventdata, handles)
% hObject    handle to open_pushbutton (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
global im;
[filename,pathname]=...
uigetfile({'*.*';'*.bmp';'*.tif';'*.png';'*.gif';'*.jpg'},'select picture');  %选择图片路径
str = [pathname filename];  %合成路径+文件名
im = imread(str);   %读取图片
axes(handles.axes1);  %使用第一个axes
imshow(im);  %显示图片
title('原图像');



% --- Executes on button press in blpf_pushbutton.
function blpf_pushbutton_Callback(hObject, eventdata, handles)
% hObject    handle to blpf_pushbutton (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
global im;
button_state = get(hObject,'Value');
if button_state == get(hObject,'Max')

        I0 = im;
        % I0 = imread('beauty.jpg');  
        % subplot(2,2,1),imshow(I0);title('原图');
        str={'半径D0大小：','阶数n：'};                             %输入相关参数：高斯模板大小、方差、二值化阈值
        sData=InputDlg(str,'输入参数',1);
        
        if ~isempty(sData)
            d0=str2double(sData(1));%模板大小
            n=str2double(sData(2));%次数
        end     
       
        I1 = imnoise(I0,'gaussian');         %对原图像加噪声
%         subplot(2,2,2),imshow(I1);title('加入噪声后图像')
         %将灰度图像的二维不连续Fourier 变换的零频率成分移到频谱的中心
        s=fftshift(fft2(I1));
        [M,N]=size(s);                      %分别返回s的行数到M中，列数到N中
        n1=floor(M/2);                       %对M/2进行取整
        n2=floor(N/2);                       %对N/2进行取整
        %IHPF滤波
%         d0=200;                               %初始化d0
        for i=1:M
            for j=1:N
                  d=sqrt((i-n1)^2+(j-n2)^2);        %点（i,j）到傅立叶变换中心的距离
                   h(i,j)=1/((1+(d/d0))^(2*n));                %butterrworth lowpass滤波函数,取n=4
                  s(i,j)=h(i,j)*s(i,j);                  %butterworth低通滤波后的频域表示
            end
        end
        s=ifftshift(s);                       %对s进行反FFT移动
        %对s进行二维反离散的Fourier变换后，取复数的实部转化为无符号8位整数
        s=uint8(real(ifft2(s)));  
        % subplot(2,2,3),imshow(h);title('传递函数');         %butterworth低通滤波的传递函数                           
%         subplot(2,2,3),imshow(s); title('butterworth低通滤波（D0=200，n=4)');  %显示butterworth高通滤波后的图像
        axes(handles.axes2);  %使用第一个axes
        imshow(s);
        title('butterworth低通滤波');
end


% --- Executes on button press in canny_pushbutton.
function canny_pushbutton_Callback(hObject, eventdata, handles)
% hObject    handle to canny_pushbutton (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
global im;
button_state = get(hObject,'Value');
if button_state == get(hObject,'Max')
        I = im;
        I1=edge(I,'canny');
%         subplot(2,3,2),imshow(I1),title('canny算子');
        axes(handles.axes2);  %使用第一个axes
        imshow(I1);
        title('canny算子演示');
end
% 
% I=imread('diaohua.jpg');
% %第一步，用3x3高斯滤波器进行滤波，消除噪声；
% k1=filter2(fspecial('gaussian',3),I); 
% % k=rgb2gray(I);
% %第二步，针对每一个像素，计算横向与纵向两方向的微分近似，以得到像素的梯度大小和方向；
% [m,n]=size(k1);
% for i=2:m-1
%     for j=2:n-1
%         zz1=k1(m,n-1)+k1(m+1,n-1);
%         zz2=k1(m,n)+k1(m+1,n);
%         zz3=k1(m,n-1)+k1(m,n);
%         zz4=k1(m+1,n-1)+k1(m+1,n);
%         kp(m,n)=0.5*(zz2-zz1);
%         kq(m,n)=0.5*(zz3-zz4);
%         kfu(m,n)=sqrt((kp(m,n)^2)+(kq(m,n)^2));%梯度大小
%         angle(m,n)=atan(kq(m,n)/(kp(m,n)+0.001));%梯度方向
%     end
% end
% %第三步，对梯度进行"非极大抑制"
% for i=2:m-1
%     for j=2:n-1
%            if angle(i,j)>=3/8*pi
%             angle(m,n)=2;
%         else if angle(i,j)>=1/8*pi
%             angle(m,n)=1;
%         else if  angle(i,j)>=-1/8*pi
%              angle(i,j)=0;
%         else if  angle(m,n)>=-3/8*pi
%              angle(i,j)=3;
%             else
%                 angle(i,j)=2;
%             end
%             end
%             end
%         end
%     end
% end
% %判断该点是否是8邻域的局部最大值（梯度方向），比如，梯度方向为左右方向的点，
% %判断其是否比左右两点的值来的大，如果不是，使该点的值为0.
% %按照各个方向分别判断
% k2=k1;
% for i=2:m-1
%     for j=2:n-1
%         if angle(i,j)==0
%            if k1(i,j)>k1(i,j-1)&&k1(ij)>k1(i,j+1);
%            else k2(i,j)=0;
%            end
%        end
%        if angle(i,j)==1
%            if k1(i,j)>k1(i+1,j-1)&&k1(i,j)>k1(i-1,j+1);
%            else k2(i,j)=0;
%            end
%        end
%        if angle(i,j)==2
%            if k1(i,j)>k1(i-1,j)&&k1(i,j)>k1(i+1,j);
%            else k2(m,n)=0;
%            end
%        end
%        if angle(i,j)==3
%            if k1(i,j)>k1(i-1,j-1)&&k1(i,j)>k1(i+1,j+1);
%            else k2(i,j)=0;
%            end
%        end
%     end
% end
% 
% %第四步，对两次梯度取阈值，t2=2*t1，若f<t1,则f=0
% %两次阈值分割
% k3=k2;%以t1为阈值分割后的矩阵
% k4=k2;%以t2为阈值分割后的矩阵
% t1=50;
% t2=2*t1;
% for i=2:m-1
%     for j=2:n-1
%        if kfu(i,j)<t1
%            k3(i,j)=0;
%        end
%        if kfu(i,j)<t2
%            k4(i,j)=0;
%        end
%     end
% end

三、运行结果

四、备注

版本：2014a
完整代码或代写添加QQ 1564658423