After the time series in the preprocessed functional data were transformed onto the six laminar surface reconstructions, we performed vertex by vertex (a vertex is a point on a triangle surface as explained before, see Fig. 1c) blind deconvolution [58] to get each vertex’s latent neural response and HRF.
