SYSTEM FOR MEASURING COMPONENTS OF SOLAR RADIATION

Abstract

Disclosed is a system for measuring solar radiation, with a camera having a hemispherical objective and light sensor, and a processor to: perform a geometric calibration of the camera establishing correspondence between coordinate systems of the image pixels and the camera; calculate the solid angle occupied by each pixel; perform a second calibration by comparing the theoretical position of the sun and its position in the image, establishing correspondence between the camera coordinate system and the cardinal points; calculate the angles between: each pixel and the zenith; each pixel and the azimuth; the sun and the zenith; and the sun and the azimuth; then each pixel and the sun. The next steps are: obtain a high-dynamic-range image of the sky; calculate the global horizontal irradiance, the direct normal irradiance, and the diffuse horizontal irradiance; and convert, into global horizontal luminance, direct normal irradiance and diffuse horizontal irradiance, respectively.

Claims

1. System of measurement of the solar radiation comprising a camera provided with a hemispherical objective and comprising a sensor adapted to receive light to generate an image and a processor able to: perform a calibration geometry of the camera in order to obtain the correspondence between the system of coordinates of pixels of the image and the system of coordinates of said camera; calculate the solid angle subtended by each pixel of said image, in order to weight the value of the luminance in subsequent calculations; carry out a second calibration of the camera by comparing the theoretical position of the Sun and its position on the image, of how to obtain the correspondence between the system of coordinates of said camera and the cardinal points; calculate the angle between each pixel and the zenith PZA, the angle between each pixel and the azimuth PAA, the angle between the Sun and the zenith SZ, and the angle between the Sun and the azimuth SAA, then the angle between each pixel and the Sun SPA, by the relationship:
cos(SPA)=cos(SZA).Math.cos(PZA)+sin(SZA).Math.sin(PZA).Math.cos(|SAA−PAA|) obtain a high dynamic image of the sky. calculate the Global Horizontal Irradiance index σ.sub.GHI, the Direct Normal Irradiance index σ.sub.DNI, and the Diffuse Horizontal Irradiance index σ.sub.DHI, defined as the sum of the luminances measured by the pixels of said image belonging to a first region of interest Λ.sub.1, to a second region of interest Λ.sub.2 and to a third region of interest Λ.sub.3, respectively, where Ω.sub.p, L.sub.p and PZA are the solid angle subtended by a pixel p, the luminance measured by this pixel and the angle between this pixel and the zenith, respectively: ${\sigma }_{\mathrm{GHI}}=\frac{{.Math.}_{p\in {\Lambda }_1}{L}_p{\Omega }_p\cos \left(\mathrm{PZA}\right)}{{.Math.}_{p\in {\Lambda }_1}{\Omega }_p}$ ${\sigma }_{\mathrm{DNI}}=\frac{{.Math.}_{p\in {\Lambda }_2}{L}_p{\Omega }_p}{{.Math.}_{p\in {\Lambda }_2}{\Omega }_p}$ ${\sigma }_{\mathrm{DHI}}=\frac{{.Math.}_{p\in {\Lambda }_3}{L}_p{\Omega }_p}{{.Math.}_{p\in {\Lambda }_3}{\Omega }_p}$ convert the Global Horizontal Irradiance index σ.sub.GHI, the Direct Normal Irradiance index σ.sub.DNI, and the Diffuse Horizontal Irradiance index σ.sub.DHI to Global Horizontal Irradiance, Direct Normal Irradiance and Diffuse Horizontal Irradiance, respectively.

2. The solar radiation measurement system according to claim 1, wherein the said high dynamic range image is obtained by the acquisition of several low dynamic range images taken with different exposure times and then the combination of these images in one high dynamic range image, all while maintaining the linearity of the sensor's camera.

3. The system of measurement of the solar radiation according to claim 1, wherein said first region of interest Λ.sub.1 is the set of pixels of the image verifying PZA<90°, said second region of interest Λ.sub.2 is the set of pixels of the image verifying SPA<β, β being the angle defined as the limit between the circumsolar zone and the rest of the sky, and said third region of interest Λ.sub.3 is the set of pixels of the image verifying {SPA<β}∩{PZA<90°}.

