Travel modes are generally derived from Global Positioning System (GPS) data on the basis of either a rule-based or machine learning classification method. The rule-based classification approach is generally easy to understand, whereas the machine learning classification method has better generalization. However, studies that jointly explore both methods are limited. The present research proposes a two-stage method that aims to impute travel modes from GPS trajectory data. In the first stage, rules are employed to detect subway modes. In the second stage, a Gaussian process classifier based on sequential forward selection methods is utilized to derive the remaining travel modes. On the basis of the four selected features constituting the feature set (i.e., average speed, average acceleration, heading change, and low-speed point rate), over 97% of the samples with subway modes are correctly identified and 93.04% of segments in the walk-based balanced test subset are accurately detected. Over 92% of the car and bus samples are correctly detected for the training and test datasets. Results provide a new perspective in selecting classification methods for the detection of travel modes and other travel characteristics from GPS trajectory data. Furthermore, high differentiation is achieved between the bus and car modes without the bus transit geographic information system sources of bus networks. Therefore, reasonable extracted features contribute to the detection of travel modes, particularly between bus and car modes. INDEX TERMS Gaussian processes, classification algorithms, global positioning system.