Background: The development of Brain-Computer Interface (BCI) technology has brought tremendous potential to various fields. In recent years, prominent research has focused on enhancing the accuracy of BCI decoding algorithms by effectively utilizing meaningful features extracted from electroencephalographic (EEG) signals. Objective: This paper proposes a method for extracting brain functional network features based on directed transfer function (DTF) and graph theory. The method incorporates the extracted brain network features with common spatial pattern (CSP) to enhance the performance of motor imagery (MI) classification task. Methods: The signals from each electrode of the EEG, utilizing a total of 32 channels, are used as input signals for the network nodes. In this study, 26 healthy participants were recruited to provide EEG data. The brain functional network is constructed in Alpha and Beta bands using the DTF method. The node degree (ND), clustering coefficient (CC), and global efficiency (GE) of the brain functional network are obtained using graph theory. The DTF network features and graph theory are combined with the traditional signal processing method, the CSP algorithm. The redundant network features are filtered out using the Lasso method, and finally, the fused features are classified using a support vector machine (SVM), culminating in a novel approach we have termed CDGL. Results: For Beta frequency band, with 8 electrodes, the proposed CDGL method achieved an accuracy of 89.13%, a sensitivity of 90.15%, and a specificity of 88.10%, which are 14.10, 16.69, and 11.50% percentage higher than the traditional CSP method (75.03, 73.46, and 76.60%), respectively. Furthermore, the results obtained with 8 channels were superior to those with 4 channels (82.31, 83.35, and 81.74%), and the result for the Beta frequency band were better than those for the Alpha frequency band (87.42, 87.48, and 87.36%). Similar results were also obtained on two public datasets, where the CDGL algorithm's performance was found to be optimal. Conclusion: The feature fusion of DTF network and graph theory features enhanced CSP algorithm's performance in MI task classification. Increasing the number of channels allows for more EEG signal feature information, enhancing the model's sensitivity and discriminative ability toward specific activities in brain regions. It should be noted that the functional brain network features in the Beta band exhibit superior performance improvement for the algorithm compared to those in the Alpha band. [ABSTRACT FROM AUTHOR]