The farinograph characteristics of wheat flour determine the quality and the end use of wheat flour. The farinograph characteristics of wheat flour are influenced by wheat variety, origin, and milling process technology. There are four important farinograph parameters: water absorption, development time, stability time and degree of softening. Near-infrared spectroscopy (NIR) is widely used to determine wheat flour composition parameters, such as moisture, protein, ash and wet gluten content. Most of them directly use linear regression algorithms to establish models, which has low prediction accuracy, and there are few studies on detecting farinograph characteristics, and the results are also affected by the lack of sample richness. In this study, 968 samples of wheat flour from different countries and regions were collected, and an ensemble method of classification model and a regression model was proposed to improve the prediction accuracy of farinograph characteristics. Spectral preprocessing methods, including standard normal variation (SNV), linear detrending, multiplicative scatter correction (MSC) and Savitzky-Golay first-order derivative, were applied to the spectral data, and the best preprocessing method was selected with cross-validation. As for the modeling methods, the classical linear regression methods, i.e., partial least squares regression (PLSR) and principal component regression (PCR), were explored. The accuracies of the two methods are approximately equivalent. The root mean squared error of calibration (RMSEC) on farinograph parameters (i.e. water absorption, development time, stability time, and degree of softening) of the PCA model were 2.186, 1.838, 4.037, 21.693 and 2.039, 1.837, 3.968, 21.252 for PLSR correspondingly. The PLSR model requires fewer factors than PCR. Secondly, the two-stage regression model proposed in this paper was explored. Gaussian process regression (GPR) results were used as the classifier to cluster the samples, PLSR models were established in different clusters to predict the farinograph characteristics, and the sigmoid function was used to fuse the PLSR models. This modeling method can significantly improve the prediction accuracy of farinograph characteristics. The RMSEC on the predictions of farinograph parameters is 1.876, 1.160, 2.459 and 14.449 correspondingly.