Spatial statistics for gaze patterns in scene viewing: Effects of repeated viewing

TL;DR

This study applies spatial point process methods, specifically the pair correlation function, to analyze eye movement patterns during scene viewing, revealing increased fixation aggregation upon repeated image viewing and suggesting attentional span as a key factor.

Contribution

It introduces the use of the pair correlation function to study fixation dependencies and demonstrates how repeated viewing affects fixation aggregation in scene viewing.

Findings

Fixation locations are more aggregated within four degrees than chance.

Aggregation of fixations increases upon second viewing of the same image.

A narrower attentional span may explain differences in fixation patterns between viewings.

Abstract

Scene viewing is used to study attentional selection in complex but still controlled environments. One of the main observations on eye movements during scene viewing is the inhomogeneous distribution of fixation locations: While some parts of an image are fixated by almost all observers and are inspected repeatedly by the same observer, other image parts remain unfixated by observers even after long exploration intervals. Here, we apply spatial point process methods to investigate the relationship between pairs of fixations. More precisely, we use the pair correlation function (PCF), a powerful statistical tool, to evaluate dependencies between fixation locations along individual scanpaths. We demonstrate that aggregation of fixation locations within four degrees is stronger than expected from chance. Furthermore, the PCF reveals stronger aggregation of fixations when the same image is presented a second time. We use simulations of a dynamical model to show that a narrower spatial attentional span may explain differences in pair correlations between the first and the second inspection of the same image.

Tables

Table 1. Table 1: Model parameters. Standard errors were calculated from five parameter estimations of the experimental data recorded during the first presentation of natural images.

		Natural Images		Texture Images
	Parameter	First	Second	First	Second	Error
Attention span	${\sigma }_1$	4.88	4.20	4.72	4.20	$\pm$0.11
Attention decay	${\log }_{10}\rho$	-1.18	-1.18	-1.36	-1.36	$\pm$0.28
Inhibition span	${\sigma }_0$	2.16	2.18	1.79	1.76	$\pm$0.25
Inhibition decay	${\log }_{10}\omega$	-4.03	-4.03	-3.80	-3.80	$\pm$0.08
Noise	${\log }_{10}\eta$	-4.04	-4.04	-4.51	-4.51	$\pm$0.07

Table 2. Table 2: Fixed effects of linear mixed effect models. For each point process we estimated separate models for saccade amplitudes ( log -transformed) and the summed deviation from complete spatial randomness of the PCF ( log -transformed). The table reports estimates of fixed effects ( β 𝛽 \beta ) with standard errors (SE) and t 𝑡 t values. | t | > 2 𝑡 2 |t|>2 are interpreted as significant effects.

		Saccade amplitudes				PCF deviation
		$\beta$	SE	$t$		$\beta$	SE	$t$
Experiment
(Intercept)		1.599	0.034	46.58		5.668	0.095	59.43
Presentation		-0.062	0.012	-5.33		0.508	0.099	5.13
Image Type		0.035	0.023	1.52		-0.073	0.123	-0.59
Presentation $\times$ Image Type		0.024	0.016	1.51		-0.144	0.187	-0.77
SceneWalk
(Intercept)		1.459	0.006	232.72		5.519	0.037	147.80
Presentation		-0.090	0.006	-16.18		0.586	0.062	9.52
Image Type		-0.010	0.013	-0.83		-0.135	0.075	-1.81
Presentation $\times$ Image Type		0.058	0.011	5.25		-0.211	0.123	-1.71
Joint Probability Control
(Intercept)		1.805	0.011	165.32		5.123	0.030	168.32
Presentation		-0.033	0.005	-6.39		0.219	0.055	4.00
Image Type		0.040	0.022	1.81		0.168	0.061	2.75
Presentation $\times$ Image Type		0.010	0.010	0.94		-0.114	0.110	-1.04

Equations

$$\lambda (x)=\underset{|dx|\to 0}{\lim }\left\{\frac{E[N(dx)]}{|dx|}\right\}.$$

$$\rho (x,y)=\underset{|dx|,|dy|\to 0}{\lim }\left\{\frac{E[N(dx)N(dy)]}{|dx||dy|}\right\}.$$

$$\hat{\rho }(r)=\sum _{x,y\in W}^{\ne }\frac{k(\Vert x-y\Vert -r)}{2\pi r{A}_{\Vert x-y\Vert }}$$

$$g_{hom}(r)=\frac{1}{{\lambda }^2}\rho (r).$$

$$g_{inhom}(r)=\sum _{x,y\in W}^{\ne }\frac{1}{\lambda (x)\lambda (y)}\rho (r)$$

$$\chi ={\int }_0^{\infty }(g(r)-1{)}^{2}dr.$$

$$a_{ij}(t+1)=\frac{{\varphi }_{ij}{A}_{ij}(t)}{\sum _{kl}{\varphi }_{kl}{A}_{kl}(t)}+(1-\rho )a_{ij}(t)$$

$$f_{ij}(t+1)=F_{ij}(t)+(1-\omega )f_{ij}(t)$$

$$u_{ij}(t)=-\frac{[a_{ij}(t){]}^{\lambda }}{\sum _{kl}[a_{kl}(t){]}^{\lambda }}+\frac{[f_{ij}(t){]}^{\gamma }}{\sum _{kl}[f_{kl}(t){]}^{\gamma }}$$

$${\pi }_{ij}(t)=\max \left(\frac{u_{ij}(t)}{\sum _{(k,l)\in S}{u}_{kl}(t)},\eta \right)$$

$$\mbox{{log(dv) $\sim$ }}\underbrace{\mbox{{1 + A + B + A:B }}}_{\mbox{fixed effects}}\underbrace{\mbox{{+ ( 1 + A + B + A:B | id )}}}_{\mbox{random effects id}}\underbrace{\mbox{{+ ( 1 + A | img )}}}_{\mbox{random effects image}}$$

Full text

**Spatial statistics for gaze patterns in

scene viewing: Effects of repeated viewing 

**

Hans A. Trukenbrod1 & Simon Barthelmé2 &

Felix A. Wichmann3,4,5 & Ralf Engbert1

1University of Potsdam, Germany

2CNRS, Gipsa-lab, Grenoble INP, France

3Eberhard Karls University of Tübingen, Germany

4Bernstein Center for Computational Neuroscience Tübingen, Germany

5Max Planck Institute for Intelligent Systems, Tübingen, Germany

(Version: )

Address correspondence to:

Hans Trukenbrod

Department of Psychology

University of Potsdam

Am Neuen Palais 10

D–14469 Potsdam

Germany

E-mail: Hans.Trukenbrod@uni-potsdam.de

