Building upon our recent study [G. S. Bali et al., Eur. Phys. J. C 78, 217 (2018)], we investigate the feasibility of calculating the pion distribution amplitude from suitably chosen Euclidean correlation functions at large momentum. We demonstrate in this work the advantage of analyzing several correlation functions simultaneously and extracting the pion distribution amplitude from a global fit. This approach also allows us to study higher-twist corrections, which are a major source of systematic error. Our result for the higher-twist parameter delta(pi)(2) is in good agreement with estimates from QCD sum rules. Another novel element is the use of all-to-all propagators, calculated using stochastic estimators, which enables an additional volume average of the correlation functions, thereby reducing statistical errors.