FR 3071057 A1 20190315 FR 1758549 A 20170914 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 10 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 15 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 kE, celui d'un suiveur solaire étant proche de 10 kE. De plus, pour fonctionner de façon optimale, cet instrument doit 20 ê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 25 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 30 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 35 rayonnement solaire). Par exemple, une bande rotative vient occulter la photodiode périodiquement afin de mesurer le DHI. 2 - 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 »). 5 Le GHI est exprimé en W.rn-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.rn-2. Dans certains cas, on cherche en outre à connaître le trouble 20 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 25 (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 : bL3 in( 30 TA=11,1 'cc) +1 avec b=0,664+0,163.exp (80001) où h est l'altitude du site considéré, Io est l'éclairement solaire extraterrestre (W.m-2), lu est le DNI par ciel clair (W.rn-2) et m est la masse d'air optique relative. 3 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 5 du DNI. En pratique, la zone circumsolaire est donc définie par la condition 0,26° 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 30 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. 4 - 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 5 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(ISAA-PAAI) - Obtenir une image à haute dynamique (HDR) du ciel. - Calculer l'indice d'éclairement global horizontal (GGFii), l'indice 10 d'éclairement normal direct (6DNi) et l'indice d'éclairement diffus horizontal (0DHI), 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 (A1), à une deuxième région d'intérêt (A2) et à une troisième région d'intérêt (A3), respectivement, par les formules ci-après, où Qp, Lp et PZA sont l'angle 15 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 : EpEA1 LpQp cos(PZA) aGHT - EpEA1 Qp LpEA2 Q EpEA2 1-pç2p aDNI I)- v EpEA3 1-pç2p apHi- v LpEA3 Qp - Convertir l'indice d'éclairement global horizontal (oGiii), l'indice d'éclairement normal direct (6DNi) et l'indice d'éclairement diffus horizontal (GDHI) en éclairement global horizontal (GHI), éclairement normal direct 20 (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 25 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 5 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 5 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 (A1) est l'ensemble des pixels de l'image vérifiant PZA 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 35 processeur, et ce procédé comprend les étapes suivantes : 6 (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 5 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. 10 (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(ISAA-PAA1) (e) On obtient une image à haute dynamique (HDR) du ciel. 15 (f) On calcule l'indice d'éclairement global horizontal (0GHI), l'indice d'éclairement normal direct (6DNi) et l'indice d'éclairement diffus horizontal (6DHI), où Qp, Lp 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 : EpEAiLpQp cos(PZA) aGHT EpEA1Qp E pEA2 1-pç2p aDNI- v LpEA2Qp E pEA3 1-pç2p apHi- v LpEA3Qp 20 (g) On effectue une conversion entre l'indice d'éclairement global horizontal (aGHI), l'indice d'éclairement normal direct (aDNI) et l'indice d'éclairement diffus horizontal (GDHI), et l'éclairement global horizontal (GHI), l'éclairement normal direct (DNI) et l'éclairement diffus horizontal (DHI), respectivement. 25 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 7 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 5 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 : 10 - 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 ; 15 - 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-à- 20 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 Metal Oxide Semiconductor ») ou un capteur CCD (« Charge Coupled Device »). Le capteur 12 comporte sur sa surface un 25 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 30 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. 8 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. 5 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. 10 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. 15 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 20 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). 25 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. 30 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, 35 par exemple à l'aide d'un damier et d'un algorithme. 9 - Calculer, de façon connue, l'angle solide Qp 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 Qp et montre la correspondance entre 5 cet angle solide Qp 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 10 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 15 (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(ISAA-PAA1) 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). 20 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 Lmin de la zone la plus sombre à la luminance Lmax de la zone la plus lumineuse (Soleil). La prise en compte de la totalité de la luminance par une caméra implique 25 que le capteur de la caméra ait une plage dynamique DRcam au moins égale à la plage dynamique du ciel DRs4 avec : DRcam=20 1og10 (1) et DRsiw=20 1og10 L,', où L1 est la luminance minimale mesurable par le capteur de la caméra et L2 est la luminance maximale mesurable par le capteur de la 30 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 35 ces images en une image HDR, tout en préservant la linéarité du capteur, 10 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 5 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 trnin et maximal tmax de telle sorte que la plage dynamique 10 pouvant être théoriquement balayée par l'image HDR vérifie la relation : tmx DRH =DRL+20 log10 Lm >DRsain iw où DRH est la plage dynamique de l'image HDR et DRS 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 15 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 DRs4. Selon l'invention, on calcule ensuite l'indice d'éclairement global horizontal (oGiii), l'indice d'éclairement normal direct (aDNI) et l'indice 20 d'éclairement diffus horizontal (6DFil), 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 (A1), à une deuxième région d'intérêt (A2) et à une troisième région d'intérêt (A3), respectivement, ces sommes étant pondérées par l'angle solide de la région considérée. 25 Ainsi, ces indices d'éclairement sont obtenus comme suit : EpEA1 LpQp cos(PZA) aGHT EpEA1 Qp E pEA21-pQp aDNI- v L pEA2 Qp E pEA3 1-pç2p apHi- v L pEA3 Qp OÙ Qp est l'angle solide sous-tendu par un pixel p (voir figure 2), Lp est la luminance mesurée par ce pixel et PZA est l'angle entre ce pixel et le zénith (voir figure 3). 11 Avantageusement, la première région d'intérêt (Ai) est l'ensemble des pixels (125) de l'image vérifiant PZA REVENDICATIONS1. Système de mesure du rayonnement solaire caractérisé en ce qu'il comporte une caméra (10) munie d'un objectif (20) hémisphérique et 5 comprenant un capteur (12) apte à capter la lumière pour générer une image, et un processeur (30) apte à : - Effectuer une calibration géométrique de la caméra (10) afin d'obtenir la correspondance entre le système de coordonnées des pixels (125) de l'image et le système de coordonnées de ladite caméra (10). 10 - Calculer l'angle solide sous-tendu par chaque pixel (125, p) de ladite image, afin de pondérer la valeur de la luminance dans les calculs ultérieurs. - Effectuer une seconde calibration de la caméra (10) en comparant la position théorique du Soleil et sa position sur l'image, de façon à 15 obtenir la correspondance entre le système de coordonnées de ladite caméra (10) et les points cardinaux. - Calculer l'angle entre chaque pixel (125) et le zénith (PZA), l'angle entre chaque pixel (125) 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 20 entre chaque pixel (125) et le Soleil (SPA), par la relation : cos(SPA) = cos(SZA) - cos(PZA) +sin(SZA)-sin(PZA)-cos(ISAA-PAAI) - Obtenir une image à haute dynamique (HDR) du ciel. - Calculer l'indice d'éclairement global horizontal (oGiii), l'indice d'éclairement normal direct (6DNi) et l'indice d'éclairement diffus horizontal (6DHI), définis comme la somme des luminances mesurées 25 par les pixels (125) de ladite image appartenant à une première région d'intérêt (A1), à une deuxième région d'intérêt (A2) et à une troisième région d'intérêt (A3), respectivement, où Qp, Lp 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 : ODNI- v LpEA2Q1) EpEA3 1-pç2p ODHI- LpEA3Qp EpEA1 LpQp cos(PZA) EpEA1Qp E pEA2 1-pç2p 16 - Convertir l'indice d'éclairement global horizontal (aGFii), l'indice d'éclairement normal direct (aDNI) et l'indice d'éclairement diffus horizontal (aDFii) en éclairement global horizontal (GHI), éclairement normal direct (DNI) et éclairement diffus horizontal (DHI), respectivement. 2. Système de mesure du rayonnement solaire selon la revendication 1, caractérisé en ce que ladite image à haute dynamique 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 à haute dynamique, tout en conservant la linéarité du capteur de la caméra. 3. Système de mesure du rayonnement solaire selon la revendication 1 ou 2, caractérisé en ce que ladite première région d'intérêt (A1) est l'ensemble des pixels (125) de l'image vérifiant PZA 4. Système de mesure du rayonnement solaire selon l'une quelconque des revendications 1 à 3, caractérisé en ce qu'il comporte un caisson de protection (40) étanche à l'eau et à la poussière dans lequel est logée ladite caméra (10) et qui comporte une fenêtre transparente (42) située devant l'objectif (20) de ladite caméra (10). 5. Système de mesure du rayonnement solaire selon la revendication 4, caractérisé en ce que ledit caisson (40) est équipé d'un régulateur de température (50) apte à réguler la température à l'intérieur dudit caisson (40). 6. Système de mesure du rayonnement solaire selon l'une quelconque des revendications 1 à 5, caractérisé en ce qu'il comprend au moins un filtre (60) qui est apte à atténuer ou moduler la lumière entrant par l'objectif (20) dans ladite caméra (10). 7. Procédé de mesure du rayonnement solaire caractérisé en ce qu'il utilise une caméra (10) munie d'un objectif (20) hémisphérique et comportant un capteur (12) apte à capter la lumière pour générer une image, et un processeur (30), et en ce qu'il comporte les étapes suivantes : 17 (a) On calibre géométriquement la caméra (10) afin d'obtenir la correspondance entre le système de coordonnées des pixels (125) de l'image fournie par ladite caméra (10) et le système de coordonnées de la caméra (10). 5 (b) On calcule l'angle solide sous-tendu par chaque pixel (125, p) de ladite image, afin de pondérer la valeur de la luminance dans les calculs ultérieurs. (c) On calibre ladite camera (10) une seconde fois en comparant la position théorique du Soleil et sa position obtenue sur l'image, de façon à 10 obtenir la correspondance entre le système de coordonnées de ladite caméra (10) et les points cardinaux. (d) On calcule l'angle entre chaque pixel (125, p) et le zénith (PZA), l'angle entre chaque pixel (125, p) 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 15 entre chaque pixel (125, p) et le Soleil (SPA), grâce à la relation : cos(SPA) = cos(SZA) - cos(PZA) +sin(SZA)-sin(PZA)-cos(ISAA-PAAI) (e) On obtient une image à haute dynamique (HDR) du ciel. (f) On calcule l'indice d'éclairement global horizontal (0GHI), l'indice d'éclairement normal direct (aDNI) et l'indice d'éclairement diffus horizontal (0DHI), définis comme la somme des luminances mesurées par les pixels 20 (125) de ladite image appartenant à une première région d'intérêt (A1), à une deuxième région d'intérêt (A2) et à une troisième région d'intérêt (A3), respectivement, où Qp, Lp 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 : EpEA1 LpQp cos(PZA) aGHT = ODHI- LpEA3Qp 25 (g) On effectue une conversion entre l'indice d'éclairement global horizontal (oGiii), l'indice d'éclairement normal direct (aDNI) et l'indice d'éclairement diffus horizontal (aDFii), et l'éclairement global horizontal (GHI), l'éclairement normal direct (DNI) et l'éclairement diffus horizontal (DHI), respectivement. ODNI- v LpEA2Qp E pEA3 1-pç2p EpEA1Qp E pEA2 I-PQP 18 8. Procédé de mesure du rayonnement solaire selon la revendication précédente caractérisé en ce qu'il comprend en outre l'étape suivante : (h) On calcule le profil solaire en mesurant l'éclairement normal direct 5 (DNI) le long d'une ligne de pixels passant par le centre du Soleil sur l'image, puis on répète cette mesure pour une pluralité d'angles, de manière à obtenir une distribution radiale de la luminance autour du Soleil.