www.gusucode.com > IPCV_Eval_Kit_R2019a_0ce6858工具箱源码程序matlab > IPCV_Eval_Kit_R2019a_0ce6858/code/demo_files/I2_07_2_3a_gaborFilt.m

    %% 僈儃乕儖僼傿儖僞
clc;clear;close all;imtool close all;

%% 攇挿:5pix/cycle丄妏搙:0搙 (墶)偺僈儃乕儖僼傿儖僞偺僆僽僕僃僋僩傪惗惉\帵
g1 = gabor(5, 0)
figure; surf(real(g1.SpatialKernel));shading interp; xlabel('X'); ylabel('Y');
             axis ij; box on; colorbar; title('攇挿:5pix/cycle丄妏搙:0搙');
             
%% 夋憸偺撉崬傒\帵
G = imread('testpat1.png');
figure; imshow(G);

%% 夋憸偵僈儃乕儖僼傿儖僞傪揔墳
mag1 = imgaborfilt(G, g1);
figure; imshow(mag1, []);

%% 攇挿:5,10pix/cycle丄妏搙:0,45搙 偺僈儃乕儖僼傿儖僞僶儞僋偺嶌惉\帵
%      攇挿偵傛傝丄僇乕僱儖偺僒僀僘偑曄壔
g2 = gabor([5,10], [0,45])

sizeMax = size(g2(end).SpatialKernel, 1);
figure;
for p = 1:4
    subplot(2,2,p);
    size1 = size(g2(p).SpatialKernel,1);
    a = ones(sizeMax)*0.9;
    a(1:size1, 1:size1) = real(g2(p).SpatialKernel);
    imshow(a, []);
    title(sprintf('攇挿= %d, 妏搙 = %d, Kernel size = %d', ...
            g2(p).Wavelength, g2(p).Orientation, size(g2(p).SpatialKernel,1)) );
end

%% 夋憸偵僈儃乕儖僼傿儖僞僶儞僋傪揔墳\帵
mag2 = imgaborfilt(G, g2);

figure;
for p = 1:4
    subplot(2,2,p)
    imshow(mag2(:,:,p),[]);
    title(sprintf('攇挿= %d, 妏搙 = %d, Kernel size = %d', ...
            g2(p).Wavelength, g2(p).Orientation, size(g2(p).SpatialKernel,1)) );
end

%%
% Copyright 2014 The MathWorks, Inc.