La présente invention concerne le domaine du rayonnement solaire et des propriétés de l'atmosphère.

L'invention concerne plus particulièrement un système de mesure du rayonnement solaire.

Pour certaines applications, il est nécessaire de connaître les différentes composantes du rayonnement solaire. En particulier, on souhaite connaître les composantes suivantes :

- La composante directe du rayonnement solaire.

Lorsque le rayonnement solaire qui n'interagit pas avec l'atmosphère terrestre lors de sa traversée est mesuré sous incidence normale, on parle d'éclairement normal direct ou DNI (« Direct Normal Irradiance »). Le DNI est exprimé en W.m ^"2 .

Pour mesurer le DNI, on utilise généralement un pyrhéliomètre, qui se présente sous la forme d'un tube dont les dimensions définissent l'angle d'ouverture et au fond duquel est placée une thermopile. Le pyrhéliomètre doit être aligné avec l'axe du Soleil, c'est pourquoi un suiveur solaire est nécessaire si l'on souhaite automatiser l'acquisition du DNI. Le coût d'un pyrhéliomètre est d'environ 2,5 k€, celui d'un suiveur solaire étant proche de 10 k€. De plus, pour fonctionner de façon optimale, cet instrument doit être entretenu quotidiennement. Face à ces contraintes de coût et de maintenance, des instruments moins onéreux et plus robustes ont été développés. Parmi eux, l'instrument le plus utilisé est le pyranomètre à bande rotative (voir ci-dessous). Le pyranomètre présente toutefois le désavantage de nécessiter une calibration spectrale précise de la photodiode. Cette calibration est complexe, produisant une incertitude sur le DNI plus élevée que pour un pyrhéliomètre.

- La composante diffuse du rayonnement solaire.

Lorsque le rayonnement solaire ayant interagi avec l'atmosphère terrestre lors de sa traversée est mesuré sur une surface horizontale, on parle d'éclairement diffus horizontal ou DHI (« Diffuse Horizontal Irradiance »). Le DHI est exprimé en W.m ^"2 .

Le plus souvent, pour mesurer le DHI, on utilise conjointement un pyranomètre et un suiveur solaire équipé d'un élément qui masque le Soleil (permettant ainsi la suppression de la composante directe du rayonnement solaire). Par exemple, une bande rotative vient occulter la photodiode périodiquement afin de mesurer le DHI.

- Le rayonnement solaire global, défini comme la somme des rayonnements direct et indirect reçus au sol. Si ce rayonnement est mesuré horizontalement, on parle d'éclairement global horizontal ou GHI (« Global Horizontal Irradiance »).

Le GHI est exprimé en W.m ^"2 .

Pour mesurer le GHI, on utilise généralement un pyranomètre, qui comprend une photodiode montée horizontalement.

À partir des mesures du GHI et du DHI et de la connaissance de la position du Soleil, le DNI peut être déterminé.

En effet, le GHI, le DNI et le DHI sont liés par la relation suivante :

GHI= cos(SZA) -DNI+DHI

où SZA est l'angle entre le Soleil et le zénith.

Dans certains cas, on cherche en outre à connaître la durée d'ensoleillement.

La durée d'ensoleillement d'une zone géographique est la durée pendant laquelle cette zone géographique est éclairée par le Soleil.

L'Organisation Météorologique Mondiale (OMM) définit la durée d'ensoleillement comme la durée pendant laquelle le rayonnement solaire direct dépasse le seuil de 120 W.m ^"2 .

Dans certains cas, on cherche en outre à connaître le trouble atmosphérique.

