
2
Projector
Object
Reected light
Human
Laptop
EEG headset
Gray wall
White wall
FIG. 1. The setup used for adaptive ghost imaging. A light
projector illuminates an object cut out from a cardboard sup-
port. Transmitted light is diffused by a ground glass that is in
contact with the cardboard support and illuminates a white,
observation wall. This part of the setup is obscured from
the observer by a wall. The distance of both the object and
the observer from this secondary wall is „0.5´1 m. The
EEG signal from the observer is recorded and processed on a
computer.
signal) and explicit conscious processing (by asking the
participant to either verbally communicate or type on a
keyboard the perceived light intensity).
Imaging protocol. Computational ghost imaging
relies on a light source that can project a series of
typically binary (black and white) patterns, Pn. These
light patterns are then reflected (or transmitted) from
the object or scene we wish to image and collected with
a bucket detector (i.e. sensitive only to total energy),
an. Then, summation of all the bucket value-weighted
patterns will produce an image, O“řanPn. A very
common choice of patterns are the Hadamard set, H,
that can be recursively defined.
Over the years, researchers have optimised ghost
imaging by using different light sources, detectors or
computational algorithms. Recent attempts have also
used the human visual system as a detector where the
visual persistence time of the retina is used to directly
perform the summation operation described above, i.e.
a series of pre-weighted patterns are pre-calculated and
are then visualised at sufficiently high rates that they
are effectively perceived by the eye as an accumulated
sum [39–41]. Conversely, here we implement a form of
computational ghost imaging in which the human visual
cortex processes visual data and also provides feedback
that allows to adapt the projected patterns in real-time
so as to minimise measurement time.
Ghost imaging with the brain. A schematic
overview of the experiments is shown in Fig. 1. We
project a series of binary Hadamard patterns using a
standard digital light projector (DLP) onto an object.
The light transmitted past the object is then observed in
reflection from a secondary white surface (white wall).
Each binary pattern is periodically switched on/off for
several periods with a frame rate that chosen in the
3-30 Hz region. We detect the SSVEP generated by
the visual cortex activity from a single electrode placed
at Oz, the medial visual cortex region (see SI). This
SSVEP is then analysed in the spectral domain and the
corresponding fundamental (i.e. at the same frequency
of the light modulation) and higher harmonic (due to
neuron nonlinearity) amplitudes are extracted. These
are then used to reconstruct an image of the object,
which as shown in the schematic overview, is hidden
behind a wall.
Linearity. The first step for any form of imaging
requires calibration of the detection system and iden-
tification of linear regions or at least regions in which
the system response is monotonic with increasing input
intensity. In this case, the ‘system’ is the visual system
and SSVEP read-out, which is known to exhibit signif-
icant nonlinearity. We characterised the (non)linearity
of the SSVEP readout with a standard LCD screen that
displayed a flickering uniform intensity with frequency
between 3 and 30 Hz and that was varied across the
full 8 bit range of the screen, i.e. in values from 0
to 255, corresponding to completely black (no light)
and very bright (corresponding to 125 Lumens). The
EEG signal is then Fourier transformed [42, 43]. Clear
harmonic peaks are observed as expected [13] and we
then consider the maximum values of the individual
harmonics (up to the fourth) as well as the total sum
of these values (the total SSVEP energy). The SSVEP
energy heatmap for each individual harmonic shows a
complicated and typically non-monotonic dependence
for varying screen intensity and flicker frequency (see
Supplementary Information).
Figure 2a shows the total SSVEP energy. Here we
can identify two ideal flicker frequency regions at 6
and 15 Hz, shown in Fig. 2b. The region around 15
Hz shows a clear monotonic increase of SSVEP energy
with increasing illumination and a similar behaviour
occurs also at 6 Hz, albeit only for a more limited screen
intensity range (between 0 and „125 bits, i.e. between 0
and „75 Lumens). The same calibration measurements
performed across three different people resulted in a
similar behaviour (see Supplementary Information). We
therefore perform most of our experiments at either 15
Hz (using the full 0-125 Lumens intensity range) or 6 Hz
(using a limited intensity range).
Ghost Imaging results. Using the setup shown in
Fig. 1, objects are illuminated with Hadamard patterns
that are each periodically flickered (see SI for full
details).
Figures 3(a) and (b) show results for the standard
ghost imaging approach for a 4 ˆ4 pixel object with a
6 Hz flicker frequency and for 4 s and 2 s illumination
time for each of the first 16 Hadamard patterns. The
columns show the ghost image reconstruction obtained
using each individual harmonic SSVEP energy and then
for the total energy (sum over all harmonics). Only the
total SSVEP energy allows to reconstruct a clear image,