Phone: +49 331 9772703, Fax: +49 331 9772794

**Scene viewing is used to study attentional selection in complex but still controlled environments. One of the main observations on eye movements during scene viewing is the inhomogeneous distribution of fixation locations: While some parts of an image are fixated by almost all observers and are inspected repeatedly by the same observer, other image parts remain unfixated by observers even after long exploration intervals. Here, we apply spatial point process methods to investigate the relationship between pairs of fixations. More precisely, we use the pair correlation function (PCF), a powerful statistical tool, to evaluate dependencies between fixation locations along individual scanpaths. We demonstrate that aggregation of fixation locations within four degrees is stronger than expected from chance. Furthermore, the PCF reveals stronger aggregation of fixations when the same image is presented a second time. We use simulations of a dynamical model to show that a narrower spatial attentional span may explain differences in pair correlations between the first and the second inspection of the same image. **

Introduction

When we move our eyes during scene viewing, fixation locations are not selected uniformly.  Instead fixations cluster in specific areas while other areas remain unfixated by observers, even after long exploration intervals. Most research on this inhomogeneity has tried to identify factors that contribute to the placement of fixations across trials, while statistical correlations within trials have mostly been ignored. Here we describe the pair correlation function (PCF), a method from spatial statistics, to investigate the relationship between fixation locations of individual scanpaths. In a first part, we provide a step-by-step tutorial how the PCF can be applied to eye movement data and demonstrate  that fixations in a scanpath are more aggregated than expected by the distribution of fixation locations of all subjects Engbert \BOthers. (\APACyear2015).  In a second part, we show that a long-term memory manipulation, i.e. the second inspection of an image, leads to even more aggregation. We discuss the results in the light of simulations of the SceneWalk model Engbert \BOthers. (\APACyear2015), a dynamical model for the generation of saccadic sequences in scene viewing. Furthermore, we demonstrate with simulated fixation locations that the PCF can be used to test specific hypotheses and that seemingly similar distributions of fixation locations, may lead to very different PCFs.

Eye movements during scene perception.

Fixation locations during scene perception are influenced by bottom-up and top-down as well as low-level and high-level factors <see¿[for an extensive discussion]Schutt.arXiv.2018. Bottom-up factors refer to parts of the image that attract gaze independent of the internal state of an observer and might differ in complexity. Simple low-level features are extracted early in the visual hierarchy while high-level features are complex shapes and extracted late in the visual hierarchy. Examples of bottom-up, low-level features that predict fixation locations are luminance contrast and edge density Mannan \BOthers. (\APACyear1997); Reinagel \BBA Zador (\APACyear1999); Tatler \BOthers. (\APACyear2005). The strength of the relationship for different image features depends on the type of image viewed Parkhurst \BOthers. (\APACyear2002). Examples of bottom-up, high-level features that predict fixation locations are objects <e.g., faces, persons, cars;¿Einhauser.JVis.2008,Cerf.AdvNeuralInfoProcSyst.2007,Judd.CompVision.2009. This interpretation is supported by the existence of a preferred viewing location close to the object center in scene viewing Nuthmann \BBA Henderson (\APACyear2010).

The influence of bottom-up factors has led to the development of computational saliency models (\citeNPItti.NatRevNeurosci.2001; cf., \citeNPKoch.HumNeurobiol.1985) and a variety of models has been put forward over the years <e.g.,¿Bruce.JVis.2009,Kienzle.JVis.2009,Kummerer.CoRR.2016,Parkhurst.VisionRes.2002,Vig.CVPR.2014,Zhang.JVis.2008. In particular since the rise of sophisticated machine-learning algorithms, these models perform well in predicting fixation locations when evaluated with a data set obtained under unconstrained (“free”) viewing <¿Bylinskii.SaliencyBenchmark.

Top-down factors refer to cognitive influences on fixation locations and depend on the internal state of an observer (e.g., aims, memory load, knowledge). As for bottom-up factors, top-down factors differ in respect to the complexity of features. Examples of top-down, low-level factors can be found in visual search, where observers search for a specific color or a line of a specific orientation. Examples of top-down, high-level factors are task instructions where observers need to judge the age of people or their wealth. In his seminal work, \citeAYarbus.Book.1967 reported anecdotal evidence for top-down control. Scanpaths were influenced by the instruction given before viewing an image <see also¿Castelhano.JVis.2009. These top-down effects strengthen when observers are engaged in natural tasks like preparation of a sandwich or during tea-making Hayhoe \BBA Ballard (\APACyear2005); Land \BBA Hayhoe (\APACyear2001). In natural tasks, fixations generally support the smooth execution of a task and occur on objects just-in-time Ballard \BOthers. (\APACyear1997) or as look-ahead fixations to inspect objects needed later during the task Pelz \BBA Canosa (\APACyear2001). Thus, the eyes do not necessarily fixate the most salient location (bottom-up) but are rather directed towards informative locations that are important for task execution. Another source of top-down control comes from memory representations due to reinspection of previously seen images Kaspar \BBA König (\APACyear2011\APACexlab\BCnt1, \APACyear2011\APACexlab\BCnt2) and due to the acquired scene or world knowledge Torralba \BOthers. (\APACyear2006). Incorporating such knowledge by contextual priors improves the predictions of saliency models Torralba \BOthers. (\APACyear2006); Judd \BOthers. (\APACyear2009).

In addition, fixation locations depend on systematic tendencies that are common to eye movements in general Tatler \BBA Vincent (\APACyear2008). These systematic tendencies lead to spatial and directional biases in the selection of fixation locations.  A well-known example is the central fixation bias during scene viewing Tatler (\APACyear2007). On average, participants prefer to fixate near the center of an image rather than towards the periphery. The central fixation bias is strongest in the beginning of a trial and reaches an asymptotic level after a few fixations Rothkegel \BOthers. (\APACyear2017). Other systematic tendencies include the preference to generate positively skewed, long-tailed distributions of saccade amplitudes or the preference to execute saccades in the cardinal directions during scene viewing Tatler \BBA Vincent (\APACyear2008).  The later effect is shaped by image features and varies systematically with the perceived horizon. Tilting an image results in an equally tilted distribution of saccade directions Foulsham \BOthers. (\APACyear2008); Foulsham \BBA Kingstone (\APACyear2010). In addition, observers tend to make saccades in the same direction as the preceding saccade (cf., saccadic momentum) and a large number of fixations bring the eyes back to the the last or penultimate fixation location <cf., facilitation of return;¿Smith.VisCogn.2009,Wilming.PLoSCompBiol.2013. Finally, inhibition of return is believed to facilitate exploration during visual search Klein \BBA MacInnes (\APACyear1999) and has been suggested as a mechanism to drive attention during scene perception Itti \BBA Koch (\APACyear2001). Inhibition of return during scene perception seems to primarily prolong fixation durations before return saccades to previously fixated locations Hooge \BOthers. (\APACyear2005); Smith \BBA Henderson (\APACyear2009) and seems to shape spatial dynamics of scanpaths in the long run Rothkegel \BOthers. (\APACyear2016). Adding systematic tendencies to saliency models further improves their predictions Le Meur \BBA Liu (\APACyear2015).

