FR 3004841 A1 20141024 FR 1353661 A 20130422 L'invention concerne un procédé et un dispositif pour la gestion dynamique de la mobilité urbaine. L'invention est plus particulièrement, mais non exclusivement, destinée à la supervision et à la prévision du stationnement des véhicules automobiles dans une agglomération urbaine. L'invention est plus particulièrement adaptée à la supervision du 5 stationnement en voirie et s'adresse ainsi à une ville, ou une agglomération urbaine, dont la part de stationnement en voirie est majoritaire ou significative, ce qui, en France, correspond à la plupart des villes de moins de 100 000 habitants. En effet, dans une ville de plus grande importance le stationnement est majoritairement organisé en ouvrages. Cependant l'invention est également adaptée à la supervision du stationnement en ouvrage 10 et plus précisément aux situations de stationnement mixtes qui en constituent l'essentiel des applications. Dans tout le texte, les termes « stationnement en ouvrage » désignent le stationnement organisé en emplacements délimités dans un espace clos dont l'accès est contrôlé, qu'il s'agisse d'un parking aérien, souterrain ou de surface. Les termes 15 « stationnement en voirie » désignent un stationnement dans un espace ouvert, payant ou gratuit, sur des emplacements délimités ou non. Les termes « espace de stationnement » désignent l'ensemble de l'offre de stationnement dans un domaine géographique donné, combinant stationnement en ouvrage et stationnement en voirie, payant ou non. 20 Bien qu'il représente pour la collectivité une recette d'exploitation significative et qu'il influe directement sur les flux de véhicules dans une ville, le stationnement, et plus particulièrement le stationnement en voirie, est mal supervisé car techniquement complexe à prévoir et à identifier. La capacité de supervision du stationnement dans une ville vise au moins trois objectifs : 25 - permettre aux usagers de se diriger rapidement vers des places disponibles et ainsi limiter les flux redondants de véhicules ; - permettre aux autorités publiques d'adapter l'offre et la demande en fonction de diverses contraintes ; - optimiser les recettes d'exploitation et l'efficacité des équipes de surveillance de la 30 voirie publique. Le document EP 1 128 350 décrit un procédé, basé sur les indications des bornes de stationnement et sur des statistiques historiques, pour prévoir le taux d'occupation d'un secteur urbain en termes de stationnement, et ainsi fournir des indications aux usagers sur la disponibilité de places dans ce secteur. Ce procédé de d'art antérieur présente plusieurs limitations. Tout d'abord il n'y a pas d'adéquation directe entre le nombre de titres de stationnement délivrés par une borne et le taux d'occupation effectif du secteur concerné par ladite borne, ni entre le temps de stationnement payé et le temps de stationnement effectif d'un véhicule. De plus, ce système est limité au cas d'une voirie comprenant des emplacements matérialisés, lesquels emplacements ne sont pas toujours respectés. Finalement, rien n'empêche un usager de retirer son titre de stationnement dans une borne correspondant à un autre secteur. Ces difficultés sont partiellement résolues par ce procédé de l'art antérieur en corrigeant les données objectives, délivrées par les bornes, par des données statistiques issues de l'observation des comportements. Par suite, ce système est inadapté aux cas engendrés par des situations inhabituelles où son utilité serait la plus grande. En regard des trois besoins énoncés supra ce procédé de l'art antérieur ne résout que le premier, celui de l'usager. Le document EP 2 082 386 décrit un procédé de supervision du stationnement en voirie, plus particulièrement adapté au cas d'un stationnement sans emplacements délimités. Ce dispositif repose sur une interaction entre l'usager et un système de supervision, l'usager déclarant audit système, d'une part, la longueur de son véhicule et d'autre part, son choix de secteur de stationnement lorsqu'il y est arrivé. Ce procédé de l'art antérieur ne résout, en regard du précédent, que la problématique de la non-délimitation des places de stationnement et présente globalement les mêmes limitations. Ces procédés de l'art antérieur, qui tirent leurs informations de l'acquittement d'une taxe de stationnement par l'usager, sont délicats de mise en oeuvre en présence d'une pluralité de moyens de paiement. En effet il n'est pas rare qu'au sein d'une même ville, plusieurs moyens d'acquittement existent, tels que bornes, badges, carte de résidence, paiement par téléphone etc. De plus, si ces procédés de l'art antérieur présentent une certaine efficacité vis-à-vis du guidage ou de l'orientation de l'usager au sein d'un même secteur, ou de plusieurs secteurs proches les uns des autres, ils s'avèrent moins efficaces à l'échelle d'une ville entière. En effet, l'intérêt de l'usager de rejoindre un secteur à stationnement disponible, dépend du temps mis par cet usager à rejoindre ce secteur. Ainsi, un secteur présentant une disponibilité importante de places, risque d'être complet au moment où l'usager rejoint ledit secteur, si celui-ci a été recommandé à un nombre 5 important d'usagers capables de s'y rendre plus rapidement. À l'inverse, il est parfois plus judicieux d'orienter un usager vers un secteur à disponibilité plus faible mais dont la rotation est à même d'assurer une meilleure probabilité de stationnement. L'invention vise à résoudre les inconvénients de l'art antérieur et concerne à cette fin un procédé pour la supervision du stationnement des véhicules dans un domaine 10 géographique, lequel procédé comprend les étapes consistant à: a. obtenir l'état d'occupation instantané de l'espace de stationnement correspondant audit domaine géographique à un instant t; b. obtenir le taux de rotation instantané du stationnement sur une période T dans ledit espace de stationnement à l'instant t; 15 c. calculer à partir des données recueillies aux étapes a) et b) le taux d'occupation futur à un instant (t+ At) et sa probabilité ; d. délivrer une information lorsque la probabilité ou le taux d'occupation déterminés à l'étape c) sont inférieurs ou supérieurs à un seuil déterminé. Dans tout le texte les termes « taux de rotation » désignent, pour un domaine 20 géographique donné, le nombre moyen de véhicules occupant successivement un emplacement sur une période considérée. Ainsi, le taux d'occupation de l'espace de stationnement visé est déterminé de manière dynamique, ainsi que sa probabilité d'occurrence. Ces informations permettent notamment de contextualiser et de pondérer l'information d'occupation en fonction de 25 données géographiques, telles que le temps de trajet jusqu'au domaine géographique considéré, ou en fonction d'une condition événementielle particulière, ceci à des fins tant de supervision en temps réel que de simulation. L'information délivrée à l'étape d) est utilisable, après un formatage adapté, tant par les usagers que par les agents de surveillance de la voirie publique ou les autorités publiques dans le but d'optimiser 30 l'occupation du domaine géographique ou d'optimiser le résultat d'exploitation de l'espace de stationnement associé. L'invention est avantageusement mise en oeuvre selon les modes de réalisation exposés ci-après, lesquels sont à considérer individuellement ou selon toute combinaison techniquement opérante. Avantageusement, l'étape a) comprend les étapes consistant à: ai. obtenir le nombre de véhicules instantanément présents dans le domaine géographique ; aii, obtenir le nombre instantané de véhicules en transit dans ledit domaine géographique ; aiii. déterminer le nombre de véhicules stationnés dans le domaine géographique en soustrayant le nombre défini à l'étape au) du nombre défini à l'étape ai). Ainsi le procédé objet de l'invention permet de déterminer le taux d'occupation de l'espace de stationnement indépendamment de l'acquittement du paiement et du respect des espaces de stationnement. De plus, ce mode de réalisation est mis en oeuvre sans 15 moyen de détection individuel de stationnement à la place, ne nécessitant que des balises de comptage séquentiel. Ainsi, ce mode de réalisation du procédé objet de l'invention est aisément mis en oeuvre dans un espace urbain sans travaux de voirie important et avec un investissement réduit. Avantageusement l'étape b) comprend les étapes consistant à : 20 bi obtenir le taux d'occupation spatial du domaine géographique par le nombre de véhicules stationnés dans ledit domaine à l'instant t; bii, déterminer le taux d'occupation temporel dudit domaine sur une période T, par la moyenne du taux d'occupation spatial sur cette période ; biii. déterminer le taux de rotation dudit domaine géographique en moyennant, 25 sur la période T, la variation du taux d'occupation temporel déterminé à l'étape bii). Ainsi, le procédé objet de l'invention se fonde sur un taux de rotation constamment mis à jour, lequel est calculé, si besoin pour différentes périodes d'intégration. Avantageusement, l'étape c) est déterminée en fonction d'un temps de trajet At vers 30 le domaine géographique au moyen d'une métaheuristique d'optimisation. Ainsi l'information délivrée correspond, pour son destinataire, à un optimum à l'instant t. Le choix d'une métaheuristique d'optimisation permet de traiter rapidement le problème d'optimisation complexe. Selon un mode de réalisation particulier, la métaheuristique utilise une règle aléatoire 5 de transition proportionnelle, issue d'un algorithme de colonie de fourmis, et l'information délivrée à l'étape d) est construite à partir de la visibilité du domaine géographique au temps (t+At). Ce choix de l'algorithme d'optimisation s'avère particulièrement adapté au traitement du problème. L'utilisation de la visibilité du domaine géographique comme paramètre conditionnel d'affichage permet de n'afficher que des informations pertinentes 10 notamment à destination des usagers. Avantageusement l'étape au) comprend la détermination de la vitesse moyenne des véhicules sur une partie du domaine géographique et la détermination de la fréquence temporelle dominante d'entrée dans ladite partie du domaine à l'instant t. La demanderesse a en effet déterminé que la dérivée de la fonction de comptage des véhicules dans une 15 zone géographique donnée est une fonction à variance finie pouvant être décomposée en une somme de fonctions sinusoïdales. Ainsi, la détermination du nombre de véhicules simultanément en transit dans ladite zone géographique est grandement facilitée, en s'appuyant sur les fréquences de passage ayant la plus forte densité spectrale tout en améliorant la robustesse de l'algorithme. 20 Selon un mode de réalisation avantageux, le procédé objet de l'invention comporte les étapes consistant à : e. obtenir une grille tarifaire du stationnement dans le domaine géographique ; f. déterminer la recette théorique à l'instant t et pour une période T en fonction du taux d'occupation du domaine déterminé à l'étape bii). 25 Ainsi, le procédé objet de l'invention permet, selon ce mode de réalisation, de simuler les recettes générées en fonction des conditions d'occupation de l'espace selon la grille tarifaire sélectionnée. Le résultat de l'étape f) est obtenu indépendamment de l'acquittement effectif des taxes de stationnement dans le domaine géographique considéré. 30 Avantageusement, le procédé selon ce dernier mode de réalisation de l'invention, comprend les étapes consistant à : g. obtenir la recette effective du domaine correspondant à l'instant t et sur une période d'observation T; h. déterminer le taux d'acquittement du domaine géographique considéré sur une période donnée en comparant sur cette période le rapport entre les résultats instantanés des étapes g) et f). Ainsi, le procédé objet de l'invention permet de superviser et de contrôler le respect de l'acquittement des taxes de stationnement, notamment afin de maximiser les recettes et d'accroître l'efficacité du contrôle. Avantageusement, le procédé objet de l'invention comprend, selon ce dernier mode de réalisation, les étapes consistant à : j. combiner l'information obtenue à l'étape c) avec l'information obtenue à l'étape h) et déterminer le taux d'acquittement futur à l'instant (t+At) et sa probabilité ; k. délivrer une information si la probabilité ou le taux d'acquittement déterminés à l'étape j) sont inférieurs ou supérieurs à un seuil donné. Ainsi, le procédé objet de l'invention délivre une information géographique contextualisée du taux d'acquittement dans l'espace de stationnement pour anticiper les mesures à prendre. Avantageusement, l'étape j) est déterminée en fonction d'un temps de trajet At d'un agent de surveillance de la voirie publique vers le domaine géographique, au moyen d'une métaheuristique d'optimisation, utilisant une règle aléatoire de transition proportionnelle, issue d'un algorithme de colonie de fourmis, et où l'information délivrée à l'étape k) est construite à partir de la visibilité du domaine géographique au temps (t+At). Ainsi, le procédé objet de l'invention permet d'optimiser l'intervention des agents de contrôle en prenant en compte des problématiques complexes. Selon un mode de réalisation particulier, le procédé objet de l'invention comprend une étape consistant à : I. modifier la grille tarifaire en fonction du taux d'occupation déterminé à l'étape a) ou de la recette théorique déterminée à l'étape f). Ainsi, le procédé objet de l'invention permet de proposer une offre de gestion dynamique de la tarification afin, notamment, d'optimiser l'occupation de l'espace de stationnement par une offre incitative ou pénalisante. L'invention concerne également un dispositif pour la mise en oeuvre du procédé objet 5 de l'invention selon l'un quelconque de ses modes de réalisation, lequel dispositif comprend : i. un ensemble de capteurs apte à déterminer la présence et le stationnement d'un véhicule dans le domaine géographique, notamment par comptage ; ii. un agrégateur de données apte à récupérer les informations générées par 10 l'ensemble des capteurs et à mettre en oeuvre les étapes a) à c) du procédé ; iii. un référentiel comportant : un référentiel dit spatial comportant une description géographique du domaine, et de l'espace de stationnement ; un référentiel tarifaire comprenant les gammes tarifaires et leur 15 association avec l'espace de stationnement ; un référentiel réglementaire comprenant les dispositions réglementaires associées à l'espace de stationnement ; iv. un module, dit applicatif, apte à mettre en oeuvre l'une quelconque des étapes d) à I) du procédé. 20 Le dispositif objet de l'invention est ainsi adaptable et paramétrable autour d'un coeur commun, à toute situation urbaine et à tout objectif d'optimisation ou de prévision de la gestion de la mobilité. L'utilisation d'une architecture commune permet d'enrichir et d'améliorer l'ensemble des dispositifs installés en fonction des retours d'expérience particuliers relatifs à chaque installation. 25 Avantageusement, le dispositif objet de l'invention comporte : y. un module de mise à disposition de données formatées de l'agrégateur de données sur internet. Ainsi, une partie des applicatifs est réalisable par des parties tierces en interaction avec d'autres données, par exemple en matière de guidage des usagers au moyen de leur 30 terminal mobile de positionnement géographique. Selon un mode de réalisation avantageux, le référentiel du dispositif objet de l'invention comporte : - un référentiel de règles définissant les modalités d'intervention des agents de surveillance de la voirie. Ainsi, lesdites règles sont utilisées dans le cadre de l'optimisation sous contrainte de la tournée desdits agents. L'invention est exposée ci-après selon ses modes de réalisation préférés, nullement limitatifs, et en référence aux figures 1 à 9, dans lesquelles : - la figure 1 illustre selon un schéma de principe un exemple de réalisation de l'architecture du dispositif objet de l'invention ; - la figure 2 montre selon une vue de dessus un exemple de carte d'une zone de stationnement comprenant un tronçon élémentaire dans une agglomération ; - la figure 3 est un exemple de relevé de l'occupation d'un espace de stationnement en fonction du temps, comparé à la prévision délivrée par le module de comptage selon un exemple de réalisation du dispositif et du procédé objets de l'invention ; - la figure 4 représente la dérivée temporelle de la fonction d'interpolation de la figure 3; - la figure 5 est un exemple, en vue de dessus, d'une zone de stationnement comprenant deux tronçons élémentaires ; - la figure 6 représente un exemple d'évolution temporelle du taux d'occupation d'une zone de stationnement pour différentes périodes d'intégration ; - la figure 7 est un exemple d'évolution temporelle du taux de rotation d'une zone de stationnement pour différentes périodes d'intégration ; - la figure 8 montre un exemple d'évolution temporelle de la probabilité de trouver une place dans une zone de stationnement donnée ; - et la figure 9 est un logigramme d'un exemple de réalisation du procédé objet de l'invention. Figure 1, le procédé objet de l'invention est mis en oeuvre par ordinateur au moyen d'un dispositif comportant un agrégateur (110) de données. Ledit agrégateur comporte un 30 module (112), dit DIC, acronyme de « Détection, Identification, Comptage », apte à collecter et à interpréter les données issues d'un premier ensemble de capteurs installés dans le domaine géographique, pour réaliser un comptage des véhicules dans ledit domaine. Ces capteurs sont connus de l'art antérieur et consistent principalement en : - des capteurs (121) à la place, mécaniques ou magnétiques ; - des caméras (123) vidéo comportant, le cas échéant, des algorithmes embarqués de reconnaissance de forme ou de lecture des plaques minéralogiques ; - des capteurs (122) de comptage et de détection séquentielle ; - des applications déclaratives basées sur des terminaux mobiles (124), tels que les téléphones cellulaires des usagers. La liste précédente n'est pas exhaustive. Le module (112) DIC reçoit en entrée les informations issues de ces capteurs par l'intermédiaire de concentrateurs (131, 132, 133, 134), agrégeant l'état desdits capteurs et une information géographique de leur localisation dans l'espace de stationnement avec une granularité adaptée. Le module (112) DIC reçoit également en entrée des informations relatives à l'entrée et à la sortie d'un véhicule dans le domaine géographique visé. Selon le mode de réalisation, cette information est délivrée par un ensemble dédié de capteurs ou par certains des capteurs du premier ensemble. Selon un exemple de réalisation, le domaine géographique est agrégé en 7 niveaux : - la ville est l'agrégat de plus haut niveau en matière de gestion du stationnement ; - le quartier est un sous-ensemble homogène de la ville ; - le secteur est un sous-ensemble homogène du quartier correspondant de manière univoque à une liste finie d'horodateurs ; - la zone est un espace géographique continu présentant une unicité tarifaire et réglementaire, à titre d'exemple non limitatif, un stationnement en ouvrage est assimilé à une zone ; - le tronçon élémentaire est une portion de voie limitée par l'axe de la voie, les façades qui limitent un des côtés de la voie, et par deux extrémités situées à des intersections entre la voie et une autre voie ; - la place est un emplacement de stationnement qui, selon le cas, est délimité par un marquage au sol, ou, en cas de stationnement libre un découpage un segments arbitraires de l'espace de stationnement. Préférentiellement le comptage d'entrée et de sortie des véhicules est réalisé à l'échelle du tronçon élémentaire. Ainsi, un tronçon élémentaire ne comporte aucune intersection avec une autre voie, publique ou privée, ou avec un ouvrage de stationnement, public ou privé. L'agrégateur (110)de données comprend un module (111), dit transactionnel, destiné à déterminer l'état des encaissements en regard l'occupation de l'espace de stationnement. Ce module (111) transactionnel reçoit en entrée, d'une part, des informations de référence, permettant de déterminer les zones tarifaires de manière géographique et temporelle, et d'autre part des informations de paiement issues à la fois des horodateurs (141) et des moyens (142) de paiement dématérialisé. Le référentiel comprend 3 types d'informations : - un référentiel (161) spatial, REF SPA, qui permet le paramétrage de l'offre géographique de stationnement ; - un référentiel (162) tarifaire, REF TAR, qui permet de définir les gammes tarifaires associées à des zones géographiques et des plages horaires ; - un référentiel (163) réglementaire, REF REG, qui permet le paramétrage de l'ensemble des dispositions réglementaires associées aux différentes zones considérées. Le dispositif comporte des moyens (non représentés) de saisie permettant de paramétrer chacun de ces référentiels de sorte à créer, supprimer ou modifier tout ou partie de ces référentiels et d'y inclure des situations particulières, tels que des travaux, des manifestations ou tout type d'exception aux conditions nominales d'exploitation. Le module (111) tarifaire reçoit en entrée et de manière continue, les informations de paiement issues des horodateurs (141) et des moyens (142) de paiement dématérialisé, par des concentrateurs (151, 152) agrégeant les données d'acquittement, les données temporelles et les données géographiques. Le module (112) DIC et le module (111) transactionnel transmettent les informations recueillies des différents capteurs, après les avoir interprétées et formatées, à un calculateur (113) qui effectue le traitement de ces données relatif aux applications visées. 30 Un centre (190) de contrôle et de décision, permet la supervision et le pilotage du dispositif. Ainsi, le module (111) transactionnel et le module (112) DIC, construisent deux informations, indépendantes, relatives à l'occupation de l'espace de stationnement : l'une relative à l'occupation effective et l'autre relative à l'occupation payée. La comparaison de ces deux informations en temps réel, selon des protocoles spécifiques, programmés dans le calculateur (113), permet d'obtenir des informations sur l'état d'occupation effectif et l'état d'occupation licite de l'espace de stationnement, lesquelles informations sont utilisables par le centre (190) de contrôle pour différentes applications. Figure 2, selon un exemple de mise en oeuvre, un tronçon élémentaire (200), ici constitué par une rue, comporte une pluralité (220) d'emplacements de stationnement. Pour simplifier l'exposé, ladite rue est considérée en sens (201) unique. Lesdits emplacements sont définis dans le référentiel REF SPA du dispositif objet de l'invention. Selon un exemple de réalisation ledit référentiel REF SPA est modifié pour tenir compte de la neutralisation d'un ensemble (227) de places, neutralisation consécutive, par exemple, à des travaux. Ainsi, l'offre de stationnement dans le tronçon élémentaire (200) est mise à jour en temps réel dans le dispositif objet de l'invention. Selon un exemple de réalisation, les emplacements (220) sont équipés chacun d'un capteur de présence (non représenté) permettant d'en détecter l'occupation. Alternativement, l'occupation d'un emplacement est déterminée par une ou plusieurs caméras de vidéosurveillance (non représentées). Selon un mode de réalisation alternatif, le tronçon élémentaire (200) comporte un capteur (222) d'entrée et un capteur (223) de sortie distants d'une distance D. Ce mode de réalisation, utilisant un comptage séquentiel, est particulièrement avantageux car simple d'installation, ne nécessitant que des compteurs en entrée et en sortie du tronçon élémentaire. Selon ce mode de réalisation deux capteurs renseignent un compteur, programmé dans le module DIC, lequel compteur est incrémenté d'une unité à chaque entrée de véhicule et décrémenté à chaque sortie. Dans le cas où la rue n'est pas en sens unique, les capteurs permettent de détecter le sens de passage du véhicule ou sont installés par demi-chaussée. Ce type de capteur est connu de l'art antérieur et consiste par exemple en une barrière optique, un dispositif mécanique sensible au poids du véhicule ou un radar à effet Doppler, sans que cette liste ne soit limitative. Un véhicule entrant dans le tronçon élémentaire (200) pour y stationner ne sort dudit tronçon qu'après un temps T STAT très supérieur au temps de transit dans ledit tronçon. La fonction de comptage s'écrit : CoiJritclge _ Eiu`;.ées_TEff, orties_ REFTEj Où T est un temps inférieur ou égal à T STAT déterminé par l'expérience ou par 5 observation. Typiquement T STAT est égal à 30 minutes et T est égal à 10 minutes. REF TE] correspond au nombre moyen de véhicules en transit simultané dans le tronçon élémentaire considéré. Figure 3, l'observation du nombre (302) de véhicules stationnés en fonction du temps (301) sur une durée de 12 heures met en évidence la similitude de la forme de la fonction 10 (312) Comptage TE] (t) comparée à l'occupation (311) réelle d'un tronçon élémentaire comportant 20 emplacements de stationnement. Ainsi, même en l'absence de capteurs de présence à la place, le paramétrage de la fonction Comptage TEj (t), par exemple au moyen d'observations de terrain, permet d'obtenir une image fidèle de l'occupation de l'espace de stationnement. Ainsi, selon ce mode de réalisation, le procédé objet de l'invention est 15 avantageusement mis en oeuvre avec un volume réduit d'installations sur le domaine public, ce qui limite le coût de mise en oeuvre. Le modèle présenté ci-avant est réalisé avec une valeur fixe du paramètre REF TEj . Avantageusement la valeur de REF TE] est déterminée dynamiquement. Ainsi, le temps de passage T LECT (i) d'un véhicule i devant l'un des capteurs (222, 223) à l'entrée 20 ou à la sortie de la zone permet de déterminer la vitesse du véhicule connaissant la longueur dudit véhicule. En prenant la longueur, L, statistique moyenne des véhicules, par exemple 5,3 mètres pour un passage comprenant des automobiles, des véhicules utilitaires et des poids lourds, la moyenne de ce temps sur N véhicules donne une bonne approximation de la vitesse moyenne de transit. vitesse = Figure 4, l'observation de la dérivée temporelle (412), Comptage TE" (t), de la fonction de comptage (312 figure 3) montre qu'il s'agit d'une fonction à variance finie qui, par suite, est décomposable en une somme de fonctions sinusoïdales. Coiî?ptage_ a . COS(.0 ç La transformée de Fourrier de la fonction d'auto covariance de cette fonction permet 5 d'identifier, à partir de la densité spectrale, la fréquence f0 qui correspond à la fréquence dominante d'entrée dans la zone considérée. Ainsi il vient : R_EF TE: D vitesse où D est la longueur du tronçon de voirie concerné. Ainsi, la fonction de comptage est auto-ajustée en fonction de la mesure instantanée du trafic. Comparativement aux solutions de l'art antérieur, la détermination du nombre de véhicules stationnés dans un 10 tronçon élémentaire et, par suite, dans un domaine géographique donné est indépendante de toute information relative à l'acquittement d'une taxe de stationnement. Calculée à partir de ces données, la quantité de véhicules en stationnement dans la zone tient compte à la fois des véhicules en stationnement régulier mais aussi des véhicules en stationnement n'ayant pas acquitté leur titre ou en ayant dépassé l'échéance, ainsi que des véhicules 15 stationnés en double file, attendant ou non la libération d'une place, et des véhicules, dits ventouse, qui occupent un emplacement pendant un temps exagérément long selon la législation en vigueur définie dans le référentiel REF REG. Ainsi, le procédé et le dispositif objets de l'invention, permettent, en comparant le taux d'occupation réel avec la capacité théorique du domaine géographique considéré de détecter, par exemple, la probabilité de 20 la présence d'un stationnement en double file. Figure 9, selon ce mode de réalisation, le procédé objet de l'invention comporte une première étape (910) visant la détermination de l'état d'occupation de l'espace de stationnement considéré, laquelle étape (910) comporte une première sous-étape (911) consistant à déterminer le nombre de véhicules présents dans ledit espace de stationnement, par exemple, au moyen d'un comptage séquentiel incrémentant les entrées et décrémentant les sorties. Selon une deuxième sous-étape (912) de correction, ledit comptage séquentiel est corrigé du nombre de véhicules en transit, déterminé de manière statique ou de manière dynamique. Les informations issues de ces sous-étapes (911,912) sont combinées dans une troisième sous-étape (913) pour déterminer le nombre de véhicules dans l'espace de stationnement considéré. Cette étape (910) de détermination de l'état d'occupation de l'espace de stationnement est mise en oeuvre par le calculateur (113) de l'agrégateur (110) de données, à partir des informations tirées du module (112) DIC. La fonction de comptage exposée ci-avant permet de calculer le taux d'occupation spatial instantané de l'espace de stationnement considéré. Ainsi le taux d'occupation spatial instantané, G T (t) , d'un tronçon élémentaire est donné, par exemple, par la relation : Où Ni correspond au nombre de places dans le tronçon élémentaire considéré. Ce taux d'occupation spatial compare le nombre de places disponibles dans le tronçon 15 considéré avec le nombre de véhicules en stationnement dans ledit tronçon. Le taux d'occupation temporel, OET (TEi, t), ou périodique, de cet espace à l'instant t et sur une période T, est donné par : -f Le taux d'occupation, temporel correspond ainsi au taux d'occupation moyen sur une période T. 20 Le taux de rotation pT (TEi,t), d'un tronçon élémentaire TEi, à l'instant t et sur une période de référence T, est égal à la moitié du nombre de changements d'état du taux d'occupation temporel dudit tronçon élémentaire sur la période (t-T, t) rapporté au nombre de places Ni du tronçon, soit : Cotage c 15 où G 'TE (x) est la dérivée temporelle du taux d'occupation spatial considéré et donne ainsi la variation de ce taux d'occupation spatial en fonction du temps. Le taux de rotation donne sur la période Tle nombre moyen de véhicules différents, occupant successivement une même place de stationnement dans le tronçon élémentaire considéré. Figure 9, selon un exemple de réalisation, la deuxième étape (920) du procédé objet de l'invention est essentiellement une étape de calcul, mise en oeuvre par le calculateur (113) de l'agrégateur (110) de données, à partir du résultat de l'étape (910) précédente et des données du référentiel (161) spatial REF SPA en ce qui concerne le nombre de places de l'espace de stationnement considéré. Cette deuxième (920) étape de calcul comporte, selon cet exemple de réalisation, une première (921) sous-étape consistant à déterminer le taux d'occupation spatial de l'espace de stationnement considéré. La deuxième (922) sous-étape de l'étape (920) de calcul, consiste à déterminer le taux d'occupation temporel de l'espace de stationnement, cette sous-étape comprend avantageusement la création d'une table de résultats correspondant à différentes périodes d'observation, ou d'intégration d'un point de vue mathématique. Une troisième (923) sous-étape de l'étape (920) de calcul permet de déterminer le taux de rotation de l'espace de stationnement considéré. Avantageusement, cette sous-étape (923) comporte le calcul d'une table de taux de rotation pour différentes périodes d'observation. Figure 5, une zone de stationnement est constituée de plusieurs (501, 502) tronçons élémentaires. Selon un exemple de réalisation, chaque tronçon (501, 502) élémentaire, est délimité par des balises (521,522, 523, 524) aptes détecter l'entrée-sortie d'un véhicule dans chacun des tronçons. Une même zone comprend, selon le mode de mise en oeuvre, des tronçons élémentaires équipés de capteurs de détection à la place et des tronçons élémentaires équipés d'un comptage séquentiel. Ainsi le taux d'occupation instantané, 04 (Zi,t), d'une zone Zi comportant plusieurs tronçons élémentaires est donné par : 3004 84 1 16 où Nia désigne le nombre de places instrumentées avec des capteurs de détection et Nib le nombre de places non instrumentées et dont l'occupation est déterminée par un comptage séquentiel comme explicité supra. Ainsi, le taux de rotation instantané de la zone Zi pour une période de référence Test 5 donné par : Ainsi la durée moyenne de stationnement (0, (Zi,t) dans une zone Zi à l'instant t et sur une période de référence T est donnée par : et le temps moyen, AT (Zi,t), de disponibilité d'une place dans une zone Zi à l'instant t et sur la période de référence T est donné par : 10 Figure 6, selon un exemple de réalisation, la courbe d'évolution du taux d'occupation (602) temporel en fonction du temps (301) d'une zone sur une durée de 24 heures, dépend de la période de référence Td'observation. Ainsi, le taux (611) d'occupation, estimé sur une période T de 24 heures, donne un taux sensiblement constant voisin de 60 % (0,6) selon cet exemple de réalisation. L'examen sur une période Tplus courte, par exemple 12 heures 15 (612) permet de détecter, par exemple, les fluctuations théoriques dans la zone considérée entre la plage horaire de stationnement payant et la plage horaire de stationnement gratuit. L'examen sur des périodes T encore plus courtes, par exemple, 1 heure (613) ou 30 minutes (614) met en évidence les fluctuations de remplissage de la zone en fonction de l'instant. Le choix de la période T d'observation dépend de l'utilisation de ces données. Ainsi, pour des données macroscopiques visant à supporter une décision sur l'implantation d'une zone de stationnement, le suivi du taux d'occupation sur 24 heures et son évolution selon les jours de la semaine ou du mois constitue une indication pertinente. En revanche, pour évaluer la disponibilité instantanée d'une place dans un domaine géographique dans un délai AT court, l'analyse sur une période T plus courte est plus pertinente. Figure 7, de manière similaire l'évolution du taux de rotation (702) d'un domaine géographique en fonction du temps (301), procure une information de pertinence différente selon la période T considérée. Ainsi, le taux de rotation (711) selon une période de 24 heures montre qu'en moyenne, sur une telle période, 7 changements de véhicules interviennent sur une même place. Lorsque le taux de rotation est examiné sur une période plus courte, par exemple 12 heures (712) ou 1 heure (713) le taux de rotation baisse. Ainsi le taux de rotation instantané, en tant que tel, est utilisable par exemple, pour déterminer la probabilité instantanée de trouver une place à un instant (t+At), At étant de l'ordre de grandeur de la période d'observation T, mais aussi, l'évolution de ce taux de rotation en fonction de la période d'observation, comme représenté figure 7, donne des renseignements sur la typologie de la zone, par exemple, s'il s'agit d'une zone d'occupation essentiellement résidentielle, ou d'une zone dont l'occupation est plutôt déterminée par une chalandise commerciale. La demanderesse a ainsi déterminé qu'il existe une corrélation directe entre le taux de rotation d'une zone Zi, sa capacité Ci en nombre de places, et la probabilité PT (Zi,t) instantanée pour un usager d'y trouver une place disponible sur une période de référence T. Figure 8, la demanderesse a ainsi déterminé que plus le taux de rotation pT (Zi,t), d'une zone est élevé et plus les chances d'y trouver une place libre se réduisent en fonction de l'éloignement du véhicule par rapport à ladite zone. De plus, le temps de trajet de l'usager jusqu'à ladite zone est sensiblement proportionnel à son éloignement de cette zone, sauf circonstances exceptionnelles. Ces données de probabilité (802) en fonction contexte, sont ainsi modélisées par une fonction mathématique (800) univoque et sont ainsi prévisibles en fonction du temps (301). À titre d'exemple non limitatif, la fonction suivante a été testée sur des cas réels et donne de bons résultats : Ci ( K- n où KT désigne une constante de la zone Zi déterminée, par exemple, de manière empirique, par des observations sur le terrain. Toutes ces informations sont enregistrées en temps réel dans les moyens de mémoire de l'agrégateur de données. Ainsi, ces données sont utilisables ultérieurement à des fins statistiques ou pour réaliser des simulations. Lesdites informations sont également affichées en temps réel sur les moyens (191 figure 1) d'affichage des moyens de commande et de supervision, afin de délivrer, en temps réel, une information géo- contextualisée et horodatée de l'occupation de l'espace de stationnement et de l'occupation licite de cet espace. Cette information, combinée aux capacités de simulation à partir de l'historique des informations mémorisées, font du dispositif et du procédé objets de l'invention un observatoire dynamique du stationnement et de la mobilité dans l'espace urbain supervisé. À titre d'exemple, la comparaison du taux de rotation d'une zone sur une période de 24 heures avec le taux de rotation moyen de cette zone, établi à partir de l'historique de ladite zone, permet de détecter la présence d'un véhicule, dit ventouse, dans ladite zone. D'un point de vue dynamique, le taux de rotation définit pour un véhicule distant d'une zone Zi, la probabilité de trouver une place de stationnement disponible dans la zone 20 considérée en fonction du nombre d'emplacements libres à un instant t. Ainsi, à un instant t donné, une place n'est disponible que pendant une durée At qui dépend du taux de rotation instantané. Il est par conséquent inutile d'orienter un véhicule à la recherche d'un stationnement vers la zone considérée, si le temps de trajet vers cette zone est supérieur à At et qu'il y a peu de places disponibles dans ladite zone. Par suite, à l'instant t 25 l'information pertinente pour un usager est : Ci.OET gi,t+At) où Ci est la capacité totale de la zone Zi et al, (Zi, t+At) est de taux d'occupation projeté à (t+At). L'affichage de la disponibilité d'une zone de stationnement, mais corrigée par le temps de trajet moyen, At, jusqu'à ladite zone, permet de renseigner l'usager qui, connaissant la ville optimisera lui-même son choix. Ledit temps de trajet moyen est déterminé soit par l'historique des déplacements observés en ville, ou également de manière dynamique. Cette information est avantageusement affichée sur les panneaux à affichage variable installés dans l'agglomération. Le problème se pose cependant pour un usager ne connaissant pas la ville et pour lequel la simple information de disponibilité, même corrigée du temps de trajet, n'est pas 10 suffisante pour réaliser un choix pertinent. À cette fin, le procédé et le dispositif objets de l'invention permettent de réaliser un choix optimisé et ainsi, par exemple, de présenter à l'usager le choix de la zone de stationnement la plus pertinente en fonction du contexte. À cette fin, les informations délivrées par l'agrégateur de données sont 15 avantageusement utilisées pour optimiser le guidage en temps réel d'un usager à la recherche d'un emplacement de stationnement. En revenant à la figure 1, les informations issues de l'agrégateur (110) de données, et plus particulièrement les informations AT (Zi,t) et pT (Zi,t) sont mises à disposition sur internet (199), pour différentes périodes T d'intégration, par l'intermédiaire d'un module 20 (195) spécifique. À titre d'exemple non limitatif, ces informations sont mises à disposition pour T=15 minutes, T=30 minutes, T=1 heure, T=2 heures. Ces données permettent de réaliser une prévision instantanée de l'état de stationnement dans le domaine géographique au temps (t+At), à partir des données correspondant aux informations délivrée par l'agrégateur de données, telles que At T, At étant de l'ordre de grandeur de 25 T. Ainsi, pour At = 10 minutes les données correspondant à une période T = 15 minutes sont utilisées. Ainsi, lesdites données sont récupérées gratuitement, ou par abonnement, par des parties tierces, aptes à mettre en oeuvre et à partager des données de guidage dynamique des usagers, notamment à travers les outils et terminaux de positionnement par satellite. De telles applications sont notamment en mesure de récupérer, par divers canaux, des informations en temps réel sur l'état du trafic. Le guidage des usagers vers la zone de stationnement la mieux adaptée relève d'un problème d'optimisation similaire au problème dit « du voyageur de commerce » lequel problème est résolu en mettant en oeuvre une métaheuristique d'optimisation. Ainsi, le procédé et le dispositif objets de l'invention permettent le recueil des données adaptées pour la mise en oeuvre d'une telle méthode d'optimisation. Selon un exemple de réalisation, la métaheuristique utilisée est celle de la colonie de fourmis. Ainsi, pour un usager se trouvant à une distance d'une zone de stationnement telle qu'un temps de trajet At est nécessaire pour se rendre dans cette zone, ladite zone est définie par sa visibilité et par son intensité. Son intensité représente l'attractivité de la zone. Elle est d'autant plus élevée que la zone comprend un nombre de places disponibles élevé. La visibilité est d'autant plus faible, que le taux de rotation de la zone, pour une période d'intégration T proche de At, est élevé. Le choix est déterminé par le produit pondéré de la visibilité et de l'intensité des différentes zones à l'instant t. Les méthodes d'optimisation par métaheuristiques et plus particulièrement par les algorithmes de colonie de fourmis sont connues de l'art antérieur et l'homme du métier les met en oeuvre à partir des indications ci-avant, à savoir en utilisant l'information Ci.OET(Zi,t+At) pour la détermination de l'intensité instantanée de la zone Zi et le paramètre PT(Zi,t) pour la détermination de la visibilité instantanée de ladite zone Zi. Le dispositif et le procédé objets de l'invention sont également avantageusement utilisés pour organiser le contrôle de l'espace de stationnement par les agents de surveillance de la voirie publique et par les forces de police. Il est en effet peu utile d'orienter vers une zone Zi des agents de surveillance de la 25 voirie publique, même si le nombre de places régulièrement acquittées est inférieur au nombre de véhicules stationnés dans ladite zone, si, dans le temps nécessaire auxdits agents pour rejoindre la zone, les véhicules stationnés ont changé. L'agrégateur de données permet de calculer à partir des données récoltées par le DIC, et les informations du référentiel REF TAR, la recette théorique, 9 (Zi,t), d'une zone 30 à stationnement payant à l'instant t, laquelle recette instantanée est donné par : Zi. TcIf(v).dv.du où Zi.Tarif (v) est le tarif horaire à l'heure y défini dans le référentiel REF TAR, et au(Zi,t) le taux d'occupation de la zone Zi à l'instant t sur la période de référence u. Figure 9, les données issues de l'étape de calcul (920) sont utilisées dans une étape (930) de caractérisation et de prévision, mise en oeuvre, selon cet exemple de réalisation, par le module (113) de calcul de l'agrégateur (110) de données, à partir de modèles, pour caractériser de manière quantitative l'espace de stationnement considéré et notamment prévoir son taux d'occupation futur et des probabilités associées. Les données issues de l'étape (930) de caractérisation et de prévision, lesquelles données sont fondées sur les informations collectées par le DIC (112), et sont utilisées, selon un exemple de réalisation, dans une étape (950) de calcul de la recette théorique de l'espace de stationnement. Avantageusement cette étape (950) de calcul est réalisée pour une table de valeur de la période d'observation. Une étape (960) de collecte et de traitement des informations issues des horodateurs et des concentrateurs de moyens de paiement, mise en oeuvre par le module (111) transactionnel, permet de déterminer la recette ou chiffre d'affaires effectif de l'espace de stationnement pour différentes périodes d'observation. Ainsi, le chiffre d'affaires instantané CA HOi (t) d'un horodateur Hoi est donné par : CA HOi (t)= transaction 0).montant La sommation étant réalisée sur le nombre de transactions. Le chiffre d'affaires instantané sur une zone Zi, CA Zi (t), est obtenu en sommant les chiffres d'affaires de tous les horodateurs contenus dans ladite zone et, le cas échéant, le chiffre d'affaires correspondant aux autres moyens de paiement, sur une période de référence u. Ces données sont agrégées par le module transactionnel. L'enregistrement de ces informations en temps réel dans les moyens de mémoire de l'agrégateur de données permet de disposer de statistiques et, par exemple, d'intégrer une 25 part prévisionnelle dans ce chiffre d'affaires, notamment en présence d'horodateurs autonomes dont l'état des recettes n'est interrogé que de manière périodique. Figure 9, le chiffre d'affaires de l'espace de stationnement considéré est comparé 3004 84 1 22 (970) à la recette théorique du même espace et pour la même période. La comparaison de la recette théorique instantanée, élaborée à partir des données du module DIC avec le chiffre d'affaires effectif, établi par le module transactionnel dans cette même zone, permet de définir un taux d'acquittement °Zi (t,T) instantané sur une 5 période de référence u=T : (t où CAL, (t) est le chiffre d'affaires établi par le module transactionnel sur la zone considérée. Ainsi, les données issues de l'étape (930) de caractérisation et de prévision, ou les données issues de la comparaison (970) de la recette théorique et de la recette réelle, sont 10 comparées (980) avec des références pour un espace de stationnement et une période d'observation donnés. Selon le résultat de cette comparaison une information est affichée (940) à l'attention des usagers, ou des agents de surveillance de la voirie publique ou encore à l'attention d'un superviseur. En revenant à la figure 1, selon un exemple de réalisation, un agent de surveillance 15 de la voirie publique dispose d'un terminal (192) mobile connecté par liaison radio au centre (190) de contrôle et de supervision du dispositif objet de l'invention. Ainsi, l'affichage, sur le terminal (192) mobile d'un agent, de l'information du taux d'acquittement instantané d'un ensemble de zones, permet à l'agent d'organiser sa tournée. Selon un autre mode de réalisation, l'équipe d'agents de surveillance est pilotée par 20 un superviseur. Ledit superviseur dispose, au moyen du dispositif objet de l'invention, d'une information en temps réel sur l'occupation de l'espace de stationnement et sur le taux d'acquittement des différentes zones Zi de cet espace, par exemple au moyen d'un terminal (191) du centre (190) de contrôle. Ainsi, ledit superviseur utilise cette information pour piloter l'équipe d'agents de surveillance en communiquant avec lesdits agents au moyen de leur terminal (192) mobile ou par tout autre moyen. Toutefois, à l'instant t, ledit superviseur ne dispose que d'un nombre fini d'agents de surveillance et est par ailleurs contraint par des dispositions légales, notamment liées à la législation du travail. Par exemple, une équipe d'agent de surveillance se rendant sur une zone, doit pouvoir verbaliser dans un temps compatible avec ses horaires de travail en comptant, le cas échéant, le retour au vestiaire. Ainsi, le centre (190) de contrôle objet de l'invention dispose avantageusement d'un applicatif permettant d'optimiser la tournée des agents de 5 surveillance ou d'assister le superviseur de ces équipes dans l'organisation desdites tournées, en fonction de l'occupation de l'espace de stationnement et des autres contraintes. Ainsi, le dispositif objet de l'invention comporte un référentiel (164), REF ASVP, comportant les règles d'intervention des agents de surveillance de la voirie publique. Ledit référentiel comporte les contraintes liées à l'intervention des agents de surveillance par 10 exemple : - des contraintes légales et réglementaires issues des dispositions de la législation du travail et des conventions collectives ; - des informations sur la disponibilité des agents tels que leurs horaires de travail, les périodes de congés, les arrêts ; 15 - le cumul des heures déjà réalisées ; - les zones qui peuvent être surveillées par un agent ; - le nombre d'agents minimum pour contrôler une zone Zi donnée ; sans que cette liste ne soit exhaustive. Ainsi, l'applicatif vise à définir la meilleure affectation des agents en fonction de 20 l'occupation de l'espace de stationnement. Il s'agit d'un problème d'optimisation sous contrainte, qui, selon un exemple de réalisation est résolu par la mise en oeuvre d'une métaheuristique d'optimisation telle qu'un algorithme de colonie de fourmis. Ainsi la tournée d'un agent de surveillance de la voirie publique est modélisable comme un chaînage fini de déplacements 8.j au sein d'une même zone et di entre deux 25 zones. Chaque déplacement di est caractérisé par un temps de déplacement Gdi : Distance( + où Distance (Zi,Zi+1) est la distance entre deux zones successives de la tournée et V la vitesse moyenne de déplacement de l'agent. Chaque déplacement 8.j est caractérisé par un temps de déplacement CY8j : C 1 -OH?' T)) où e est le temps nécessaire à l'établissement d'un procès-verbal de contravention, il) le temps de parcours entre deux places et Cj la capacité de la zone Zj. Ainsi, à chaque déplacement di et 8.j correspond une efficacité économique où (di) est fonction du nombre de véhicules en stationnement illicite et : Le rendement du déplacement de l'agent est le rapport entre l'efficacité économique du déplacement et le temps de déplacement. La connaissance en temps réel et l'archivage de ces données dans les moyens de mémoire de l'agrégateur de données du dispositif objet de l'invention permettent d'enrichir les modèles par l'analyse des données et de réaliser des simulations, notamment par l'analyse statistique de ces données. Ainsi, la mise en oeuvre de la métaheuristique d'optimisation de la tournée des agents utilise l'inverse du temps de trajet de l'agent vers la zone à contrôler comme la visibilité de ladite zone et le rendement du déplacement vers cette zone comme l'intensité de ladite 15 zone. En revenant à la figure 2, le domaine géographique comprend des emplacements de stationnement critiques, par exemple, un emplacement (221) réservé aux personnes à mobilité réduite ou un passage pour piéton. Ces emplacements sont identifiés à la fois dans le référentiel spatial, REF SPA, et dans le référentiel réglementaire REF REG. À titre 20 d'exemple non limitatif, les emplacements particuliers comprennent : - les places réservées aux personnes à mobilité réduite ; - les places réservées aux livraisons ; - les places réservées aux véhicules électriques ; - les places réservées aux véhicules de transport de fonds ; 25 - les places réservées à une catégorie particulière d'usager (police, corps diplomatique ...). Avantageusement la zone (200) de stationnement comprend des capteurs spécifiques, par exemple, une ou plusieurs caméras de vidéosurveillance pour détecter la présence d'un véhicule à ces emplacements critiques et le cas échéant, vérifier que ledit véhicule est autorisé à utiliser un tel emplacement. À titre d'exemple, la caméra de vidéosurveillance de l'emplacement (221) pour personne à mobilité réduite, dispose de moyens pour lire la plaque minéralogique du véhicule ou identifier un macaron spécifique sur le pare-brise dudit véhicule. Ainsi, le dispositif objet de l'invention est en mesure de déclencher une alerte en cas d'utilisation abusive d'un emplacement (221, 225) réservé ou interdit. De plus, la connaissance à tout instant pour une zone Zi de l'information Ci.ocT(Zi,t+At) permet, lorsque cette valeur devient inférieure à un seuil défini, d'identifier un risque d'occurrence de stationnement en double file dans une zone donnée et dans un temps At. Ainsi, le superviseur du stationnement est en mesure d'anticiper l'envoi d'agents de police dans ladite zone afin de prévenir ce risque. Par ces possibilités d'anticipation, le dispositif et le procédé objets de l'invention permettent d'accroître l'efficacité de tous les préposés au contrôle du stationnement. En revenant à la figure 1, le centre (190) de contrôle et de décision comporte une pluralité de modules applicatifs pour mettre en oeuvre les modes de réalisation exposés ci-avant. Plus particulièrement, selon un exemple de réalisation avantageux, ledit centre de contrôle comporte un module applicatif apte à modifier le référentiel (162) tarifaire en fonction des informations délivrées par l'agrégateur (110) de données. Ainsi, l'offre tarifaire est modifiée en fonction du taux d'occupation de l'espace de stationnement afin de modifier le taux de rotation dans certaines zones selon un ajustement dynamique incitatif. La description ci-avant et les exemples de réalisation montrent que l'invention atteint les objectifs visés. En particulier elle permet un pilotage dynamique de l'offre de stationnement dans un domaine géographique, en anticipant le taux d'occupation et le taux de rotation dans un espace de stationnement à partir de données acquises en temps réel. En séparant les fonctions de comptage et de mesure du taux d'occupation de l'espace de stationnement, des fonctions de collecte des données d'acquittement des taxes de stationnement dans cet espace, le dispositif et le procédé objets de l'invention permettent d'optimiser la surveillance dudit espace de stationnement. Le procédé objet de l'invention est aisément mis en oeuvre à partir de balises de comptage, facilement installées en voirie sans travaux onéreux. REVENDICATIONS1. Dispositif pour pour la supervision du stationnement des véhicules dans un domaine géographique caractérisé en ce qu'il comprend : un ensemble de capteurs apte à déterminer la présence et le stationnement d'un véhicule dans le domaine géographique, notamment par comptage ; un agrégateur (110) de données apte à récupérer les informations générées par l'ensemble des capteurs et comprtant un caculateur (113) apte à réaliser des calculs statistiques sur la base de ces données ; un référentiel, sous forme numérique, comportant : un référentiel (161), dit spatial, comportant une description géographique du domaine, et de l'espace de stationnement ; un référentiel (162), dit tarifaire, comprenant les gammes tarifaires et leur association avec l'espace de stationnement ; un référentiel (163), dit réglementaire, comprenant les dispositions réglementaires associées à l'espace de stationnement ; iv. un module (190), dit applicatif, apte à délivrer des informations à des terminaux (161,162, 191). 2. Dispositif selon la revendication 1, comportant : un module (195) de mise à disposition de données formatées de l'agrégateur (110) de données sur internet (199). 3. Dispositif selon la revendication 1, dans lequel le référentiel comporte : un référentiel (164) de règles définissant les modalités d'intervention des agents de surveillance de la voirie. 4. Procédé pour la supervision du stationnement des véhicules dans un domaine géographique au moyen du dispositf selon la revendication 1, caractérisé en ce qu'il comprend les étapes consistant à: a. obtenir (910), au moyen de l'agrégateur (110), l'état d'occupation instantané de l'espace de stationnement correspondant audit domaine géographique à un instant t; b. obtenir (920), de l'agrégateur (110), le taux de rotation instantané du stationnement sur une période Tdans ledit espace de stationnement à l'instant t; c. calculer (930) à partir des données recueillies aux étapes a) et b) le taux d'occupation futur à un instant (t+At) et sa probabilité ; d. délivrer (940), par le module (190) applicatif, une information lorsque la probabilité ou le taux d'occupation déterminés à l'étape c) sont inférieurs ou supérieurs à un seuil déterminé. 5 Procédé selon la revendication 4, dans lequel l'étape a) comprend les étapes consistant à : ai. obtenir (911) le nombre de véhicules instantanément présents dans le domaine géographique ; aii. obtenir (912) le nombre instantané de véhicules en transit dans ledit domaine géographique ; aiii. déterminer (913) le nombre de véhicules stationnés dans le domaine géographique en soustrayant le nombre défini à l'étape au) du nombre défini à l'étape ai). 6. Procédé selon la revendication 4, dans lequel l'étape b) comprend les étapes consistant à: bi. obtenir (921) le taux d'occupation spatial du domaine géographique par le nombre de véhicules stationnés dans ledit domaine à l'instant t; 7. 8. 9. 10. 11. bii. déterminer (922) le taux d'occupation temporel dudit domaine sur une période T par la moyenne du taux d'occupation spatial sur cette période ; biii. déterminer (923) le taux de rotation dudit domaine géographique en moyennant, sur la période T, la variation du taux d'occupation temporel déterminé à l'étape bii) Procédé selon la revendication 4, dans lequel l'étape c) est déterminée en fonction d'un temps de trajet At vers le domaine géographique au moyen d'une métaheuristique d'optimisation. Procédé selon la revendication 7, dans lequel la métaheuristique utilise une règle aléatoire de transition proportionnelle issue d'un algorithme de colonie de fourmis et où l'information délivrée à l'étape d) est construite à partir de la visibilité du domaine géographique au temps (t+At). Procédé selon la revendication 5, dans lequel l'étape au) comprend la détermination de la vitesse moyenne des véhicules sur une partie du domaine géographique, et la détermination de la fréquence temporelle dominante d'entrée dans ladite partie du domaine à l'instant t. Procédé selon la revendication 5, comportant les étapes (950) consistant à: e. obtenir une grille tarifaire du stationnement dans le domaine géographique ; déterminer la recette théorique à l'instant t et pour une période T en fonction du taux d'occupation du domaine déterminé à l'étape bii). Procédé selon la revendication 10, comportant les étapes consistant à: g- obtenir (960) la recette effective du domaine correspondant à l'instant t et sur une période d'observation T; h. déterminer (970) le taux d'acquittement du domaine géographique considéré sur une période donnée en comparant sur cette période lerapport entre les résultats instantanés des étapes g) et f). 12. Procédé selon la revendication 11, comportant les étapes consistant à: j. combiner l'information obtenue à l'étape c) avec l'information obtenue à l'étape h) et déterminer le taux d'acquittement futur à l'instant (t+At) et sa probabilité ; k. délivrer une information si la probabilité ou le taux d'acquittement déterminés à l'étape j) sont inférieurs ou supérieurs à un seuil donné. 13. Procédé selon la revendication 12, dans lequel l'étape j) est déterminée en fonction d'un temps de trajet At d'un agent de surveillance de la voirie publique vers le domaine géographique au moyen d'une métaheuristique d'optimisation utilisant une règle aléatoire de transition proportionnelle issue d'un algorithme de colonie de fourmis et où l'information délivrée à l'étape k) est construite à partir de la visibilité du domaine géographique au temps (t+At). 14. Procédé selon la revendication 10, comportant une étape consistant à: modifier la grille tarifaire en fonction du taux d'occupation déterminé à l'étape a) ou de la recette théorique déterminée à l'étape f).