Hello, I'm trying to write the Gabor ternary pattern, Please help me with this!

Question

Mehran Kh il 22 Ott 2014

1
Link

Link diretto a questa domanda

https://it.mathworks.com/matlabcentral/answers/159600-hello-i-m-trying-to-write-the-gabor-ternary-pattern-please-help-me-with-this

Modificato: NM il 23 Ago 2015

Risposta accettata: Mehran Kh

Apri in MATLAB Online

Hello, I'm trying to write the GAbor ternary pattern:

The answer should be like below:

Please tell me where did I wrong in code!

MATLAB Code:

clc
close all
clear all
[fn fp]=uigetfile('*.*');
I=(rgb2gray(imread([fp,fn])));
[m n]=size(I);
kv=pi/2;
GTP=zeros(m,n);
for x=1:m
      for y=1:n
          for ii=0:3
              phi=ii*pi/4;
              z=[x y]';
              k=kv*[cos(phi) sin(phi)]';
              psi(x,y)=((norm(k))^2)*exp(-(((norm(k))^2)*((norm(z))^2))/2)*...
                  (exp(sqrt(-1)*k'*z)-exp(-1/2));
              G(x,y)=imag(psi(x,y));
              f(x,y)=conv(G(x,y),I(x,y));
              if f(x,y)<-.03
                  pp=-1;
              else pp=0;
              end
              GTP(x,y)=3^ii*(pp+2*(f(x,y)>.03))+GTP(x,y);
          end
      end
  end
figure,imshow(GTP)

Thanks for your considration

2 Commenti
Mostra NessunoNascondi Nessuno

Geoff Hayes il 22 Ott 2014

Mehran - perhaps describe what problems you are having with your code. Is the problem that the output image is all zeros, or something else?

Mehran Kh il 25 Ott 2014

Yes the output image is zeros but it should be like the Posted Image! Don't know what is wrong with my code!

Accedi per commentare.

Accedi per rispondere a questa domanda.

Answer 1

Mehran Kh il 25 Ott 2014

0
Link

Link diretto a questa risposta

https://it.mathworks.com/matlabcentral/answers/159600-hello-i-m-trying-to-write-the-gabor-ternary-pattern-please-help-me-with-this#answer_156575

Modificato: Mehran Kh il 25 Ott 2014

Apri in MATLAB Online

I found why I get Zeros in output image!Convolution (f and GTP) should be out of loop!

clc
close all
clear all
[fn fp]=uigetfile('*.*');
I=im2double(rgb2gray(imread([fp,fn])));
I=zscore(I);
[m n]=size(I);
kv=pi/2;
t=0.03;
GTP=zeros(m,n);
for x=1:m
    for y=1:n
              z=[x y]';
              k0=kv*[cos(0) sin(0)]';
              psi0(x,y)=norm(k0)^2*exp(-norm(k0)^2*norm(z)^2/2)*...
                  (exp(sqrt(-1)*k0'*z)-exp(-1/2));
               k1=kv*[cos(pi/4) sin(pi/4)]';
              psi1(x,y)=norm(k1)^2*exp(-norm(k1)^2*norm(z)^2/2)*...
                  (exp(sqrt(-1)*k1'*z)-exp(-1/2));
              k2=kv*[cos(pi/2) sin(pi/2)]';
              psi2(x,y)=norm(k2)^2*exp(-norm(k2)^2*norm(z)^2/2)*...
                  (exp(sqrt(-1)*k2'*z)-exp(-1/2));
              k3=kv*[cos(3*pi/4) sin(3*pi/4)]';
              psi3(x,y)=norm(k3)^2*exp(-norm(k3)^2*norm(z)^2/2)*...
                  (exp(sqrt(-1)*k3'*z)-exp(-1/2));
          end
  end
G0=imag(psi0);
G1=imag(psi1);
G2=imag(psi2);
G3=imag(psi3);
f0=conv2(G0,I);f0=f0(1:m,1:n);
f1=conv2(G1,I);f1=f1(1:m,1:n);
f2=conv2(G2,I);f2=f2(1:m,1:n);
f3=conv2(G3,I);f3=f3(1:m,1:n);
GTP0=3^0*(-(f0<-t)+2*(f0>t));
GTP1=3^1*(-(f1<-t)+2*(f1>t));
GTP2=3^2*(-(f2<-t)+2*(f2>t));
GTP3=3^3*(-(f3<-t)+2*(f3>t));
figure,
subplot(2,2,1),imshow(GTP0)
subplot(2,2,2),imshow(GTP1)
subplot(2,2,3),imshow(GTP2)
subplot(2,2,4),imshow(GTP3)
GTP=GTP0+GTP1+GTP2+GTP3;
figure, imshow(GTP)