Mengubah Spasial Belah ketupat Menjadi Frekuensi di Matlab

script berikut memanggil fungsi yang telah di buat di postingan sebelumnya

clear all;
clc;
[x,y] = meshgrid(1:256, 1:256)
b = (x+y < 329) & (x+y > 182) & (x-y > -67) & (x-y <73)
imshow(b)
bf = fftshift(fft2(b));
figure, fftshow(bf)

hasil output :