4. The solar radiation measuring system according to claim 1, comprising a waterproof and dustproof protective box in which is housed said camera and which comprises a transparent window located in front of the objective of said camera.

5. The solar radiation measurement system according to claim 4, wherein said box is equipped with a temperature controller capable of regulating the temperature inside said box.

6. The solar radiation measurement system according to claim 1, at least one filter which is adapted to mitigate or modulate the light entering through the objective of said camera.

7. A method of measuring the solar radiation wherein the method uses a camera provided with a hemispherical objective and having a sensor adapted to receive light to generate an image and a processor, and comprises the following stages: (a) geometrically calibrating the camera in order to obtain the correspondence between the pixel coordinate system of the image supplied by said camera and the coordinate system of the camera; (b) calculating the solid angle subtended by each pixel of said image, in order to weight the value of the luminance in the subsequent calculations; (c) calibrating said camera a second time by comparing the theoretical position of the Sun and its position obtained on the image, in order to obtain the correspondence between the coordinates of said camera system and the cardinal points; (d) calculating the angle between each pixel and the zenith PZA, the angle between each pixel and the azimuth PAA, the angle between the Sun and the zenith SZA, and the angle between the Sun and the azimuth SAA, then the angle between each pixel and the Sun SPA, thanks to the relationship:
cos(SPA)=cos(SZA).Math.cos(PZA)+sin(SZA).Math.sin(PZA).Math.cos(|SAA−PAA|) (e) obtaining a High Dynamic Range (HDR) image of the sky; (f) calculating the value of the Global Horizontal Irradiance index σ.sub.GHI, the value of the Direct Normal Irradiance index σ.sub.DNI, and the value of the Diffuse Horizontal Irradiance index σ.sub.DHI, defined as the sum of the luminances measured by the pixels of said image belonging to a first region of interest Λ.sub.1, to a second region of interest Λ.sub.2, and to a third region of interest Λ.sub.3, respectively, where Ω.sub.p, L.sub.p and PZA are the solid angle subtended by a pixel p, the luminance measured by this pixel and the angle between this pixel and the zenith, respectively: ${\sigma }_{\mathrm{GHI}}=\frac{{.Math.}_{p\in {\Lambda }_1}{L}_p{\Omega }_p\cos \left(\mathrm{PZA}\right)}{{.Math.}_{p\in {\Lambda }_1}{\Omega }_p}$ ${\sigma }_{\mathrm{DNI}}=\frac{{.Math.}_{p\in {\Lambda }_2}{L}_p{\Omega }_p}{{.Math.}_{p\in {\Lambda }_2}{\Omega }_p}$ ${\sigma }_{\mathrm{DHI}}=\frac{{.Math.}_{p\in {\Lambda }_3}{L}_p{\Omega }_p}{{.Math.}_{p\in {\Lambda }_3}{\Omega }_p}$ (g) carrying out a conversion between the Global Horizontal Irradiance index σ.sub.GHI, the Direct Normal Irradiance index σ.sub.DNI, and the Diffuse Horizontal Irradiance index σ.sub.DHI, and the Global Horizontal Irradiance, Direct Normal Irradiance and Diffuse Horizontal Irradiance, respectively.

8. The method for measuring solar radiation according to the claim 7, further comprising the following step: (h) calculating the solar profile by measuring the Direct Normal Irradiance along a row of pixels passing through the center of the Sun in the image, then repeating this measurement for a plurality of angles, in order to obtain a radial distribution of the luminance around the Sun.

9. The system of measurement of the solar radiation according to claim 2, wherein said first region of interest Λ.sub.1 is the set of pixels of the image verifying PZA<90°, said second region of interest Λ.sub.2 is the set of pixels of the image verifying SPA<β, β being the angle defined as the limit between the circumsolar zone and the rest of the sky, and said third region of interest Λ.sub.3 is the set of pixels of the image verifying {SPA<β}∩{PZA<90°}.

10. The solar radiation measuring system according to claim 2, comprising a waterproof and dustproof protective box in which is housed said camera and which comprises a transparent window located in front of the objective of said camera.

11. The solar radiation measuring system according to claim 3, comprising a waterproof and dustproof protective box in which is housed said camera and which comprises a transparent window located in front of the objective of said camera.

