Brevet FR3065529A1 — publié le 20181026

OPTIMISATION DE LA DISTRIBUTION SPATIALE DE MOYENS DE MESURE DE LA QUALITE DE L'AIR Domaine technique L'invention concerne un système de mesure d'au moins une grandeur physique représentative de la qualité de l'air dans une zone d'observation. L'invention trouve notamment son application dans la surveillance et le contrôle de la qualité de l'air dans des zones urbaines. Etat de la technique antérieure Un système de mesure connu de l'état de la technique, notamment du document « Guidelines for Ambient Air Quake Monitoring» établi par le CPCB (Central Pollution Control Board) du ministère de l'environnement du gouvernement de l'Inde, daté d'avril 2003, (ci-après Dl) comporte un ensemble de stations de mesure de la qualité de l'air. Dl mentionne (cf. §4.2.2) des critères pour sélectionner les emplacements des stations de mesure dans la zone d'observation. Toutefois, Dl précise (cf. §4.2.2.a) qu'il est difficile de savoir à l'avance si les emplacements sélectionnés vont refléter et rendre compte de la qualité de l'air dans la zone d'observation. D1 propose (cf. §4.2.2.a) alors d'effectuer des mesures préalables afin de vérifier la pertinence des emplacements. Un tel système de mesure de l'état de la technique n'est pas entièrement satisfaisant, en particulier lorsque l'on souhaite augmenter la résolution spatiale des mesures. La solution proposée par Dl est très fastidieuse et coûteuse lorsque les stations de mesure présentent une densité spatiale élevée dans la zone d'observation. Exposé de l'invention L'invention vise à remédier en tout ou partie aux inconvénients précités. A cet effet, l'invention a pour objet un système de mesure d'au moins une grandeur physique représentative de la qualité de l'air dans une zone d'observation, le système comportant : - une cartographie de la zone d'observation, comprenant un ensemble, noté V, de valeurs 30 modélisées représentatives de la grandeur physique ; - des moyens de mesure de la grandeur physique, possédant un nombre N de positions ou un nombre N de trajectoires dans la zone d'observation, les N positions ou les N trajectoires étant destinées à présenter une distribution spatiale, notée Sopt , dans la zone d'observation ; - des moyens de calcul de la distribution spatiale Sopt, configurés pour : 3065529 2 - effectuer un maillage de la zone d'observation, le maillage comprenant un nombre G de points ; - calculer, pour une distribution spatiale donnée des N positions ou des N trajectoires, notée S, un estimateur de l'ensemble V, noté g, pour chacun des G 5 points du maillage ; - calculer une fonction-coût, notée y(S), représentative de la différence ou de la vraisemblance entre 11 et les valeurs modélisées de l'ensemble V extraites aux G points du maillage, notées V; - extraire la distribution spatiale Sop, qui minimise ou maximise la fonction-coût 10 selon que la fonction-coût est représentative de la différence ou de la vraisemblance entre V et V. Ainsi, un tel système de mesure selon l'invention permet d'obtenir des mesures reflétant et rendant compte de la qualité de l'air dans la zone d'observation, autorisant une résolution 15 spatiale supérieure à l'état de la technique, tout en s'affranchissant de mesures préalables locales. En effet, une telle cartographie et de tels moyens de calcul permettent, en conjugaison, l'extraction d'une distribution spatiale optimale des moyens de mesure à partir de valeurs modélisées. Définitions 20 - Par « cartographie », on entend une représentation spatiale des valeurs modélisées, pouvant être par exemple en deux dimensions (2D) ou en trois dimensions (3D). - Par «représentatives de la grandeur physique », on entend que les valeurs modélisées représentent directement la grandeur physique, ou représentent indirectement la grandeur physique par corrélation. 25 - N est un nombre entier naturel. - Les «N positions » des moyens de mesure sont définies par des coordonnées spatiales, par exemple des coordonnées cartésiennes (x, y, z). - Les «N trajectoires » des moyens de mesure sont N ensembles de T positions successives au cours du temps. Les N trajectoires sont définies par des coordonnées spatiales, 30 par exemple des coordonnées cartésiennes (x, y, z), et par une composante temporelle (t), i étant compris entre 1 et T. - Par « distribution spatiale », on entend la répartition des N positions (ou des N trajectoires) dans l'espace de la zone d'observation. On peut parler de distribution volumique ou de densité volumique dans le cas 3D. 3065529 3 - Par « maillage », on entend la discrétisation spatiale de la zone d'observation. - Par « estimateur », on entend un estimateur statistique où les N positions (ou les N trajectoires) sont des variables aléatoires. - Par « fonction-coût », on entend toute fonction capable de quantifier la cohérence entre 5 17 et P. On parle également de fonction-objectif. Le système de mesure selon l'invention peut comporter une ou plusieurs des caractéristiques suivantes. 10 Selon une caractéristique de l'invention, les moyens de calcul sont configurés pour calculer l'estimateur de l'ensemble V selon la formule : Wi- Vi = j E où : - V, sont les valeurs modélisées de l'ensemble V extraites pour les positions ou pour les 15 trajectoires des moyens de mesure présentant la distribution spatiale donnée S, - w est une fonction de transfert entre si et mi, où si sont les positions ou les trajectoires des moyens de mesure présentant la distribution spatiale donnée S, et où mi sont les positions des points du maillage dans la zone d'observation. Ainsi, un avantage procuré par la fonction de transfert est de pondérer les valeurs V, en 20 tenant compte des positions relatives si et mi, par exemple selon une distance ou un angle issu des positions relatives si et mi. Selon une caractéristique de l'invention, les moyens de calcul sont configurés pour calculer l'estimateur de l'ensemble V selon la formule : 25 -2 vs [d(mps3] j -2 E [[1, G]] où : - mi sont les positions des points du maillage dans la zone d'observation, - si sont les positions ou les trajectoires des moyens de mesure présentant la distribution spatiale donnée S, - V, sont les valeurs modélisées de l'ensemble V extraites pour les positions ou pour les trajectoires si, 3065529 4 - d est une distance entre une position d'un point du maillage et une position ou une trajectoire des moyens de mesure dans la zone d'observation. Ainsi, un avantage procuré par un tel estimateur est sa simplicité de mise en oeuvre pour les calculs. 5 Selon une caractéristique de l'invention, la fonction-coût est représentative de la différence entre V et V, et la fonction-coût est une norme de la différence entre V et V. Ainsi, la minimisation d'une telle fonction-coût permet d'envisager l'obtention d'une distribution spatiale Sopt reflétant au mieux la qualité de l'air dans la zone d'observation en 10 référence aux valeurs modélisées. Selon une caractéristique de l'invention, la fonction-coût est représentative de la différence entre 17 et ri, et les moyens de calcul sont configurés pour calculer la fonction-coût selon la formule : 15 G 1 q(S)= UI(Vi ril)2 i=1 où vi sont les valeurs modélisées de l'ensemble V extraites pour les points du maillage. Selon une caractéristique de l'invention, la fonction-coût est représentative de la différence entre 11 et V, et dans lequel les moyens de calcul sont configurés pour effectuer 20 une modification aléatoire de la distribution spatiale donnée S selon une probabilité, notée p)' vérifiant : Pji = 1 si (p(SW) cp(S(i)) PJL = exp (cp(S(i)) - cp(SU))) si «,(s-n) > cio(S(0) où : 25 - 5(i) est une distribution spatiale donnée initiale, - SO) est une distribution spatiale donnée modifiée aléatoirement, - c est un paramètre ; et les moyens de calcul sont configurés pour itérer la modification aléatoire de la distribution spatiale donnée S jusqu'à l'extraction de St. 30 Ainsi, un avantage procuré est d'obtenir un algorithme robuste vis-à-vis du choix de la distribution spatiale donnée initiale, c'est-à-dire que la distribution spatiale Sopt ne dépend pas de la distribution spatiale donnée initiale. 3065529 5 Selon une caractéristique de l'invention, la fonction-coût est représentative de la vraisemblance entre I7 et V, et la fonction-coût comprend un produit scalaire entre I7 et V, de préférence normalisé. 5 Ainsi, la maximisation d'une telle fonction-coût permet d'envisager l'obtention d'une distribution spatiale Sopt reflétant au mieux la qualité de l'air dans la zone d'observation en référence aux valeurs modélisées. Selon une caractéristique de l'invention, la fonction-coût est représentative de la 10 vraisemblance entre I7 et V, et dans lequel les moyens de calcul sont configurés pour effectuer une modification aléatoire de la distribution spatiale donnée S selon une probabilité, notée pi, vérifiant : Pji = 1 si cp(SW) (p(S(i)) = exp (q)(.5(i)) - v(S(0)) Si v(sw) N. En effet, il est alors possible de désactiver les moyens de mesure dont les positions ou les trajectoires dans la zone d'observation n'appartiennent pas aux N positions ou N trajectoires présentant la distribution spatiale la plus proche de S. 5 Selon une caractéristique de l'invention, la zone d'observation comporte un nombre M de positions autorisées ou un nombre M de trajectoires autorisées, vérifiant M>N, dans lesquelles les moyens de mesure sont respectivement autorisés à posséder des positions ou des trajectoires. 10 Selon une caractéristique de l'invention, les moyens de mesure comportent au moins N dispositifs de mesure équipant chacun un véhicule, le véhicule étant de préférence routier ou aérien. Ainsi, un avantage procuré est de pouvoir utiliser la circulation de véhicules dans la zone 15 d'observation comme déploiement de dispositifs de mesure selon différentes distributions spatiales de positions ou de trajectoires. Selon une caractéristique de l'invention, les moyens de mesure comportent au moins N dispositifs de mesure équipant chacun un mobilier urbain, le mobilier urbain étant de 20 préférence sélectionné dans le groupe comportant des feux de circulation et des abribus. Ainsi, un avantage procuré est d'éviter de créer des emplacements dédiés pour les dispositifs de mesure. Selon une caractéristique de l'invention, les moyens de mesure comportent des capteurs 25 spectroscopiques, de préférence des capteurs infrarouges non-dispersifs. Ainsi, un avantage procuré est leur compacité par rapport à des stations de mesure conventionnelles, permettant une densité spatiale plus élevée et à moindre coût. L'invention a également pour objet un procédé de mesure d'au moins une grandeur 30 physique représentative de la qualité de l'air dans une zone d'observation, le procédé comportant les étapes : a) prévoir une cartographie de la zone d'observation, comprenant un ensemble, noté V, de valeurs modélisées représentatives de la grandeur physique ; 3065529 7 b) prévoir des moyens de mesure de la grandeur physique, possédant un nombre N de positions ou un nombre N de trajectoires dans la zone d'observation, les N positions ou les N trajectoires étant destinées à présenter une distribution spatiale, notée S'p, ; c) effectuer un maillage de la zone d'observation, le maillage comprenant un nombre G de 5 points ; d) calculer, pour une distribution spatiale donnée des N positions ou des N trajectoires, notée S, un estimateur de l'ensemble V, noté 11, pour chacun des G points du maillage ; e) calculer une fonction-coût, notée y(S), représentative de la différence ou de la vraisemblance entre V et les valeurs modélisées de l'ensemble V extraites aux G points du 10 maillage, notées V; f) extraire la distribution spatiale S.pt qui minimise ou maximise la fonction-coût selon que la fonction-coût est représentative de la différence ou de la vraisemblance entre 12 et V; g) agencer les moyens de mesure dans la zone d'observation de sorte que les N positions ou les N trajectoires présentent la distribution spatiale Scpt. 15 Brève description des dessins D'autres caractéristiques et avantages apparaîtront dans l'exposé détaillé de différents modes de réalisation de l'invention, l'exposé étant assorti d'exemples et de référence aux dessins joints. 20 Figure 1 est une vue schématique partielle d'un système de mesure selon l'invention, illustrant des moyens de mesure, agencés dans une zone d'observation (un quartier de la ville de Paris), et possédant N positions présentant la distribution spatiale Sopt (ronds blancs, N=10), sélectionnées parmi M positions autorisées (ronds noirs). Figure 2 est une vue schématique partielle d'un système de mesure selon l'invention, 25 illustrant un maillage d'une zone d'observation (un quartier de la ville de Paris). Figure 3 est une vue schématique partielle d'un système de mesure selon l'invention, illustrant une cartographie d'une zone d'observation (un quartier de la ville de Paris), la cartographie comprenant un ensemble de valeurs modélisées représentatives de la grandeur physique à mesurer. 30 Figure 4 est un graphique illustrant en ordonnées la fonction-coût et en abscisses le nombre d'itérations (i.e. le nombre de modifications de la distribution spatiale des N positions ou N trajectoires). Figures 5 à 7 sont des vues analogues à la figure 1, illustrant respectivement une distribution spatiale initiale (ronds blancs) des N positions (N=10), une distribution spatiale 3065529 8 intermédiaire (ronds blancs) des N positions, et la distribution spatiale Scpt (ronds blancs) des N positions, sélectionnées parmi M positions autorisées (ronds noirs). Figure 8 est un graphique illustrant en ordonnées la fonction-coût et en abscisses le nombre N de positions ou de trajectoires. 5 Figure 9 est une vue schématique partielle d'un système de mesure selon l'invention, illustrant des moyens de mesure, agencés dans une zone d'observation (un quartier de la ville de Paris), et possédant N positions présentant la distribution spatiale Sopt (ronds blancs, N=10). Figure 10 est une vue schématique partielle d'un système de mesure selon l'invention, 10 comportant des moyens d'activation configurés pour activer des moyens de mesure fixes dans la zone d'observation, et possédant N positions présentant la distribution spatiale la plus proche de S.p, (ronds noirs encerclés, N=10). La distribution spatiale S.pt extraite est illustrée par des ronds blancs. L'ensemble des positions fixes des moyens de mesure est illustré par des ronds noirs. 15 Figure 11 est une vue schématique partielle d'un système de mesure selon l'invention, comportant des moyens d'activation configurés pour activer des moyens de mesure mobiles dans la zone d'observation, et possédant N positions, à un instant donné, présentant la distribution spatiale la plus proche de Sopt (ronds noirs encerclés, N=10). La distribution spatiale Sop, extraite est illustrée par des ronds blancs. L'ensemble des positions des moyens de 20 mesure, à un instant donné, est illustré par des ronds noirs. Figure 12 est une vue schématique partielle d'un système de mesure selon l'invention, illustrant des moyens de mesure, agencés dans une zone d'observation, et possédant N positions présentant la distribution spatiale S.pt (points encerclés, N=7), sélectionnées parmi M positions autorisées (points). 25 Figure 13 est une vue schématique partielle d'un système de mesure selon l'invention, illustrant des moyens de mesure, agencés dans une zone d'observation, et possédant N trajectoires présentant la distribution spatiale S.pt (pointées par les flèches, N=2), sélectionnées parmi M trajectoires autorisées (M=4). 30 Exposé détaillé des modes de réalisation Les éléments identiques ou assurant la même fonction porteront les mêmes références pour les différents modes de réalisation, par souci de simplification. 3065529 9 Un objet de l'invention est un système de mesure d'au moins une grandeur physique représentative de la qualité de l'air dans une zone d'observation 1, le système comportant : - une cartographie 2 de la zone d'observation 1, comprenant un ensemble, noté V, de valeurs modélisées représentatives de la grandeur physique ; 5 - des moyens de mesure de la grandeur physique, possédant un nombre N de positions ou un nombre N de trajectoires dans la zone d'observation 1, les N positions ou les N trajectoires étant destinées à présenter une distribution spatiale, notée S.pt , dans la zone d'observation 1; - des moyens de calcul de la distribution spatiale Sopt, configurés pour : 10 - effectuer un maillage 3 de la zone d'observation 1, le maillage 3 comprenant un nombre G de points ; - calculer, pour une distribution spatiale donnée des N positions ou des N trajectoires, notée S, un estimateur de l'ensemble V, noté I7, pour chacun des G points du maillage 3; 15 - calculer une fonction-coût, notée p(S), représentative de la différence ou de la vraisemblance entre 17 et les valeurs modélisées de l'ensemble V extraites aux G points du maillage 3, notées V; - extraire la distribution spatiale Sopt qui minimise ou maximise la fonction-coût selon que la fonction-coût est représentative de la différence ou de la 20 vraisemblance entre 12 et V. Grandeur physique A titre d'exemples non limitatifs, la ou les grandeurs physiques peuvent être des concentrations de molécules polluantes ou des concentrations de particules nocives pour l'environnement. On peut citer le monoxyde de carbone CO, le dioxyde de soufre S02, le 25 dioxyde d'azote NO2, les matières particulaires en suspension (SPM pour « SueendedParticulate Malter»), les matières particulaires en suspension respirables (RSPM pour «Respirable Sueended Particulate Matter»), les hydrocarbures aromatiques polycycliques (PAHs pour «Polygclic AromaticHydrocarbons»), l'ozone 03. Zone d'observation 30 A titre d'exemples non limitatifs, la zone d'observation 1 peut être une zone urbaine, une zone périurbaine, une zone industrielle. La zone d'observation 1 peut comporter un nombre M de positions autorisées ou un nombre M de trajectoires autorisées, vérifiant M>N, dans lesquelles les moyens de mesure 3065529 10 sont respectivement autorisés à posséder des positions ou des trajectoires. M est un entier naturel. Cartographie Il existe dans l'état de la technique des cartographies 2 de zones d'observation 1, 5 généralement urbaines, pouvant comprendre des valeurs modélisées de la distribution spatiale de concentrations en molécules polluantes ou en particules nocives. Ces cartographies 2 sont classiquement issues de modèles développés pour étudier la pollution atmosphérique et son évolution temporelle. A titre d'exemple non limitatif, la publication de R. Berkowicz et al., «Modeling traffic pollution in streets », National Environmental Research Institute, datée de 10 janvier 1997, décrit un modèle de dispersion spatiale de polluants adapté aux spécificités des zones urbaines. En milieu urbain, la topographie particulière formée de rues séparées par des bâtiments justifie une approche spécifique, tenant compte de la formation de vortex de circulation d'air dans les rues, ces vortex jouant un rôle déterminant dans la dispersion de la pollution atmosphérique. De tels modèles sont désignés par les termes « Street Canyon Model » 15 ou « Street Model ». La publication précitée décrit un modèle d'estimation de la pollution en milieu urbain désigné par l'acronyme OSPM (« Operational Street Pollution Model »), c'est-à-dire un modèle opérationnel de pollution urbaine. Selon ce modèle, à partir de l'émission d'un polluant dans une rue, dépendant du nombre de véhicules et d'une émission moyenne par véhicule, le modèle prend en compte le vortex de recirculation formé dans la rue, la 20 turbulence aérologique résultant du trafic routier, la pollution ambiante, provenant d'autres rues, ainsi que le vent circulant dans la canopée, au-dessus du milieu urbain. La publication de R. Berkowicz « OSPM - A Parameterised Street Pollution Model », Environmental Monitoring and Assessment, 65(1), pp. 323-331, 2000, présente également les hypothèses sur lesquelles est fondé le modèle OSPM, ainsi qu'une validation expérimentale de ce modèle. 25 La cartographie 2 de l'ensemble V des valeurs modélisées peut représenter directement la grandeur physique mesurée. Selon une variante, la cartographie 2 de l'ensemble V des valeurs modélisées peut représenter indirectement, par corrélation, la grandeur physique mesurée. A titre d'exemples non limitatifs, on peut citer des valeurs modélisant les nuisances sonores, l'ensoleillement, le vent, la température, le trafic routier, les ondes électromagnétiques (radio, 30 téléphone portable). La cartographie 2 présente une résolution spatiale supérieure à celle du maillage 3. 3065529 11 Moyens de mesure de la grandeur physique Les moyens de mesure comportent des dispositifs de mesure. Les dispositifs de mesure comportent avantageusement des capteurs spectroscopiques, de préférence des capteurs infrarouges non-dispersifs. 5 Selon un mode de réalisation, les moyens de mesure sont agencés dans la zone d'observation 1 de manière à posséder N positions ou N trajectoires présentant la distribution spatiale Sept. Selon une variante, le système de mesure comporte des moyens d'activation configurés pour activer les moyens de mesure possédant N positions ou N trajectoires présentant la distribution spatiale la plus proche de Sept. A titre d'exemple, les moyens 10 d'activation peuvent être commandés depuis un poste central de régulation. Selon un mode de réalisation, les moyens de mesure comportent au moins N dispositifs de mesure équipant chacun un véhicule, le véhicule étant de préférence routier ou aérien. A titre d'exemples non limitatifs, le véhicule routier peut être une voiture ou un bus ; le véhicule 15 aérien peut un être un aéronef sans pilote tel qu'un drone. Selon une variante, les moyens de mesure comportent au moins N dispositifs de mesure équipant chacun un mobilier urbain, le mobilier urbain étant de préférence sélectionné dans le groupe comportant des feux de circulation et des abribus. 20 Exemple n°1 Dans cet exemple illustré à la figure 9, les moyens de mesure comportent N dispositifs de mesure fixes dans la zone d'observation 1. Les N dispositifs de mesure sont agencés dans la zone d'observation 1 de manière à posséder N positions présentant la distribution spatiale Sept (ronds blancs). 25 Exemple n°2 Dans cet exemple illustré à la figure 10, les moyens de mesure comportent un nombre total Nt, de dispositifs de mesure (Nt., > N, ronds noirs) fixes dans la zone d'observation 1. Le système de mesure comporte des moyens d'activation configurés pour activer, sur une 30 période donnée, les dispositifs de mesure possédant N positions présentant la distribution spatiale la plus proche de Sept (ronds noirs encerclés). La distribution spatiale Sept extraite est illustrée par des ronds blancs. Exemple n°3 3065529 12 Dans cet exemple illustré à la figure 11, les moyens de mesure comportent un nombre total N,0, de dispositifs de mesure (N,Q, > N, ronds noirs) mobiles dans la zone d'observation 1. Le système de mesure comporte des moyens d'activation configurés pour activer, à instant donné, les dispositifs de mesure possédant N positions présentant la 5 distribution spatiale la plus proche de S.p, (ronds noirs encerclés). La distribution spatiale Scipt extraite est illustrée par des ronds blancs. Exemple n°4 Dans cet exemple illustré à la figure 12, la zone d'observation 1 comporte M positions 10 autorisées (points). Les moyens de mesure comportent N dispositifs de mesure agencés dans la zone d'observation 1 de manière à posséder N positions (parmi les M positions autorisées) présentant la distribution spatiale S'p, (points encerclés). Exemple n°5 15 Dans cet exemple illustré à la figure 13, la zone d'observation 1 comporte M trajectoires autorisées. Les moyens de mesure comportent N dispositifs de mesure agencés dans la zone d'observation 1 de manière à posséder N trajectoires (parmi les M trajectoires autorisées) présentant la distribution spatiale S'p, (pointées par les flèches). Moyens de calcul de la distribution spatiale Sao 20 Les moyens de calcul sont avantageusement mis en oeuvre par ordinateur ou par tout dispositif programmable comportant les instructions de calcul adaptées. Estimateur de l'ensemble V Les moyens de calcul sont avantageusement configurés pour calculer l'estimateur de l'ensemble V selon la formule : 25 Wi- Vi = E LAÏ=1 Wii où : - V, sont les valeurs modélisées de l'ensemble V extraites pour les positions ou pour les trajectoires des moyens de mesure présentant la distribution spatiale donnée S, - w est une fonction de transfert entre si et mi, où si sont les positions ou les trajectoires 30 des moyens de mesure présentant la distribution spatiale donnée S, et où mi sont les positions des points du maillage 3 dans la zone d'observation 1. 3065529 13 La fonction de transfert permet de pondérer les valeurs V, en tenant compte des positions relatives si et mi, par exemple selon une distance ou un angle issu des positions relatives si et m. ). Les moyens de calcul sont avantageusement configurés pour calculer l'estimateur de 5 l'ensemble V selon la formule : Ele_1 y.[d (mi, si)] -2 _2 ,j E 111, G]] ri:1=1[d (mi , si)] où : - mi sont les positions des points du maillage 3 dans la zone d'observation 1, - s, sont les positions ou les trajectoires des moyens de mesure présentant la distribution 10 spatiale donnée S, - V, sont les valeurs modélisées de l'ensemble V extraites pour les positions ou pour les trajectoires si, - d est une distance entre une position d'un point du maillage 3 et une position ou une trajectoire des moyens de mesure dans la zone d'observation 1. 15 La distance d peut être une distance routière ou une distance euclidienne (i.e. à vol d'oiseau). Fonction-coût Lorsque la fonction-coût est représentative de la différence entre 17 et V, la fonction-coût est avantageusement une norme de la différence entre V et V. Lorsque la fonction-coût est 20 représentative de la différence entre V et V, les moyens de calcul sont avantageusement configurés pour calculer la fonction-coût selon la formule : G (S) = ril)2 j=1 où g sont les valeurs modélisées de l'ensemble V extraites pour les points du maillage 3. 25 Lorsque la fonction-coût est représentative de la vraisemblance entre I7 et V, la fonction- coût comprend avantageusement un produit scalaire entre V et V, de préférence normalisé. Algorithme d'extraction de Sopt Lorsque la fonction-coût est représentative de la différence entre I7 et V, les moyens de calcul sont avantageusement configurés pour effectuer une modification aléatoire de la 30 distribution spatiale donnée S selon une probabilité, notée p, vérifiant : 3065529 14 pii = 1 Si «S(i)) (p(S(0) où : - S(i) est une distribution spatiale donnée initiale, 5 - S(i) est une distribution spatiale donnée modifiée aléatoirement, - c est un paramètre. Si l'on considère une itération, notée t, le paramètre c est de préférence choisi de sorte que : c(t 1) = a c(t) ; a1. Les moyens de calcul sont configurés pour itérer la modification aléatoire de la 10 distribution spatiale donnée S jusqu'à l'extraction de Sol, Les moyens de calcul sont configurés pour appliquer un critère d'arrêt aux itérations. A titre d'exemple, le critère d'arrêt peut être défini lorsque la différence frp(SU+1)) - (p(SU))1 est inférieure à un seuil prédéterminé. 15 Lorsque la fonction-coût est représentative de la vraisemblance entre ri et V, les moyens de calcul sont avantageusement configurés pour effectuer une modification aléatoire de la distribution spatiale donnée S selon une probabilité, notée Pi, vérifiant : Pji = 1 si p(S(1)) ye(i)) = exp (q)(SW) - (p(S(0)) si (p(sw) REVENDICATIONS1. Système de mesure d'au moins une grandeur physique représentative de la qualité de l'air dans une zone d'observation (1), le système comportant : - une cartographie (2) de la zone d'observation (1), comprenant un ensemble, noté V, de valeurs modélisées représentatives de la grandeur physique ; - des moyens de mesure de la grandeur physique, possédant un nombre N de positions ou un nombre N de trajectoires dans la zone d'observation (1), les N positions ou les N trajectoires étant destinées à présenter une distribution spatiale, notée S.pt , dans la zone 10 d'observation (1) ; - des moyens de calcul de la distribution spatiale Sopt, configurés pour : - effectuer un maillage (3) de la zone d'observation (1), le maillage (3) comprenant un nombre G de points ; - calculer, pour une distribution spatiale donnée des N positions ou des N 15 trajectoires, notée S, un estimateur de l'ensemble V, noté g, pour chacun des G points du maillage (3) ; - calculer une fonction-coût, notée p(S), représentative de la différence ou de la vraisemblance entre 17 et les valeurs modélisées de l'ensemble V extraites aux G points du maillage, notées V; 20 - extraire la distribution spatiale Sot, qui minimise ou maximise la fonction-coût selon que la fonction-coût est représentative de la différence ou de la vraisemblance entre V et V. 2. Système selon la revendication 1, dans lequel les moyens de calcul sont configurés pour 25 calculer l'estimateur de l'ensemble V selon la formule : - - ,_,-N E où : - Vi sont les valeurs modélisées de l'ensemble V extraites pour les positions ou pour les trajectoires des moyens de mesure présentant la distribution spatiale donnée S, 30 - w est une fonction de transfert entre si et mi, où si sont les positions ou les trajectoires des moyens de mesure présentant la distribution spatiale donnée S, et où mi sont les positions des points du maillage (3) dans la zone d'observation (1). 3065529 17 3. Système selon la revendication 1 ou 2, dans lequel les moyens de calcul sont configurés pour calculer l'estimateur de l'ensemble V selon la formule : Eiv=i [c/ (mi, si)] -2 _2 E [[1, G]] où : 5 - mi sont les positions des points du maillage (3) dans la zone d'observation (1), - s, sont les positions ou les trajectoires des moyens de mesure présentant la distribution spatiale donnée S, - V, sont les valeurs modélisées de l'ensemble V extraites pour les positions ou pour les trajectoires si, - d est une distance entre une position d'un point du maillage (3) et une position ou une trajectoire des moyens de mesure dans la zone d'observation (1). 4. Système selon l'une des revendications 1 à 3, dans lequel la fonction-coût est représentative de la différence entre V et V, et la fonction-coût est une norme de la différence entre 1-7 et V. 5. Système selon l'une des revendications 1 à 4, dans lequel la fonction-coût est représentative de la différence entre V et V, et dans lequel les moyens de calcul sont configurés pour calculer la fonction-coût selon la formule : G 1 où ri sont les valeurs modélisées de l'ensemble V extraites pour les points du maillage (3). 6. Système selon l'une des revendications 1 à 5, dans lequel la fonction-coût est représentative de la différence entre V et V, et dans lequel les moyens de calcul sont 25 configurés pour effectuer une modification aléatoire de la distribution spatiale donnée S selon une probabilité, notée ni, vérifiant : pii = 1 si (p(SW) cp(S(i)) pii = exp( s - q)((0) cp(SU))) si (feu)) > (tem) où : j [cl (mi, si)] 30 - S(i) est une distribution spatiale donnée initiale, 3065529 18 - 50) est une distribution spatiale donnée modifiée aléatoirement, - c est un paramètre ; et dans lequel les moyens de calcul sont configurés pour itérer la modification aléatoire de la distribution spatiale donnée S jusqu'à l'extraction de Sept. 5 7. Système selon l'une des revendications 1 à 3, dans lequel la fonction-coût est représentative de la vraisemblance entre 11 et 7, et la fonction-coût comprend un produit scalaire entre f1 et V, de préférence normalisé. 10 8. Système selon l'une des revendications 1 à 3 ou selon la revendication 7, dans lequel la fonction-coût est représentative de la vraisemblance entre 11 et V, et dans lequel les moyens de calcul sont configurés pour effectuer une modification aléatoire de la distribution spatiale donnée S selon une probabilité, notée pi' vérifiant : Pji = 1 si cp(SW) > ye(i)) 15 = exp (q)(S(D) - v(S(0)) si v(sw) 9. Système selon l'une des revendications 1 à 8, dans lequel les moyens de mesure sont agencés dans la zone d'observation (1) de manière à posséder N positions ou N trajectoires 25 présentant la distribution spatiale Sept. 10. Système selon l'une des revendications 1 à 8, comportant des moyens d'activation configurés pour activer les moyens de mesure possédant N positions ou N trajectoires présentant la distribution spatiale la plus proche de Sopt- 30 11. Système selon l'une des revendications 1 à 10, dans lequel la zone d'observation (1) comporte un nombre M de positions autorisées ou un nombre M de trajectoires autorisées, 3065529 19 vérifiant M>N, dans lesquelles les moyens de mesure sont respectivement autorisés à posséder des positions ou des trajectoires. 12. Système selon l'une des revendications 1 à 11, dans lequel les moyens de mesure 5 comportent au moins N dispositifs de mesure équipant chacun un véhicule, le véhicule étant de préférence routier ou aérien. 13. Système selon l'une des revendications 1 à 11, dans lequel les moyens de mesure comportent au moins N dispositifs de mesure équipant chacun un mobilier urbain, le 10 mobilier urbain étant de préférence sélectionné dans le groupe comportant des feux de circulation et des abribus. 14. Système selon l'une des revendications 1 à 13, dans lequel les moyens de mesure comportent des capteurs spectroscopiques, de préférence des capteurs infrarouges non-15 dispersifs. 15. Procédé de mesure d'au moins une grandeur physique représentative de la qualité de l'air dans une zone d'observation (1), le procédé comportant les étapes : a) prévoir une cartographie (2) de la zone d'observation (1), comprenant un ensemble, 20 noté V, de valeurs modélisées représentatives de la grandeur physique ; b) prévoir des moyens de mesure de la grandeur physique, possédant un nombre N de positions ou un nombre N de trajectoires dans la zone d'observation (1), les N positions ou les N trajectoires étant destinées à présenter une distribution spatiale, notée S.pt ; c) effectuer un maillage (3) de la zone d'observation (1), le maillage (3) comprenant un 25 nombre G de points ; d) calculer, pour une distribution spatiale donnée des N positions ou des N trajectoires, notée S, un estimateur de l'ensemble V, noté V, pour chacun des G points du maillage (3) ; e) calculer une fonction-coût, notée p(S), représentative de la différence ou de la vraisemblance entre I7 et les valeurs modélisées de l'ensemble V extraites aux G points du 30 maillage (3), notées V; f) extraire la distribution spatiale Sopt qui minimise ou maximise la fonction-coût selon que la fonction-coût est représentative de la différence ou de la vraisemblance entre 17 et V; g) agencer les moyens de mesure dans la zone d'observation (1) de sorte que les N positions ou les N trajectoires présentent la distribution spatiale Sopt.