As we drive along a road, a continuous optic flow pattern projected onto our retina is available to our visual system, and at any given point in time it can be used to compute the time remaining until we pass an object of interest. This time to passage (TTP), together with cues about motion trajectory, allows us to anticipate and judge oncoming objects regarding their path of movement and to prepare time-critical motor actions. Despite ample research on the topic, it still remains unresolved that how TTP is computed and which other optical (such as object velocity and expansion cues) are being exploited when observers are asked to provide judgments regarding the time to passage of an oncoming object. One reason why these cues have been eluding identification may lie in the adaptive nature of the visual system. In a recent psychophysical study, we have shown that the visual system appears to employ an adaptive strategy that changes with the task at hand  . We presented a moving cloud of randomly placed dots viewed through a square aperture, consistent with forward observer motion but devoid of local expansion cues. Therefore, the dots remained of constant retinal size throughout the motion display. In each trial, two dots were colored red, and upon occluding the display, observers had to indicate which of the two dots would pass first the observer’s eye plane. Combined with two coherence levels (0% and 75%), that is, the proportion of background dots that could not deviate from the motion they should perform when thought of as stationary points in the 3D flow-field approaching the observer. Our results showed that, when no coherent optic flow was available (coherence 0%), observers resorted to the use of a relative velocity strategy and picked the dot with the faster screen velocity. However, when optic flow was highly coherent (75%), observers used a more complex strategy involving the global flow-field information.
Global tau is a property of coherent optic flow that relies on the systematic change of distances on the retina and between-objects. Global tau makes the assumption that 3D-distances between the objects in the world remain constant thus producing coherent optic flow, and it can be computed from the relative rate of change of the angular displacement of the target from the observer’s line of sight  . This flow-field analysis is exploiting properties of coherent optic flow which cannot be reduced to, but sometimes are correlated with local expansion cues (expanding retinal object size). For instance, in simulated forward motion through a cloud of fixed spherical objects, the systematic change of retinal distances among these objects specifies the direction of the observer’s motion through the cloud, provided that the objects remain static. The angular subtense between the observer’s path (track vector) and a given object, or more precisely the relative rate of change of this angular subtense, gives away the object’s TTP. Far-away objects typically produce less centrifugal retinal motion than close-by objects  . Interestingly, TTP is specified for expansionless objects as long as they do not coincide with the tracking direction. Thus, the optic flowfield provides TTP information even when it is devoid of local expansion cues   .
In naturalistic scenarios, in which the retinal size of the targets does expand, both the local expansion cues (local tau) and the global tau cues can be exploited by the visual system to predict TTP of the target. Studies on the utilization of local tau information report that the neural substrate involved in the extraction of local tau expansion cues is the locus rotundus in pigeons , and more recently fMRI studies in humans point to the superior colliculus, the pulvinar nucleus of the thalamus, and cortical regions associated with motor preparation   . fMRI studies of time-to-contact estimation tasks (TTC) have demonstrated significant cortical activation in left inferior parietal regions  , superior parietal, motor and cortical regions around the central sulcus , the insula, and inferior and middle frontal areas  . Bilateral hMT areas play important roles in optic flow processing    . All of these regions show more or less specific activation in response to the local retinal expansion of looming objects that move toward the observer. However, the specific neural substrate associated with global tau has not yet been identified.
In the present study we used fMRI to identify such regions in human observers. Given that at the behavioral level, observers differ in the strategies they use between the case of local expansion scenarios and expansionless global optic flow, one would expect some shared but also some specific cortical areas to be involved in TTP judgment.
Imagine a cloud of fixed expansionless objects (dots) through which an observer is moving. Two dots are marked red while all the others are white. If these dots are at equal lateral distance on opposite sides of the track vector, then the dot that is sagittally farther away from the observer will project closer to the focus of expansion in the retinal flow pattern. If the observer is asked to judge which of the two marked dots is closer, she/he could base the decision on this fact. In other words, in the case of such symmetrical lateral spacing, observers might use an image-based strategy once they have discerned the track vector from the optic flow. Reducing the coherence of the optic flow makes it harder to determine the track vector, and performance should break down or resort to some other strategy. For instance, subjects may merely base their judgments on how far a target is from the center of the screen. We have previously found that observers employ flexible strategies that can use a combination of global flow analysis and image-based cues  . Thus, we created stimuli that provide information about the direction of self-motion (track vector and track velocity) and others that do not. The former provides global information containing a certain amount of noise (75% coherence), the latter preserves the local motion magnitude but removes all global information (0% coherence). Note that local tau information was absent at all times.
We have collected task-based fMRI data while observers were making TTP judgments in the absence of local tau information. By manipulating the initial positions of the target objects relative to the observers’ track vector and by changing the coherence of optic flow (75% or 0%), we used a limited number of specific information sources that observers could exploit when making TTP judgments. Based on our previous psychophysical study , we hypothesized that observers use global flow information in the case of 75% coherent flow, but resort to guesswork when global flow is incoherent.
Our behavioral results are consistent with previous psychophysics findings: TTP judgments reflected the differential use and integration of multiple sources of information, including global optic flow, object retinal velocities, and other depth cues    . In the experiment detailed below, we will focus on and interpret the cortical and subcortical activations in light of the likely response strategies.
Seven subjects (5 females, 2 males, mean age = 24.42 years, SD = 4.82 years) participated in the study. They were graduate students at Boston University, recruited from our pool of subjects. All of them had normal or corrected-to-normal vision. All underwent a psychophysical testing session prior to the scan, to make sure that their performance was at least 70% correct for sec for the symmetric configuration regardless of the initial x-offset and background motion coherence (0% or 75%)  . All participants signed an informed consent form before the start of the experimental sessions in accordance with the requirements on research involving human subjects, as approved by the Massachusetts General and Boston University Institutional Research Boards. All subjects fully satisfied the inclusion criteria for participating in MRI/fMRI studies and none of the exclusion criteria were met. They participated previously in other psychophysical and functional imaging tasks conducted by our research team. Those studies had no similarity with the task reported here. All subjects reassured us they could pay attention throughout an experimental task, maintain fixation, and stay still during the imaging experiments.
2.2. Apparatus and Data Acquisition
The stimuli were generated on an Intel-based Macintosh laptop and displayed at a resolution of 1024 × 768 pixels and a refresh rate of 75 Hz. Two of the dots, referred to as target dots, were red (51.20 cd/m2) and the rest of them were white (79.55 cd/m2), all displayed against a gray background (10.22 cd/m2). They were back-projected onto a translucent screen (27.3 cm × 36.5 cm) using a LCD projector. Subjects viewed the translucent screen through a mirror mounted on the head coil of a whole-body scanner. The distance between the eyes of the subject and the mirror was approximately 4 cm and the distance between the mirror and the screen was approximately 81 cm, therefore, the total viewing distance was about 85 cm. This setup provided a square viewing aperture subtending 17˚ × 17˚. fMRI data were acquired at Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, using a 3T Siemens whole body scanner and a standard 8-channel head coil. Structural images were obtained as T1 weighted magnetization prepared rapid acquisition gradient echo images (MPRAGE) (128 slices with slice thickness of 1.33 mm, voxel size: 1.00 × 1.00 × 1.33, FOV = 256, TR = 2.53 sec, TE = 3.39 msec, flip angle = 90˚). Two T1-weighted images were collected for each subject. Functional images were obtained with gradient echo, echo planar (EPI) interleaved sequence (33 slices with slices oriented along the AC-PC line, slice thickness of 3 mm with 20% distance factor, FOV = 200, TR = 2.00 sec, TE = 30 msec, flip angle = 90˚) for measurement of BOLD contrast.
2.3. Stimuli and Experimental Procedure
A field of moving white dots simulated the observer’s forward self-motion in 3D and was presented through a square viewing aperture. The dots remained stationary with respect to one-another in the simulated space. Subjects were asked to indicate which one of two red dots would pass their eye plane first, mimicking they were moving forward through the field. All of the dots subtended 2 pixels × 2 pixels (4 arcmin × 4 arcmin) throughout the simulated approach and were placed such that they maintained a density of 2 dots/deg2. The screen size of all dots, including the targets, remained constant, thus eliminating all local tau cues. The motion of the dots simulated the subjects’ forward self-motion along a straight-line trajectory at a speed of 150 cm/s. In each trial, the direction of simulated self-motion was toward the center of the aperture. Dots that moved out of the volume behind the observer’s eye plane were randomly assigned to new locations such that the density of the dots remained constant (Figure 1)  .
Figure 1. Schematic of the virtual trapezoidal volume. White dots simulating forward self-motion were randomly distributed between 260 cm and 2060 cm from the observer. The two red target dots were embedded in the flow field and moved with the same speed as the flow field according to their instantaneous position within the volume. The direction of self-motion matched the center of the aperture  . The dots maintained constant screen size.
The target dots were placed at different depths such that the difference between their passage times (tau difference Δtau) was set to be 0.5 sec for all trials. This value was chosen because it should be just above detection threshold, rendering the task of identifying the leading target meaningful but still allowing for errors. The initial depth of the reference target was 1200 cm. Thus, the possible TTP values from stimulus onset until both targets would have passed the observer, ranged from 7.5 s to 8.5 s. This left ample time for the 2AFC responses, which had to be made before the leading target passed the vertical eye-plane of the observer. Observers were not required to respond as quickly as possible.
The two target dots were placed such that they were always on opposite sides of the track vector. Their lateral distances from the track vector (x-offsets) were either 10 cm or 50 cm in the simulated space. This resulted in four different lateral target offset combinations: Two equidistant symmetric placements (leading target 10 cm to one side-trailing target 10 cm to the other side, or 50 cm - 50 cm). In the other two combinations, the targets were placed with asymmetric x-offsets (10 cm - 50 cm, and 50 cm - 10 cm). These 4 stimuli were paired with two coherence levels of the background dots. Remember that an entirely incoherent flow-field no longer specifies observer motion. Our previous psychophysical data  showed that 50% coherence of the dots begins to provide prominent global flow information relative to the condition of 0% coherence dots. This design led to eight unique stimuli, each of which was repeated 16 times within a run, using a randomized event-related design paradigm.
The visual stimulus was occluded after 3 seconds. The next trial would not be presented until a decision had been made. The timing and order were randomized using optseq2 (http://www.freesurfer.net/optseq/). Inter-stimulus Intervals (ISIs) between trials varied from 1 - 7 s. Frames with static dots were presented within the ISIs, serving as a baseline condition. During the whole scanning period, subjects were required to fixate a small central cross (40 × 40 arcmin). Stimuli were presented binocularly in a two-alternative forced choice (2AFC) paradigm without feedback. The subjects’ task was to determine which of the two targets would arrive at their eye plane first. Subjects entered their responses by pressing a designated key on a magnet-compatible button box.
A separate block design employing a MT localizer task was performed by all subjects in two runs. The human middle temporal complex (hMT) has been shown to be highly involved in motion processing, including optic flow. Accordingly, area hMT was functionally localized by utilizing moving and static dot patterns , so that we could locate the exact position of hMT for further analysis  . Other anatomical regions were defined with normalized functional images, using the Automated Anatomical Labeling (AAL) atlas.
2.4. Data Analysis
Imaging data analysis was performed using the Statistical Parametric Mapping software package (SPM12, Wellcome Department of Cognitive Neurology, London, UK) and utilizing MATLAB (The MathWorks, Natick, Massachusetts, USA). The preprocessing steps were as follows: 1) format conversion by converting data from original DICOM files to Nifti files, 2) slice timing eliminated the time shift for all voxels of functional images, 3) realignment by motion correction, 4) co-registration of the anatomical image with the mean functional image, 5) spatial normalization normalized all images to a standard space (Montreal Neurological Institute, MNI) with a voxel size of 3 × 3 × 3 mm, and 6) smoothing: all functional data were smoothed with a FWHM (Full Width at Half Maximum) kernel of 4 mm.
For each subject, the onset and duration of each condition was modeled by a general linear model (GLM). The motion parameters from realignment were also used as multiple regressors in generating the design matrix. Trials with correct responses and incorrect responses were separated as different conditions. With the contrast images of each subject, a group level randomeffect analysis was performed for each condition. The resulting t-value maps were set as uncorrected for multiple comparisons, p < 0.05. Clusters with less than 10 contiguous voxels were excluded.
Based on group-level activation maps in normalized space, we defined several functional regions of interest (ROI) for each single subject. Every functional ROI was defined as a sphere, with its center at the respective local maximum of the activation cluster and with a 5-mm radius. Subsequently, we calculated the percent BOLD signal change for each functional ROI using Marsbar  .
Bilateral hMT areas were defined using localizer tasks  . We set the minimum overlapped proportion on individual activation maps as 0.5, based on the group-constrained subject-specific (GSS) method  .
3.1. Behavioral Performance
Response accuracy for each condition was first calculated per subject and then averaged across subjects (Figure 2). As expected, subjects could do the task and performed well above chance when the two target dots were symmetric around the track vector of the simulated motion. This was the case for both 0% coherence and 75% coherence conditions. Thus, subjects exploited the simple image-base cue of eccentricity to guide their answer choices. In contrast, when the two target dots were placed asymmetrically, such that the correlation of target eccentricity and proximity to the observer was severely reduced, the response accuracy was around chance level. Thus, the behavioral performance showed that in the absence of local expansion cues, TTP judgments were based on retinal eccentricities.
Interestingly, when the flow field contained additional information about the track vector, performance improved. When the target dots were symmetric, subjects performed better under the 75% coherence condition compared to 0% coherence. Thus, global motion information provided by background dots
Figure 2. Behavioral performance accuracy in % as a function of lateral target-offset pairing and coherence of optic flow. Offset pairings indicate the lateral distance of the leading and trailing targets from the track vector. They were always placed on opposite sides. Thus, a pairing of 10 vs 50 signifies that the leading target was placed at a lateral offset of 10 cm in virtual space, whereas the trailing target was offset by 50 cm to the opposite side. Error bars indicate standard derivations.
enhanced subjects’ performance only if the targets were spaced symmetrically around the track vector. However, when the target dots were asymmetric, there were no significant differences between 0% coherence and 75% coherence (p > 0.05 in paired t-test). These results replicate previous behavioral results we collected with a similar experimental design  .
3.2. Functional Imaging
We contrasted activation during the trials against the activation within the static dots presentation (baseline) to obtain significance activation maps. The analysis was performed separately for 0% and 75% coherence levels in the optic flow field Table 1 & Table 2. The runs with the same coherence value were grouped together and the contrasts were done separately for each experimental condition. Based on behavioral performance, we separated correct responses and incorrect responses for each condition. Therefore, the experimental conditions were the cases where initial target x-offsets were symmetric and 10 cm from the center of the aperture (10 vs 10), symmetric and 50 cm from the center of the aperture (50 vs 50), and the leading target’s initial x-offset was 10 cm from the center of the aperture (10 vs 50) or the leading target’s initial x-offset was 50 cm from the center of the aperture (50 vs 10), either with correct responses or incorrect responses (Figure 3).
Table 1. Activation position of the local maxima in MNI space for 0% coherence trials by lateral target-offset pairing.
Table 2. Activation position of the local maxima in MNI space for 75% coherence trials by lateral target-offset pairing.
AG: angular gyrus, IFG: inferior frontal gyrus, IPG: inferior parietal gyrus, ITG: inferior temporal gyrus, MFG: middle frontal gyrus, MOG: middle occipital gyrus, MTG: middle temporal gyrus, ParaCL: paracentral lobule, PostCS: postcentral cortex, PreCS: precentral cortex, SFG: superior frontal gyrus, SMG: supramarginal gyrus, SOG: superior occipital gyrus, SPG: superior parietal gyrus.
Figure 3. Cortical activation maps as a function of lateral target-offset pairing and coherence of optic flow. The threshold t-value was at 1.943, uncorrected p < 0.05, extent 10 voxels. Panel (a) condition 10 vs 10; (b) 10 vs 50; (c) 50 vs 10; (d) 50 vs 50. Area hMT showed activation in all conditions based on the findings from the localizer task.
Figure 3 ((a) condition 10 cm vs 10 cm; (b) 10 cm vs 50 cm; (c) 50 cm vs 10 cm; (d) 50 cm vs 50 cm) shows the activation maps of trials with correct responses with a threshold t-value 1.943, extent of 10 voxels (p < 0.05, uncorrected), across all the subjects. Table 1 and Table 2 show the corresponding coordinates of the local maxima of these clusters, in MNI space.
In general, across all the subjects, the activation areas were distributed along the motion processing pathway  . Consistent with previous research , the activation was not only showing along the ventral pathway (also known as “what” pathway), it also showed along the dorsal pathway (also known as “where” pathway). Activation during the stimulus motion was significantly elevated compared to baseline in the occipital lobe, including the superior occipital gyrus (SOG), temporal lobe including middle temporal gyrus (MTG), inferior temporal gyrus (ITG) and superior temporal sulcus (STS), parietal lobe including superior parietal gyrus (SPG) and intraparietal sulcus (IPS), frontal lobe including middle frontal gyrus (MFG) and inferior frontal gyrus (IFG). In addition, the pre-central sulcus and post-central sulcus were activated bilaterally across all subjects.
For all x-offset target conditions, when we compared the activation for 0% coherence and 75% coherence conditions, the latter resulted in more distributed activation aroundthe bilateral pre-central and post-central temporal areas, MTG, IFG and IPS.
Based on previous functional imaging data related to TTP and TTC processing (e.g. Field & Wann, 2005), we defined several functional regions of interest based on the Automated Anatomical Labeling (AAL) atlas. Figure 4 shows the BOLD percentsignal changes in these functional ROIs. By and large the BOLD percentsignal changes increased from 0% to 75% coherence levels in bilateral precentral, postcentral, and middle temporal areas, whereas they decreased in superior and middle frontal areas. In the inferior frontal and parietal cortical regions, including intra-parietal sulcus, there was not much difference between the two coherence levels. Specifically, in the left precentralsulcus (PreCS: p = 0.045), left postcentralsulcus (PostCS: p = 0.027) and bilaterally in the middle temporal (MTG: p = 0.025) cortical regions, percentsignal changes at 75% coherence were significantly larger than those at the 0% coherence level, for the symmetric condition (10 vs 10). In the right hemisphere MTG: p = 0.024), significant activations were also found in the asymmetric condition (10 vs 50) at both coherence levels. Bilaterally in the middle frontalgyrus (MFG: left: p = 0.012, right: p = 0.005) and right PreCS (p = 0.046) areas, there were significant differences between the 0% and the 75% coherence level in the asymmetric condition (50 vs 10). In the left inferior frontal gyrus (IFG: 0.034), percentsignal changes were also found significantly higher under 75% coherence than 0% coherence, in the symmetric condition (50 vs 50). In MFG in the left hemisphere (p = 0.030), percent signal changes were found to be significantly higher at 0% coherence than at 75% coherence, in the symmetric condition (10 vs 10). Interestingly, bilateral significant differences in percent signal change were found in hMT areas for asymmetric conditions, but not for symmetric conditions (left hMT: p = 0.013 in 10 vs 50, right hMT: p = 0.032 in 10 vs 50 and p = 0.036 in 50 vs 10).
The percent signal changes suggest that bilateral precentral and postcentral sulci as well as a MTG are highly involved in the processing of global optic flow. The activation in hMT bilaterally suggests that more complicated visual processing is performed when there is more than one cue that subjects might use (e.g. global
Figure 4. BOLD percent signal changes plotted by functional region of interest (ROI) and lateral target-offset condition (*indicates p < 0.05; **indicates p < 0.025).
optic flow, object velocities, and symmetry heuristics) for their TTP judgments. The activation of IFG and MFG may be underlying the process of decision making for solving the task.
3.3. Correlation of Behavioral Responses with BOLD Percent Signal Change
For each condition described above, we also compared the trials in which subjects gave the correct responses to those with incorrect responses. For the symmetric conditions (offsets of 10 cm vs 10 cm and 50 cm vs 50 cm), the activations were more extended for trials with correct responses than for those with incorrect responses. The activations in precentral postcentral, IPS, MTG and MFG regions, in both correct and incorrect trials, suggest the involvement of these areas in solving the underlying task. For asymmetric conditions with the leading target at 10 cm x-offset (10 vs 50), invalid image velocity information was provided to the subjects, resulting in more distributed activation in incorrect trials than in correct trials, at 0% coherence. When coherent background dots were presented (75% coherence), thereby providing global motion information, subjects exploited multiple cues, resulting in more distributed activation in correct trials than in incorrect trials. Conversely, for the asymmetric conditions with the leading target at 50 cm x-offset (50 vs 10), valid image velocity information was provided, and activations were more distributed in correct trials than incorrect ones in both 0% and 75% coherence level.
In this study, we have used fMRI to record observers’ TTP judgments in the absence of local expansion information. During simulated forward motion, the observer had to judge, which of two red dots would pass him/her first. We have presented the information indicative of forward motion of the observer (global flow information) by manipulating the coherence of the flow field (no coherence vs. 75% coherence). We also manipulated the lateral offsets of the targets from the track vector and the initial target depths from the observer. Since local expansion information was not present in the optic flow, only global flow information could be used for the task. From this global flow information observers could, in principle, utilize simple image cues and/or the complex pattern provided by the entire dot field, but they could do so only when the dot field (RDK) moved coherently. In the case of incoherent motion of the dots, observers could have used other cues, such as relative screen velocity of the targets, which remained available in the display. However, the latter image cues were only valid when the targets were spaced symmetrically around the track vector. In cases of asymmetric spacing and incoherent optic flow, no useful information about relative TTP remained and, as expected, subjects’ behavioral performance was at chance.
For scenarios with asymmetric target spacing and coherent optic flow, we expected above chance performance if and only if the global flow information could be fully exploited. However, this was not the case. Coherence of optic flow had a small positive effect, only when the targets were spaced symmetrically around the track vector. In the presence of coherent optic flow in the background, this produced stronger and temporally extended cortical activation along the middle temporal gyrus, the precentral, and postcentral sulcus regions.
When the two targets were spaced asymmetrically, coherent optic-flow could in principle produce good performance if the relative rates of change of the respective angle between the target object and the track vector are considered. However, our results showed that this was not the case suggesting that in the presence of asymmetric targets, observers failed to exploit the global flow information for judging TTP. Simpler image cues, such as the targets’ positions and relative velocities would only provide valid information if the targets are spaced symmetrically to either side of the track vector (or the center of the aperture). Thus, only in the presence of symmetric targets, could above-chance performance be reached with incoherent flow. This did in fact improve performance but failed to approach perfection. Thus, with multiple sources of information, when judging TTP, subjects appear to integrate several cues through an economic strategy that mostly rely on image cues. This strategy becomes clearly noticed when local tau information is missing and the symmetry assumption holds.
The cortical activities during TTP judgments reflect this economic strategy. In general, subjects showed higher and broader activations on trials with 75% coherence than on those with 0% coherence. This suggests that they did processoptic flow information when making TTP judgments, which is consistent with previous studies  . When the two targets were located symmetrically around the observer’s track vector, the percent-signal changes on bilateral PreCS, PostCS and MTG proved significantly higher when global optic flow was available, as compared to when it was absent. This is where coherent global flow provided a behavioral advantage.
Previous retinotopic mapping and fMRI studies in humans have established a continuum of several motion-selective regions, including cortical areas hMT and superior parietal gyrus   . In our study, bilateral hMT cortical regions were activated across all subjects, regardless of their performance on judging TTP.
Significant differences in percent signal change were found bilaterally in hMT when the two targets were asymmetric. Activity during stimuli with 0% coherence was higher than during stimuli with 75% coherence in the 10 vs 50 condition, whereas activity during stimuli with 75% coherence was higher than during stimuli with 0% coherence in the 50 vs 10 condition. This points to lateralized differences that reflect the complex reaction of hMT to changes in global and local information. Remember that only when global cues were unavailable, subjects based their TTP judgments on the velocity discrepancy between the targets   .
Consistent with previous research, we also found activation in bilateral superior colliculus (SC), which is an area involved in motor preparation and attention  . The percent signal change was not significantly different in the conditions of 0% coherence and 75% coherence, which is expected as there should not be a difference in motor (response) preparation between different coherence levels.
In summary, in this study we investigated the neural substrate of the mechanisms involved in TTP judgments in the absence of local expansion cues. Previous behavioral results suggested that the subjects base their TTP judgments on the integration of multiple sources of information, with emphasis on image cues, such as target velocity, which are supplemented by global optic flow information, if the latter is coherent. Accordingly, and consistent with previous studies    , out fMRI results show a broad range of activation along the visual motion processing pathway, which reflects the complex information processing strategy. Unlike in pigeons, there does not seem to be one area dedicated to TTP processing    . Instead, the BOLD percent-signal changes show that PreCS, PostCS and MTG are involved in global information processing. Strong activation has also been found in bilateral hMT areas. Further investigation of the cortical involvement in TTP judgments will contribute towards a better understanding of how temporal and spatial perceptual mechanisms are integrated.
 Beardsley, S.A., Sikoglu, E.M., Hecht, H. and Vaina, L.M. (2011) Global Flow Impacts Time-to-Passage Judgments Based on Local Motion Cues. Vision Research, 51, 1880-1887.
 Kaiser, M.K. and Mowafy, L. (1993) Optical Specification of Time-to-Passage: Observers’ Sensitivity to Global Tau. Journal of Experimental Psychology: Human Perception and Performance, 19, 1028-1040.
 Tresilian, J.R. (1995) Perceptual and Cognitive Processes in Time-to-Contact Estimation: Analysis of Prediction-Motion and Relative Judgment Tasks. Perception & Psychophysics, 57, 231-45.
 Billington, J., Wilkie, R.M., Field, D.T. and Wann, J.P. (2011) Neural Processing of Imminent Collision in Humans. Proceedings of the Royal Society B-Biological Sciences, 278, 1476-1481.
 Uesaki, M. and Ashida, H. (2015) Optic-Flow Selective Cortical Sensory Regions Associated with Self-Reported States of Vection. Frontiers in Psychology, 6, 775.
 Assmus, A., Marshall, J.C., Ritzl, A., Noth, J., Zilles, K. and Fink, G.R. (2003) Left Inferior Parietal Cortex Integrates Time and Space during Collision Judgments. Neuroimage, 20, S82-S88.
 Assmus, A., Marshall, J.C., Noth, J., Zilles, K. and Fink, G.R. (2005) Difficulty of Perceptual Spatiotemporal Integration Modulates the Neural Activity of Left Inferior Parietal Cortex. Neuroscience, 132, 923-927.
 Delucia, P.R. and Warren, R. (1994) Pictorial and Motion-Based Depth Information during Active Control of Self-Motion—Size Arrival Effects on Collision-Avoidance. Journal of Experimental Psychology-Human Perception and Performance, 20, 783-798.
 DeLucia, P.R. (2004) Multiple Sources of Information Influence Time-to-Contact Judgments: Do Heuristics Accommodate Limits in Sensory and Cognitive Processes? Advances in Psychology, 135, 243-285.
 Tootell, R.B.H., Reppas, J.B., Dale, A.M., Look, R.B., Sereno, M.I., Malach, R., Brady, T.J. and Rosen, B.R. (1995) Visual-Motion after Effect in Human Cortical Area Mt Revealed by Functional Magnetic-Resonance-Imaging. Nature, 375, 139-141.
 Vaina, L.M., Soloviev, S., Calabro, F.J., Buonanno, F., Passingham, R. and Cowey, A. (2014) Reorganization of Retinotopic Maps after Occipital Lobe Infarction. Journal of Cognitive Neuroscience, 26, 1266-1282.
 Fedorenko, E., Hsieh, P.J., Nieto-Castanon, A., Whitfield-Gabrieli, S. and Kanwisher, N. (2010) New Method for fMRI Investigations of Language: Defining ROIs Functionally in Individual Subjects. Journal of Neurophysiology, 104, 1177-1194.
 Pitzalis, S., Sereno, M.I., Committeri, G., Fattori, P., Galati, G., Patria, F. and Galletti, C. (2010) Human v6: The Medial Motion Area. Cerebral Cortex, 20, 411-424.
 Tootell, R.B., Mendola, J.D., Hadjikhani, N.K., Ledden, P.J., Liu, A.K., Reppas, J.B., Sereno, M.I. and Dale, A.M. (1997) Functional Analysis of V3A and Related Areas in Human Visual Cortex. Journal of Neuroscience, 17, 7060-7078.
 Kerzel, D., Hecht, H. and Kim, N.-G. (1999) Image Velocity, Not Tau, Explains Arrival-Time Judgments from Global Optical Flow. Journal of Experimental Psychology: Human Perception and Performance, 25, 1540-1555.
 Wang, Y.-C., Jiang, S. and Frost, B.J. (2009) Visual Processing in Pigeon Nucleus Rotundus: Luminance, Color, Motion, and Looming Subdivisions. Visual Neuroscience, 10, 21-30.