12. The solar radiation measurement system according to claim 2, at least one filter which is adapted to mitigate or modulate the light entering through the objective of said camera.

13. The solar radiation measurement system according to claim 3, at least one filter which is adapted to mitigate or modulate the light entering through the objective of said camera.

14. The solar radiation measurement system according to claim 4, at least one filter which is adapted to mitigate or modulate the light entering through the objective of said camera.

15. The solar radiation measurement system according to claim 5, at least one filter which is adapted to mitigate or modulate the light entering through the objective of said camera.

16. The solar radiation measuring system according to claim 9, comprising a waterproof and dustproof protective box in which is housed said camera and which comprises a transparent window located in front of the objective of said camera.

17. The solar radiation measurement system according to claim 9, at least one filter which is adapted to mitigate or modulate the light entering through the objective of said camera.

18. The solar radiation measurement system according to claim 10, at least one filter which is adapted to mitigate or modulate the light entering through the objective of said camera.

19. The solar radiation measurement system according to claim 11, at least one filter which is adapted to mitigate or modulate the light entering through the objective of said camera.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0054] The invention will be better understood and its advantages will become more apparent upon reading the detailed description which follows, of an embodiment shown by way of non-limiting example. The description refers to the accompanying drawings in which:

[0055] FIG. 1 is a schematic view of a measurement system according to the invention;

[0056] FIG. 2 is a diagram illustrating the correspondence between a camera sensor of the pixel of the system according to the invention and a solid angle;

[0057] FIG. 3 is a diagram showing the different angles used to define the position of the Sun and of a pixel of the camera of the system according to the invention with respect to a reference.

DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0058] The system according to the invention comprises a camera 10, equipped with a hemispherical objective 20 (also called a “fisheye” objective), i.e. a lens with a large angle of field, with a field of vision greater than 180°. The camera 10 includes a sensor 12 capable of capturing light to generate an image, for example a CMOS (Complementary Metal Oxide Semiconductor) sensor or a CCD (Charge Coupled Device) sensor. The sensor 12 comprises on its surface a certain number of pixels (125) on which the light terminates to form an image.

[0059] Advantageously, the camera 10 is located on a horizontal platform, and is located in a place where there is no shading (for example, high up and far from buildings and vegetation). Thus, one can obtain an image of all, or almost all of the sky, without obstacles hiding parts of the sky.

[0060] The camera 10 is oriented along a Z axis which points from the ground towards the zenith.

[0061] Advantageously, the measurement system comprises a protective case 40 that is waterproof and dustproof, and in which the camera 10 is housed.

[0062] The box 40 has a transparent window 42, for example a transparent dome as shown in FIG. 1, which is located in front of the objective 20 of the camera 10 and which protects this objective 20. Thus, the transparent window 42 is located between the objective 20 and the sun.

[0063] Thus, the camera 10 is protected from water and from dust, while letting sunlight reach the objective 20.

[0064] Advantageously, the protective box 40 is equipped with a temperature regulator 50 capable of regulating the temperature inside the box 40.

[0065] Thus, condensation inside the box 40 is reduced or prevented. Consequently, the images taken by the camera 10 are not polluted by humidity.

[0066] Advantageously, the system comprises one or more filters 60, which are capable of attenuating or modulating the light entering the camera 10.

[0067] The filter(s) 60 are positioned between the Sun and the camera 10. Thus, the filter(s) 60 are located either between the camera 10 and the objective 20, or in the objective 20, or in front of the objective 20, or a combination of these positions.

[0068] For example, the internal face of the transparent window 42 is covered with a filter 60, as shown in FIG. 1 in the case where the window 42 is a dome.

[0069] Advantageously, the transparent window 42 is tinted in order to attenuate the stray reflections (which could, mistakenly, be taken for clouds by the algorithms developed).

[0070] Alternatively, the transparent window 42 is not tinted. In this case, the processor 30 (see below) includes an algorithm that suppresses these parasitic reflections.

[0071] The system according to the invention comprises a processor 30 connected to the camera 10.

