V1 is ideally suited to study the spatial organization of neurovascular coupling at the level of synapses, neurons, individual blood vessels and laminar-resolution fMRI. This is because at least in layers 1 and 2/3 of V1, the functional micro-architecture for neurons, synapses and blood vessels has been determined using 2-photon imaging (Ohki et al 2005 Nature; Kara and Boyd 2009 Nature; O'Herron et al 2016 Nature). Hence, feature selectivity, e.g., orientation and direction selectivity of spiking, synaptic and hemodynamic activity in layer 2/3 is known. However, the micro-architecture of layer 4 neural activity (spiking and synaptic) along with individual blood vessel responses is unknown because conventional 2-photon imaging cannot access deeper cortical layers. The organizing principles of neural maps and the selectivity of hemodynamic responses is of paramount importance for laminar processing because the thalamic inputs arriving into layer 4 are untuned. 3-photon imaging triples the imaging depth compared to 2-photon imaging. Using this optical technique and high-resolution fMRI, we have determined the extent to which different types of neural (spiking, synaptic) and vascular signals (blood flow from individual vessels and fMRI voxels) are coupled across cortical layers in the cat V1. Our data show systematic changes in selectivity of hemodynamic signals across cortical layers that have clear underpinnings in neural circuitry and the propagation of hemodynamic signals.