Eye movements and long-term memory: Repeated presentation

Humans have a remarkable capacity to store images in long-term memory Brady \BOthers. (\APACyear2008); Konkle \BOthers. (\APACyear2010); Standing \BOthers. (\APACyear1970). These representations are not limited to the gist of a scene but include abstract representations of objects, in particular of previously fixated objects Hollingworth \BOthers. (\APACyear2001); Hollingworth \BBA Henderson (\APACyear2002). By viewing the same image multiple times scene memory accumulates Melcher (\APACyear2001); Melcher \BBA Kowler (\APACyear2001). Due to this scene-specific memories repeated presentations can be used to investigate effects of long-term memory on eye movements.

In two studies, Kaspar \BBA König; Kaspar \BBA König studied the interaction of bottom-up and top-down influences and long-term memories by presenting images five times. In general, fixation locations tend to be similar across different inspections of the same image by the same participant Kaspar \BBA König (\APACyear2011\APACexlab\BCnt2). Similarity was strongest for successive presentations and decreased with increasing distance between presentations. The correlation between fixation locations and low-level features, however, remained rather constant. In contrast, the number of fixated regions decreased after multiple presentations of images, as did the average number of fixations. Furthermore, saccade amplitudes were largest during the first presentation and decreased on subsequent presentations.

In a second study \citeAKaspar.PLoSOne.2011 explored the effects of repeated presentations on top-down influences on saccade selection. Motivation (as measured via the reported interestingness of the viewed images) and a personality trait of participants (action orientation) influenced repeated viewing of images. In addition, fixation durations and variability between participants’ fixation locations increased, whereas saccade frequency, saccade length, and entropy of fixation locations decreased. The authors concluded that the locus of attention became increasingly local with repeated presentations of images. Thus, participants scrutinized individual regions during later presentations. This interpretation was supported by participants’ self-reports and was augmented in participants who found the images more interesting.

Research Questions

Much progress has been made to understand where observers fixate in an image. This research primarily focused on fixation locations across observers while neglecting the fixation history during a trial of a single observer. Fixation locations, however, exhibit strong spatial correlations during a trial and are not independent of one another Engbert \BOthers. (\APACyear2015) and adding mechanisms that generate more realistic scanpaths improves performance of saliency models <¿LeMeur.VisionRes.2015. \citeABarthelme.JVis.2013 introduced spatial point processes as a theoretical framework for the study of gaze patterns, and demonstrated how this helps to turn qualitative into quantitative questions. Here, we present a method from spatial statistics, i.e., the pair correlation function (PCF), to estimate spatial correlations between fixation locations during a trial in the presence of spatial inhomogeneity Engbert \BOthers. (\APACyear2015). Before applying the PCF to eye movement data during the first and second inspection of an image, we briefly describe the theoretical details of the PCF. We demonstrate (i) that the PCF provides rigorous statistical evidence for aggregation of fixation locations in single trials, (ii) that this effect is not well explained by the tendency of participants to generate short saccade amplitudes,  (iii) that the PCF is differentially affected by a memory manipulation (first vs. second inspection of an image), and (iv) that these differences can be explained by modulations of the attentional span within our SceneWalk model Engbert \BOthers. (\APACyear2015).

Pair Correlation Function

We refer the reader to \citeADiggle.Book.2013 for an introduction to the statistical analysis of point patterns and to \citeALaw.JEcol.2009 for a detailed description of the PCF and the application to point patterns in plant ecology.

Analyzing spatial point patterns: Pair correlation function

The density estimation of a point pattern, i.e., the probability of observing a point at a given location, is a first-order statistic for a spatial point process, which, therefore, plays the role of the mean value in classical statistics. In the upcoming sections we denote this first-order statistic as the intensity $\lambda(x)$. In the case of eye movements the intensity represents the local average spatial density of fixations at a location $x$. Point patterns that are generated by a homogeneous point process are uniformly distributed and the underlying intensity $\lambda(x)=\lambda$ is constant for all $x$. For inhomogeneous point processes, where the 2D density of fixation locations is non-uniformly distributed, the intensity $\lambda(x)$ is estimated for each location $x$ separately, i.e.,

[TABLE]

Here, the density is given by the expected number of fixations $E[N]$  falling into a disc of infinitesimal size $|dx|$.

While first-order statistics are concerned with locations of single points (while ignoring spatial correlations between points), second-order statistics describe the relation between pairs of points. A crucial second-order statistic for the computation of PCFs is the pair density $\rho(r)$. The pair density (or second-order intensity function) describes the probability of simultaneously observing points generated by a point process  in two disjoint discs with centers $x$ and $y$ of infinitesimal size $dx$ and $dy$,

[TABLE]

For a stationary, isotropic point process the pair density $\rho(x,y)=\rho(r)$ depends on the distances of pairs of points only, where $r$ corresponds to the distance between pairs of points $\|x-y\|$ Diggle (\APACyear2013). Mathematically we estimate the pair density $\hat{\rho}(r)$ at distance $r$ by

[TABLE]

where $k$ represents an Epanechnikov kernel111The Epanechikov kernel $\epsilon(x)=\frac{3}{4w}(1-\frac{x^{2}}{w^{2}})_{+}$ with $(x)=max(0,x)$ is a quadratic function that is truncated to the interval [-w,w]. Baddeley \BOthers. (\APACyear2015)  and $A_{\|x-y\|}$ is an edge-correction factor to counter the loss of pairs of points near the boundary of the inspection window, which is particularly important at large distances $\|x-y\|$. As edge correction we chose to use the translation correction, that weights each pair of points $(x_{i},x_{j}$) with the reciprocal of the fraction of the window area, in which a first point $x_{i}$ could be placed, so that both points $x_{i},x_{j}$ would be observable Baddeley \BOthers. (\APACyear2015)

Since the pair density depends on the number of points, the PCF is computed as a normalized version of the pair density (i.e., the PCF is given by the intensity-weighted pair density). In the case of a homogeneous PCF the pair density is weighted by a constant at all distances,

[TABLE]

