Cet article présente une implémentation efficace de l'approximation de Vecchia sur GPU afin de réduire la complexité de calcul des processus gaussiens. Contrairement aux logiciels multicœurs existants, nous nous concentrons sur l'amélioration de la vitesse et de la précision de l'analyse de données statistiques spatiales à grande échelle grâce aux GPU. Nous comparons et analysons trois implémentations GPU, étudions l'influence du type de mémoire et proposons une nouvelle méthode optimisée. Nous développons le package R GpGpU incluant la nouvelle méthode et comparons ses performances à celles des logiciels existants. Nous démontrons également la supériorité de GpGpU par des expériences sur divers jeux de données, notamment des données satellitaires d'observation de la Terre à grande échelle d'une taille de $n>10^6$.
