Neurovascular coupling represents the relationship between changes in neuronal activity and cerebral hemodynamics. Concurrent Electroencephalography (EEG) and functional near-infrared spectroscopy (fNIRS) recording and integration analysis has emerged as a promising multi-modal neuroimaging approach to study the neurovascular coupling as it provides complementary properties with regard to high temporal and moderate spatial resolution of brain activity. In this study we developed an EEG-informed-fNIRS analysis framework to investigate the neuro-correlate between neuronal activity and cerebral hemodynamics by identifying specific EEG rhythmic modulations which contribute to the improvement of the fNIRS-based general linear model (GLM) analysis. Specifically, frequency-specific regressors derived from EEG were used to construct design matrices to guide the GLM analysis of the fNIRS signals collected during a hand grasp task. Our results showed that the EEG-informed fNIRS GLM analysis, especially the alpha and beta band, revealed significantly higher sensitivity and specificity in localizing the task-evoked regions compared to the canonical boxcar model, demonstrating the strong correlations between hemodynamic response and EEG rhythmic modulations. Results also indicated that analysis based on the deoxygenated hemoglobin (HbR) signal slightly outperformed the oxygenated hemoglobin (HbO)-based analysis. The findings in our study not only validate the feasibility of enhancing fNIRS GLM analysis using simultaneously recorded EEG signals, but also provide a new perspective to study the neurovascular coupling of brain activity.