AbstractIn this paper, we used time-domain functional near infrared spectroscopy (TD-fNIRS) to evaluate the haemodynamic response function (HRF) in the occipital cortex following visual stimulation in glaucomatous eyes as compared to healthy eyes. A total of 98 subjects were enrolled in the study and clinically classified as healthy subjects, glaucoma patients (primary open-angle glaucoma) and mixed subjects (i.e. with a different classification for the two eyes). After quality check data were used from HRF of 73 healthy and 62 glaucomatous eyes. The amplitudes of the oxygenated and deoxygenated haemoglobin concentrations, together with their latencies with respect to the stimulus onset, were estimated by fitting their time course with a canonical HRF. Statistical analysis showed that the amplitudes of both haemodynamic parameters show a significant association with the pathology and a significant discriminating ability, while no significant result was found for latencies. Overall, our findings together with the ease of use and noninvasiveness of TD-NIRS, make this technique a promising candidate as a supporting tool for a better evaluation of the glaucoma pathology.