Optical spatial modulation (OSM) is a multiple-transmitter technique that can provide higher data rates with low system complexity as compared with single-input single-output systems. Orthogonal frequency division multiplexing (OFDM) is widely implemented to achieve better spectral efficiency in wireless channels. Asymmetrically clipped optical OFDM (O-OFDM) and DC-biased O-OFDM are two well-known O-OFDM techniques suitable for intensity-modulation direct-detection optical systems. In this work, sample indexed spatial OFDM (SIS-OFDM) is proposed to combine OSM and O-OFDM in a novel way and achieve significant performance gain. By assigning time-domain samples of the O-OFDM transmit symbol to different transmitters, SIS-OFDM achieves much better spectral efficiency and reduces computational complexity at the transmitter as compared with previous work that combines OSM with O-OFDM in the frequency domain. We also consider the impact of optical source biasing on overall performance, and the relative performance of imaging receiver (ImR) versus non-imaging receiver (NImR) design for our proposed SIS-OFDM technique. Results indicate that for an Ntx × Nrx multiple-input multiple-output configuration where Ntx = Nrx = 4, SIS-OFDM using ImR can achieve up to 135 dB of signal-to-noise ratio gain over comparable system using a NImR. Also, using Nsc number of O-OFDM subcarriers provides up to Nsc × log2(Ntx) additional bits per symbol of spectral efficiency over techniques that combine OSM and O-OFDM in the frequency domain.