FR 3125197 A1 20230113 FR 2107437 A 20210708 Le présent document un procédé de détermination d’un ensemble de trajectoires d’un utilisateur portant un appareil mobile ( ) dans un réseau de téléphonie mobile, sur un territoire donné, comprenant une pluralité d’antennes ( ), le procédé consistant essentiellement à collecter des données de signalisation ( ) dudit appareil mobile ( ), chaque donnée de signalisation comprenant un horodatage ( ) associé à des coordonnées ( ) d’une antenne communiquant avec ledit appareil mobile, et déterminer au moins une série de données de signalisation ( ) successives associées à une seule antenne commune, et calculer un temps cumulé pour chaque série de données de signalisation, et déterminer l’ensemble de trajectoires à partir dudit ensemble de séries mobiles ( ). Figure de l’abrégé : Figure 2 Détermination de trajectoires à partir de données de téléphonie mobile Le présent document se rapporte au domaine de détermination de trajectoires à partir de données de téléphonie mobile. Plus précisément, la présente invention se rapporte à un procédé de détermination d’un ensemble de trajectoires d’un utilisateur portant un appareil de téléphonie mobile, à un dispositif de traitement de données permettant la mise en œuvre de ce procédé, à un programme d’ordinateur comprenant des instructions qui permettent la mise en œuvre du procédé et à un support d’enregistrement lisible par ordinateur comprenant des instructions qui permettent à un ordinateur de mettre en œuvre le procédé. Afin de comprendre et d’analyser la mobilité de la population en zone urbaine, les enquêtes sur les ménages ont été la seule source de données pour déterminer les trajectoires de déplacement des individus de la population. D’autres sources de données existent aujourd’hui telles que les cartes à puce, le système de positionnement global (GPS), les réseaux sociaux géolocalisés ou les enregistrements de téléphones mobiles. Parmi les nouvelles sources de données sur la mobilité, les données de signalisation des appareils mobiles émis par les opérateurs de réseaux mobiles constituent une option privilégiée pour suivre la trajectoire des appareils mobiles connectés à ces réseaux mobiles. En effet, ces données de signalisation sont disponibles pour l'ensemble des abonnés d'un fournisseur de réseau, ce qui couvre généralement un nombre important d’utilisateurs d’environ des dizaines ou des centaines de millions d'utilisateurs. Ces données de signalisation sont collectées en continu sur de longues périodes, ce qui permet des études de mobilité sur des mois ou des années. En outre, ces données de signalisation sont collectées de façon passives dans le cadre de la facturation des utilisateurs, et ne nécessitent donc pas un dispositif particulier pour les collecter. Cependant, ces données de signalisation manquent de précision en termes de localisation géographique et temporelle. En effet, l'emplacement de l’appareil mobile à partir des données de signalisation ne peut être mis en correspondance qu'avec la zone de couverture dudit appareil mobile ou la position des antennes avec lesquelles il communique. En outre, la fréquence d’acquisition des données est conditionnée par la fréquence de communication de l’appareil mobile soit des appels téléphoniques, l’échange de message texte ou de données internet, etc. qui restent des évènements sporadiques. Enfin, le nombre de données de signalisation collecté est important et nécessite des ressources de calcul importantes. En conséquence, les données de signalisation présentent une source prometteuse pour suivre la mobilité des personnes en zone urbaine mais il existe un besoin de méthodes pour déduire la localisation des personnes de façon précise et fiable à partir de ces données de signalisation. Résumé A cet effet, selon un premier aspect, le présent document propose un procédé de détermination d’un ensemble de trajectoires d’un utilisateur portant un appareil mobile dans un réseau de téléphonie mobile, sur un territoire donné, comprenant une pluralité d’antennes, le procédé consistant essentiellement à : (a) collecter des données de signalisation dudit appareil mobile, chaque donnée de signalisation comprenant un horodatage associé à des coordonnées d’une antenne communiquant avec ledit appareil mobile, (b) déterminer au moins une série de données de signalisation successives associées à une seule antenne commune, et calculer un temps cumulé pour chaque série de données de signalisation, (c) pour chaque série de données de signalisation, lorsque le temps cumulé est supérieur à un premier seuil temporel, labeliser ladite série de données de signalisation comme étant une série statique, (d) comparer deux séries statiques successives, et lorsqu’une première série statique et une seconde série statique successives sont associées à une même antenne, et lorsque des données de signalisation intermédiaires situées entre la première série statique et la seconde série statique sont associées à d’autres antennes, différentes de l’antenne en commun ; et que le nombre des antennes associées aux données de signalisation intermédiaires est inférieur à un seuil d’antennes, fusionner la première série statique et la seconde série statique en une seule série statique, (f) construire un ensemble de séries mobiles comprenant les données de signalisation exclues des séries statiques, et (g) déterminer l’ensemble de trajectoires à partir dudit ensemble de séries mobiles. Les étapes (c) et (d) permet de détecter les sessions statiques et mobiles de l’utilisateur. En particulier, l’étape (d) permet de détecter et de lisser les données de signalisation issues des phénomènes d’oscillation, ces phénomènes sont inhérents aux communications de téléphonie mobile. En effet, l’affectation d’un appareil mobile à une antenne est régie selon des stratégies propres au fournisseur de téléphonie mobile et se base sur des mesures dynamiques telles que la puissance du signal émis par l’appareil mobile ou le niveau de disponibilité des antennes environnantes. Ceci peut résulter dans la bascule de l’appareil entre deux ou plusieurs antennes et induit donc les phénomènes d’oscillation. Ainsi, les inventeurs ont constaté le procédé permet de détecter les périodes de mobilités où l’utilisateur se déplace réellement, avec une grande fiabilité. En effet, les inventeurs ont constaté la détermination des périodes de mobilités avec une précision de 80% et un rappel de 96%. Les données de signalisation de l’appareil mobile peuvent être collectées pendant une durée de temps d’un ou plusieurs jours. A l’étape (d), les données de signalisation intermédiaires peuvent être supprimées ou stockées dans une série distincte des séries statiques et des séries mobiles. Les données de signalisation de l’appareil mobile peuvent être collectées pour différents évènements de communication de l’appareil mobile, telle que : - les communications vocales ou l’échange de messages textes SMS ( « Short Message Service » en anglais), - les transferts intercellulaires, - les mises à jour des zones de localisation (LA, « localisation area » en anglais) et des zones de suivi (TA, « tracking area » en anglais), c'est-à-dire les mises à jours lors de changements de cellule ou d’antenne d’une zone LA, généralement couverte par un réseau 2G/3G et une zone TA, généralement couverte par un réseau 4G, - les demandes périodiques émises par le réseau de téléphonie mobile pour mettre à jour l'emplacement de l’appareil mobile, - les connexions et/ou les déconnexions du réseau de téléphonie mobile, c'est-à-dire lorsqu’un appareil mobile qui rejoint ou quitte le réseau de téléphonie mobile à son activation ou sa désactivation, - les échanges de données internet, qui comprennent entre autre les requêtes d'affectation de ressources internet à des applications mobiles en cours d’exécution par l’appareil mobile. Selon un mode de réalisation, l’étape (g) peut comprendre les sous-étapes suivantes : (g1) construire un ensemble de trajectoires récurrentes comprenant les séries mobiles correspondant à une même trajectoire, (g2) déterminer une durée moyenne de ladite trajectoire à partir de l’ensemble de trajectoires récurrentes, (g3) pour chaque série mobile déterminée à l’étape (g1), mettre temporellement à l’échelle ladite série mobile pour avoir une durée totale égale à la durée moyenne, (g4) pour chaque série mobile mise à l’échelle à l’étape (g3), échantillonner ladite série mobile selon un pas de temps déterminé. La détermination des trajectoires récurrentes permet d’enrichir temporellement les données d’une même trajectoire et de corriger les biais spatiaux. Les inventeurs ont constaté que le procédé permet de déduire les trajectoires de l’utilisateur avec une précision de l’ordre de 190 m. Le pas de temps déterminé peut être compris entre 0.5 et 2 minutes. L’étape (g1) peut comprendre les sous-étapes suivantes : (g11) pour chaque paire de séries mobiles, déterminer la distance de Hausdorff, entre lesdites séries mobiles, (g12) déterminer l’ensemble de trajectoires récurrentes par partitionnement des distances de Hausdorff obtenues à l’étape (g11). L’étape (g) peut comprendre en outre après la sous-étape (g2), les sous-étapes suivantes : (g21) calculer la médiane des durées de trajectoires de l’ensemble des trajectoires récurrentes, (g22) filtrer de l’ensemble de trajectoires récurrentes, les séries mobiles présentant une durée de trajectoire supérieure à 50% de ladite médiane calculée. L’étape (g12) peut être réalisée par l’algorithme DBSCAN (« density-based spatial clustering of applications with noise » en anglais). Chaque trajectoire de l’ensemble de trajectoires peut comprendre les moyennes des coordonnées de l’antenne des données de signalisation à chaque pas de temps déterminé. Le procédé peut comprendre après l’étape (d), pour chaque série statique, une étape (e) consistant à retirer le label statique de ladite série statique lorsque le temps cumulé pour cette série statique est inférieur à un second seuil temporel. Le procédé peut comprendre après l’étape (e), une étape consistant à déterminer une position statique de l’utilisateur, ladite étape comprenant : - calculer, pour chaque série statique, un centroïde des coordonnées des antennes associées à des données de signalisation compris dans ladite série statique, - partitionnement des centroïdes obtenus en plusieurs partitions comprenant les centroïdes tels que partitionnés et les séries statiques associées auxdits centroïdes, - calculer, pour chaque partition, un barycentre des centroïdes, - affecter à chaque série statique une position statique, ladite position statique étant le barycentre de la partition comprenant ladite série statique. Le premier seuil temporel peut être compris entre 5 et 60 minutes, en particulier égal à 20 minutes. Le second seuil temporel peut être compris entre 5 et 60 minutes, en particulier égal à 20 minutes. La valeur du second seuil temporel peut être égale à la valeur de premier seuil temporel. Le seuil d’antennes peut être compris entre 1 et 10, en particulier égal à 2. Selon un autre aspect, le présent document concerne un procédé de caractérisation de la mobilité d’une population d’une agglomération, le procédé comprenant les étapes suivantes : (1) déterminer des ensembles de trajectoires selon l’une quelconque des revendications précédentes pour une pluralité d’utilisateurs, (2) déterminer des trajectoires communes pour un nombre déterminé d’utilisateurs. Le procédé peut comprendre la superposition d’une ou plusieurs trajectoires communes à une carte du territoire donné. Selon un autre aspect de l’invention, il est proposé un dispositif de traitement de données comprenant un contrôleur centralisé comportant une mémoire et permettant la mise en œuvre du procédé selon l’invention. Selon un autre aspect de l’invention, il est proposé un programme d’ordinateur comprenant des instructions qui, dès lors qu’elles sont chargées et implémentées sur un ordinateur, permettent la mise en œuvre du procédé selon l’invention. Selon un autre aspect de l’invention, il est proposé un support d’enregistrement lisible par ordinateur comprenant des instructions qui, lorsqu’elles sont exécutées par un ordinateur, conduisent celui-ci à mettre en œuvre le procédé selon l’invention. D’autres caractéristiques, détails et avantages apparaîtront à la lecture de la description détaillée ci-après, et à l’analyse des dessins annexés, sur lesquels : La illustre schématiquement une trajectoire, obtenue par géolocalisation, d’un utilisateur portant un appareil mobile connecté à un réseau de téléphonie mobile ; La montre un schéma bloc d’un exemple des étapes de traitement des données de signalisation de l’appareil mobile pour déterminer une trajectoire de l’utilisateur portant l’appareil mobile ; La montre un schéma bloc d’un exemple des étapes de détermination de la trajectoire de l’utilisateur à partir des données de signal après leur traitement par les étapes de la ; La figure 4a illustre les données de signalisation brutes pour un utilisateur et ces les données de signalisation après leur traitement par les étapes de la , La montre un ensemble de trajectoires récurrentes obtenues par les étapes de la . La montre une trajectoire de l’utilisateur déterminée selon les étapes de la en comparaison avec la trajectoire réelle de l’utilisateur obtenue par géolocalisation. Procédé de détermination d’un ensemble de trajectoires d’un utilisateur portant un appareil mobile ( ) dans un réseau de téléphonie mobile, sur un territoire donné, comprenant une pluralité d’antennes ( ), le procédé consistant essentiellement à : (a) collecter des données de signalisation ( ) dudit appareil mobile ( ), chaque donnée de signalisation comprenant un horodatage ( ) associé à des coordonnées ( ) d’une antenne communiquant avec ledit appareil mobile, (b) déterminer au moins une série de données de signalisation ( ) successives associées à une seule antenne commune, et calculer un temps cumulé pour chaque série de données de signalisation, (c) pour chaque série de données de signalisation ( ), lorsque le temps cumulé est supérieur à un premier seuil temporel ( ), labeliser ladite série de données de signalisation comme étant une série statique ( ), (d) comparer deux séries statiques successives, et lorsqu’une première série statique ( ) et une seconde série statique ( ) successives sont associées à une même antenne ( ), et lorsque des données de signalisation ( ) intermédiaires situées entre la première série statique ( ) et la seconde série statique ( ) sont associées à d’autres antennes, différentes de l’antenne en commun ; et que le nombre des antennes associées aux données de signalisation intermédiaires est inférieur à un seuil d’antennes ( ), fusionner la première série statique ( ) et la seconde série statique ( ) en une seule série statique ( ), (f) construire un ensemble ( ) de séries mobiles ( ) comprenant les données de signalisation ( ) exclues des séries statiques ( ), et (g) déterminer l’ensemble ( ) de trajectoires à partir dudit ensemble ( ) de séries mobiles ( ). Procédé selon la revendication 1, dans lequel l’étape (g) comprend les sous-étapes suivantes : (g1) construire un ensemble ( ) de trajectoires récurrentes comprenant les séries mobiles correspondant à une même trajectoire, (g2) déterminer une durée moyenne de ladite trajectoire à partir de l’ensemble ( ) de trajectoires récurrentes, (g3) pour chaque série mobile déterminée à l’étape (g1), mettre temporellement à l’échelle ladite série mobile pour avoir une durée totale égale à la durée moyenne, (g4) pour chaque série mobile mise à l’échelle à l’étape (g3), échantillonner ladite série mobile selon un pas de temps déterminé. Procédé selon la revendication 2, dans lequel l’étape (g1) comprend les sous-étapes suivantes : (g11) pour chaque paire de séries mobiles ( , ), déterminer la distance de Hausdorff ( ( , )), entre lesdites séries mobiles ( , ), (g12) déterminer l’ensemble ( ) de trajectoires récurrentes par partitionnement des distances de Hausdorff ( ( , )) obtenues à l’étape (g11). Procédé selon la revendication 2 ou 3, dans lequel l’étape (g) comprend en outre après la sous-étape (g2), les sous-étapes suivantes : (g21) calculer la médiane des durées de trajectoires de l’ensemble ( ) des trajectoires récurrentes, (g22) filtrer de l’ensemble ( ) de trajectoires récurrentes, les séries mobiles présentant une durée de trajectoire supérieure à 50% de ladite médiane calculée. Procédé selon la revendication 2 à 4, dans lequel chaque trajectoire de l’ensemble ( ) de trajectoires comprend les moyennes des coordonnées ( ) de l’antenne des données de signalisation à chaque pas de temps déterminé. Procédé selon l’une des revendications précédentes, comprenant, après l’étape (d), pour chaque série statique ( ), une étape (e) consistant à retirer le label statique de ladite série statique lorsque le temps cumulé pour cette série statique est inférieur à un second seuil temporel ( ). Procédé selon la revendication précédente, comprenant, après l’étape (e), une étape consistant à déterminer une position statique de l’utilisateur (i), ladite étape comprenant : - calculer, pour chaque série statique ( ), un centroïde des coordonnées ( ) des antennes associées à des données de signalisation ( ) compris dans ladite série statique, - partitionnement des centroïdes obtenus en plusieurs partitions comprenant les centroïdes tels que partitionnés et les séries statiques associées auxdits centroïdes, - calculer, pour chaque partition, un barycentre des centroïdes, - affecter à chaque série statique une position statique, ladite position statique étant le barycentre de la partition comprenant ladite série statique. Procédé selon l’une des revendications précédentes en combinaison avec la revendication 6, dans lequel le premier seuil temporel ( ) et/ou le second seuil temporel ( ) est compris entre 5 et 60 minutes. Procédé selon l’une des revendications précédentes, dans lequel le seuil d’antennes ( ) est compris entre 1 et 10. Procédé selon la revendication 3, dans lequel l’étape (g12) est réalisée par l’algorithme DBSCAN (« density-based spatial clustering of applications with noise » en anglais). Procédé de caractérisation de la mobilité d’une population d’une agglomération, le procédé comprenant les étapes suivantes : (1) déterminer des ensembles ( ) de trajectoires selon l’une quelconque des revendications précédentes pour une pluralité d’utilisateurs, (2) déterminer des trajectoires communes pour un nombre déterminé d’utilisateurs ( ). Procédé selon la revendication précédente, comprenant la superposition d’une ou plusieurs trajectoires communes à une carte du territoire donné. Dispositif de traitement de données comprenant un contrôleur centralisé comportant une mémoire et permettant la mise en œuvre du procédé selon l’une des revendications 1 à 11 ou du procédé 11 ou 12. Programme d’ordinateur comprenant des instructions qui, dès lors qu’elles sont chargées et implémentées sur un ordinateur, permettent la mise en œuvre du procédé selon l’une au moins des revendications 1 à 11 ou du procédé 11 ou 12. Support d’enregistrement lisible par ordinateur comprenant des instructions qui, lorsqu’elles sont exécutées par un ordinateur, conduisent celui-ci à mettre en œuvre le procédé selon l’une au moins des revendications 1 à 11 ou du procédé 11 ou 12.