clear all close all S = load('matlab.mat'); b= cell2mat(struct2cell(S)); [l m n]=size(b); mm =mean(b,3); figure(1) subplot(1,3,1) imshow(mm,'DisplayRange',[]);title('Original Image'); for i=1:n spe1(:,:,i)= fftshift(fftn(b(:,:,i)-mm)); spe(:,:,i)=(abs(spe1(:,:,i)).^2); end w = (mean(spe,3)); subplot(1,3,2) imshow(w,'DisplayRange',[]); title('DFT'); %Radial average [XX, YY] = meshgrid( (1:m)-m/2, (1:m)-m/2); R = sqrt(XX.^2 + YY.^2); profile = []; for i = 1:1:(m+1) mask = (i-1