FR 3037701 A1 20161223 FR 1555518 A 20150616 1 Procédé d'estimation et de prévision d'indicateur de production d'énergie d'un système solaire. [001] Domaine technique [002] L'invention concerne le domaine de l'estimation et de la prévision de production d'énergie d'un ou plusieurs systèmes de conversion d'énergie solaire (nommés systèmes solaires par la suite) comme, par exemple, des panneaux solaires ou des capteurs solaires thermiques. [003] Etat de la technique [4] Les moyens de conversion des systèmes solaires utilisés pour la génération d'énergie telle que de l'énergie électrique ou thermique, possèdent une production qui dépend de l'irradiation solaire incidente à laquelle ils sont exposés (nommée irradiation ci-après). Typiquement, un panneau solaire produit une énergie électrique qui varie en fonction des conditions d'ensoleillement et de luminosité de l'environnement proche du panneau. [5] Une prévision de l'irradiation du système solaire permet d'anticiper la quantité d'énergie qui sera produite sur un horizon de temps futur. Il existe différentes approches pour prévoir l'irradiation d'un tel système dont notamment : les prévisions utilisant des modèles météorologiques (horizon de prévision allant de quelques heures à plusieurs jours) ; les prévisions basées sur des mesures d'irradiation prises par un instrument de mesure au sol (tel qu'un capteur pyranométrique), des modèles mathématiques étant utilisés pour prévoir les valeurs des prochaines mesures sur la base d'un historique de mesures (horizon de prévision allant de quelques minutes à plusieurs heures) ; 3037701 2 les prévisions par traitement d'images satellitaires, estimant le rayonnement solaire reçu au sol à partir de données satellitaires (horizon de prévision de l'ordre de quelques minutes à quelques heures) ; les prévisions utilisant les images du ciel prises à partir du sol pour déterminer 5 un déplacement de couvertures nuageuses et en déduire l'irradiation au sol (horizon de prévision pouvant aller jusqu'à quelques minutes). [6] Avec ces dernières, i.e. les prévisions par traitement d'images du ciel prises à partir du sol, l'estimation et la prévision du rayonnement solaire au sol peut être effectuée avec une bonne précision, que ce soit en temps réel ou sur la base d'une 10 prévision à horizon de temps proche (5 à 15 minutes par exemple). Il est entendu par horizon de prévision l'intervalle de temps entre le moment où la prévision est effectuée et le moment auquel se passe l'évènement prévu. [7] Généralement, les prévisions basées sur des images du ciel prises à partir du sol comprennent au moins trois phases qui sont : 15 - la localisation d'un masque nuageux, - l'estimation du mouvement du masque nuageux et - la détermination du rayonnement solaire au sol en fonction du ou des masques nuageux détectés et de leurs mouvements estimés. [8] Lors de la première phase, la présence d'un masque nuageux est classiquement 20 reconnue à partir d'une image prise depuis une zone d'intérêt (à proximité de panneaux solaires typiquement). A partir de la localisation du masque nuageux sur l'image, un indice de couverture nuageuse au niveau de la zone d'intérêt peut être calculé. [9] Dans la deuxième phase, le déplacement du ou des masques détectés est étudié. [10] Dans un troisième temps, les images sont traduites en valeurs de rayonnement 25 au sol en fonction des couvertures nuageuses détectées et de leurs déplacements. A cet effet, des méthodes de traitement connues en soi peuvent être utilisées pour convertir le niveau de chaque pixel d'images en une valeur brute d'irradiation au sol. A titre purement illustratif, une des méthodes de traitement possibles est décrite dans le document US20100309330. 3037701 3 [011] Toutefois, il s'avère que les couvertures nuageuses sont habituellement considérées comme impactant le rayonnement au sol indépendamment de leur position par rapport à la zone d'intérêt. Or, une couverture nuageuse qui est à l'horizon terrestre par rapport à la zone d'intérêt n'a pas une même influence sur l'irradiation perçue au 5 sol qu'une autre couverture nuageuse qui serait quant à elle au zénith par exemple. Il s'avère qu'une couverture nuageuse au zénith a des pixels qui impactent plus directement le rayonnement perçu au sol. Des divergences importantes peuvent ainsi apparaître entre la valeur de rayonnement au sol estimée ou prévue et l'irradiation réellement perçue sur la zone d'intérêt. 10 [012] Ces divergences créent un biais dans l'estimation et la prévision de la production d'énergie des systèmes solaires et, par répercussion, nuit à l'optimisation de la gestion d'énergie (avec en l'occurrence une évaluation erronée de la quantité d'énergie réellement produite et/ou à produire en compensation par d'autres dispositifs de production). 15 [013] Il existe ainsi un besoin de mettre en place un procédé plus précis d'estimation et de prévision de production d'énergie d'un système solaire. [14] L'invention vient améliorer la situation en ce sens. [15] Résumé de l'invention 20 [016] L'objet de l'invention est de remédier aux inconvénients précités en proposant notamment de prendre en compte la localisation des nuages sur une image du ciel prise à partir du sol en vue d'estimer la production d'énergie d'un système solaire. [017] A cet effet, un premier aspect de l'invention concerne un procédé d'estimation d'indicateur de production d'énergie d'un système solaire comprenant les étapes de : 25 a) réception d'images du ciel prises au sol depuis une zone terrestre comportant le système solaire, chaque image reçue comprenant au moins un repère optique ; b) extraction des valeurs de pixels et localisation de pixels d'une image issue des images reçues ; 3037701 4 c) détermination des distances d'éloignement entre les pixels localisés et l'au moins un repère optique de l'image issue des images reçues ; d) génération d'un vecteur de paramètres à partir d'une 5 classification des pixels selon les valeurs extraites et les distances d'éloignement déterminées ; e) comparaison des paramètres du vecteur généré avec des paramètres prédéterminés, les paramètres prédéterminés étant respectivement associés à des indicateurs de production d'énergie ; 10 f) estimation de l'indicateur de production d'énergie du système solaire à partir des indicateurs de production d'énergie des paramètres prédéterminés comparés au vecteur généré. [018] Le repère optique des images reçues peut être par exemple : - le centre optique de la caméra ou de l'appareil photo dont est issu l'image, 15 notamment lorsque cette caméra ou cet appareil est muni d'un objectif ayant un angle de champ supérieur à 170 degrés (appareil de prise d'image à objectif de type « fish-eye », qui signifie oeil de poisson en anglais), - le centre de l'image lorsqu'elle est prise par une caméra ou un appareil photo muni d'un objectif ayant un angle de champ inférieur à 170 degrés, 20 - la localisation du soleil sur l'image issue des images reçues, - ou autre. [019] Selon une réalisation préférée, lorsque le repère optique est la localisation du soleil, la distance au soleil est déterminée entre les pixels et le centre du soleil sur l'image du ciel. 25 [020] L'indicateur de production d'énergie est une donnée qui peut être directement relative à la production d'énergie du système solaire ou, être une donnée correspondant à une grandeur qui impacte la production d'énergie. L'indicateur d'énergie peut en particulier être : 3037701 5 - un rayonnement solaire au niveau du système solaire ou de la zone terrestre depuis laquelle les images reçues sont prises ; - une production d'énergie électrique du système solaire ; - une production d'énergie thermique du système solaire ; 5 - ou autre donnée relative à une production d'énergie s'expliquant au premier ordre à l'aide du rayonnement solaire. [021] Selon des réalisations préférées mais non limitatives, les indicateurs de production d'énergie sont des données de rayonnement solaire et/ou de production d'énergie électrique. 10 [022] La classification des pixels en fonction de leurs valeurs respectives et distances d'éloignement par rapport au repère optique permet d'identifier les éléments du ciel auxquels correspondent les pixels de l'image (ciel dégagé, nuage ou soleil typiquement) et la localisation de ces éléments dans le ciel par rapport à la zone dans laquelle le système solaire est installé. 15 [023] Les distances déterminées entre les pixels et le repère optique permettent notamment de qualifier les contributions respectives des éléments du ciel identifiés à la production d'énergie du système solaire. A titre d'exemple, plus la distance d'éloignement entre le pixel correspondant à un nuage et le repère optique est grande, moins la contribution du pixel pour l'estimation de la production d'énergie peut être 20 importante, et inversement. [024] L'estimation d'indicateur de production d'énergie peut ainsi être basée sur des indicateurs de production d'énergie prédéterminés (irradiation au niveau de panneaux solaires, production d'énergie électrique effective de panneaux solaires ou autre) dont l'importance peut être modérée par des coefficients relatifs, voire proportionnels, aux 25 distances d'éloignement des pixels par rapport au centre optique de l'image ou à un ou plusieurs points d'intérêt de cette image, comme la localisation du soleil. L'estimation est alors plus précise qu'avec l'attribution d'une même contribution pour tous les pixels de l'image comme dans l'état de la technique, ou tout au moins l'attribution d'une même contribution indépendamment des positions dans l'image des éléments de ciel détectés. 3037701 6 [025] Un poids particulier peut notamment être attribué aux indicateurs de production d'énergie issus de paramètres correspondant à des nuages clairs. Ce poids particulier permet de tenir compte dans l'estimation de la contribution supplémentaire de nuages clairs. En effet, de tels nuages clairs sont des sources secondaires importantes de 5 rayonnement solaire, d'intensité d'autant plus importante que leur distance au soleil est proche sur l'image. À ce titre, le poids des indicateurs associés aux paramètres prédéterminés correspondant à des nuages clairs peuvent être augmentés dans l'estimation. Cette pondération complémentaire rend l'estimation d'indicateur de production d'énergie d'autant plus précise. 10 [026] On entend par système solaire tout système configuré pour transformer l'énergie solaire en une autre énergie, par exemple une énergie électrique ou thermique. Un système solaire peut être par exemple un ou plusieurs panneaux photovoltaïques, panneaux thermiques, chauffe-eau solaires ou autre. [27] Selon une réalisation avantageuse, à l'étape f), l'estimation d'indicateur de 15 production d'énergie est en outre issue d'une pondération des indicateurs de production associés aux paramètres prédéterminés, la pondération étant fonction d'écarts de valeurs entre les paramètres du vecteur et les paramètres prédéterminés. [28] Il est alors possible d'attribuer des poids plus ou moins forts aux paramètres prédéterminés en fonction de leur pertinence par rapport au vecteur de paramètres, de 20 sorte notamment à cibler plus particulièrement les paramètres prédéterminés qui sont égaux ou proches des paramètres du vecteur. A titre d'exemple purement illustratif, le poids des indicateurs associés aux paramètres prédéterminés les plus proches du vecteur peuvent être augmentés dans l'estimation, et le poids des indicateurs associés aux autres paramètres prédéterminés peut être diminué. Une telle pondération permet 25 d'améliorer encore la qualité de l'estimation de production d'énergie. [29] Selon différentes réalisations possibles, la pondération peut être appliquée aux indicateurs associés à : - l'ensemble de paramètres prédéterminés, ou - des paramètres prédéterminés avec des valeurs proches du vecteur de 30 paramètres généré comme par exemple des paramètres prédéterminés ayant des 3037701 7 valeurs comprises dans une plage prédéfinie autour des valeurs du vecteur, telle qu'une plage comprise entre +30% et -30% des valeurs du vecteur généré. [030] Avantageusement, à l'étape d), le vecteur de paramètres dans lequel sont classifiés les pixels comprend au moins : 5 ^ une première catégorie de paramètres correspondant à un ciel clair, ^ une deuxième catégorie de paramètres correspondant au soleil, ^ une troisième catégorie de paramètres correspondant aux 10 nuages. [31] Dans ce mode de réalisation, les pixels sont ainsi classifiés selon trois catégories permettant de définir à quel élément du ciel ils correspondent. [32] Selon une réalisation possible, seuls les paramètres correspondant aux nuages sont comparés aux paramètres prédéterminés. Ceci permet de limiter le nombre de 15 calculs (par rapport à une comparaison avec tous les paramètres à l'étape e)) tout en conservant une estimation précise de la production d'énergie en se basant sur l'une des catégories d'élément du ciel qui impacte le plus directement les variations de production d'un système solaire. [33] Avantageusement, la troisième catégorie de paramètres comprend : 20 ^ une première sous-catégorie de paramètres correspondant aux nuages clairs, dans laquelle sont classifiés des pixels de valeur supérieure ou égale à un seuil, ^ une deuxième sous-catégorie de paramètres correspondant aux nuages sombres, dans laquelle sont classifiés des 25 pixels de valeur inférieure au seuil. [034] A titre d'exemple purement illustratif, le seuil correspond à une valeur de luminance du pixel dont le signal d'intensité I est défini par la combinaison linéaire suivante (où les coefficients sont des pondérations universelles standardisées ITURBT709) : 30 / = 0.2126 Ir + 0.7152 Iv + 0.0722 lb 3037701 8 [35] Les première et deuxième sous-catégories permettent de faire une distinction entre les nuages clairs et les nuages sombres, lesquels n'ont pas le même impact sur la production d'énergie d'un système solaire. En effet, un nuage clair peut se comporter comme une source de rayonnement solaire secondaire alors qu'un nuage sombre peut 5 quant à lui absorber complètement les rayonnements solaires. [36] Selon un mode de réalisation de l'invention, le procédé comprend en outre une étape de détermination d'une proportion de couverture nuageuse sur l'image issue des images reçues en fonction des pixels classifiés dans la troisième catégorie de paramètres. 10 [037] La proportion de couverture nuageuse peut être utilisée pour des évaluations des conditions météorologiques, en vue notamment de vérifier ou prévoir la présence dans le ciel de masque nuageux. Par exemple, si la proportion de pixels classifiés dans la troisième catégorie (i.e. catégorie « nuage ») est importante (nombre de pixels de l'image appartenant la catégorie « nuage » est égal ou supérieur à un seuil prédéfini), 15 l'image issue des images reçues peut être considérée comme une image d'un ciel nuageux. [038] Par ailleurs, cette évaluation peut être affinée en segmentant l'image en plusieurs parties (huit en l'occurrence) et en déterminant la proportion de couverture nuageuse dans chaque partie de l'image. 20 [039] Selon une réalisation avantageuse, à l'étape b), l'image issue des images reçues est une image générée selon les étapes suivantes : - récupération parmi les images reçues d'images prises à au moins un premier instant et un deuxième instant, le premier instant étant antérieur au deuxième instant ; 25 - localisation des valeurs de pixels sur les images prises respectivement au premier et au deuxième instant ; - comparaison des localisations des valeurs de pixels au premier instant et au deuxième instant ; - détermination des déplacements, entre les premier et deuxième 30 instants, d'au moins une partie des valeurs de pixels localisées ; 3037701 9 - génération de l'image du ciel pour un troisième instant en fonction des déplacements déterminés, le troisième instant étant postérieur au deuxième instant. [040] Le troisième instant peut correspondre notamment à un horizon de temps donné 5 (i.e. instant futur pour lequel l'estimation d'indicateur de production d'énergie du système solaire est souhaitée). En appliquant les étapes b) à e) du procédé d'estimation à une image générée pour le troisième instant permet de prévoir la production d'énergie à un instant futur. Le procédé proposé peut ainsi être un procédé d'estimation et de prévision d'indicateur de production d'énergie. 10 [041] Avantageusement, le procédé comprend en outre une étape de détermination d'une plage d'indicateur de production d'énergie qui est fonction d'un écart type calculé à partir de la racine d'une moyenne pondérée d'écarts quadratiques entre des indicateurs de production des paramètres prédéterminés et la production estimée, la moyenne étant pondérée en fonction d'écarts de valeurs entre les paramètres du 15 vecteur généré et les paramètres prédéterminés. [42] Cet écart type permet de définir un intervalle de confiance autour de l'indicateur de production d'énergie estimé, conférant à l'estimation une borne supérieure et une borne inférieure de production d'énergie attendue. [43] Le procédé peut en outre comprendre une étape de réception d'une donnée de : 20 ^ rayonnement solaire effectif au sol dans la zone terrestre, ou ^ production d'énergie effective du système solaire. La donnée reçue peut en particulier être utilisée pour étalonner ou corriger les indicateurs de production d'énergie associés aux paramètres prédéterminés. [044] De cette manière, les paramètres prédéterminés peuvent être étalonnés, 25 réinitialisés et/ou ajustés selon des mesures d'irradiation ou de production effective au sol, permettant d'affiner les données utilisées pour l'estimation et de fait, rendre la prévision de production d'énergie du système solaire encore plus précise et robuste. Par exemple, la donnée reçue peut être issue d'un moyen de mesure pyranométrique ou d'une mesure de la production d'énergie du système solaire lui-même. 3037701 10 [045] Selon un deuxième aspect, l'invention propose un dispositif d'estimation d'indicateur de production d'énergie d'un système solaire, le dispositif d'estimation comprenant au moins : - une interface d'entrée pour la réception d'images du ciel prises au sol 5 depuis une zone terrestre comportant le système solaire ; - une mémoire configurée pour stocker : ^ les images du ciel reçues, et ^ des paramètres prédéterminés respectivement associés à des indicateurs de production d'énergie ; 10 - une unité de traitement de données connectée à l'interface d'entrée et à la mémoire. [46] L'unité de traitement de données étant configurée pour mettre en oeuvre le procédé précité pour estimer l'indicateur de production d'énergie du système solaire. [47] Selon un mode de réalisation du dispositif, l'interface d'entrée est : 15 o connectée à un appareil de prise d'image orienté vers le ciel au niveau de la zone terrestre, l'appareil de prise d'image étant muni d'un objectif ayant un angle de champ supérieur à 170 degrés ; et/ou o configurée pour recevoir des données de : ^ rayonnement solaire effectif au sol dans la zone terrestre, ou 20 ^ production d'énergie effective du système solaire. [48] La donnée reçue peut être plus particulièrement une donnée émise par un pyranomètre placé au voisinage du système solaire. [49] La présente invention vise également un programme informatique comportant des instructions pour la mise en oeuvre du procédé précédemment décrit, lorsque ce 25 programme est exécuté par une unité de traitement de données telle qu'un processeur. [50] Ce programme peut utiliser n'importe quel langage de programmation (par exemple, un langage objet ou autre), et être sous la forme d'un code source interprétable, d'un code partiellement compilé ou d'un code totalement compilé. 3037701 11 [51] Les figures 3, 5 et 7 décrites en détails ci-après, peuvent former des organigrammes de l'algorithme général d'un tel programme informatique. [52] Brève description des figures 5 [053] D'autres caractéristiques et avantages de l'invention apparaîtront à l'examen de la description détaillée ci-après, et des figures annexées sur lesquelles : - la figure 1 illustre un exemple de dispositif d'estimation d'un indicateur de production d'énergie d'un système solaire selon un exemple de réalisation de l'invention ; 10 - les figures 2 (figures 2A à 2D) illustrent des exemples dans lesquels des images du ciel sont prises au sol depuis une zone terrestre où se trouve des panneaux photovoltaïques ; - la figure 3 représente un organigramme composé d'un exemple de succession d'étapes du procédé d'estimation selon l'invention ; 15 - la figure 4 est une représentation graphique d'une base de données comportant des paramètres prédéterminés ; - la figure 5 représente un organigramme composé d'un exemple de succession d'étapes pour la construction d'une base de données de paramètres prédéterminés ; 20 - la figure 6 représente un organigramme composé d'un exemple de succession d'étapes permettant de générer une image de ciel à partir d'images de ciel prises au sol depuis la zone terrestre ; - les figures 7 (figures 7A à 7C) illustrent des exemples d'images prises à partir du sol de la zone de terrestre dans laquelle est installé le système solaire et 25 d'une image généré à partir de ces images ; et - la figure 8 représente des exemples de courbes de rayonnements solaires mesurés, estimés et prévus selon le procédé d'estimation proposé. 3037701 12 [054] Pour des raisons de clarté, les dimensions des différents éléments représentés sur ces figures ne sont pas en proportion avec leurs dimensions réelles. Sur les figures, des références identiques correspondent à des éléments identiques pour les différents modes de réalisation exposés. 5 [55] Description détaillée [56] On se réfère tout d'abord à la figure 1 qui illustre un exemple de dispositif d'estimation d'un indicateur de production d'énergie d'un système solaire selon un exemple de réalisation de l'invention. 10 [057] Dans ce mode de réalisation, le dispositif comporte un ordinateur 100, comprenant une mémoire 102 pour stocker : - des instructions permettant la mise en oeuvre du procédé, - d'images du ciel, - de paramètres prédéterminés associés à des indicateurs de production 15 d'énergie, et - de données temporaires pour réaliser les différentes étapes du procédé tel que décrit précédemment et détaillé plus loin. [058] L'ordinateur comporte en outre une unité de traitement de données 104. Cette unité de traitement de données peut être un circuit comme par exemple : 20 - un processeur apte à interpréter des instructions sous la forme de programme informatique, ou - une carte électronique dont les étapes du procédé de l'invention sont décrites dans le silicium, ou encore - une puce électronique programmable comme une puce FPGA (pour « Field- 25 Programmable Gate Array » en anglais). [059] Cet ordinateur 100 comporte aussi : o une interface d'entrée 106 pour : 3037701 13 ^ la réception d'images du ciel prises au sol depuis une zone terrestre comportant le système solaire, par exemple par un appareil de prise d'images telle qu'une caméra 109 ou un appareil photo, 5 ^ d'une donnée de rayonnement solaire effectif au sol dans la zone terrestre ou de production d'énergie effective du système solaire, o une interface de sortie 108 pour la fourniture à un dispositif distant 116 de commandes de contrôle ou de résultats d'estimation selon le procédé. 10 [060] Le système solaire peut par exemple être un ou plusieurs panneaux solaires 110 installés dans la zone terrestre depuis laquelle les images du ciel sont prises par la caméra 109. Le système solaire peut être configuré pour mesurer une production d'énergie effective et pour la communiquer à l'ordinateur 100 via une liaison de communication avec l'interface d'entrée 106. 15 [061] La caméra 109 peut être située à une même localisation que le système solaire ou dans une zone autour du système solaire, par exemple à une localisation située dans un rayon inférieur ou égal à quelques kilomètres (1 à 10 km par exemple) par rapport aux panneaux solaires 110. Une calibration de la caméra au préalable peut être effectuée de sorte à mettre en correspondance chaque pixel de l'image avec des angles 20 de visée (élévation, azimut) déterminés pour que les images prises soient représentatives du ciel au-dessus du système solaire. [62] L'ordinateur peut par ailleurs comporter, pour permettre une interaction aisée avec un utilisateur, un écran 112 et un clavier 114. Bien entendu, le clavier est facultatif, notamment dans le cadre d'un ordinateur ayant la forme d'une tablette tactile 25 ou d'un Smartphone, par exemple. [63] Typiquement, le dispositif distant 116 peut être une plateforme de service ou un terminal de client apte à recevoir les estimations envoyées par le dispositif via l'interface de sortie 108. A cet effet, le dispositif 116 peut comporter une interface de communication apte à recevoir des données communiquées par le dispositif et une 30 unité de traitement de données pour les interpréter. Ainsi, les données d'estimation 3037701 14 d'un indicateur de production d'énergie reçues par le dispositif distant 116 peuvent être utilisées par ce dernier pour gérer et anticiper la production du système solaire. A titre d'exemple purement illustratif, un utilisateur peut recevoir sur son téléphone les estimations de production des panneaux photovoltaïques installés sur le toit de son 5 domicile. Ainsi, l'utilisateur peut optimiser sa consommation en anticipant la production d'énergie des panneaux selon la production estimée. Des services variés peuvent être proposés à des utilisateurs en fonction de différentes interprétations possibles des estimations reçues par le dispositif distant 116. [64] L'interface 106 peut en outre recevoir des données d'un pyranomètre relatives 10 au rayonnement solaire effectif au sol ou des données de production d'énergie du panneau solaire lui-même. Ces données peuvent par exemple être une valeur brute d'irradiation mesurée par le pyranomètre à proximité des panneaux solaires 110 ou une déduction d'une telle valeur brute d'irradiation en fonction de la production d'énergie effective des panneaux 110. Cette donnée peut notamment permettre de corriger les 15 indicateurs de production d'énergie associés aux paramètres prédéterminés, comme décrit plus loin. Une telle correction peut être effectuée une seule fois pour un étalonnage des données ou bien au fur et à mesure pour une correction rétroactive. [65] Aux figures 2 (figures 2A à 2D), il est illustré des exemples dans lesquels des 20 images du ciel sont prises au sol depuis une zone terrestre où se trouvent des panneaux photovoltaïques tels que les panneaux 110, et qui sont reçues par l'ordinateur 100 via l'interface 106. [66] La figure 2A représente une situation dans laquelle les panneaux solaires 110 sont par exemple installés sur le toit d'un bâtiment 201 (le domicile d'un utilisateur 25 typiquement). Dans cet exemple, la caméra 109 est placée à proximité des panneaux solaires 110, sur le toit du bâtiment 201, et est orientée vers le ciel afin de réaliser des captures d'image à partir du sol de la zone terrestre dans laquelle se trouve les panneaux 110. [67] Toutefois, selon d'autres réalisations possibles, la caméra 109 peut être placée : 30 à une position confondue avec celle des panneaux solaires 110 ou 3037701 15 - dans un périmètre approprié autour des panneaux. [68] La zone terrestre dans laquelle se trouve un système solaire tel que les panneaux 110 peut être définie par un espace autour de ce système dans lequel une image du ciel prise depuis le sol correspond au ciel qui est réellement au-dessus du 5 système solaire. En effet, en vue de l'estimation de production du système solaire, il convient que l'image prise depuis la zone terrestre soit représentative du ciel effectivement au-dessus du système. [69] Il convient néanmoins de noter que plus la prise d'images du ciel est effectuée près du système solaire, plus les images prises seront représentatives du ciel au-dessus 10 du système solaire, permettant d'obtenir une meilleure qualité d'estimation de l'indicateur de production d'énergie. [70] La caméra 109 est ici une caméra comportant un objectif de type « fish-eye » (ayant un angle de champ supérieur à 170 degrés). Selon d'autres alternatives possibles, la caméra peut être un appareil de prise de vue (un appareil photo ou une 15 webcam par exemple) ayant un objectif de type « fish-eye » ou un objectif ayant un angle de champ inférieur à 170 degrés. La caméra 109 est orientée vers le ciel pour prendre des images de celui-ci avec des éléments de ciel tels que le soleil 205, des nuages clairs 207 et des nuages sombres 208. [71] A la figure 2B, une image 200 prise par la caméra 109 est illustrée. 20 L'image 200 comprend un premier repère optique 210 correspondant au centre optique de la caméra 109 et un deuxième repère optique 212 correspondant à la localisation du centre du soleil 205 sur l'image. L'image 200 est représentative des conditions de ciel au-dessus des panneaux 110. A partir de l'image 200, il est possible d'analyser les valeurs des pixels pour déterminer à quel type d'élément du ciel (ciel dégagé, soleil 25 205, nuages clairs 207, nuages sombres 208) correspond chaque pixel de l'image. [72] La figure 2C illustre un autre exemple d'image du ciel prise par la caméra 109 et la figure 2D représente cette image suite à un traitement des valeurs des pixels de l'image. Ce traitement peut notamment consister à : - identifier les valeurs des pixels de l'image, 3037701 16 - déterminer à quel élément de ciel les valeurs de pixels correspondent (ciel clair, nuage clair, nuage sombre, disque solaire, région circumsolaire), - classer dans un vecteur de paramètres les pixels selon l'élément de ciel déterminé pour leur valeur, 5 - afficher les pixels de l'image selon une couleur correspondant à l'élément de ciel auquel ils ont été attribués (par exemple rouge pour le disque du soleil, jaune pour la région circumsolaire, bleu pour le ciel clair, cyan pour les nuages clairs et violet pour les nuages sombres). [073] Ce traitement de l'image prise par la caméra correspond à une segmentation des 10 pixels selon des catégories d'élément de ciel correspondant au moins au soleil, à des nuages et à un ciel clair. Elles peuvent être complétées par d'autres éléments comme ceux précités : nuages clairs ou sombres, centre du soleil, disque solaire, région circumsolaire, etc. 15 [074] On se réfère maintenant à la figure 3 représentant un organigramme composé d'un exemple de succession d'étapes du procédé d'estimation d'indicateur de production d'énergie de système solaire. [075] L'étape 300 correspond à la réception d'images du ciel prises par la caméra 109. Les images reçues comprennent au moins un repère optique qui peut être : 20 - le centre optique de la caméra lorsque l'objectif est de type « fish-eye » - le centre de l'image lorsque l'objectif de la caméra à un angle de champ inférieur à 170 degrés, - la localisation du soleil sur l'image. [076] Selon une étape 310, les images reçues peuvent être traitées de sorte à extraire 25 des valeurs de pixels et à localiser des pixels des images. L'image issue des images reçues peut être directement l'une des images reçues ou une image générée à partir d'images reçues, comme détaillé plus loin. Il est alors possible de déterminer par rapport au(x) repère(s) optique(s) de l'image, la distance d'éloignement des pixels 3037701 17 localisés. Cette distance est typiquement déterminée en nombre de pixels qui séparent le pixel localisé et le repère optique, ou alors selon un angle d'élévation ou autre. [77] Comme précité, le repère optique sur l'image peut être le soleil. La distance au soleil des pixels correspondant aux nuages clairs et sombres peut aussi être déterminée. 5 Cette distance peut notamment être déterminée entre le centre du soleil et la position du pixel sur l'image. [78] Ensuite, une étape 320 de génération d'un vecteur de paramètres peut consister à classer des pixels de l'image en fonction des valeurs extraites et des distances d'éloignement déterminées. Les paramètres contenus dans le vecteur peuvent 10 respectivement correspondre à des catégories d'éléments de ciel telles que celles précitées : ciel clair, nuage clair, nuage sombre, disque solaire, région circumsolaire. [79] Le vecteur de paramètres dans lequel sont classifiés les pixels peut notamment comprendre : ^ une première catégorie de paramètres correspondant à un 15 ciel clair, dans laquelle sont classés les pixels dont la valeur extraite : o est comprise dans une première plage de valeurs de pixels ou o dépasse un premier seuil de valeur, 20 ^ une deuxième catégorie de paramètres correspondant au soleil, dans laquelle sont classés les pixels qui sont proches du soleil (par exemple dans un rayon déterminé) dont la valeur extraite : o est comprise dans une deuxième plage de 25 valeurs de pixels ou o dépasse un deuxième seuil de valeur, ^ une troisième catégorie de paramètres correspondant aux nuages, dans laquelle sont classés les pixels dont la valeur extraite : 30 o est comprise dans une troisième plage de valeurs de pixels ou 3037701 18 o dépasse un troisième seuil de valeur. [80] Par ailleurs, la troisième catégorie de paramètres correspondant aux nuages peut comprendre au moins deux sous-catégories permettant de distinguer les nuages sombres des nuages clairs. Les pixels sont affectés à la première ou la deuxième sous- 5 catégorie si sa valeur est inférieure ou supérieure à un seuil. [81] Par exemple, le seuil qui limite la première et la deuxième sous-catégorie peut être par exemple une valeur de luminance dont le signal d'intensité I est défini par la combinaison linéaire suivante (où les coefficients sont des pondérations universelles standardisées ITU-RBT709) : 10 I = 0.2126 Ir + 0.7152 lv + 0.0722 lb [82] Le seuil qui limite la deuxième et la troisième catégorie de paramètres peut quant à lui être une valeur de luminance dont le signal d'intensité I est défini par la combinaison linéaire suivante : I = 0.2989 Ir + 0.5870 Iv + 0.1140 lb 15 [083] Toutefois, il convient de noter que les seuils et plages qui caractérisent les différentes catégories et sous-catégories de paramètres peuvent évoluer, par exemple en fonction de l'appareil de prise d'image utilisé. [84] Lors de la génération du vecteur de paramètres, des coefficients peuvent être associés aux pixels. Ces coefficients prennent en compte notamment la catégorie 20 d'élément de ciel correspondant au pixel et la distance de pixel par rapport au repère optique. [85] Par exemple, les coefficients Cl et C2 associés aux pixels sont définis selon les équations suivantes : Cl = 1 ( ndg(pixel) * sin (y) ) /I sin (y) 25 C2 = 1 ( luminance(pixel) * distS ) /I distS où « ndg » est la valeur des pixels mise en niveau de gris et y est l'élévation pour chaque pixel et où « luminance » est la valeur des pixels mise en luminance et distS correspond à la distance du pixel au soleil sur l'image. 3037701 19 [86] A l'étape 330, les paramètres du vecteur sont comparés avec des paramètres prédéterminés enregistrés dans la mémoire 102. Dans la mémoire 102, les paramètres prédéterminés en question sont respectivement associés à des indicateurs de production d'énergie. A l'issue de cette comparaison, il est possible de situer le vecteur de 5 paramètres par rapport à l'historique constitué par les paramètres prédéterminés. Des opérations complémentaires vont pouvoir ainsi être procédées pour estimer plus finement la production d'énergie du système solaire en fonction des paramètres prédéterminés auquel a été comparé le vecteur généré. [87] A l'étape 340, une moyenne pondérée des indicateurs de production peut être 10 effectuée, la pondération étant fonction de la distance des paramètres prédéterminés associés aux indicateurs par rapport au vecteur de paramètres. Typiquement, le poids des indicateurs peut être : - augmenté pour les paramètres prédéterminés proches du vecteur, et - diminué pour les autres paramètres prédéterminés. 15 [088] A titre d'exemple purement illustratif, la moyenne pondérée est nommée Kmoy et est défini pour k vecteurs de paramètres selon l'équation : E poids(k) * K (k) Kmoy = E poids (k) La variable poids(k) peut être calculée à partir d'un noyau gaussien donné par : poids(k) = e-ciponci(k)z )hz où h est une variable relative à la largeur du noyau et dpond correspond à la distance pondérée entre les paramètres prédéterminés X et le vecteur de paramètres x : dpond (k) = px (i) * (x(k, i) - x(i))2 EL Px (i) 20 [089] Les pondérations px permettent d'accorder plus ou moins d'importance aux différents paramètres selon le type d'irradiation à estimer (globale, directe, ou diffuse). D'autres noyaux peuvent être utilisés comme par exemple une décroissance exponentielle ou un noyau uniforme (i.e. constant sur une certaine distance puis nul au-delà de cette distance). 3037701 20 [90] La moyenne pondérée Kmoy peut être effectuée : sur l'ensemble des paramètres prédéterminés, ou seulement pour une partie des paramètres prédéterminés, par exemple en ne prenant en compte que : 5 o les paramètres prédéterminés situés dans un rayon donné autour du vecteur de paramètres, ou o les paramètres prédéterminés les plus proches du vecteur à partir d'une distance donnée autour du vecteur de paramètres. [91] Par exemple, la figure 4 représente graphiquement, sous forme d'un 10 hypercube, une base de données comportant des paramètres prédéterminés associés à des indicateurs de production photovoltaïque. Ces paramètres prédéterminés sont représentés dans un repère comportant les axes Pl, P2 et P3, chaque axe correspondant à une catégorie de paramètres telle que ciel dégagé, nuages clair, nuages sombres et soleil. Le vecteur de paramètres généré est représenté selon un cube 410 et les 15 paramètres prédéterminés sont illustrés par des cercles 420, étoiles 430 et croix 440. [92] Dans cet exemple, la moyenne pondérée Kmoy est calculée à partir des paramètres prédéterminés qui sont dans un rayon réduit 450 autour du vecteur 410. [93] Par ailleurs, à l'étape 340, un écart type de l'erreur d'estimation de l'indicateur de production d'énergie peut être déterminé. Cet écart type peut permettre de définir un 20 intervalle de confiance (i.e. plage d'indicateur de production d'énergie attendue) autour de l'estimation de l'indicateur de production avec des limites supérieure et inférieure. Cet écart type peut notamment être calculé à partir de la racine d'une moyenne pondérée d'écarts quadratiques entre des indicateurs de production d'énergie des paramètres prédéterminés et l'indicateur de production d'énergie estimé, la 25 moyenne étant pondérée en fonction d'écarts de valeurs entre les paramètres du vecteur généré et les paramètres prédéterminés. [94] Selon une réalisation possible, l'équation pour calculer l'écart type u est la suivante : 3037701 21 ,\IE poids (k) * (Kmoy - K(k))2 5 = E poids (k) [95] Il convient de noter que l'intervalle de confiance est d'autant plus précis que les paramètres prédéterminés considérés pour le calcul des moyennes Kmoy sont nombreux et proches du vecteur de paramètres 410. [96] Selon une étape 350, le vecteur de paramètres, les moyennes Kmoy et/ou l'écart 5 type u sont normalisés par l'intermédiaire d'un modèle ciel clair (modèle de valeurs de rayonnement pour un ciel dégagé, prenant en compte l'évolution zénithale du soleil). [97] Ainsi à une étape 360, l'indicateur de production d'énergie du système solaire peut être estimé à partir des indicateurs de production d'énergie des paramètres prédéterminés comparés au vecteur. 10 [98] A la figure 5, un organigramme composé d'un exemple de succession d'étapes est représenté pour la construction d'une base de données de paramètres prédéterminés. [99] Aux étapes 500 et 502, des mesures d'irradiations à un instant donné INST T 15 et des images du ciel concomitantes prises au niveau de la zone terrestre sont reçues. [0100] A l'étape 510, un modèle ciel clair est appliqué aux mesures reçues pour normaliser les mesures et tenir compte de l'évolution zénithale du soleil. [0101] Selon une étape 512, les valeurs de pixels sont extraites et localisées dans les images reçues. Des paramètres sont alors calculés en fonction des valeurs extraites et 20 des distances d'éloignement déterminées selon des catégories d'éléments de ciel telles que le ciel clair, nuage clair, nuage sombre, disque solaire, région circumsolaire. [0102] A l'étape 520, les mesures d'irradiations normalisées sont associées aux paramètres prédéterminés. La mesure d'irradiation associée au paramètre prédéterminé est un indicateur de la production d'énergie du système solaire. Les paramètres 25 prédéterminés ainsi générés sont ensuite stockés dans la mémoire 102. [0103] En variante ou en complément, l'indicateur de production d'énergie peut être une production électrique effective mesurée au niveau du système solaire ou autre. 3037701 22 [0104] On se réfère maintenant à la figure 6 qui représente un organigramme composé d'un exemple de succession d'étapes permettant de générer une image de ciel à partir d'images de ciel prises au sol depuis de la zone terrestre. 5 [0105] A l'étape 600, les données DAT_IMG d'images prises à un premier instant INST1 et un deuxième instant INST2 sont reçues par l'interface 106 du dispositif. Les données DAT_IMG comportent des informations sur les valeurs et localisations de pixels des images prises par la caméra 109 aux instants INST1 et INST2. Il est considéré que l'instant INST1 est antérieur à l'instant INST2. 10 [0106] En fonction des données DAT IMG, n valeurs de pixels aux instants INST1, INST2 peuvent être obtenues, n étant un entier naturel (supérieur ou égal à zéro) : n.VALINsTi, n.VALINsT2. Par exemple, n peut correspondre au nombre de pixels des images prises aux instants INST1 et INST2, une valeur VAL brute de luminance étant associée à chacun des n pixels de l'instant INST1 et de l'instant INST2. 15 [0107] Selon une étape 610, les valeurs n.VALINsTiet n.VALINST2 obtenues sont localisées. Cette localisation peut être effectuée par identification du positionnement de chacune des n valeurs mesurées aux instants INST1 et INST2. [0108] Selon une étape 620, les localisations déterminées pour les valeurs des instants INST1 et INST2 sont comparées entre elles (via une fonction nommée COMP LOC 20 par exemple). Les valeurs localisées de l'instant INST1 sont comparées aux valeurs localisées de l'instant INST2 en vue de déterminer les déplacements DEP des n valeurs entre l'instant INST1 et INST2. [0109] Selon une réalisation possible, n déplacements DEP des n valeurs localisées sont déterminées par identification des déplacements suivis par des valeurs de même 25 niveau aux instants INST1 et INST2. [0110] En alternative, seulement p déplacements DEP peuvent être déterminés à partir des valeurs localisées, p étant un entier positif strictement inférieur à n. Cette alternative est à utiliser de préférence pour sélectionner certaines valeurs dont le déplacement est à déterminer, ce qui évite d'une part de calculer le déplacement de 30 toutes les valeurs de l'image et, d'autre part, des erreurs de calcul lorsqu'une valeur à 3037701 23 l'instant INST1 n'apparait pas dans les valeurs de l'instant INST2. Par exemple, il est possible de fixer un pas r permettant de ne calculer des vecteurs que tous les r pixels, réduisant ainsi le nombre de vecteurs calculés et donc le temps et les ressources de calcul. 5 [0111] A une étape 630, les déplacements DEP des mêmes valeurs de pixels entre les instants INST1 et INST2 sont déterminés pour les n valeurs de pixels. [0112] Bien entendu, dans des conditions de ciel clair (pas de nuage), les déplacements de valeurs susmentionnés peuvent ne pas être détectés (ils sont nuls) car les valeurs seront identiques aux premier et deuxième instants. 10 [0113] Selon une étape 640, une image pour un instant INST3 à venir (i.e. instant postérieur au deuxième instant INST2) est générée en appliquant les n déplacements n.DEP aux n pixels de l'image prise à l'instant INST2. [0114] L'horizon de prévision est défini selon le temps qui sépare l'instant INST2 (auquel la prévision pour l'instant 1NST3 est effectuée) de l'instant INST3. L'horizon 15 de prévision peut être de plusieurs minutes (5, 10, 15 voire 60 minutes par exemple). [0115] Aux figures 7, des exemples d'images prises à partir du sol de la zone de terrestre dans laquelle est installé le système solaire et d'une image générée à partir de ces images. En l'occurrence, les figures 7A et 7B sont des prises de vues prises 20 respectivement à l'instant INST1 et à l'instant INST2. Entre l'instant INST1 et l'instant INST2 se sont par exemple écoulés 10 secondes et le masque nuageux s'est déplacé dans le ciel entre ces deux instants. [0116] Les valeurs de l'instant INST1 sont comparées aux valeurs de l'instant INST2. Des déplacements des valeurs entre l'instant INST1 et l'instant INST2 peuvent être 25 déterminés selon des vecteurs. Les vecteurs correspondent au déplacement du masque nuageux sur l'image, c'est-à-dire à la translation suivie par des mêmes valeurs de pixel entre les premier et deuxième instants. [0117] A la figure 7C, les vecteurs déterminés entre les instants INST1 et INST2 ont été appliqués à partir des pixels de l'image prise à l'instant INST2 pour générer une 3037701 24 image à venir à l'instant INST3, par exemple à un horizon de prévision de 5 minutes. Cette image générée est une prévision des conditions de ciel à venir. L'application du procédé d'estimation à l'image ainsi générée permet de prévoir une production d'énergie électrique à venir, dans 5 minutes en l'occurrence. 5 [0118] On se réfère maintenant à la figure 8 sur laquelle sont représentés des exemples de courbes de rayonnement solaire global normalisé (W/m2) évoluant dans le temps (échelle horaire). Parmi ces courbes sont illustrées : - une courbe 810 de rayonnement solaire prévue à un horizon de 5 minutes selon 10 le procédé d'estimation proposé (i.e. images utilisées pour l'estimation sont ici des images prévues à un horizon de 5 minutes, lesquelles images prévues sont générées à partir d'images du ciel prises par la caméra 109, selon les étapes susmentionnées) ; - des courbes 812 et 814 correspondant aux bornes supérieure et inférieure de 15 l'intervalle de confiance qui peut être calculé à partir de l'écart type u suivant : jE poids (k) * (Kmoy - K(k))2 5 = E poids (k) - une courbe 820 de mesures effectives et normalisées du rayonnement au niveau du système solaire (mesures prises par un instrument de type capteur pyranométrique par exemple, placé sur les panneaux solaires 110). [0119] Dans cet exemple, l'indicateur de production d'énergie est donc un 20 rayonnement solaire. Il convient de constater que les mesures effectives 820 de rayonnement sont très proches des valeurs de rayonnement prévues 810 selon le procédé d'estimation proposé. [0120] Cette qualité d'estimation et de prévision autorisée par le procédé peut permettre de connaître et anticiper précisément l'évolution d'indicateurs de production 25 d'énergie du système solaire. Ainsi, la capacité de production du système solaire peut être déterminée avec justesse, même à des horizons proches (5, 10, 15 voire 60 minutes). 3037701 25 [0121] Par ailleurs, l'intervalle de confiance (espace compris entre les courbes 812 et 814) permet de prévoir une plage de valeurs autour de l'estimation de l'indicateur de production d'énergie, afin notamment de prendre en compte d'éventuelles légères variations liées par exemple à des changements brutaux des conditions climatiques. 5 [0122] L'invention a été décrite en référence à des modes de réalisations particuliers qui ne sont pas limitatifs. Bien entendu, la présente invention ne se limite pas à la forme de réalisation décrite à titre d'exemple et elle s'étend à d'autres variantes. Par exemple, les images reçues peuvent être prises par plusieurs caméras distinctes situées 10 à différentes positions dans la zone terrestre. L'utilisation de multiples caméras (qui peuvent être reliées entre elles par un réseau de communication maillé par exemple) peut permettre de prévoir plus précisément les déplacements de masques nuageux et ainsi affiner les estimations de production d'énergie à venir. Une telle réalisation permet en outre d'établir une carte de rayonnement estimé autour du système solaire 15 ou de l'ensemble de systèmes solaires considéré. REVENDICATIONS1. Procédé d'estimation d'indicateur de production d'énergie d'un système solaire, ledit procédé comprenant les étapes de : a) réception (300) d'images du ciel prises au sol depuis une zone terrestre comportant le système solaire (110), chaque image reçue comprenant au moins un repère optique (210, 212) ; b) extraction des valeurs de pixels et localisation de pixels d'une image issue des images reçues ; c) détermination (310) des distances d'éloignement entre les pixels localisés et l'au moins un repère optique de l'image issue des images reçues ; d) génération (320) d'un vecteur de paramètres à partir d'une classification des pixels selon les valeurs extraites et les distances d'éloignement déterminées ; e) comparaison (330) des paramètres du vecteur généré avec des paramètres prédéterminés, les paramètres prédéterminés étant respectivement associés à des indicateurs de production d'énergie ; f) estimation (360) de l'indicateur de production d'énergie du système solaire à partir des indicateurs de production d'énergie des paramètres prédéterminés comparés au vecteur généré. 2. Procédé selon la revendication 1, dans lequel à l'étape f), l'estimation de l'indicateur de production d'énergie du système solaire est en outre issue d'une pondération des indicateurs de production d'énergie associés aux paramètres prédéterminées, la pondération étant fonction d'écarts de valeurs entre les paramètres du vecteur et les paramètres prédéterminés. 3. Procédé selon la revendication 1 ou 2, comprenant en outre une étape de détermination d'une plage d'indicateur de production d'énergie qui est fonction d'un écart type calculé à partir de la racine d'une moyenne pondérée d'écarts 3037701 27 quadratiques entre des indicateurs de production d'énergie des paramètres prédéterminés et l'indicateur de production d'énergie estimé, la moyenne étant pondérée en fonction d'écarts de valeurs entre les paramètres du vecteur généré et les paramètres prédéterminés. 5 4. Procédé selon l'une des revendications précédentes, dans lequel l'image comprend un premier repère optique (210) correspondant au centre optique d'un appareil de prise d'image utilisé pour capturer les images reçues. 10 5. Procédé selon l'une des revendications précédentes, dans lequel l'image comprend un deuxième repère optique (212) correspondant à la localisation du soleil sur les images reçues. 6. Procédé selon l'une des revendications précédentes, dans lequel à l'étape d), le 15 vecteur de paramètres dans lequel sont classifiés les pixels comprend au moins : ^ une première catégorie de paramètres correspondant à un ciel clair, ^ une deuxième catégorie de paramètres correspondant au 20 soleil, ^ une troisième catégorie de paramètres correspondant aux nuages. 7. Procédé selon la revendication 6, dans lequel la troisième catégorie de 25 paramètres comprend : ^ une première sous-catégorie de paramètres correspondant aux nuages clairs, dans laquelle sont classifiés des pixels de valeur supérieure ou égale à un seuil, ^ une deuxième sous-catégorie de paramètres correspondant 30 aux nuages sombres, dans laquelle sont classifiés des pixels de valeur inférieure au seuil. 3037701 28 8. Procédé selon la revendication 6 ou 7, comprenant en outre une étape de détermination d'une proportion de couverture nuageuse sur l'image issue des images reçues en fonction des pixels classifiés dans la troisième catégorie de 5 paramètres. 9. Procédé selon l'une des revendications précédentes, dans lequel à l'étape b), l'image issue des images reçues est une image générée selon les étapes suivantes : 10 - récupération (600) parmi les images reçues d'images prises à au moins un premier instant et un deuxième instant, le premier instant étant antérieur au deuxième instant ; - localisation (610) des valeurs de pixels sur les images prises respectivement au premier et au deuxième instant ; 15 - comparaison (620) des localisations des valeurs de pixels au premier instant et au deuxième instant ; - détermination (630) des déplacements, entre les premier et deuxième instants, d'au moins une partie des valeurs de pixels localisées ; - génération (640) de l'image du ciel pour un troisième instant en 20 fonction des déplacements déterminés, le troisième instant étant postérieur au deuxième instant. 10. Procédé selon l'une des revendications précédentes, dans lequel les images reçues sont des images prises au sol au niveau de la zone terrestre par un 25 appareil de prise d'image (109) muni d'un objectif ayant un angle de champ supérieur à 170 degrés. 11. Procédé selon l'une des revendications précédentes, dans lequel une donnée de : 30 ^ rayonnement solaire effectif au sol dans la zone terrestre, ou ^ production d'énergie effective du système solaire, 3037701 29 est reçue, la donnée reçue étant utilisée pour étalonner ou corriger les indicateurs de production d'énergie associés aux paramètres prédéterminés. 12. Produit programme informatique comportant des instructions pour la mise en 5 oeuvre du procédé selon l'une des revendications 1 à 11, lorsque ce programme est exécuté par une unité de traitement de données telle qu'un processeur. 13. Dispositif (100) d'estimation d'indicateur de production d'énergie d'un système solaire (110), ledit dispositif d'estimation comprenant au moins : 10 - une interface d'entrée (106) pour la réception d'images du ciel prises au sol depuis une zone terrestre comportant le système solaire ; - une mémoire (102) configurée pour stocker : ^ les images du ciel reçues, et ^ des paramètres prédéterminés respectivement associés à 15 des indicateurs de production d'énergie ; - une unité de traitement de données (104) connectée à l'interface d'entrée et à la mémoire ; ladite unité de traitement de données étant configurée pour mettre en oeuvre un procédé selon l'une des revendications 1 à 11. 20 14. Dispositif selon la revendication 13, dans lequel l'interface d'entrée est connectée à un appareil de prise d'image (109) orienté vers le ciel au niveau de la zone terrestre, l'appareil de prise d'image étant muni d'un objectif ayant un angle de champ supérieur à 170 degrés. 25 15. Dispositif selon la revendication 13 ou 14, dans lequel l'interface d'entrée est en outre configurée pour recevoir des données de : ^ rayonnement solaire effectif au sol dans la zone terrestre, ou ^ production d'énergie effective du système solaire.