1  import HoloProcessing: load_holo, plan_fft, reconst, FFTW, ldr, smpi
 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(similar(lo, size(lo) .÷ N))
 9  Dx, Dy = 50, 100
10  scale = 1200
11  
12  ldr_img = ldr(lo, N, (Dx, Dy), P, scale, shift=true)
13  
14  println(smpi(noised=lo, filtered=ldr_img))