However, fixation locations are not distributed homogeneously and show aggregation due to bottom-up, top-down, and systematic oculomotor factors. Fortunately, the PCF can be computed for inhomogeneous processes by weighting the pair density with the intensities $\lambda(x_{i})$ of an inhomogeneous point process at each fixation location,

[TABLE]

The resulting PCFs will be non-negative, $g(r)\geq 0$, at all distances $r$. Values of the PCF close to one, $g(r)\approx 1$, indicate that pairs of points at distance $r$ are independent. Points at distance $r$ occur solely due to the underlying intensity $\lambda(x)$. For larger values, i.e., $g(r)>1$, point patterns are more abundant at distance $r$ than expected by the intensity $\lambda(x)$. Thus, pairs of points at distance $r$ interact and observing a point $x$ increases the probability of observing a point $y$ at distance $r$. The probability of observing point $y$ is higher  than predicted by the local intensity $\lambda(y)$. Conversely, smaller values, $g(r)<1$, reveal that points are less abundant than the spatial average at distance $r$. Observing a point reduces the probability of observing a second point at distance $r$.

Figure 1 shows the PCF of three different point patterns. In all examples we computed the PCF assuming a homogeneous point process with constant intensity $\lambda(x)=\lambda$ (cf., Eq. 4) since deviations from uniformity are easier to interpret visually. The same interpretation, however, can be applied to inhomogeneous PCFs. The first example shows a regular point pattern (left column). Visual inspection of the points indicate a grid-like arrangement. The distance between neighboring points is relatively constant. The resulting PCF (bottom row) summarizes this behavior. At short distances, $r<4$, the PCF reveals a strong inhibitory effect, ${g}(r)\approx 0$. The existence of a point impedes the occurrence of other points within this radius. At medium distances, $4<r<6$, the PCF reveals aggregation of points, ${g}(r)>1$. Observing a point boosts the occurrence of points at this distance. Hence, the grid-like appearance. At larger distances, $r>6$, the PCF lends support to the hypothesis of independence of points, since ${g}(r)\approx 1$. We observe no long-range interaction of pairs of points and the distribution of distances can be explained by the density distribution, $\lambda(x)$.

The second example (Fig. 1, central panels) shows the realization of a point pattern with complete spatial randomness (CSR) . Points are distributed uniformly. The resulting PCF reveals the independence of points at all distances, ${g}(r)\approx 1$. Note, the aggregation at short distances $r$ is an artifact generated by the estimation process. Finally, the third example illustrates an aggregated point process (right panels). The PCF at short distances, $r<2$, reveals aggregation, ${g}(r)>1$, while the PCF at longer distances reveals independence, ${g}(r)\approx 1$. Thus, observing a point increases the likelihood of observing other points in close proximity. The occurrence of distant points can be explained by the uniform distribution.

Finally, deviations from complete spatial randomness , i.e., $g(r)\neq 1$, can be summed up and serve as a useful summary statistic of the overall behavior:

[TABLE]

In practical applications, the deviation from complete spatial randomness  $\chi$ is computed over a finite interval, e.g., $\int_{a}^{b}(g(r)-1)^{2}dr$ with $a<b$.

Application of the PCF to fixation locations

In this section we demonstrate how to compute the pair correlation function (PCFs) for eye movement data in three steps.  The PCF reveals whether the distribution of fixation locations during a single scanpath can be explained by the overall inhomogeneity observed across all observers or whether fixation locations of a single scanpath contain additional spatial correlations . All analyses and graphs reported have been implemented in R using the spatstat Baddeley \BBA Turner (\APACyear2005); Baddeley \BOthers. (\APACyear2015) and ggplot2 packages Wickham (\APACyear2009). We provide R-code at \urlhttp://www.rpubs.com/hans/PCF.

0.0.1 1. Simulate inhomogeneous and homogeneous control processes

To evaluate our PCF computations, we simulate two control point processes, namely a homogeneous and an inhomogeneous point process. Points (fixation locations) are sampled independently from each other in both control processes and due to the independence of points, we do not expect to observe any correlations between points at distance $r$. Any observed correlations would be spurious and depend on the data structure (e.g., length of fixation sequences) or a wrong parameterization of the method. Hence, both control processes ensure that correlations in the PCF arise from the empirical data and not by the method itself. In addition, the inhomogeneous point process is used in the second step to estimate an optimal bandwidth for the intensity estimation of the PCF.

Figure 2 (left panel) shows fixations on an image from participants viewing the same scene for 10 seconds (see Methods for details). As expected fixation locations are not uniformly distributed and indicate inhomogeneity. The estimated intensity $\hat{\lambda}_{s}(x)$ of all fixation locations is depicted by gray shading where darker areas represent higher intensities. We used Scott’s rule of thumb to compute the smoothing bandwidth for the intensity estimation (R-function: bw.scott) from the spatstat package Baddeley \BBA Turner (\APACyear2005) and  estimated an optimal bandwidth for each image. The estimated intensity $\hat{\lambda}_{s}(x)$ is used to simulate an inhomogeneous control point process. The inhomogeneous point process (central panel) samples points proportionally to the intensity $\hat{\lambda}_{s}(x)$. Hence, the resemblance of the experimental and simulated distributions. The homogeneous point process (right panel) samples from a uniform distribution across the entire image area, therefore the subsequently estimated intensity (gray shading) is approximately constant, $\hat{\lambda}_{s}(x)\approx\lambda$. For every empirical scanpath we simulated one scanpath of equal length (same number of fixations) for the inhomogeneous point process and for the homogeneous point process obtained from Monte-Carlo simulations.

Examples of the empirical and simulated scanpaths are visualized in Figure 3. Each row shows matching scanpaths of the empirical, inhomogeneous and homogeneous point process. The estimated overall intensity of each point process on an image is displayed via gray shading. Fixations are likely to be located in areas of high average intensity $\hat{\lambda}(x)$. However, each sequence consists of a unique set of points where some scanpaths explore otherwise ignored locations or “missed” locations of high intensity (Fig. 3). Overall, scanpaths of the inhomogeneous and homogeneous point processes reveal less systematic exploration behavior than the empirical data. Hence, saccade amplitudes increase considerably.

0.0.2 2. Choose optimal bandwidth for intensity estimation of PCF

