1  import HoloProcessing: load_holo, plan_fft, reconst, FFTW, sdm
 2  import Makie: save
 3  
 4  lo = load_holo("set", "111.jpg"; convert=true)
 5  FFTW.set_num_threads(Sys.CPU_THREADS)
 6  
 7  N = 2
 8  P = plan_fft(lo)
 9  Nx, Ny = size(lo) .÷ N
10  Dx, Dy = 50, 100
11  scale = 600
12  
13  sdm_img = sdm(lo, (Nx, Ny), (Dx, Dy), P, scale)
14  
15  save("out_sdm.png", sdm_img)