use dft para el procesamiento de imágenes en matlab sin función matlab - matlab, dft

Tengo 2 funciones en matlab

%%%%%%% test_dft
clc;clear all;close all;
%%


img=(imread("images/test_image.jpg"));
if isrgb(img)
a=double(rgb2gray(img));
else
a=double(img);
end
temp = zeros(size(a));
f = double(zeros(size(a)));
[m n]=size(a)
for i = 1:m
temp(:,i)=dft(a(:,i));
end
for j = 1:n
f(j,:) = dft(a(j,:));
end


%%%%%%%% dft.m
function X=dft(x1)
N=length(x1);
X=zeros(size(x1));
i=sqrt(-1);
for k=0:N-1
for n=0:N-1
X(k+1)=(X(k+1)+x1(n+1)*exp(-2*pi*i*k*n/N)) ;
end
end
end

Quiero calcular el dft de la imagen y luego usarlo para enmascarar en el dominio de la frecuencia (sin la función matlab). este código se ocupa de mucho tiempo, pero no funciona. No sé por qué.

Respuestas

imread devuelve valores enteros. Convertirlos al doble primero.

a = double(a);