Next we need to choose an optimal bandwidth for estimation of the intensity $\hat{\lambda}(x)$ used to calculate the inhomogeneous PCF $\hat{g}(r)$ (see Eq. 5). Note, this is different from the bandwidth estimated in Step 1 to simulate the inhomogeneous point process.  Since fixation locations in scanpaths of both control point processes are sampled independent from the preceding fixation history, average PCFs of both point processes are expected to reveal no spatial correlations, i.e. we expect $\hat{g}(r)\approx 1$ at all distances $r$. We computed the deviation from complete spatial randomness  (Eq. 6) for PCFs computed with different bandwidths. We varied bandwidths from $0\degree$ and $10\degree$ in steps of $0.1\degree$ and computed the deviation from complete spatial randomness  for each scanpath. The average deviation at each bandwidth is plotted in Figure 4. Lines represent individual images. For all images, the deviation increases for small bandwidths and large bandwidths with an optimal bandwidth between $1.5\degree$ and $5\degree$. The bandwidth yielding the smallest deviation was chosen for the intensity estimation of the PCF.

0.0.3 3. Compute PCF for each trial

In the last step we compute the PCF of our empirical data. For estimation, we use the intensity $\hat{\lambda}(x)$ that resulted in the smallest deviation from complete spatial randomness  of the PCF $\hat{g}(r)$ for the inhomogeneous point process (see previous step). PCFs of individual scanpaths on an image are displayed in Figure 5 (gray lines). The three example scanpaths from Figure 3 are plotted in black. PCFs vary strongly between individual trials for all point processes. The average empirical PCF across all scanpaths on an image (red line) deviates from complete spatial randomness  $\hat{g}(r)\neq 1$ for distances smaller than $4\degree$. At distances beyond 4° the average PCF suggests independence of points, i.e., $\hat{g}(r)\approx 1$. Thus, fixations co-occur in close proximity during individual trials. Conversely, areas further away are fixated as predicted by chance, i.e., the overall inhomogeneity $\hat{\lambda}(x)$ observed across all participants. Inspection of the control point processes demonstrates the absence of spatial correlations. The average PCF of the inhomogeneous and homogeneous point process are constant with $\hat{g}(r)\approx 1$. The artifact of the estimation procedure at short distances is present in all estimates.

The same procedure can be repeated for each image. Figure 6 shows PCFs of each image (gray lines) as well as the average across all images (colored lines). While inhomogeneous and homogeneous point processes reveal no spatial correlations, empirical PCFs show spatial aggregation at short distances, $r<4\degree$ in all conditions. For a more detailed discussion see the Results section.

Methods

For our experiment we had participants view two types of images twice. The repeated presentation of images can be understood as a form of visual long-term memory manipulation. From previous work it can be expected that this leads to similar fixation densities but shortens saccade amplitudes during the second inspection Kaspar \BBA König (\APACyear2011\APACexlab\BCnt1, \APACyear2011\APACexlab\BCnt2). The resulting point patterns are similar but differ slightly in the overall inhomogeneity, which makes a direct comparison of the eye movement behavior difficult. The same problem is true for the comparison of different image types. The PCF takes differences in the underlying inhomogeneity into account and allows a direct comparison of the spatial correlations under different viewing conditions (first vs. second presentation) and for different image types (natural vs. texture images).

Participants

We recorded eye movements of 35 participants (15 male, 20 female) aged 17–36 years (mean: 24.0). Participants received study credits or 8€ for participation and were recruited at the University of Postdam and from a local school (32 students, 3 pupils). All participants had normal or corrected-to-normal vision as assessed by the Freiburg Vision Test <FrACT;¿Bach.OptomVisSci.1996.222\urlhttp://michaelbach.de/fract/ The study conformed to the Declaration of Helsinki as well as national ethics guidelines. We obtained written informed consent from all participants.

Apparatus

Stimuli were presented on a 20" CRT monitor (Mitsubishi Diamond Pro 2070; refresh rate 120 Hz; resolution: 1280$\times$1024 pixels). Eye movements were recorded binocularly using the video-based Eyelink 1000 system (SR Research, Osgoode, ON, Canada) with a sampling rate of 1000 Hz. In order to reduce head movements we asked participants to position their head on a chin-rest in front of the computer screen (viewing distance: 70 cm). Stimulus presentation and response collection were implemented in MATLAB (The MathWorks, Natick, MA, USA) using the Psychophysics Toolbox <PTB–3;¿Brainard.SpatialVision.1997,Kleiner.Perception.2007,Pelli.SpatialVision.1997 and Eyelink Toolbox Cornelissen \BOthers. (\APACyear2002).

Stimulus material

We used two sets of colored photographs in our experiment. In the first part of the experiment each set consisted of 15 images (Fig. 7). Image Set 1 contained photographs of natural landscapes and rural scenes. Image Set 2 contained photographs of textures. During the memory test in the second part of the experiment we added 15 novel images of the same category to each image set.

Task and procedure

Each trial began with the presentation of a black fixation cross in front of a uniform gray background. The fixation cross was placed randomly within the image boundaries. After successful fixation, the fixation cross was replaced by the image for 10 s. Participants were instructed to explore each image for a subsequent memory test. The first block of the experiment consisted of 60 trials where each image was presented twice. Results from the first inspection of images have been described previously Engbert \BOthers. (\APACyear2015). The second inspection has not been published earlier. In a second block participants completed a memory test with 60 trials. The memory test contained all presented images and thirty new images with natural scenes and textures. Participants had to judge whether the image had been presented during Block 1.

To minimize the potential influence of the monitor frame and since accuracy of eye trackers falls off towards the edges of a monitor, images were presented centrally with gray borders extending 32 pixels to the top/bottom and 40 pixels to the left/right of the image. The resulting size of the image was 1200$\times$960 pixels (31.1°$\times$24.9°).

Data preprocessing

We detected saccades by using a velocity-based algorithm Engbert \BBA Kliegl (\APACyear2003); Engbert \BBA Mergenthaler (\APACyear2006). Saccades were defined as fast movements of both eyes that exceeded the average velocity during a trial by 6 standard deviations for at least 6 ms with a minimum amplitude of 0.5∘. Eye traces between successive saccades were tagged as fixations. Fixation positions were computed by averaging the mean eye position of both eyes. Since trials started with a fixation check, first fixations on images were removed from the data set ($N=2100$). In addition, fixations containing a blink or with a blink during an adjacent saccade were excluded from subsequent analyses ($N=2214$). Overall, 55.526 fixations remained for further analyses.

0.1 SceneWalk model

