Reposity for Course 图形绘制技术
As dct can't do perfect FFT because of downsampling, the frequency image looks a bit wierd. As a compromise, I write simulate functions to simulate the real frequency images.
Due to the algorithm it runs very slowly while calculating isotropic situation...