Procédé de détermination d’un mode de transport (MDP) d’un trajet (traj) parcouru par au moins un utilisateur, au moyen d’un réseau (Res) de téléphonie mobile, du téléphone mobile (Tel) de l’utilisateur apte à se connecter à des antennes dudit réseau (Res) de téléphonie mobile et au moyen de graphes de transport (G_MDP) pour chaque mode de transport ou combinaison de mode de transport, chaque graphe de transport (G_MDP) comprenant des nœuds, des portions de routes reliant les différents nœuds et les vitesses moyennes de déplacement sur chaque portion de route, chaque portion de route du graphe de transport (G_MDP) représentant les routes possibles par le mode de transport associé, dans lequel on met en œuvre au moins les étapes suivantes : a) on acquiert (Acq) des données spatiales et temporelles (DSP) de connexion du téléphone (Tel) au réseau (Res) de téléphone mobile, les données spatiales et temporelles (DSP) comprenant les positions des antennes auxquels le téléphone (Tel) s’est connecté et les instants de début et de fin auxquels ont lieu ces connexions du téléphone aux antennes, on détermine une succession de couples emplacements-instants définis par des emplacements correspondant à certaines desdites positions des antennes issues des données spatiales et temporelles (DSP) acquises et par les instants acquis de début et/ou de fin correspondant à ces emplacements, et on forme (For) au moins une sous-route reliant deux couples emplacements-instants successifs dudit trajet (traj), une succession des sous-routes déterminant le trajet (traj) parcouru ; b) Pour chaque graphe de transport (G_MDP), on détermine (Det) les positions successives du téléphone de l’utilisateur sur des nœuds identifiés successifs (NIS) du graphe de transport (G_MDP) à partir des couples emplacements-instants, les nœuds identifiés successifs (NIS) étant positionnés le plus près possible de chaque emplacement ; c) Puis pour chaque graphe de transport (G_MDP) et pour chaque sous-route, on détermine (Opt), par des optimisations de plus court chemin, un nombre prédéterminé de chemins possibles (CP) permettant de relier les différents nœuds identifiés successifs (NIS) par des portions de route, chaque chemin possible (CP) comprenant les nœuds de passage reliant toutes les portions de route du chemin possible (CP), et on calcule l’instant de passage de chaque nœud de passage à partir des vitesses moyennes associées à chaque portion de route, chaque chemin possible (CP) étant déterminé par la liste de tous les couples nœuds de passage - instants de passages dudit chemin possible (CP) ; d) on détermine (Corr) un indice de corrélation (Ind) entre chacun des chemins possibles (CP) déterminés et le trajet (traj) parcouru déterminé, l’indice de corrélation (Ind) étant représentatif de la proximité spatiale et/ou de la proximité temporelle de chaque nœud identifié successif (NIS) avec les emplacements du trajet (traj) ; et e) on détermine (Det_MDP) le mode de transport (MDP) de chaque sous-route du trajet, par détermination du mode de transport qui optimise l’indice de corrélation (Ind) pour la sous-route du trajet. Procédé selon la revendication 1, pour lequel les données spatiales et temporelles (DSP) comprennent, de préférence exclusivement, des données d’enregistrement de transactions téléphoniques constituées d’un identifiant du téléphone (Tel) de l’utilisateur sur le réseau (Res), du type d’évènements de connexion du téléphone (Tel) au réseau (Res), un instant de début et un instant de chaque évènement de connexion du téléphone au réseau de téléphonie mobile et la position de l’antenne du réseau de téléphonie mobile à laquelle s’est connecté le téléphone. Procédé selon l’une des revendications précédentes, pour lequel on extrait (Det_SG) un sous-graphe (SG) de chaque graphe de transport (G_MDP), le sous-graphe (SG) étant une partie du graphe de transport limité à une largeur prédéterminée (lg) autour de chaque sous-route déterminée à l’étape a) et on utilise le sous-graphe (SG) au lieu du graphe de transport (G_MDP) pour chacune des étapes b) à e). Procédé selon l’une des revendications précédentes, dans lequel on identifie (S/D) des phases statiques de l’utilisateur à partir des données spatiales et temporelles (DSP) et on détermine des phases de mobilité (Mob) comme étant l’ensemble des phases qui ne sont pas des phases statiques, et à chaque phase de mobilité (Mob) on détermine un trajet (traj) selon l’étape a), et de préférence, une phase statique étant identifiée par le fait que pendant une durée prédéterminée, les positions des antennes qui ont été connectées au téléphone sont situées dans un espace restreint autour d’une antenne prédéfinie parmi ces antennes, l’espace restreint étant de préférence déterminé par l’ensemble des antennes voisines de l’antenne prédéfinie à un ordre prédéfini,. Procédé selon l’une des revendications précédentes, pour lequel à l’étape a), on forme le trajet (traj) depuis un point de départ à un instant de départ vers un point d’arrivée à un instant d’arrivée en passant de préférence par des points-étapes définis par des emplacements des couples emplacements-instants à des instants-étapes entre l’instant de départ et l’instant d’arrivée, les points de départ et d’arrivée correspondant aux emplacements des couples emplacements-positions, respectivement à l’instant de départ et à l’instant d’arrivée, l’instant de départ correspondant à l’instant de fin de la connexion de l’antenne au point de départ et l’instant d’arrivée correspondant à l’instant de début de la connexion de l’antenne au point d’arrivée. Procédé selon la revendication 5, dans lequel on détermine successivement les points-étapes : - en calculant la distance entre l’emplacement du point-étape précédent, le premier point-étape correspondant au point de départ, et à l’instant correspondant à cet emplacement et chaque position d’antenne des données spatiales et temporelles à laquelle le téléphone s’est connecté après l’instant du point-étape précédent, et - en calculant l’accroissement de cette distance au cours du temps, le point-étape étant déterminé à la première position d’antenne au cours du temps pour laquelle l’accroissement de la distance est inférieur à un seuil prédéfini. Procédé selon l’une des revendications précédentes, dans lequel, à l’étape b), on identifie (Det) les nœuds successifs (NIS) de la manière suivante : i) on identifie le nœud de départ le plus proche possible de l’emplacement du couple emplacement-instant du premier instant du trajet ; ii) puis pour chaque couple emplacement-instant suivant, on identifie le nœud le plus proche de l’emplacement de ce couple emplacement-instant qui soit relié par des portions de route au nœud identifié du couple emplacement-instant précédent et s’il n’y a pas de solution possible : a1) on réitère l’étape ii) pour le couple emplacement-instant précédent pour identifier le nœud suivant le plus proche ; a2) si l’étape a1) ne permet pas de trouver une solution, on revient à l’étape i) où on identifie le nœud suivant le plus proche l’emplacement du couple emplacement-instant du premier instant du trajet Procédé selon l’une des revendications précédentes, dans lequel, à l’étape c), on détermine le nombre prédéterminé de chemins possibles (CP) par une optimisation au plus court chemin de manière séquentielle et en affectant un coefficient de pénalité à la vitesse des portions de route utilisées par les chemins possibles (CP) issue d’une optimisation au plus court chemin précédente. Procédé selon l’une des revendications précédentes, dans lequel à l’étape d), on détermine (Corr) l’indice de corrélation (Ind) à partir de la composante de similarité spatiale J spat suivante : Où ) correspond aux emplacements des couples emplacements-instants à chaque instant ti et où correspond aux nœuds de passage du chemin possible à chaque instant tj, la fonction d mesurant la distance entre ) et et N sample correspondant aux nombres de nœuds de passage. Procédé selon l’une des revendications précédentes, dans lequel à l’étape d), on détermine (Corr) l’indice de corrélation (Ind) à partir de la composante de similarité de mouvement J mouv suivante : Où et Où ) correspond aux emplacements des couples emplacements-instants à chaque instant ti et où correspond aux nœuds de passage du chemin possible à chaque instant tj, la fonction d mesurant la distance entre ) et et F med étant la fonction de filtre médian. Procédé selon les revendications 9 et 10, dans lequel, l’indice de corrélation (Ind) est déterminé par la formule suivante : Où λ est compris entre 0 et 1. Procédé selon la revendication 11, dans lequel, à l’étape e), on détermine un score S par la formule suivante : Avec Où et sont respectivement l’espérance et l’écart-type de et on détermine le mode de transport maximise le score S. Procédé selon l’une des revendications précédentes, pour lequel on prétraite (Pre) les données spatiales et temporelles (DSP) acquises pour ne retenir que les données correspondant à celles du téléphone de l’utilisateur (DSP_t). Produit programme d’ordinateur téléchargeable depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou exécutable par un processeur ou un serveur, comprenant des instructions de code de programme pour la mise en œuvre du procédé selon l’une des revendications précédentes, lorsque ledit programme est exécuté sur un ordinateur ou un téléphone portable.