[0072] This processor 30 is able to perform the following operations: [0073] Perform a first calibration of the camera 10 to obtain the correspondence between the coordinate system of the pixel 125 of the image produced by the sensor 12 and the coordinate system of the camera 10. This calibration geometry is performed, in a known manner, for example, with the assistance of a checkerboard and an algorithm. [0074] Calculate, in a known manner, the solid angle Ωp subtended by each pixel (125, p) of the sensor 12 of the camera 10 recording the images, in order to weight the value of the luminance L in the subsequent calculations.

[0075] FIG. 2 illustrates a solid angle Ωp and shows the correspondence between this solid angle Ωp and a pixel (125, p) of the sensor 12. [0076] Perform a second calibration of the camera 10 by comparing the theoretical position of the Sun with its position obtained on the image. In this way, the correspondence between the coordinate system of the camera 10 and the cardinal points (mark X, Y, Z) is obtained. The theoretical position of the Sun is obtained in a known manner, for example, with the aid of an algorithm such as SPA (Solar Position Algorithm) or SG (Solar Geometry). [0077] Calculate the angle between each pixel (125, p) and the zenith (PZA), the angle between each pixel and the azimuth (PAA), the angle between the Sun and the zenith (SZA), and the angle between the Sun and the azimuth (SAA), then the angle between each pixel (125, p) and the Sun (SPA), by the relationship:

cos(SPA)=cos(SZA).Math.cos(PZA)+sin(SZA).Math.sin(PZA).Math.cos(|SAA−PAA|)

[0078] FIG. 3 shows the relationship between these different angles, the Sun S, and the zenith Z, as well as the north reference direction (N, cardinal points).

[0079] A High Dynamic Range (HDR) image is defined as an image on which is faithfully restored (i.e. without underexposure or overexposure) all the luminance of the sky, from the luminance L.sub.min of the darkest zone at the luminance L.sub.max of the brightest area (Sun). Taking account of all of the luminance by a camera means that the camera sensor has a dynamic range DR.sub.cam at least equal to the dynamic range of the sky DR.sub.sky with:

[00004]${\mathrm{DR}}_{\mathrm{cam}}=20{\log }_{10}\left(\frac{L_2}{L_1}\right)\mathrm{and}{\mathrm{DR}}_{\mathrm{sky}}=20{\log }_{10}\left(\frac{L_{\max }}{L_{\mathrm{mi}n}}\right)$

[0080] where L.sub.1 is the luminance minimum measurable by the sensor of the camera and L.sub.2 is the maximum luminance measurable by the camera sensor.

[0081] If this condition is not verified, i.e. if the sensor 12 of the camera 10 is capable of taking only Low Dynamic Range (LDR) images, the High Dynamic Range image is obtained by the acquisition of several LDR images and then the combination of these images into one HDR image, while preserving the linearity of the sensor, in a known manner. An example of a method for obtaining an HDR image is described below.

[0082] Either several HDR images taken with different exposure times, wherein the total acquisition time of the LDR image sequence must be short enough to guarantee that the movement of the clouds between each LDR image is negligible, in other words that the luminance of the scene is constant during the acquisition. In addition, it is important to properly configure the desired number of LDR images as well as the minimum exposure times t.sub.min and maximum exposure t.sub.max so that the dynamic range that can theoretically be scanned by the HDR image verifies the relationship:

[00005]${\mathrm{DR}}_H={\mathrm{DR}}_L+20{\log }_{10}\left(\frac{t_{\max }}{t_{\min }}\right)>{\mathrm{DR}}_{\mathrm{sky}}$

[0083] where DR.sub.H is the dynamic range of the HDR image and DR.sub.L is the dynamic range of the LDR image. The estimated value of each pixel of the generated HDR image is obtained from a weighted average of the pixel value of each LDR image (considering only the pixels within the linearity range of the sensor).

[0084] Alternatively, the HDR image is obtained using a sensor 12 with a dynamic range greater than the dynamic range of the DR.sub.sky.

[0085] According to the invention, one then calculates the Global Horizontal Irradiance (σ.sub.GHI) index, the value of the Direct Normal Irradiance (σ.sub.DNI) index, and the Diffuse Horizontal Irradiance (σ.sub.DHI) index, defined as the sum of the luminances measured by the pixels (125) of the image belonging to a first region of interest (Λ.sub.1), to a second region of interest (Λ.sub.2), and to a third region of interest (Λ.sub.3), respectively, these sums being weighted by the solid angle of the region considered.