For the interpretation of our results, we simulated fixation sequences with the SceneWalk model (\citeNPEngbert.JVis.2015; cf., \citeNPSchutt.PsycholRev.2017). Fixation sequences are generated in the computational model by two competing activation maps: (i) an excitatory attention map that provides potential saccade targets and (ii) an inhibitory fixation map that tags previously fixated locations. Maps in the model are discretized with dimensions $k\times l=128\times 128$ Stensola \BOthers. (\APACyear2012). Activations in the attention map $a_{ij}$ at coordinates $(i,j)$ evolve over time. To approximate potential saccade targets we compute the empirical density of fixation locations on a given image. The empirical density contains all effects generated by bottom-up and top-down processing that contribute to the inhomogeneous distribution of fixation locations. In our model simulations, the empirical density feeds into the attention map. Extraction from the empirical density is highest at fixation and decreases with increasing eccentricity. This corresponds to the attentional window of our model. Mathematically, the empirical density of fixation locations is weighted by a Gaussian envelope of size $\sigma_{1}$. Position of the attentional window $a_{ij}$ changes after each saccade and remains constant otherwise. In addition leakage leads to a temporal decay of activations. The updating rule of the attention map $a_{ij}$ is given by

[TABLE]

with the 2D-Gaussian $A_{ij}(t)$ centered upon fixation at time $t$, the distribution of fixation locations $\phi_{ij}$, and the rate of decay $\rho$.

Temporal evolution of activations in the inhibitory map $f_{ij}$ is very similar to the dynamics in the attention map $a_{ij}$. Temporal evolution consists of activation accumulation centered at fixation and a proportional temporal decay across the map. The updating rule for the fixation map $f_{ij}$ is given by

[TABLE]

with a 2D-Gaussian $F_{ij}(t)$ with standard deviation $\sigma_{0}$ centered at fixation at time $t$ and a decay rate of the fixation map $\omega$. The fixation map tracks fixated areas and is motivated by inhibition of return Klein \BBA MacInnes (\APACyear1999), which has been suggested as a mechanism to drive exploration in scenes Itti \BBA Koch (\APACyear2001). Although the role of inhibition of return has been questioned Smith \BBA Henderson (\APACyear2009), model simulations support inhibitory tagging as an important mechansim during scene perception (\citeNPRothkegel.VisionRes.2016; cf. \citeNPBays.JVis.2012).

While attention and fixation map evolve independently over time, both maps are subsequently normalized and combined into a single map for target selection. The potential $u_{ij}$ is given by

[TABLE]

where the exponents $\lambda$ and $\gamma$ are free parameters. \citeAEngbert.JVis.2015 fixed these parameters to $\lambda=1$ to reproduce the densities of gaze positions and $\gamma=0.3$ to reproduce spatial correlations between fixation locations. We kept these values in our simulations. The probability of a location $(i,j)$ to be chosen as the next saccade target can be extracted from the potential, i.e.

[TABLE]

where $S$ contains all positions on the grid with $u_{ij}(t)\leq 0$ and a free parameter $\eta$ that adds noise to the selection process so that every position has at least a minimal probability to be chosen as the next saccade target. Target selection in the SceneWalk model occurs at the end of fixation where the eyes move instantaneously. The intervals between successive saccades were drawn from a Gamma distribution with a shape parameter of $9$ and a scale parameter of $\sim\!0.031$, which corresponds to a mean fixation duration of 275 ms.

0.2 Parameter estimation

Simulations were based on different parameters for each of the four experimental conditions (presentation $\times$ image type). As a starting point we used the parameters reported in \citeAEngbert.JVis.2015. These parameters were estimated for fixations during the first presentation of natural images. We hypothesized that image type (natural scenes vs. textures) might affect target selection and decided to estimate all parameters for the first presentation of texture images anew. Previous work suggested that reinspection of images leads to a decreased attentional span Kaspar \BBA König (\APACyear2011\APACexlab\BCnt1, \APACyear2011\APACexlab\BCnt2). Hence, we decided to fix all parameters except for the sizes of the attention span $\sigma_{1}$ and inhibition span $\sigma_{0}$ for simulation of the second inspection of images. We used a genetic algorithm approach Mitchell (\APACyear1998) to estimate model parameters. Parameter estimation was based on the first five images for each image type. The remaining ten images in each image set were used for model evaluations. Limiting the analysis to the predicted images did not alter effects. We used first-order statistics (2D density of fixation locations) and the distribution of saccade lengths as an objective function to evaluate parameters. A list of estimated parameter values and standard errors can be found in Table 1.

0.3 Control Model: Joint Probability of Saccade Amplitude and Fixation Density

To investigate the influence of saccade amplitudes on the aggregation of points, we simulated a control model that generates both realistic distributions of fixation positions and saccade amplitudes. Fixation sequences began at the initial fixation position of the empirically observed scanpath and subsequent fixation positions were simulated iteratively. The next fixation position was chosen proportional to the joint probability of the empirical density of saccade amplitudes and the empirical density of fixation positions on an image. For each empirical fixation sequence we simulated a scanpath with the same number of fixations to exclude effects due differences in the number of fixations or sequences. The optimal bandwidth for density estimation was chosen by applying Scott’s rule of thumb (bw.scott) for fixation densities and unbiased cross-validation (bw.ucv) for saccade amplitudes.

0.4 Statistical modeling

For statistical analyses, we computed linear mixed effect models for each dependent variable using the lme4 package Bates \BOthers. (\APACyear2015) in R R Core Team (\APACyear2018). We log-transformed both dependent variables, since they deviated considerably from normal distributions. For the statistical model of the empirical data, we used the maximal possible random effect structure Barr \BOthers. (\APACyear2013) and ensured that none of the models was degenerate Bates \BOthers. (\APACyear2015). For our results we interpret all |t|$$>2 as significant fixed effects Baayen \BOthers. (\APACyear2008).

In R the linear mixed effect model for the empirical data can be written as

[TABLE]

with a fixed effects of presentation A, image type B and their interaction A:B and the corresponding random effect structure. As each image belongs to only one image type, image type is a between factor for images and was not included in the random effect structure for images.

The random effect structure for the simulated data differed from Equation 11, since fixation sequences were only based on one participant, that is, the average participant. Therefore, we did not estimate random effects for participants. Due to convergence problems we also removed the random slope of presentation for each image. This reduced the random effect structure to a single intercept per image for the simulated data.

Results

With our first analysis we expected to replicate the results of \citeAKaspar.PLoSOne.2011,Kaspar.JVis.2011 that a second inspection of the same image leads to decreased saccade amplitudes. In addition, we investigated whether the two image types (natural scenes vs. textures) differentially affected saccade amplitudes. In a second analysis, we tested the sensitivity of the pair correlation function (PCF) to our experimental manipulations. All experimental results were compared to model simulations of the SceneWalk model Engbert \BOthers. (\APACyear2015) and the joint probability control model that reproduces fixation densities and saccade amplitudes. In addition, we display the results of the inhomogeneous and the homogeneous point processes. These processes are plotted to demonstrate that results of the PCF are not generated by the method itself.

Saccade amplitudes

