(1) To calculate the mean
Ī(x,y) of each pixel in the pre-stimulus baseline recording (
n frames):
(2) To calculate the standard deviation
σ(
x,
y) of each pixel in the pre-stimulus baseline recording (
n frames):
(3) To conduct spatiotemporal filtering of potential noises: Because blood flow changes dynamically, the variability of light intensity at the blood vessels is much larger than that at the blood-free area; that is, before the stimulus, the temporal
σ(
x,
y) of blood flow is much larger than that of photoreceptors. Upon stimulation, blood flow may increase, but within the short recording time (1 second), hemodynamic change is slow. Therefore, the temporal
σ(
x,
y) change of blood flow is insignificant compared with the fast IOSs from the photoreceptors. To reject noise attributable to blood flow, values three standard deviations above or below the mean at each pixel were used as a filtering criterion. This filter (3-
σ) allowed us to plot the vasculature profile as shown in
Figure 2B. In other words, the pixel change will be assumed to reflect noise, if
Therefore, a high threshold is used to define stimulus-evoked IOS in the retinal area superimposed by blood vessels. The signals at pixel (
x,
y) with light intensity greater than the mean above three standard deviations are positive and those with light intensity less than the mean below three standard deviations are negative. IOS images with pixels that fall into the noise range are forced to be zero, and only positive or negative IOSs are left. Therefore, after dynamic spatiotemporal filtering (
Fig. 3C), most hemodynamic-driven optical signals (
Fig. 3B) can be rejected.