[0086] Thus, these lighting indices are obtained as follows:

[00006]${\sigma }_{\mathrm{GHI}}=\frac{{.Math.}_{p\in {\Lambda }_1}{L}_p{\Omega }_p\cos \left(\mathrm{PZA}\right)}{{.Math.}_{p\in {\Lambda }_1}{\Omega }_p}$${\sigma }_{\mathrm{DNI}}=\frac{{.Math.}_{p\in {\Lambda }_2}{L}_p{\Omega }_p}{{.Math.}_{p\in {\Lambda }_2}{\Omega }_p}$${\sigma }_{\mathrm{DHI}}=\frac{{.Math.}_{p\in {\Lambda }_3}{L}_p{\Omega }_p}{{.Math.}_{p\in {\Lambda }_3}{\Omega }_p}$

[0087] where Ωp is the solid angle subtended by a pixel p (see FIG. 2), Lp is the luminance measured by this pixel and PZA is the angle between this pixel and the zenith (see FIG. 3).

[0088] Advantageously, the first region of interest (Λ.sub.1) is the set of pixels (125) of the image verifying PZA<90°, the second region of interest (Λ.sub.2) is the set of pixels (125) image verifying SPA<β, where β is the angle defined as the boundary between the circumsolar zone (as defined previously above) and the rest of the sky, and the third region of interest (Λ.sub.3) is the set of pixels (125) of the image verifying {SPA<β}∩{PZA<90° }.

[0089] This gives a more precise evaluation of the lighting indices.

[0090] Then Global Horizontal Irradiance (σ.sub.GHI) index, the Direct Normal Irradiance (σ.sub.DNI) index, and the Diffuse Horizontal Irradiance (σ.sub.DHI) index, are converted to Global Horizontal Irradiance (GHI), Direct Normal Irradiance (DNI) and Diffuse Horizontal Irradiance (DHI), respectively, by the conversion functions

[0091] These conversion functions consist of correlations between the irradiation indices (σ.sub.GHI, σ.sub.DNI and σ.sub.DHI) and the various components of the solar irradiance (GHI, DNI and DHI).

[0092] A simple or even multiple linear regression when several lighting indices and, possibly, their histories are taken into account, should make it possible to approximate each conversion function.

[0093] Alternatively, it is possible to approximate these functions through the use of artificial intelligence tools, for example a network of artificial neurones or a neuro-fuzzy system (such a system combines the structure of a connectionist neural network and writing fuzzy rules; fuzzy logic is an extension of Boolean logic). A digital learning phase allows, from a base of examples (learning is supervised here), one to determine the topology of the considered tool (structure) and to identify the parameters.

[0094] The invention also relates to a process for measuring solar radiation using a camera 10 provided with a hemispherical objective 20 and comprising a sensor adapted to capture light to generate an image, and a processor 30, which comprises the following steps, already described above:

(a) Geometric calibration to obtain the correspondence between the coordinate system of the pixels of the image provided by the camera and the coordinates system of the camera.
(b) Calculation of the solid angle subtended by each pixel of the image, in order to weight the value of the luminance in the subsequent calculations.
(c) Second calibration based on the comparison between the theoretical position of the sun and its position obtained on the image, in order to obtain the correspondence between the coordinate system of the camera and the cardinal points.
(d) Calculation of the angle between each pixel (125, p) and the zenith (PZA), of the angle between each pixel and the azimuth (PAA), of the angle between the Sun and the zenith (SZA), and the angle between the sun and the azimuth (SAA), and the angle between each pixel (125, p) and the sun (SPA), by the relationship:

cos(SPA)=cos(SZA).Math.cos(PZA)+sin(SZA).Math.sin(PZA).Math.cos(|SAA−PAA|)

(e) Obtaining an HDR image of the sky.
(f) Calculation of a Global Horizontal Irradiance (σ.sub.GHI) index, of a Direct Normal Irradiance (σ.sub.DNI) index, and of the Diffuse Horizontal Irradiance (σ.sub.DHI) index, where Ω.sub.p, L.sub.p and PZA are the solid angle subtended by a pixel p, the luminance measured by this pixel, and the angle between this pixel and the zenith, respectively:

[00007]${\sigma }_{\mathrm{GHI}}=\frac{{.Math.}_{p\in {\Lambda }_1}{L}_p{\Omega }_p\cos \left(\mathrm{PZA}\right)}{{.Math.}_{p\in {\Lambda }_1}{\Omega }_p}$${\sigma }_{\mathrm{DNI}}=\frac{{.Math.}_{p\in {\Lambda }_2}{L}_p{\Omega }_p}{{.Math.}_{p\in {\Lambda }_2}{\Omega }_p}$${\sigma }_{\mathrm{DHI}}=\frac{{.Math.}_{p\in {\Lambda }_3}{L}_p{\Omega }_p}{{.Math.}_{p\in {\Lambda }_3}{\Omega }_p}$

(g) Conversion between σ.sub.GHI, σ.sub.DNI and σ.sub.DHI and GHI, DNI and DHI.

[0095] Advantageously, the method according to the invention further comprises the following step:

(i) Calculation of the sunshine duration from the direct sunlight.

[0096] The method according to the invention then comprises providing a time counter that measures the time when the duration of sunlight is greater than the threshold previously described above, and which stops when the duration of sunlight becomes less than this threshold. We then know the duration of sunshine over a chosen period, for example, a day, a month, or a year.

[0097] Advantageously, the method according to the invention further comprises the following step:

[0098] (j) Calculation of the pheric disorder from the determination of the “clear sky” moments and the DNI.

[0099] The atmospheric disorder relates to the attenuations of the solar radiation during its passage through of the Earth's atmosphere. This attenuation is mainly due to the phenomena of diffusion by aerosols and absorption by the various atmospheric components (ozone, water vapor, oxygen . . . ).

[0100] The atmospheric disorder gives a good indication of the quality of the solar resource or the atmosphere for a given site.

[0101] From the images acquired by the camera 10, the system according to the invention is capable of determining the evolution of the atmospheric disorder on the site studied.

[0102] For this, the system detects, with the aid of an algorithm, the moments when the sun is not obscured by a cloud. This is called the “clear sky” position.

[0103] An algorithm for detecting “clear sky” instants therefore proceeds to acquire an image or several successive images taken with different exposure times, in a relatively short period of time (for example, a few milliseconds). The system then analyzes the distribution of the luminance in the circumsolar zone and, if necessary, detects a contour of the Sun, in a known manner.

[0104] From these parameters so obtained, the system uses this algorithm to determine whether the instant corresponds or not to a so-called “clear sky” situation.

[0105] As an option, time filtering, based on several successive images, can be used to ensure that this situation is correctly detected.

[0106] Each time a “clear sky” instant is determined by the algorithm, an estimate of the atmospheric disorder (TA) is made.

[0107] For example, this estimate is made using the following relationship:

[00008]$\mathrm{TA}=11.1\frac{\ln \left(b\frac{I_0}{I_{\mathrm{cc}}}\right)}{m}+1\mathrm{with}b=0.664+0.163.Math.\exp \left(\frac{-h}{8000}\right)$

[0108] where h is the altitude of the site considered, l.sub.0 is the extraterrestrial solar illumination (W.Math.m.sup.2), I.sub.cc is the DNI per clear sky (W.Math.m.sup.−2), and m is the relative optical air mass.

[0109] Advantageously, the method according to the invention further comprises the following step:

(h) Calculation of the solar profile from the analysis of the radial distribution of the luminance.

[0110] According to the invention, to obtain the solar profile, the DNI is measured along a line of pixels passing through the center of the Sun on the image. By repeating this measurement for a plurality of angles, we thus obtain a radial distribution of the luminance L. This avoids advantageously the use of the instrument SAM (see above), the cost of which is very high.

[0111] In addition, the measurement of the solar profile according to the invention is more precise than a measurement with the SAM instrument, because the measurement according to the invention takes into account the anisotropy of the luminance around the Sun (whereas the SAM instrument supposes, incorrectly, a radial symmetry of solar radiation).