Figure 8 shows the distribution of saccade amplitudes  on natural scenes (top) and texture images (bottom) during the first and second presentation (solid vs. dashed lines). As expected, empirical saccade amplitudes distributions are positively skewed and long-tailed. Distributions of saccade amplitudes of the SceneWalk model and the joint probability control model deviate less from the empirically observed saccade amplitudes than the inhomogeneous and homogeneous point processes. At closer inspection, the SceneWalk model generates slightly shorter amplitudes while the joint probability control model generates slightly longer saccade amplitudes than our participants.

A linear mixed effects model (LME) for the experimental data <see Methods section; ¿Bates.JStatSoft.2015 revealed a significant fixed effect of Presentation but no effect of Image Type and no interaction (Tab. 2, left columns). Saccades were larger during the first inspection than during the second inspection. LME models of saccade amplitudes generated by the SceneWalk model and the joint probability control model replicated this effect qualitatively. However, the interaction of Presentation and Image Type also reached significance for the SceneWalk model. The reduction of saccade amplitudes during the second inspection was stronger on natural scenes than on texture images.

Second order statistics: Pair correlation function

We computed inhomogeneous pair correlation functions (PCFs) for each condition of the five point processes (Fig. 9).  We observed spatial correlations of fixation locations during individual trials in all conditions of our experimental data (red lines). Fixations locations were more abundant than expected from the overall inhomogeneity of fixation locations at short distances $r$. The estimated PCFs deviated from complete spatial randomness, i.e., $\hat{g}(r)>1$, at distances $r<4\degree$. More importantly, the second presentation of an image (dashed lines) led to increased PCFs for both natural scenes (top row) and texture images (bottom row). Statistically, we evaluated spatial correlations by computing the deviation from complete spatial randomness of each PCF, i.e., the summed deviation of the PCF from one for distances $0.1\leq r\leq 6.5$ (Fig. 10; cf. Eq. 6). A linear mixed effect model (LME) revealed a significant deviation from complete spatial randomness (intercept) and an effect of Presentation (Tab. 2). All other fixed effects were non-significant.  Thus, deviations from complete spatial randomness  were present in all conditions of our experiment with larger deviations during the second inspection of an image irrespective of image type.

The SceneWalk model replicated this pattern of results qualitatively (Fig. 9, purple lines). All conditions showed strong spatial correlations. PCFs deviated from complete spatial randomness  for distances $r<6\degree$. The effect extended to larger distances in our model simulations than in the experimental data.  We analyzed deviations from complete spatial randomness  with another LME for the SceneWalk model (Tab. 2; cf. Fig. 10). PCFs deviated from complete spatial randomness  (intercept). The effect was larger for the second inspection. No other fixed effect was significant for the deviation score of the SceneWalk model.

In order to check whether spatial correlations between fixation locations are primarily generated by the tendency of participants to generate short saccade amplitudes, we simulated a control model based on the joint probability of saccade amplitudes and the distribution of fixations (see Methods). The PCF of the control model showed correlations at all evaluated distances. In direct comparison to the experimental data the control model generated considerably weaker correlations at small distances $r<3\degree$ and generated stronger correlations at large distances $r>4\degree$. Thus, a model that is solely based on the generation of realistic saccade amplitudes and fixation densities generates a qualitatively different correlation pattern. We analyzed deviations from complete spatial randomness with another LME for the control model (Tab. 2; cf. Fig. 10). PCFs deviated from complete spatial randomness (intercept). The effect was larger for the second inspection and on texture image. We observed no interaction of Presentation and Image Type in the deviation score.

Finally, the average PCFs of the inhomogeneous point process indicated only small spatial correlations, i.e., $\hat{g}(r)\approx 1$ at all distances (Fig. 9, green lines). The absence of spatial correlations for the average PCF was expected, since the optimal bandwidth $\lambda$ for the estimation of the PCF was chosen to minimize the deviation from complete spatial randomness of the inhomogeneous process (see Step 2 of the estimation process). PCFs of the homogeneous point process were similar to PCFs of the inhomogeneous point process, i.e., $\hat{g}(r)\approx 1$ at all distances (Fig. 9, blue lines).

1 Discussion

During scene perception, fixations are not uniformly  distributed on an image. Instead fixations cluster in parts of an image due to bottom-up factors, top-down factors, and systematic tendencies of gaze control Tatler \BBA Vincent (\APACyear2008). We propose to use the pair correlation function (PCF) to investigate the relation of fixation positions within single trials and demonstrate that the PCF reveals aggregation of points at distances $r<4\degree$, that is, it is more likely to observe fixation locations in the proximity of another fixation location than expected by the overall distribution of fixation locations. This effect cannot be explained by the tendency of participants to generate short saccade amplitudes alone, as simulations of a control model that samples fixation locations from the joint probability of the density of saccade amplitudes and the density of fixation locations led to a qualitatively different PCFs. The control model underestimated aggregation at short distances and overestimated aggregation at large distances. In addition, the PCF responded sensitively to a memory manipulation in our experiment and revealed stronger aggregation of points during the second inspection than during the first inspection of an image.  Simulations of the SceneWalk model Engbert \BOthers. (\APACyear2015) demonstrated that a reduced attentional span could lead to reduced saccade amplitudes and  explain the observed results of the PCF during the second inspection.

1.1 Pair correlation function

Research on eye movements during scene perception has focused on the distribution of fixation locations across observers, as for example during the evaluation of saliency models Bylinskii \BOthers. (\APACyear2015). However, this approach neglects dependencies between fixation locations during a trial Engbert \BOthers. (\APACyear2015). The PCF is a method from spatial statistics to evaluate the relation of pairs of points Diggle (\APACyear2013) and reveals whether points solely depend on the inhomogeneity of a point process or whether pairs of points affect each other mutually at a given distance $r$. The PCF can be applied to eye movement data (i.e., fixation locations) in three steps. In a first step, inhomogeneous and homogeneous point processes need to be simulated to evaluate the PCF estimation. Both point processes generate fixation locations that are independent at all distances $r$. Hence, PCFs of both processes are expected to show no spatial correlations at any distance $r$. During a second step an optimal bandwidth needs to be chosen for the estimation of the PCF. As a criterion we suggest to use a bandwidth for which the PCF of the simulated  inhomogeneous point process has the least deviation from complete spatial randomness, i.e. shows no spatial correlations. In a last step, the computed bandwidth is used to compute PCFs for each individual trial.