Le trouble atmosphérique est relatif aux atténuations du rayonnement solaire lors de sa traversée de l'atmosphère terrestre. Cette atténuation est principalement due aux phénomènes de diffusion par les aérosols et d'absorption par les différents composants atmosphériques (ozone, vapeur d'eau, oxygène...).

Le trouble atmosphérique donne une indication fiable de la qualité de la ressource solaire ou de l'atmosphère pour un site donné.

Par exemple, une estimation du Trouble Atmosphérique (TA) est effectuée en utilisant la relation :

TA=11,1

où h est l'altitude du site considéré, Io est l'éclairement solaire extraterrestre (W.rrf ² ), I _cc est le DNI par ciel clair (W.m ^~2 ) et m est la masse d'air optique relative.

Dans certains cas, on cherche en outre à connaître le profil solaire, c'est-à-dire la distribution de la luminance L dans la zone circumsolaire, c'est-à-dire la zone comprise entre l'extérieur du disque solaire (environ 0,26°) et le demi-angle d'ouverture classique des instruments de mesure du DNI. En pratique, la zone circumsolaire est donc définie par la condition 0,26° < SPA < β, où SPA est l'angle entre le pixel et le Soleil et β est compris entre 1° et 5°.

La luminance énergétique, notée L et exprimée en W.m ^~2 .sr ^_1 , est le flux rayonné par unité d'angle solide et de surface. La luminance énergétique L est liée à l'éclairement I par la relation :

où dQ est l'angle solide élémentaire, Θ et φ définissant la direction de propagation de la lumière.

Pour un rayonnement isotrope, on a I = π-L.

L'instrument le plus précis pour effectuer cette mesure est l'instrument SAM (« Sun Auréole Measurement »). Cet instrument présente l'inconvénient d'un coût très élevé (environ 80 k€).

Chacun des instruments mentionnés ci-avant a un coût minimum de plusieurs milliers d'euros et nécessite un entretien régulier, tant pour le suiveur solaire que pour les capteurs eux-mêmes. La mesure des différentes composantes du rayonnement solaire est donc un processus coûteux et contraignant.

La présente invention vise à remédier à ces inconvénients.

L'invention propose un système permettant de mesurer plus aisément et à moindre coût les composantes du rayonnement solaire.

Ce but est atteint grâce au fait que le système comporte une caméra munie d'un objectif hémisphérique et d'un capteur apte à capter la lumière pour générer une image, et un processeur apte à :

Effectuer une calibration géométrique de la caméra afin d'obtenir la correspondance entre le système de coordonnées des pixels de l'image et le système de coordonnées de la caméra. Calculer l'angle solide sous-tendu par chaque pixel de l'image, ce qui permettra de pondérer la valeur de la luminance dans les calculs ultérieurs.

Effectuer une calibration de la caméra donnant la position du Soleil dans le ciel de façon à obtenir la correspondance entre le système de coordonnées de la caméra et les points cardinaux.

Calculer l'angle entre chaque pixel et le zénith (PZA), l'angle entre chaque pixel et l'azimut (PAA), l'angle entre le Soleil et le zénith (SZA) et l'angle entre le Soleil et l'azimut (SAA), puis l'angle entre chaque pixel et le Soleil (SPA), en utilisant la relation :

cos(SPA) = cos(SZA) ^■ cos(PZA) +sin(SZA)-sin(PZA)-cos(|SAA-PAA|) Obtenir une image à haute dynamique (HDR) du ciel. Calculer l'indice d'éclairement global horizontal (OGHI), l'indice d'éclairement normal direct (O _D NI) et l'indice d'éclairement diffus horizontal (ODHI), définis comme la somme des luminances mesurées par les pixels de l'image appartenant à une première région d'intérêt (Ai), à une deuxième région d'intérêt (A ₂ ) et à une troisième région d'intérêt (A ₃ ), respectivement, par les formules ci-après, où Ω _ρ , L _p et PZA sont l'angle solide sous-tendu par un pixel p, la luminance mesurée par ce pixel et l'angle entre ce pixel et le zénith, respectivement :

_∑ΡΕ _Λι LpQp COs(PZA)

∑PEAJ ¾

∑ρεΛ3 '-ρΩρ

° ^" DHI ⁼

∑ρεΛ3 ¾

- Convertir l'indice d'éclairement global horizontal (OGHI), l'indice d'éclairement normal direct (ODNI) et l'indice d'éclairement diffus horizontal (ODHI) en éclairement global horizontal (GHI), éclairement normal direct (DNI) et éclairement diffus horizontal (DHI), respectivement.

Grâce à ces dispositions, on obtient une estimation des différentes composantes du rayonnement solaire (DNI, GHI, DHI) avec une seule caméra, qui remplace avantageusement plusieurs instruments de mesure. Le système selon l'invention est donc moins onéreux en achat et en maintenance et moins complexe. De plus, le système selon l'invention n'utilise pas de masque qui occulte le Soleil. Il y a donc suppression des coûts de fonctionnement et de maintenance d'un masque et des coûts associés à un système mécanique assurant le suivi de la trajectoire du Soleil. En outre, étant donné que le Soleil reste visible durant les mesures, les algorithmes d'estimation des différentes composantes sont plus performants.

Avantageusement, l'image à haute dynamique (HDR) est obtenue par l'acquisition de plusieurs images à basse dynamique prises avec différents temps de pose puis la combinaison de ces images en une image HDR, tout en conservant la linéarité du capteur de la caméra.

Avantageusement, la première région d'intérêt (Ai) est l'ensemble des pixels de l'image vérifiant PZA < 90°, la deuxième région d'intérêt (A ₂ ) est l'ensemble des pixels de l'image vérifiant SPA < β, où β est l'angle défini comme la limite entre la zone circumsolaire et le reste du ciel, et la troisième région d'intérêt (Λ ₃ ) est l'ensemble des pixels de l'image vérifiant {SPA < β} Π {PZA < 90°}.

Ainsi, on obtient des valeurs fiables pour les diverses composantes du rayonnement solaire.

Avantageusement, le système de mesure comporte un caisson de protection étanche à l'eau et à la poussière dans lequel est logée la caméra, et qui comporte une fenêtre transparente située devant l'objectif de la caméra.

Ainsi, la caméra est protégée des intempéries (pluie et vent) et des corps étrangers.

Avantageusement, le caisson de protection est équipé d'un régulateur de température apte à réguler la température à l'intérieur du caisson.

Ainsi, le système est apte à fonctionner en environnement hostile, par exemple sur le site d'implantation d'une centrale solaire, où les températures sont parfois très élevées.

Avantageusement, le système comprend au moins un filtre qui est apte à atténuer ou moduler la lumière entrant par l'objectif dans la caméra.

L'invention concerne également un procédé de mesure du rayonnement solaire. Selon l'invention, ce procédé de mesure du rayonnement solaire utilise une caméra munie d'un objectif hémisphérique et comportant un capteur apte à capter la lumière pour générer une image, et un processeur, et ce procédé comprend les étapes suivantes :

(a) On calibre géométriquement la caméra afin d'obtenir la correspondance entre le système de coordonnées des pixels de l'image fournie par la caméra et le système de coordonnées de la caméra,

(b) On calcule l'angle solide sous-tendu par chaque pixel de l'image, afin de pondérer la valeur de la luminance dans les calculs ultérieurs.

(c) On calibre une seconde fois la caméra en comparant la position théorique du Soleil avec sa position obtenue sur l'image. De cette manière, la correspondance entre le système de coordonnées de la caméra et les points cardinaux est obtenue.

(d) On calcule l'angle entre chaque pixel et le zénith (PZA), l'angle entre chaque pixel et l'azimut (PAA), l'angle entre le Soleil et le zénith (SZA) et l'angle entre le Soleil et l'azimut (SAA), puis l'angle entre chaque pixel et le Soleil (SPA), grâce à la relation :

cos(SPA) = cos(SZA) ^■ cos(PZA) +sin(SZA)-sin(PZA)-cos(|SAA-PAA|)

(e) On obtient une image à haute dynamique (HDR) du ciel.

(f) On calcule l'indice d'éclairement global horizontal (OGHI), l'indice d'éclairement normal direct (O _D NI) et l'indice d'éclairement diffus horizontal

(ODHI), OÙ Ω _ρ , L _p et PZA sont l'angle solide sous-tendu par un pixel p, la luminance mesurée par ce pixel et l'angle entre ce pixel et le zénith, respectivement :

_ Σ _ΡΕΛι LpQp COs(PZA)

° ^{GHI ~} Z vjpsA _j " ôp

(g) On effectue une conversion entre l'indice d'éclairement global horizontal (OGHI), l'indice d'éclairement normal direct (O _D NI) et l'indice d'éclairement diffus horizontal (o _D m), et l'éclairement global horizontal (GHI), l'éclairement normal direct (DNI) et l'éclairement diffus horizontal (DHI), respectivement. Avantageusement, le procédé selon l'invention comprend en outre l'étape suivante :

(h) On calcule le profil solaire en mesurant l'éclairement normal direct (DNI) le long d'une ligne de pixels passant par la centre du Soleil sur l'image, puis on répète cette mesure pour une pluralité d'angles, de manière à obtenir la distribution radiale de la luminance autour du Soleil.

Ainsi, on évite l'utilisation d'équipements supplémentaires coûteux, et on obtient une mesure du profil solaire qui est plus précise car cette mesure prend en compte l'anisotropie de la luminance autour du Soleil.

L'invention sera bien comprise et ses avantages apparaîtront mieux, à la lecture de la description détaillée qui suit, d'un mode de réalisation représenté à titre d'exemple non limitatif. La description se réfère aux dessins annexés sur lesquels :

- la figure 1 est une vue schématique d'un système de mesure selon l'invention ;

- la figure 2 est un schéma illustrant la correspondance entre un pixel du capteur de la caméra du système selon l'invention et un angle solide ;

- la figure 3 est un schéma montrant les différents angles utilisés pour définir la position du Soleil et d'un pixel de la caméra du système selon l'invention par rapport à un repère.

Le système selon l'invention comprend une caméra 10, munie d'un objectif 20 hémisphérique (également appelé objectif « fisheye »), c'est-à- dire un objectif à grand angle de champ, avec un champ de vision supérieur à 180°. La caméra 10 comprend un capteur 12 apte à capter la lumière pour générer une image, par exemple un capteur CMOS (« Complementary Métal Oxide Semiconductor ») ou un capteur CCD (« Charge Coupled Device »). Le capteur 12 comporte sur sa surface un certain nombre de pixels (125) sur lesquels aboutit la lumière pour former une image.

Avantageusement, la caméra 10 est située sur une plate-forme horizontale, et est située en un lieu où il n'y a pas d'ombrage (par exemple en hauteur et loin de bâtiments et de végétation). Ainsi, on obtient une image de la totalité ou quasi-totalité du ciel, sans obstacles cachant des parties du ciel. La caméra 10 est orientée selon un axe Z qui pointe depuis le sol vers le zénith.

Avantageusement, le système de mesure comporte un caisson de protection 40 étanche à l'eau et à la poussière dans lequel est logée la caméra 10.

Le caisson 40 comporte une fenêtre transparente 42, par exemple un dôme transparent comme représenté sur la figure 1, qui est située devant l'objectif 20 de la caméra 10 et qui protège cet objectif 20. Ainsi, la fenêtre transparente 42 est située entre l'objectif 20 et le Soleil.

Ainsi, la caméra 10 est protégée de l'eau et de la poussière, tout en laissant les rayons solaires accéder à l'objectif 20.

Avantageusement, le caisson de protection 40 est équipé d'un régulateur de température 50 apte à réguler la température à l'intérieur du caisson 40.

Ainsi, on diminue ou en empêche la condensation à l'intérieur du caisson 40. En conséquence les images prises par la caméra 10 ne sont pas polluées par l'humidité.

Avantageusement, le système comprend un ou plusieurs filtres 60 qui sont aptes à atténuer ou moduler la lumière entrant dans la caméra 10.

Le ou les filtres 60 sont positionnés entre le Soleil et la caméra 10.

Ainsi, le ou les filtres 60 sont situés soit entre la caméra 10 et l'objectif 20, soit dans l'objectif 20, soit devant l'objectif 20, soit une combinaison de ces positions.

Par exemple, la face interne de la fenêtre transparente 42 est recouverte d'un filtre 60, comme représenté sur la figure 1 dans le cas où la fenêtre 42 est un dôme.

Avantageusement, la fenêtre transparente 42 est teintée afin d'atténuer les réflexions parasites (qui pourraient, de façon erronée, être prises pour des nuages par les algorithmes développés).

Alternativement, la fenêtre transparente 42 n'est pas teintée. Dans ce cas, le processeur 30 (voir ci-dessous) comprend un algorithme qui supprime ces réflexions parasites.

Le système selon l'invention comprend un processeur 30 relié à la caméra 10.

Ce processeur 30 est apte à effectuer les opérations suivantes : - Effectuer une première calibration de la caméra 10 afin d'obtenir la correspondance entre le système de coordonnées des pixels 125 de l'image produite par le capteur 12 et le système de coordonnées de la caméra 10. Cette calibration géométrique est effectuée, de façon connue, par exemple à l'aide d'un damier et d'un algorithme.

- Calculer, de façon connue, l'angle solide Ω _ρ sous-tendu par chaque pixel (125, p) du capteur 12 de la caméra 10 enregistrant les images, afin de pondérer la valeur de la luminance L dans les calculs ultérieurs.

La figure 2 illustre un angle solide Ω _ρ et montre la correspondance entre cet angle solide Ω _ρ et un pixel (125, p) du capteur 12.

- Effectuer une seconde calibration de la caméra 10 en comparant la position théorique du Soleil avec sa position obtenue sur l'image. De cette manière, la correspondance entre le système de coordonnées de la caméra 10 et les points cardinaux (repère X, Y, Z) est obtenue. La position théorique du Soleil est obtenue de façon connue, par exemple à l'aide d'un algorithme tel que SPA (« Solar Position Algorithm ») ou SG (« Solar Geometry »).

- Calculer l'angle entre chaque pixel (125, p) et le zénith (PZA), l'angle entre chaque pixel et l'azimuth (PAA), l'angle entre le Soleil et le zénith (SZA) et l'angle entre le Soleil et l'azimut (SAA), puis l'angle entre chaque pixel (125, p) et le Soleil (SPA), par la relation :

cos(SPA) = cos(SZA) ^■ cos(PZA) +sin(SZA)-sin(PZA)-cos(|SAA-PAA|) La figure 3 montre la relation entre ces différents angles, le Soleil S, et le zénith Z, ainsi que la direction de référence nord (N, points cardinaux).

Une image à haute dynamique (HDR) est définie comme une image sur laquelle est restituée fidèlement (c'est-à-dire sans sous-exposition ou surexposition) toute la luminance du ciel, de la luminance L _min de la zone la plus sombre à la luminance L _max de la zone la plus lumineuse (Soleil). La prise en compte de la totalité de la luminance par une caméra implique que le capteur de la caméra ait une plage dynamique DR _ca m au moins égale à la plage dynamique du ciel DR _sky avec :

DR _cam =20 log ₁₀ (¾ et DR _sky =20 log ₁₀ (¾ où Li est la luminance minimale mesurable par le capteur de la caméra et L ₂ est la luminance maximale mesurable par le capteur de la caméra. Si cette condition n'est pas vérifiée, c'est-à-dire si le capteur 12 de la caméra 10 est apte à prendre uniquement des images à basse dynamique (LDR, pour « Low Dynamic Range »), l'image à haute dynamique est obtenue par l'acquisition de plusieurs images LDR puis la combinaison de ces images en une image H DR, tout en préservant la linéarité du capteur, de façon connue. On décrit ci-dessous un exemple de méthode pour l'obtention d'une image HDR.

Soit plusieurs images HDR prises avec différents temps de pose, le temps total d'acquisition de la séquence d'images LDR doit être suffisamment bref pour garantir que le mouvement des nuages entre chaque image LDR est négligeable, autrement dit que la luminance de la scène est constante au cours de l'acquisition. De plus, il est important de bien paramétrer le nombre souhaité d'images LDR ainsi que les temps de pose minimal t _min et maximal t _max de telle sorte que la plage dynamique pouvant être théoriquement balayée par l'image HDR vérifie la relation :

où DR _H est la plage dynamique de l'image HDR et DR _L la plage dynamique de l'image LDR. La valeur estimée de chaque pixel de l'image HDR générée est obtenue à partir d'une moyenne pondérée de la valeur des pixels de chaque image LDR (en considérant uniquement les pixels se trouvant dans la plage de linéarité du capteur).

Alternativement, on obtient l'image HDR en utilisant un capteur 12 avec une plage dynamique supérieure à la plage dynamique du ciel DR _sky .

Selon l'invention, on calcule ensuite l'indice d'éclairement global horizontal (o _G m), l'indice d'éclairement normal direct (O _D NI) et l'indice d'éclairement diffus horizontal (ODHI), définis comme la somme des luminances mesurées par les pixels (125) de l'image appartenant à une première région d'intérêt (Ai), à une deuxième région d'intérêt (A ₂ ) et à une troisième région d'intérêt (A ₃ ), respectivement, ces sommes étant pondérées par l'angle solide de la région considérée.

Ainsi, ces indices d'éclairement sont obtenus comme suit :

° ^" GHI _∑ρεΛ ₃ '-ρΩρ

<¾HI- - η ^~~

Ζ.ρεΛ ₃ "ρ

où Ω _ρ est l'angle solide sous-tendu par un pixel p (voir figure 2), L _p est la luminance mesurée par ce pixel et PZA est l'angle entre ce pixel et le zénith (voir figure 3).

Avantageusement, la première région d'intérêt (Ai) est l'ensemble des pixels (125) de l'image vérifiant PZA < 90°, la deuxième région d'intérêt (A ₂ ) est l'ensemble des pixels (125) de l'image vérifiant SPA < β, où β est l'angle défini comme la limite entre la zone circumsolaire (telle que définie ci-avant) et le reste du ciel, et la troisième région d'intérêt (A ₃ ) est l'ensemble des pixels (125) de l'image vérifiant {SPA < β } n {PZA < 90°}.

On obtient ainsi une évaluation plus précise des indices d'éclairement.

Ensuite, l'indice d'éclairement global horizontal (OGHI), l'indice d'éclairement normal direct (ODNI) et l'indice d'éclairement diffus horizontal (ODHI) sont convertis en éclairement global horizontal (GHI), éclairement normal direct (DNI) et éclairement diffus horizontal (DHI), respectivement, par des fonctions de conversion.

Ces fonctions de conversion consistent en des corrélations entre les indices d'éclairement (OGHI, O _D NI et o _D w) et les différentes composantes du rayonnement solaire (GHI, DNI et DHI).

Une régression linéaire simple, voire multiple lorsque plusieurs indices d'éclairement et, éventuellement, leurs historiques sont pris en compte, doit permettre d'approximer chaque fonction de conversion.

Alternativement, il est possible d'approximer ces fonctions grâce aux outils de l'intelligence artificielle, par exemple un réseau de neurones artificiels ou un système neuro-flou (un tel système combine la structure connexion niste d'un réseau de neurones et l'écriture de règles floues ; la logique floue est une extension de la logique booléenne). Une phase d'apprentissage numérique permet, à partir d'une base d'exemples (l'apprentissage est ici supervisé), de déterminer la topologie de l'outil considéré (sa structure) et d'en identifier les paramètres.

L'invention concerne également un procédé de mesure du rayonnement solaire utilisant une caméra 10 munie d'un objectif 20 hémisphérique et comportant un capteur apte à capter la lumière pour générer une image, et un processeur 30, qui comprend les étapes suivantes, déjà décrites ci -avant :

(a) Calibration géométrique afin d'obtenir la correspondance entre le système de coordonnées des pixels de l'image fournie par la caméra et le système de coordonnées de la caméra.

(b) Calcul de l'angle solide sous-tendu par chaque pixel de l'image, afin de pondérer la valeur de la luminance dans les calculs ultérieurs.

(c) Seconde calibration s'appuyant sur la comparaison entre la position théorique du Soleil et sa position obtenue sur l'image, de façon à obtenir la correspondance entre le système de coordonnées de la caméra et les points cardinaux.

(d) Calcul de l'angle entre chaque pixel (125, p) et le zénith (PZA), de l'angle entre chaque pixel et l'azimut (PAA), de l'angle entre le Soleil et le zénith (SZA) et de l'angle entre le Soleil et l'azimut (SAA), puis de l'angle entre chaque pixel (125, p) et le Soleil (SPA), par la relation :

cos(SPA) = cos(SZA) ^■ cos(PZA) +sin(SZA)-sin(PZA)-cos(|SAA-PAA|)

(e) Obtention d'une image H DR du ciel.

(f) Calcul de l'indice d'éclairement global horizontal (o _G m), de l'indice d'éclairement normal direct (O _D NI) et de l'indice d'éclairement diffus horizontal (o _D w), où Ω _ρ , L _p et PZA sont l'angle solide sous-tendu par un pixel p, la luminance mesurée par ce pixel et l'angle entre ce pixel et le zénith, respectivement :

∑ _peAl LpQp COS(PZA)

° ^" GHI - ^'

∑PEAJ Ωρ

∑ρεΛ3 '-ρΩρ

° ^" DHI ⁼

∑ρεΛ3 ¾

(g) Conversion entre o _GH i, CF _D NI et o _D m et GHI, DNI et DHI.

Avantageusement, le procédé selon l'invention comprend en outre l'étape suivante :

(i) Calcul de la durée d'ensoleillement à partir du rayonnement solaire direct.

Le procédé selon l'invention comprend alors la fourniture d'un compteur temporel qui mesure le temps lorsque la durée d'ensoleillement est supérieure au seuil précisé ci-avant, et qui s'arrête lorsque la durée d'ensoleillement est inférieure à ce seuil. On connaît alors la durée d'ensoleillement sur une période choisie, par exemple une journée, un mois, ou une année.

Avantageusement, le procédé selon l'invention comprend en outre l'étape suivante :

G) Calcul du trouble atmosphérique à partir de la détermination des instants « ciel clair » et du DNI.

Le trouble atmosphérique est relatif aux atténuations du rayonnement solaire lors de sa traversée de l'atmosphère terrestre. Cette atténuation est principalement due aux phénomènes de diffusion par les aérosols et d'absorption par les différents composants atmosphériques (ozone, vapeur d'eau, oxygène...).

Le trouble atmosphérique donne une bonne indication de la qualité de la ressource solaire ou de l'atmosphère pour un site donné.

À partir des images acquises par la caméra 10, le système selon l'invention est capable de déterminer l'évolution du trouble atmosphérique sur le site étudié.

Pour cela, le système détecte, à l'aide d'un algorithme, les instants où le Soleil n'est pas occulté par un nuage. On parle alors de situation « ciel clair ».

Un algorithme de détection des instants « ciel clair » procède donc à l'acquisition d'une image ou de plusieurs images successives prises avec différents temps de pose, dans un laps de temps relativement court (par exemple, quelques millisecondes). Le système analyse ensuite la distribution de la luminance dans la zone circumsolaire et effectue si nécessaire une détection de contour du Soleil, de façon connue.

À partir de ces paramètres obtenus, le système détermine grâce à cet algorithme si l'instant correspond ou non à une situation dite « ciel clair ».

En option, un filtrage temporel, s'appuyant sur plusieurs images successives, peut être utilisé pour s'assurer de la bonne détection de cette situation.

Chaque fois qu'un instant « ciel clair » est déterminé par l'algorithme, une estimation du Trouble Atmosphérique (TA) est réalisée.

Par exemple, cette estimation est effectuée en utilisant la relation suivante : ^|η (^) / -h \

TA=ll,l -^ +l avec b=0,664+0,163-exp

où h est l'altitude du site considéré, I ₀ est l'éclairement solaire extraterrestre (W.rrf ² ), I _cc est le DNI par ciel clair (W.m ^~2 ) et m est la masse d'air optique relative.

Avantageusement, le procédé selon l'invention comprend en outre l'étape suivante :

(h) Calcul du profil solaire à partir de l'analyse de la distribution radiale de la luminance.

Selon l'invention, pour obtenir le profil solaire, on mesure le DNI le long d'une ligne de pixels passant par le centre du Soleil sur l'image. En répétant cette mesure pour une pluralité d'angles, on obtient une distribution radiale de la luminance L. On évite ainsi avantageusement l'utilisation de l'instrument SAM (voir ci-avant) dont le coût est très élevé.

De plus, la mesure du profil solaire selon l'invention est plus précise qu'une mesure avec l'instrument SAM car le mesure selon l'invention prend en compte l'anisotropie de la luminance autour du Soleil (alors que l'instrument SAM suppose, à tort, une symétrie radiale du rayonnement solaire).