© 2017 SPIE. We present a single pixel frequency division multiplexing imaging system with two metamaterial spatial light modulators (SLMs) for THz light field imaging. One SLM is used for slicing/modulating the 4D light field from various sub-apertures, while the second one together with a single pixel detector to implement 2D multiplexing measurement. This system is in essential a programming aperture in the frequency domain of carrier signals in our categorization of the traditional light field acquisition systems. We propose a prototype adaption design with available elements for THz imaging. Besides, we present the frequency selection method and reconstruction algorithm for square wave modulation/carrier signals for at most 4096-voxel light field dataset.