In scene perception, the PCF characterizes whether fixation locations can be explained by the underlying distribution of all fixation locations (no spatial correlations between pairs of points) or whether fixation positions interact with each other at distance $r$. PCFs revealed spatial correlations of fixation locations in all conditions in our experimental data. Fixations were more abundant at distances $r<4\degree$ than we would expect from the inhomogeneity of fixation locations alone. Beyond $4\degree$ fixation locations were independent of each other. Thus, observing a fixation increased the probability of observing more fixations than expected by the overall inhomogeneity within $4\degree$. Beyond $4\degree$ fixations were as likely as predicted by the local intensity of fixation locations. As expected, neither the inhomogeneous nor the homogeneous point process revealed strong spatial correlations, since fixation locations were independent of each other for these point processes. Therefore, aggregation observed in our PCFs is generated by the empirical data and is not the result of the method itself.

Finally, the PCF provides a quantitative statistic of spatial correlations and can be used to compare data sets generated by different processes or under different conditions. The data sets may even differ in the overall inhomogeneity of the fixation location distribution, since the inhomogeneous PCF takes this inhomogeneity into account. Hence, we were able to compare point patterns of the first and second inspection of an image, point patterns on different image types as well as empirical and simulated point patterns. Even though the distribution of fixation locations and the distribution of saccade amplitudes were similar among these data sets, the PCF revealed considerable differences in the correlation patterns. For example, the experimental data, simulations of the SceneWalk model and simulations of the joint probability control model unveiled strong differences in spatial correlations. While the SceneWalk model generated stronger but qualitatively similar spatial correlations when compared to the experimental data, the joint probability control model produced a different pattern. Fixation locations aggregated at all evaluated distances, but produced much weaker correlations for short distances. Thus, short saccade amplitudes are not sufficient to produce the strong correlations at distances $r<2\degree$. There must be additional mechanisms that let participants fixate the same locations in a scanpath as through direct regressions <cf. facilitation of return; ¿Smith.VisCogn.2009 and through reinspections later during a trial.

Interestingly, image type did not affect spatial correlations. However, we only tested two types of images and a broader range of image types is needed to see the generalizability of this result. To conclude the PCF is a powerful tool to compare spatial correlations of point patterns and will help to understand the eye-movement dynamics under different experimental conditions as well as the dynamics of different models of eye-movement control during scene perception.

1.2 Repeated presentation of images

To test the sensitivity of the PCF, we recorded eye movements of participants while viewing an image twice. The repeated presentation was expected to result in shorter saccade amplitudes due to a reduced attentional span Kaspar \BBA König (\APACyear2011\APACexlab\BCnt1, \APACyear2011\APACexlab\BCnt2). We replicated shorter saccade amplitudes during the second inspection independent of image type in the experimental data.  In addition, we observed stronger spatial correlations within $4\degree$ during the second inspection.

Reduced saccade amplitudes as well as increased aggregation might be generated by a reduced attentional window Kaspar \BBA König (\APACyear2011\APACexlab\BCnt1, \APACyear2011\APACexlab\BCnt2). We tested this hypothesis with simulations of the SceneWalk model Engbert \BOthers. (\APACyear2015). Parameter estimation led to a reduced attentional span in the SceneWalk model during the second inspection of images which in turn led to shorter saccade amplitudes and stronger aggregation of fixation locations as quantified by the PCFs. Hence, our simulation results are in agreement with the interpretation of a reduced attentional span during repeated inspection of images.

In its current form the SceneWalk model overestimates aggregation in particular during the second inspection. However, our results are predictions and are not optimized to account for the experimentally observed PCFs. Beside the size of the aggregation, we observed a difference in the functional form of the PCF. PCFs of the SceneWalk model did not decrease as fast as those estimated from experimental data. This might have resulted from the Gaussian form of our attentional window with its long tails. A revision of the model using a modified attentional window might improve model fit in this respect <cf.,¿Schutt.PsycholRev.2017.

1.3 Relation to other metrics

Several measures can be used to study the dynamics of scanpaths. The measure that is most strongly related to the PCF is based on Voronoi diagrams Over \BOthers. (\APACyear2006). The suggested method provides a measure for the uniformity of a pattern of fixation locations and is normalized by the number of fixation locations. Hence, the measure can be used to compare the uniformity of different data sets. While the Voronoi method allows to estimate the inhomogeneity of a point pattern, the inhomogeneous PCF takes this inhomogeneity into account and reveals spatial correlations between points that cannot be explained by the inhomogeneity. Most importantly, the inhomogeneous PCF can be used to compare point patterns with different levels of inhomogeneity. Thus, the Voronoi method and the PCF complement each other.

A number of scanpath comparison methods have been proposed that provide metrics to describe the similarity of the dynamics of two scanpaths <for a review see¿Anderson.BehavResMeth.2015. Each metric quantifies unique aspects of the scanpaths and helps to understand which aspects resemble each other or differ between two scanpaths. However, the focus of all methods is the comparison of scanpaths and requires matching scanpaths. In contrast, the PCF does not require pairs of scanpaths as each scanpath is evaluated by its own. This allows to compare arbitrary scanpaths. In addition, PCF and scanpath comparison metrics might come to very different conclusions. Since the PCF describes spatial correlation of points, two completely different point patterns might produce the same correlations. At the same time, two seemingly similar point patterns might result in very different spatial correlations.

Other approaches have described a large variety of individual effects that influence eye movement behavior during scene perception <e.g.,¿Smith.VisCogn.2009,Tatler.JVis.2007,Tatler.JEMR.2008 and best practices have been suggested to estimate appropriate baselines for some of them, as for example for the central fixation bias Clarke \BBA Tatler (\APACyear2014); Clarke \BOthers. (\APACyear2017). Each individual effect describes parts of the dynamics present in eye movements, but none is a sufficient metric to capture the overall complexity. Since the PCF takes all fixations of a sequence into account, it uncovers the consequences of these effects that remain unnoticed on shorter time scales. Hence, the PCF acts in concert with these other methods to understand eye movements during scene perception.

Conclusions

The pair correlation function (PCF) is a powerful tool to analyze spatial correlations of fixation locations. During scene perception the PCF reveals aggregation of fixations during individual trials and reacts differentially to experimental manipulations. Simulations of a computational model demonstrate that a reduced attentional span leads to increased aggregation of fixation locations. Our work provides an example how spatial statistics and computational modeling can be combined to investigate general statistical properties of eye movement control.

Acknowledgments

This work was funded by Deutsche Forschungsgemeinschaft (DFG, German Research Foundation, grant EN 471/13–1 to R.E.) and the Collaborative Research Centre SFB 1294, project number B05 and by Bundesministerium für Bildung und Forschung (BMBF) via Bernstein Center for Computational Neuroscience Berlin (Project B3, Förderkennzeichen 01GQ1001F to R.E. & F.A.W.). Simon Barthelme was supported by an ANR grant (GenGP, ANR-16-CE23-0008).