The recently proposed permutationally invariant polynomial-neural network (PIP-NN) based ∆-machine learning (∆-ML) approach (PIP-NN ∆-ML, J. Phys. Chem. Lett. 2022, 13, 4729) is a flexible, general, and highly cost-efficient method to develop full dimensional accurate potential energy surface (PES). Only a small portion of points, which can be actively selected from the low-level (often DFT) points, with high-level energies are needed to bring a low-level PES to a high-level of quality. The hydrogen abstraction reaction between methanol and hydroxyl radical, OH + CH3OH, has been studied by theories and experiments for a long time due to its great importance in combustion, atmospheric and interstellar chemistry. However, it is not trivial to develop the full dimensional accurate PES for it. In this work, the PIP-NN ∆-ML method is successfully applied to the title reaction. The DFT PES was fitted by using 140192 points. Only 5% of the DFT dataset was needed to be calculated at the level of UCCSD(T)-F12a/AVTZ, and to improve the DFT PES to the target high-level, UCCSD(T)-F12a/AVTZ. More than 92% of the original unaffordable calculation cost were saved. The kinetics, including rate coefficients and branching ratios, were then studied by performing quasi-classical trajectory calculations on this newly fitted PES for the title reaction.