Shielding area determination method and device, storage medium, electronic device and computer program product

By constructing the first matrix and calculating the vector intersection of the photovoltaic array, the problem of low accuracy and low efficiency of the traditional photovoltaic module shading area calculation method is solved, a more accurate and efficient shadow shading area calculation is achieved, and the performance of the photovoltaic system is optimized.

CN119322911APending Publication Date: 2025-01-17HUANENG CLEAN ENERGY RES INST
View PDF 0 Cites 0 Cited by

Patent Information

Application Number
CN202411418301.6
Authority / Receiving Office
CN · China
Patent Type
Applications(China)
Current Assignee / Owner
Filing Date
2024-10-11
Publication Date
2025-01-17

AI Technical Summary

Technical Problem

The calculation accuracy and efficiency of traditional photovoltaic module occlusion area calculation methods are low, and cannot effectively solve the performance degradation caused by shadow occlusion of photovoltaic arrays.

Method used

By constructing the first matrix, recording the point coordinates in each photovoltaic array set, and calculating the vector intersection points related to the solar altitude angle and azimuth angle, the shadow blocking area is determined.

Benefits of technology

The accuracy and efficiency of photovoltaic module shading area calculation are improved, and the shadow shading area of ​​the photovoltaic array can be calculated more accurately, thereby optimizing the layout and performance of the photovoltaic system.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN119322911A_ABST
    Figure CN119322911A_ABST
Patent Text Reader

Abstract

The invention discloses a method and device for determining a shielding area, a storage medium, an electronic device and a computer program product, and relates to the field of photovoltaic power station detection.The method for determining the shielding area comprises the steps that under the condition that N photovoltaic arrays exist in a target area, a first matrix corresponding to the N photovoltaic arrays is determined; determining N * P * M vectors according to N * P * M coordinates of N * P * M points in the first matrix; n * P * M intersection points of the N * P * M vectors and N planes corresponding to the N photovoltaic arrays are determined, the number of target grids in multiple grids corresponding to the ith photovoltaic array in the N photovoltaic arrays is determined, and the ith number is obtained; and determining the shadow shielding area of the ith photovoltaic array according to the ith number and the area of each grid.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] The present application relates to the field of photovoltaic power station detection, and specifically, to a method and device for determining a shielding area, a storage medium, an electronic device, and a computer program product. Background Art

[0002] In the design and management of solar photovoltaic systems, ensuring that the photovoltaic array receives the maximum possible sunshine is one of the key factors in improving energy efficiency and system performance. However, due to environmental factors, installation location, and the influence of surrounding buildings or other obstacles, photovoltaic arrays are often subject to partial or complete shadowing. Shadowing not only reduces the amount of light received by the components, but may also cause the performance of the photovoltaic array to deteriorate, and may even cause hot spot effects due to inconsistent currents generated by different parts of the photovoltaic array, causing damage.

[0003] Traditional methods for calculating the shading area usually rely on simplified calculation models or manual measurement and adjustment of each photovoltaic array. However, these methods are often time-consuming, labor-intensive, inaccurate, and have low calculation efficiency.

[0004] In view of the problem of low calculation accuracy and efficiency of the traditional photovoltaic module shading area calculation method in the related technology, no effective solution has been proposed so far.

[0005] Therefore, it is necessary to improve the related technology to overcome the above-mentioned defects in the related technology. Summary of the invention

[0006] The embodiments of the present application provide a method and device for determining a shading area, a storage medium, an electronic device, and a computer program product to at least solve the problem of low calculation accuracy and efficiency of traditional photovoltaic module shading area calculation methods.

[0007] According to one aspect of an embodiment of the present application, a method for determining a blocked area is provided, comprising: in a case where there are N photovoltaic arrays in a target area, determining a first matrix corresponding to the N photovoltaic arrays, wherein the dimension of the first matrix is ​​N*P*M*3, the first matrix records the coordinates of M points corresponding to each photovoltaic array in the i-th array set, the i-th array set includes P photovoltaic arrays adjacent to the i-th photovoltaic array in the N photovoltaic arrays, the coordinates of each point in the first matrix are the coordinates of the photovoltaic array corresponding to the point in the spatial coordinate system, i is an integer greater than or equal to 1 and less than or equal to N, P is an integer greater than or equal to 1 and less than or equal to N-1, and N and M are both integers greater than or equal to 2; Determine N*P*M vectors according to the N*P*M coordinates of the N*P*M points in the first matrix, wherein the j-th vector among the N*P*M vectors is a vector determined according to the coordinates of the j-th point and the altitude angle and azimuth of the sun at the current moment, and j is an integer greater than or equal to 1 and less than or equal to N*P*M; determine N*P*M intersections of the N*P*M vectors and the N planes corresponding to the N photovoltaic arrays, and determine the number of target grids in the multiple grids corresponding to the i-th photovoltaic array among the N photovoltaic arrays to obtain the i-th number, wherein the target grid is a grid having the intersection; determine the shadow shielding area of ​​the i-th photovoltaic array according to the i-th number and the area of ​​each grid.

[0008] In an exemplary embodiment, determining a first matrix corresponding to the N photovoltaic arrays includes: determining a second matrix corresponding to the N photovoltaic arrays, wherein the dimension of the second matrix is ​​N*P*M*3, the second matrix records the coordinates of M points corresponding to each photovoltaic array in the i-th array set, and the coordinates of the points in the second matrix are coordinates in a geodetic coordinate system; and performing a coordinate system conversion on the coordinates of the points in the second matrix to obtain the first matrix.

[0009] In an exemplary embodiment, the coordinates of the points in the second matrix are transformed into a coordinate system to obtain the first matrix, including: determining a third matrix, wherein the third matrix records the three-dimensional vectors of the spatial coordinate system of the photovoltaic array corresponding to each point in the first matrix on the X-axis, the Y-axis and the Z-axis, the dimension of the third matrix is ​​N*P*3*3, the X-axis of the spatial coordinate system of the photovoltaic array is parallel to the line segment where the width of the photovoltaic array is located, the Y-axis of the spatial coordinate system of the photovoltaic array is parallel to the line segment where the length of the photovoltaic array is located, and the Z-axis of the spatial coordinate system of the photovoltaic array is perpendicular to the plane where the photovoltaic array is located; based on the third matrix, the coordinates of the points in the second matrix are transformed into a coordinate system to obtain the first matrix.

[0010] In an exemplary embodiment, determining the third matrix includes: obtaining the vertex coordinates of the two side end columns of each photovoltaic array according to the central axis of the long side of each photovoltaic array in the N photovoltaic arrays, obtaining a fourth matrix and a fifth matrix, wherein the fourth matrix photovoltaic has the vertex coordinates of the left end column of each photovoltaic array, and the fifth matrix has the vertex coordinates of the right end column of each photovoltaic array, and the dimensions of the fourth matrix and the fifth matrix are both N*3; determining a sixth matrix, wherein the sixth matrix has the deviation angle between each photovoltaic array in the N photovoltaic arrays and the corresponding end column; determining a seventh matrix, wherein the seventh matrix has the center point coordinates of each photovoltaic array in the N photovoltaic arrays, the dimension of the seventh matrix is ​​N*3, and the coordinates in the fourth matrix, the fifth matrix and the seventh matrix are coordinates in the geodetic coordinate system; determining the three-dimensional vectors of the spatial coordinate system corresponding to each photovoltaic array on the X-axis, the Y-axis and the Z-axis according to the fourth matrix, the fifth matrix, the sixth matrix and the seventh matrix, obtaining an eighth matrix, wherein the dimension of the eighth matrix is ​​N*3*3; determining the third matrix according to the eighth matrix.

[0011] In an exemplary embodiment, determining a second matrix corresponding to the N photovoltaic arrays includes: determining a maximum length and a maximum width of the N photovoltaic arrays, and constructing a standard rectangle according to the maximum length and the maximum width; rasterizing the standard rectangle to form M points; determining M points corresponding to each of the N photovoltaic arrays, and obtaining a ninth matrix according to the coordinates of the M points corresponding to each photovoltaic array, wherein the dimension of the ninth matrix is ​​N*M*3, the ninth matrix records the coordinates of the M points corresponding to each photovoltaic array, the coordinates of the points in the ninth matrix are coordinates in a geodetic coordinate system, and when the size of the photovoltaic array is smaller than the standard rectangle, the coordinates of the points among the M points that fall outside the photovoltaic array are determined to be (0, 0, 0); and determining the second matrix according to the ninth matrix.

[0012] In an exemplary embodiment, determining the shadow blocking area of ​​the ith photovoltaic array according to the ith number and the area of ​​each grid includes: multiplying the ith number by the area of ​​each grid to obtain the shadow blocking area of ​​the ith photovoltaic array.

[0013] According to another aspect of an embodiment of the present application, a device for determining a blocked area is also provided, including: a first determination module, for determining a first matrix corresponding to the N photovoltaic arrays when there are N photovoltaic arrays in a target area, wherein the dimension of the first matrix is ​​N*P*M*3, the first matrix records the coordinates of M points corresponding to each photovoltaic array in the i-th array set, the i-th array set includes P photovoltaic arrays adjacent to the i-th photovoltaic array in the N photovoltaic arrays, the coordinates of each point in the first matrix are the coordinates of the photovoltaic array corresponding to the point in the spatial coordinate system, i is an integer greater than or equal to 1 and less than or equal to N, P is an integer greater than or equal to 1 and less than or equal to N-1, and N and M are both integers greater than or equal to 2; a second determination module, for determining a first matrix corresponding to the N photovoltaic arrays when there are N photovoltaic arrays in a target area, wherein the dimension of the first matrix is ​​N*P*M*3, the first matrix records the coordinates of M points corresponding to each photovoltaic array in the i-th array set, the i-th array set includes P photovoltaic arrays adjacent to the i-th photovoltaic array in the N photovoltaic arrays, the coordinates of each point in the first matrix are the coordinates of the photovoltaic array corresponding to the point in the spatial coordinate system, i is an integer greater than or equal to 1 and less than or equal to N, P is an integer greater than or equal to 1 and less than or equal to N-1, and N and M are both integers greater than or equal to 2; The method comprises the following steps: determining N*P*M vectors according to the N*P*M coordinates of the N*P*M points in the first matrix, wherein the j-th vector among the N*P*M vectors is a vector determined according to the coordinates of the j-th point and the altitude angle and azimuth angle of the sun at the current moment, and j is an integer greater than or equal to 1 and less than or equal to N*P*M; a third determination module is used to determine the N*P*M intersections of the N*P*M vectors and the N planes corresponding to the N photovoltaic arrays, and determine the number of target grids among the multiple grids corresponding to the i-th photovoltaic array among the N photovoltaic arrays to obtain the i-th number, wherein the target grid is a grid having the intersection; and a fourth determination module is used to determine the shadow shielding area of ​​the i-th photovoltaic array according to the i-th number and the area of ​​each grid.

[0014] According to another aspect of the embodiments of the present application, a computer-readable storage medium is further provided, wherein the computer-readable storage medium includes a stored program, wherein the program is configured to execute the above-mentioned method for determining the occlusion area when running.

[0015] According to another aspect of the embodiments of the present application, there is further provided an electronic device, including a memory and a processor, wherein the memory stores a computer program, wherein the processor is configured to execute the above-mentioned method for determining the occlusion area through the computer program.

[0016] According to another aspect of the embodiments of the present application, a computer program product is provided, including a computer program, and when the computer program is executed by a processor, the above-mentioned method for determining the occlusion area is implemented.

[0017] In the present application, the coordinates of M points corresponding to each photovoltaic array in the ith array set are recorded through a first matrix, and the N*P*M intersections of the N*P*M vectors determined according to the first matrix and the N planes corresponding to the N photovoltaic arrays are determined, so that the number of grids with intersections can be obtained, and then the shadow blocking area of ​​the ith photovoltaic array is obtained according to the number of grids and the area of ​​each grid. Since parallel calculation is used when calculating the shadow blocking area of ​​the ith photovoltaic array, the shadow blocking area of ​​the ith photovoltaic array is calculated more accurately, thereby improving the calculation accuracy and efficiency, thereby solving the problem of low calculation accuracy and efficiency of the traditional photovoltaic module blocking area calculation method. BRIEF DESCRIPTION OF THE DRAWINGS

[0018] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the present application and, together with the description, serve to explain the principles of the present application.

[0019] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, for ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative labor.

[0020] Figure 1 It is a hardware structure block diagram of a mobile terminal of a method for determining an obstruction area according to an embodiment of the present application;

[0021] Figure 2 is a flow chart of a method for determining an occlusion area according to an embodiment of the present application;

[0022] Figure 3 is a schematic diagram of an optional plane space coordinate system for constructing an array according to an embodiment of the present application;

[0023] Figure 4 is a schematic diagram of an optional construction of a standard rectangle according to an embodiment of the present application;

[0024] Figure 5 It is a structural block diagram of a device for determining an occlusion area according to an embodiment of the present application. DETAILED DESCRIPTION

[0025] In order to enable those skilled in the art to better understand the solution of the present application, the technical solution in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of the present application.

[0026] It should be noted that the terms "first", "second", etc. in the specification and claims of the present application and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchangeable where appropriate, so that the embodiments of the present application described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any of their variations are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device comprising a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0027] The method embodiments provided in the embodiments of the present application can be executed in a mobile terminal, a computer terminal or a similar computing device. Taking running on a mobile terminal as an example, Figure 1 FIG. 1 is a hardware structure block diagram of a mobile terminal of a method for determining an obstruction area according to an embodiment of the present application. Figure 1 As shown, the mobile terminal may include one or more ( Figure 1 Only one is shown in the figure) a processor 102 (the processor 102 may include but is not limited to a processing device such as a microprocessor (MP) or a programmable logic device (FPGA)) and a memory 104 for storing data, wherein the mobile terminal may also include a transmission device 106 and an input / output device 108 for communication functions. It can be understood by those skilled in the art that Figure 1 The structure shown is only for illustration and does not limit the structure of the mobile terminal. Figure 1 More or fewer components as shown, or with Figure 1 Different configurations are shown.

[0028] The memory 104 can be used to store computer programs, for example, software programs and modules of application software, such as the computer program corresponding to the method for detecting pneumatic imbalance in the embodiment of the present application. The processor 102 executes various functional applications and data processing by running the computer program stored in the memory 104, that is, to implement the above method. The memory 104 may include a high-speed random access memory, and may also include a non-volatile memory, such as one or more magnetic storage devices, flash memory, or other non-volatile solid-state memory. In some examples, the memory 104 may further include a memory remotely arranged relative to the processor 102, and these remote memories can be connected to the mobile terminal via a network. Examples of the above-mentioned network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

[0029] The transmission device 106 is used to receive or send data via a network. The specific example of the above network may include a wireless network provided by a communication provider of the mobile terminal. In one example, the transmission device 106 includes a network adapter (Network Interface Controller, referred to as NIC), which can be connected to other network devices through a base station so as to communicate with the Internet. In one example, the transmission device 106 can be a radio frequency (RF) module, which is used to communicate with the Internet wirelessly.

[0030] In order to solve the above problem, this embodiment provides a method for determining a blocked area, including but not limited to being applied to the above mobile terminal. Figure 2 is a flow chart of a method for determining an occlusion area according to an embodiment of the present application, such as Figure 2 As shown, the process includes the following steps S202-S208:

[0031] Step S202: in the case where there are N photovoltaic arrays in the target area, determine a first matrix corresponding to the N photovoltaic arrays, wherein the dimension of the first matrix is ​​N*P*M*3, the first matrix records the coordinates of M points corresponding to each photovoltaic array in the i-th array set, the i-th array set includes P photovoltaic arrays adjacent to the i-th photovoltaic array in the N photovoltaic arrays, the coordinates of each point in the first matrix are the coordinates of the photovoltaic array corresponding to the point in the spatial coordinate system, i is an integer greater than or equal to 1 and less than or equal to N, P is an integer greater than or equal to 1 and less than or equal to N-1, and N and M are both integers greater than or equal to 2;

[0032] Optionally, the dimension of the first matrix is ​​N*P*M*3, where "3" represents coordinate information in three-dimensional space, that is, the X, Y, and Z coordinate values ​​of each point in the M points corresponding to each photovoltaic array in the i-th array set in the spatial coordinate system of the photovoltaic array corresponding to the point.

[0033] Optionally, P may be equal to 5, indicating that there are 5 adjacent photovoltaic arrays to the i-th photovoltaic array, and the dimension of the first matrix corresponding to the N photovoltaic arrays is N*5*M*3.

[0034] It should be noted that the i-th array set includes the P photovoltaic arrays adjacent to the i-th photovoltaic array among the N photovoltaic arrays in order to comprehensively and accurately analyze the shadow occlusion situation. In a photovoltaic system, the shadow of a photovoltaic array may not only be caused by obstacles directly above it, but may also be affected by other surrounding arrays, especially adjacent arrays. When the position of the sun changes, the projection of the adjacent array may block the i-th photovoltaic array, resulting in a reduction in the amount of sunlight it receives. Therefore, taking into account the information of the P arrays adjacent to the i-th photovoltaic array can more accurately simulate and predict the shadow occlusion situation. By constructing a three-dimensional matrix including the coordinates of these adjacent array points, the shadow changes of each array under different sun positions can be analyzed in more detail, which is crucial for optimizing the array layout, reducing shadow occlusion, and improving the power generation efficiency of the entire photovoltaic system.

[0035] It should be noted that the selection of the P value depends on the compactness of the array layout and the influence range of the surrounding environment. Usually, P will choose a smaller value to ensure calculation efficiency while covering adjacent arrays that may have a significant impact on the occlusion analysis.

[0036] Step S204: determining N*P*M vectors according to the N*P*M coordinates of the N*P*M points in the first matrix, wherein the j-th vector among the N*P*M vectors is a vector determined according to the coordinates of the j-th point and the altitude angle and azimuth angle of the sun at the current moment, and j is an integer greater than or equal to 1 and less than or equal to N*P*M;

[0037] Optionally, taking P equal to 5 as an example, the N*5*M points in the first matrix (with a dimension of N*5*M*3) are determined as the N*P*M vectors formed by the N*5*M vectors according to the solar altitude angle (ranging from 0 to 90°, for example, 45°) and the solar azimuth angle (ranging from -180° to 180°, for example, 5°) at a certain moment.

[0038] It should be noted that the solar altitude angle refers to the height of the sun's position in the sky relative to the horizon, ranging from 0° (horizontal) to 90° (directly above). When performing shadow shading analysis on photovoltaic modules, the solar altitude angle is one of the important parameters for calculating the incident angle and direction of sunlight; the solar azimuth angle refers to the angle of the sun on the horizontal plane relative to the north direction, ranging from -180° to 180°, where 0° points to the north, 90° points to the east, and 180° points to the south.

[0039] It should be noted that the solar altitude angle and solar azimuth angle are different at different times.

[0040] Step S206: determining N*P*M intersection points of the N*P*M vectors and the N planes corresponding to the N photovoltaic arrays, and determining the number of target grids in the multiple grids corresponding to the i-th photovoltaic array in the N photovoltaic arrays, to obtain the i-th number, wherein the target grid is a grid having the intersection points;

[0041] Step S208: determining the shadowing area of ​​the i-th photovoltaic array according to the i-th number and the area of ​​each grid.

[0042] It should be noted that if a grid has an intersection (a grid may have multiple intersections), it means that the grid is shaded, and the shadow area of ​​the i-th photovoltaic array can be determined by the number of target grids.

[0043] In the above steps, the coordinates of the M points corresponding to each photovoltaic array in the i-th array set are recorded through the first matrix, and the N*P*M intersections of the N*P*M vectors determined according to the first matrix and the N planes corresponding to the N photovoltaic arrays are determined, so that the number of grids with intersections can be obtained, and then the shadow blocking area of ​​the i-th photovoltaic array is obtained according to the number of grids and the area of ​​each grid. Since parallel calculation is used when calculating the shadow blocking area of ​​the i-th photovoltaic array, the shadow blocking area of ​​the i-th photovoltaic array is calculated more accurately, thereby improving the calculation accuracy and efficiency, thereby solving the problem of low calculation accuracy and efficiency of the traditional photovoltaic module blocking area calculation method.

[0044] In an exemplary embodiment, determining the first matrix corresponding to the N photovoltaic arrays may be implemented by the following steps S11-S12:

[0045] Step S11: determining a second matrix corresponding to the N photovoltaic arrays, wherein the dimension of the second matrix is ​​N*P*M*3, the second matrix records the coordinates of M points corresponding to each photovoltaic array in the i-th array set, and the coordinates of the points in the second matrix are coordinates in the geodetic coordinate system;

[0046] It should be noted that the geodetic coordinate system, also known as the geographic coordinate system, is a global three-dimensional coordinate system used to uniquely determine the position of any point on the earth's surface, usually including three dimensions of longitude, latitude, and altitude. In the embodiment of the present application, the geodetic coordinate system can be used to accurately describe the geographical location of the photovoltaic array and the location of surrounding obstacles. Since photovoltaic arrays may be deployed over a wide geographical range, the use of the geodetic coordinate system can ensure the accuracy and versatility of the coordinate information, allowing the algorithm to calculate and analyze shadow occlusion at any location.

[0047] Step S12: performing coordinate system transformation on the coordinates of the points in the second matrix to obtain the first matrix.

[0048] Optionally, each photovoltaic array in the i-th array set corresponds to M points, and each point has 3 coordinate values ​​(coordinate values ​​x, y, z in the geodetic coordinate system). Therefore, each photovoltaic array has M*3 data points, and for N*5 photovoltaic arrays, there are (N*5)*M*3 data points. Therefore, the coordinate information of the M points of the N*5 photovoltaic arrays constitutes a N*5*M*3 four-dimensional matrix, i.e., the second matrix.

[0049] In an exemplary embodiment, performing coordinate system transformation on the coordinates of the points in the second matrix to obtain the first matrix can be achieved by the following steps S21-S22:

[0050] Step S21: Determine a third matrix, wherein the third matrix records the three-dimensional vectors of the spatial coordinate system of the photovoltaic array corresponding to each point in the first matrix on the X-axis, the Y-axis and the Z-axis, the dimension of the third matrix is ​​N*P*3*3, the X-axis of the spatial coordinate system of the photovoltaic array is parallel to the line segment where the width of the photovoltaic array is located, the Y-axis of the spatial coordinate system of the photovoltaic array is parallel to the line segment where the length of the photovoltaic array is located, and the Z-axis of the spatial coordinate system of the photovoltaic array is perpendicular to the plane where the photovoltaic array is located;

[0051] Step S22: performing coordinate system transformation on the coordinates of the points in the second matrix based on the third matrix to obtain the first matrix.

[0052] Optionally, the coordinates of the points in the second matrix are transformed into a coordinate system based on the third matrix by the following formula to obtain the first matrix:

[0053]

[0054] in,( B P x, B P y,B P z) is the coordinate of a point in the second matrix, ( A Px, A P y, A P z) is the coordinate of the point after the coordinate system transformation, A P x is the coordinate of the point on the x-axis, A P y is the coordinate of the point on the y-axis, A P z is the coordinate of the point on the z-axis;

[0055] xA is the three-dimensional vector of the point in the space coordinate system of the corresponding photovoltaic array on the X axis, yA is the three-dimensional vector of the point in the space coordinate system of the corresponding photovoltaic array on the Y axis, and zA is the three-dimensional vector of the point in the space coordinate system of the corresponding photovoltaic array on the Z axis;

[0056] x B is the three-dimensional vector of the X axis of the geodetic coordinate system corresponding to the point, y B is the three-dimensional vector of the Y axis of the geodetic coordinate system corresponding to the point, and z B is the three-dimensional vector of the Z axis of the geodetic coordinate system corresponding to the point.

[0057] In an exemplary embodiment, determining the third matrix may be implemented by the following steps S31-S35:

[0058] Step S31: acquiring the vertex coordinates of the two side end columns of each photovoltaic array according to the central axis of the long side of each photovoltaic array in the N photovoltaic arrays, and obtaining a fourth matrix and a fifth matrix, wherein the fourth matrix photovoltaic has the vertex coordinates of the left end column of each photovoltaic array, and the fifth matrix has the vertex coordinates of the right end column of each photovoltaic array, and the dimensions of the fourth matrix and the fifth matrix are both N*3;

[0059] Optionally, refer to Figure 3 , for each photovoltaic array (a total of N) in the general plan layout, obtain the vertex coordinates of the two end columns (representing the left and right sides respectively) on the central axis of its long side. These coordinates will eventually be organized into two matrices (the fourth matrix An and the fifth matrix Bn), each with a dimension of N*3, representing N photovoltaic arrays, and the vertex coordinates of each array are represented by 3 values ​​(X, Y, Z coordinates in the geodetic coordinate system).

[0060] Step S32: determining a sixth matrix, wherein the sixth matrix includes a deviation angle between each photovoltaic array in the N photovoltaic arrays and a corresponding end pole;

[0061] Optionally, refer to Figure 3 The deviation angle can be determined jointly using three values, that is, the dimension of the sixth matrix can be N*3.

[0062] Step S33: determining a seventh matrix, wherein the seventh matrix has the coordinates of the center point of each of the N photovoltaic arrays, the dimension of the seventh matrix is ​​N*3, and the coordinates in the fourth matrix, the fifth matrix and the seventh matrix are coordinates in the geodetic coordinate system;

[0063] Optionally, refer to Figure 3 The seventh matrix Dn contains the coordinates of the center point of each of the N photovoltaic arrays, which is also expressed by the X, Y, and Z coordinates of the geodetic coordinate system, and has a dimension of N*3.

[0064] Step S34: determining the three-dimensional vectors of the spatial coordinate system corresponding to each photovoltaic array on the X-axis, the Y-axis and the Z-axis according to the fourth matrix, the fifth matrix, the sixth matrix and the seventh matrix to obtain an eighth matrix, wherein the dimension of the eighth matrix is ​​N*3*3;

[0065] Optionally, the fourth matrix (left end column coordinates), the fifth matrix (right end column coordinates), the sixth matrix (deviation angle) and the seventh matrix (center point coordinates) are used to determine the three-dimensional vectors of the spatial coordinate system corresponding to each photovoltaic array on the X, Y and Z axes. These vectors will be used to construct the corresponding spatial coordinate system of the plane where each array is located, and the results are stored in an eighth matrix with a dimension of N*3*3, which means that each of the N photovoltaic arrays has 3 basis vectors to describe its coordinate system.

[0066] Step S35: Determine the third matrix according to the eighth matrix.

[0067] Optionally, first create an empty third matrix with a dimension of N*P*3*3, which will be used to store the spatial coordinate system information of each photovoltaic array in the N photovoltaic arrays and the P photovoltaic arrays adjacent to the i-th photovoltaic array; then copy the spatial coordinate system information of each photovoltaic array in the eighth matrix to the corresponding position of the third matrix. Specifically, for the i-th photovoltaic array, copy its information to the first sub-matrix of the i-th row in the third matrix (i.e., i, 1, x, y in the third matrix); for the i-th photovoltaic array, find its P adjacent photovoltaic arrays, and then obtain the spatial coordinate system information of the P adjacent photovoltaic arrays; then fill the spatial coordinate system information of each adjacent array into the subsequent P-1 sub-matrices of the i-th row of the third matrix (if the i-th photovoltaic array has the j-th adjacent array (j from 2 to P+1), copy the spatial coordinate system information of the adjacent array to (i, j, x, y) in the third matrix), and so on, to obtain the target third matrix with a dimension of N*P*3*3, where x and y are any values ​​within [-1, 1].

[0068] Optionally, if the i-th photovoltaic array does not have enough P adjacent arrays (such as at the edge of the general plan layout), the third matrix needs to be properly filled, such as using a zero matrix or using specific default coordinate system information to maintain the integrity and consistency of the matrix.

[0069] In an exemplary embodiment, determining the second matrix corresponding to the N photovoltaic arrays may be implemented by the following steps S41-S44:

[0070] Step S41: determining the maximum length and the maximum width of the N photovoltaic arrays, and constructing a standard rectangle according to the maximum length and the maximum width;

[0071] Optionally, Figure 4 The figure shows the size of a standard rectangle in a plane rectangular coordinate system.

[0072] Step S42: rasterizing the standard rectangle to form M points;

[0073] Optionally, the standard rectangle is divided into points at a resolution of 0.2*0.2 to form M points.

[0074] It should be noted that the higher the resolution, the smaller the grid obtained by rasterizing the standard rectangle, and the closer the calculated shadow occlusion area is to the actual situation.

[0075] Step S43: determining M points corresponding to each of the N photovoltaic arrays, and obtaining a ninth matrix according to the coordinates of the M points corresponding to each photovoltaic array, wherein the dimension of the ninth matrix is ​​N*M*3, and the coordinates of the M points corresponding to each photovoltaic array are recorded in the ninth matrix, and the coordinates of the points in the ninth matrix are coordinates in the geodetic coordinate system, and when the size of the photovoltaic array is smaller than the standard rectangle, the coordinates of the points among the M points that fall outside the photovoltaic array are determined to be (0, 0, 0);

[0076] It should be noted that the above steps can accurately process the coordinate information of all photovoltaic arrays and can adapt to photovoltaic arrays of different sizes, especially when the array size is smaller than the defined standard rectangle, the robustness and continuity of the calculation can be maintained.

[0077] Specifically, M points are determined from each of the N photovoltaic arrays to form a ninth matrix of N*M*3, where the M points are determined based on the resolution set by the rasterized standard rectangle (e.g., 0.2*0.2) to ensure that a uniform dot matrix is ​​created on each photovoltaic array, and these points can be used to represent the receiving surface of the array, thereby more accurately calculating the shadow occlusion. When the size of some photovoltaic arrays is smaller than the standard rectangle, the coordinates of the points outside the photovoltaic array are set to (0, 0, 0), because these points outside the photovoltaic array do not actually participate in the calculation of shadow occlusion.

[0078] Step S44: Determine the second matrix according to the ninth matrix.

[0079] In an exemplary embodiment, determining the shadow blocking area of ​​the i-th photovoltaic array based on the i-th number and the area of ​​each grid can be achieved by the following steps: multiplying the i-th number by the area of ​​each grid to obtain the shadow blocking area of ​​the i-th photovoltaic array.

[0080] It should be noted that since there may be multiple intersections in the same grid, but these intersections still cover the area of ​​the same grid, no matter how many intersections there are in the grid, the area of ​​the grid is only calculated once. By counting the number of grids occupied within each photovoltaic array (that is, there is at least one intersection in the grid) and multiplying it by the area of ​​each grid, the total shadow area of ​​the photovoltaic array under specific conditions can be obtained.

[0081] It should be noted that with the development of computing technology, especially the widespread application of graphics processing units (GPUs), it has become possible to use the high parallel processing capabilities of GPUs to perform shadow occlusion analysis. GPU parallel computing can significantly improve computing efficiency, making real-time shadow analysis in complex environments a reality.

[0082] In addition, this algorithm can also be integrated into the monitoring and management tools of photovoltaic systems to help operators monitor system performance in real time and adjust system configuration or perform maintenance in a timely manner to deal with possible shading problems.

[0083] Obviously, the above-described embodiments are only a part of the embodiments of the present invention, rather than all the embodiments. In order to better understand the above method, the above process is described below in conjunction with the embodiments, but it is not intended to limit the technical solutions of the embodiments of the present invention. Figure 3 The overall flow chart of a method for determining a shielding area according to an embodiment of the application is shown, specifically:

[0084] 1. Assume that there are n photovoltaic arrays in the general layout diagram, and each photovoltaic array is a rectangle. Obtain the coordinates of the two end columns An (n*3 matrix) and Bn (n*3 matrix) along the central axis of the long side of the rectangle, and obtain the deviation angle Cn (n*3 matrix) between each photovoltaic array and the end column. Obtain the coordinates of the center point Dn (n*3 matrix) of the photovoltaic array rectangle, and obtain the length and width En (n*2 matrix) of each photovoltaic array;

[0085] 2. Construct the plane space coordinate system of the photovoltaic array based on An, Bn, Cn, and Dn, and obtain the three-dimensional vector (n*3*3) matrix of Xn, Yn, and Zn;

[0086] 3. Obtain the maximum length and width of all photovoltaic arrays, and use them to construct a standard rectangle;

[0087] Fourth, rasterize the standard rectangle, take points at a resolution of 0.2*0.2, and form m points;

[0088] 5. When n photovoltaic arrays are smaller than the standard rectangle, the points outside the photovoltaic arrays are adjusted to (0,0,0), forming a matrix of n*m*3;

[0089] 6. Select the information of 5 adjacent photovoltaic arrays of each photovoltaic array to form 5n photovoltaic arrays, forming a matrix of n*5*m*3;

[0090] 7. Given the vector coordinates of the spatial coordinate system xyz of 5n photovoltaic arrays, a matrix of n*5*3*3 is formed;

[0091] 8. Use the following formula to perform spatial rotation and transform the plane coordinates into a matrix of n*5*m*3 spatial coordinates:

[0092]

[0093] 9. Intersect the vector formed by the n*5*m point according to the solar altitude angle and azimuth at a certain moment with n planes to form n*5*m intersection points;

[0094] 10. Determine whether the n*5*m intersection point is within each grid within the range of each photovoltaic array. If multiple points are contained in the same grid, count them as one. Count the number of occupied grids. The grid area * the number of occupied grids is the shadow shielding area of ​​each photovoltaic array.

[0095] Through the description of the above implementation methods, those skilled in the art can clearly understand that the method according to the above embodiment can be implemented by means of software plus a necessary general hardware platform, and of course by hardware, but in many cases the former is a better implementation method. Based on such an understanding, the technical solution of the present application, or the part that contributes to the prior art, can be embodied in the form of a software product, which is stored in a storage medium (such as ROM / RAM, magnetic disk, optical disk), and includes a number of instructions for a terminal device (which can be a mobile phone, computer, server, or network device, etc.) to execute the methods of each embodiment of the present application.

[0096] In this embodiment, a device for determining a shielding area is also provided, and the device is used to implement the above-mentioned embodiments and preferred implementation modes, and the descriptions that have been made will not be repeated. As used below, the term "module" can implement a combination of software and / or hardware of a predetermined function. Although the devices described in the following embodiments are preferably implemented in software, the implementation of hardware, or a combination of software and hardware, is also possible and conceivable.

[0097] Figure 5 is a structural block diagram of a device for determining an obstruction area according to an embodiment of the present application, the device comprising:

[0098] A first determination module 502 is used to determine a first matrix corresponding to the N photovoltaic arrays when there are N photovoltaic arrays in the target area, wherein the dimension of the first matrix is ​​N*P*M*3, the first matrix records the coordinates of M points corresponding to each photovoltaic array in the i-th array set, the i-th array set includes P photovoltaic arrays adjacent to the i-th photovoltaic array in the N photovoltaic arrays, the coordinates of each point in the first matrix are the coordinates of the photovoltaic array corresponding to the point in the spatial coordinate system, i is an integer greater than or equal to 1 and less than or equal to N, P is an integer greater than or equal to 1 and less than or equal to N-1, and N and M are both integers greater than or equal to 2;

[0099] A second determining module 504 is configured to determine N*P*M vectors according to the N*P*M coordinates of the N*P*M points in the first matrix, wherein the j-th vector among the N*P*M vectors is a vector determined according to the coordinates of the j-th point and the altitude angle and azimuth angle of the sun at the current moment, and j is an integer greater than or equal to 1 and less than or equal to N*P*M;

[0100] The third determination module 506 is used to determine N*P*M intersection points of the N*P*M vectors and the N planes corresponding to the N photovoltaic arrays, and determine the number of target grids in the multiple grids corresponding to the i-th photovoltaic array in the N photovoltaic arrays to obtain the i-th number, wherein the target grid is a grid having the intersection points;

[0101] The fourth determining module 508 is configured to determine the shadowing area of ​​the i-th photovoltaic array according to the i-th number and the area of ​​each grid.

[0102] The above-mentioned device records the coordinates of M points corresponding to each photovoltaic array in the i-th array set through a first matrix, and determines the N*P*M intersections of the N*P*M vectors determined according to the first matrix and the N planes corresponding to the N photovoltaic arrays, so that the number of grids with intersections can be obtained, and then the shadow blocking area of ​​the i-th photovoltaic array is obtained according to the number of grids and the area of ​​each grid. Since parallel calculation is used when calculating the shadow blocking area of ​​the i-th photovoltaic array, the shadow blocking area of ​​the i-th photovoltaic array is calculated more accurately, thereby improving the calculation accuracy and efficiency, thereby solving the problem of low calculation accuracy and efficiency of the traditional photovoltaic module blocking area calculation method.

[0103] In an exemplary embodiment, the first determination module 502 is further used to determine a second matrix corresponding to the N photovoltaic arrays, wherein the dimension of the second matrix is ​​N*P*M*3, and the second matrix records the coordinates of M points corresponding to each photovoltaic array in the i-th array set, and the coordinates of the points in the second matrix are coordinates in the geodetic coordinate system; the coordinates of the points in the second matrix are converted into a coordinate system to obtain the first matrix.

[0104] In an exemplary embodiment, the first determination module 502 is also used to determine a third matrix, wherein the third matrix records the three-dimensional vectors of the spatial coordinate system of the photovoltaic array corresponding to each point in the first matrix on the X-axis, Y-axis and Z-axis, the dimension of the third matrix is ​​N*P*3*3, the X-axis of the spatial coordinate system of the photovoltaic array is parallel to the line segment where the width of the photovoltaic array is located, the Y-axis of the spatial coordinate system of the photovoltaic array is parallel to the line segment where the length of the photovoltaic array is located, and the Z-axis of the spatial coordinate system of the photovoltaic array is perpendicular to the plane where the photovoltaic array is located; based on the third matrix, the coordinates of the points in the second matrix are transformed into a coordinate system to obtain the first matrix.

[0105] In an exemplary embodiment, the first determination module 502 is further used to obtain the vertex coordinates of the two side end columns of each photovoltaic array according to the central axis of the long side of each photovoltaic array in the N photovoltaic arrays, and obtain a fourth matrix and a fifth matrix, wherein the fourth matrix photovoltaic has the vertex coordinates of the left end column of each photovoltaic array, and the fifth matrix has the vertex coordinates of the right end column of each photovoltaic array, and the dimensions of the fourth matrix and the fifth matrix are both N*3; determine a sixth matrix, wherein the sixth matrix has the deviation angle between each photovoltaic array in the N photovoltaic arrays and the corresponding end column; determine a seventh matrix, wherein the seventh matrix has the center point coordinates of each photovoltaic array in the N photovoltaic arrays, the dimension of the seventh matrix is ​​N*3, and the coordinates in the fourth matrix, the fifth matrix and the seventh matrix are coordinates in the geodetic coordinate system; determine the three-dimensional vectors of the spatial coordinate system corresponding to each photovoltaic array on the X axis, the Y axis and the Z axis according to the fourth matrix, the fifth matrix, the sixth matrix and the seventh matrix, and obtain an eighth matrix, wherein the dimension of the eighth matrix is ​​N*3*3; determine the third matrix according to the eighth matrix.

[0106] In an exemplary embodiment, the first determination module 502 is also used to determine the maximum length and the maximum width of the N photovoltaic arrays, and construct a standard rectangle according to the maximum length and the maximum width; rasterize the standard rectangle to form M points; determine M points corresponding to each of the N photovoltaic arrays, and obtain a ninth matrix according to the coordinates of the M points corresponding to each photovoltaic array, wherein the dimension of the ninth matrix is ​​N*M*3, and the coordinates of the M points corresponding to each photovoltaic array are recorded in the ninth matrix, and the coordinates of the points in the ninth matrix are coordinates in the geodetic coordinate system, and when the size of the photovoltaic array is smaller than the standard rectangle, the coordinates of the points among the M points that fall outside the photovoltaic array are determined to be (0, 0, 0); and the second matrix is ​​determined according to the ninth matrix.

[0107] In an exemplary embodiment, the fourth determining module 508 is further configured to multiply the ith number by the area of ​​each grid to obtain the shadow shielding area of ​​the ith photovoltaic array.

[0108] An embodiment of the present application further provides a computer-readable storage medium, in which a computer program is stored, wherein the computer program is configured to execute the steps of any of the above method embodiments when run.

[0109] Optionally, in this embodiment, the storage medium may be configured to store a computer program for performing the following steps:

[0110] S1, in the case where there are N photovoltaic arrays in the target area, determine a first matrix corresponding to the N photovoltaic arrays, wherein the dimension of the first matrix is ​​N*P*M*3, the first matrix records the coordinates of M points corresponding to each photovoltaic array in the i-th array set, the i-th array set includes P photovoltaic arrays adjacent to the i-th photovoltaic array in the N photovoltaic arrays, the coordinates of each point in the first matrix are the coordinates of the photovoltaic array corresponding to the point in the spatial coordinate system, i is an integer greater than or equal to 1 and less than or equal to N, P is an integer greater than or equal to 1 and less than or equal to N-1, and N and M are both integers greater than or equal to 2;

[0111] S2, determining N*P*M vectors according to the N*P*M coordinates of the N*P*M points in the first matrix, wherein the j-th vector among the N*P*M vectors is a vector determined according to the coordinates of the j-th point and the altitude angle and azimuth angle of the sun at the current moment, and j is an integer greater than or equal to 1 and less than or equal to N*P*M;

[0112] S3, determining N*P*M intersection points of the N*P*M vectors and the N planes corresponding to the N photovoltaic arrays, and determining the number of target grids in a plurality of grids corresponding to an i-th photovoltaic array among the N photovoltaic arrays, to obtain an i-th number, wherein the target grid is a grid having the intersection points;

[0113] S4, determining the shadowing area of ​​the i-th photovoltaic array according to the i-th number and the area of ​​each grid.

[0114] In an exemplary embodiment, the computer-readable storage medium may include, but is not limited to, various media that can store computer programs, such as a USB flash drive, a read-only memory (ROM), a random access memory (RAM), a mobile hard disk, a magnetic disk or an optical disk.

[0115] For specific examples in this embodiment, reference may be made to the examples described in the above embodiments and exemplary implementation modes, and this embodiment will not be described in detail herein.

[0116] An embodiment of the present application further provides a computer program product, including a computer program, and when the computer program is executed by a processor, the steps in any one of the above method embodiments are performed.

[0117] An embodiment of the present application further provides an electronic device, including a memory and a processor, wherein a computer program is stored in the memory, and the processor is configured to run the computer program to execute the steps in any one of the above method embodiments.

[0118] Optionally, in this embodiment, the processor may be configured to perform the following steps through a computer program:

[0119] S1, in the case where there are N photovoltaic arrays in the target area, determine a first matrix corresponding to the N photovoltaic arrays, wherein the dimension of the first matrix is ​​N*P*M*3, the first matrix records the coordinates of M points corresponding to each photovoltaic array in the i-th array set, the i-th array set includes P photovoltaic arrays adjacent to the i-th photovoltaic array in the N photovoltaic arrays, the coordinates of each point in the first matrix are the coordinates of the photovoltaic array corresponding to the point in the spatial coordinate system, i is an integer greater than or equal to 1 and less than or equal to N, P is an integer greater than or equal to 1 and less than or equal to N-1, and N and M are both integers greater than or equal to 2;

[0120] S2, determining N*P*M vectors according to the N*P*M coordinates of the N*P*M points in the first matrix, wherein the j-th vector among the N*P*M vectors is a vector determined according to the coordinates of the j-th point and the altitude angle and azimuth angle of the sun at the current moment, and j is an integer greater than or equal to 1 and less than or equal to N*P*M;

[0121] S3, determining N*P*M intersection points of the N*P*M vectors and the N planes corresponding to the N photovoltaic arrays, and determining the number of target grids in a plurality of grids corresponding to an i-th photovoltaic array among the N photovoltaic arrays, to obtain an i-th number, wherein the target grid is a grid having the intersection points;

[0122] S4, determining the shadowing area of ​​the i-th photovoltaic array according to the i-th number and the area of ​​each grid.

[0123] In an exemplary embodiment, the electronic device may further include a transmission device and an input / output device, wherein the transmission device is connected to the processor, and the input / output device is connected to the processor.

[0124] For specific examples in this embodiment, reference may be made to the examples described in the above embodiments and exemplary implementation modes, and this embodiment will not be described in detail herein.

[0125] Obviously, those skilled in the art should understand that the above modules or steps of the present application can be implemented by a general computing device, they can be concentrated on a single computing device, or distributed on a network composed of multiple computing devices, they can be implemented by a program code executable by a computing device, so that they can be stored in a storage device and executed by the computing device, and in some cases, the steps shown or described can be executed in a different order from that herein, or they can be made into individual integrated circuit modules, or multiple modules or steps therein can be made into a single integrated circuit module for implementation. Thus, the present application is not limited to any specific combination of hardware and software.

[0126] The above is only a preferred implementation of the present application. It should be pointed out that for ordinary technicians in this technical field, several improvements and modifications can be made without departing from the principles of the present application. These improvements and modifications should also be regarded as the scope of protection of the present application.

Claims

1. A method for determining an occlusion area, characterized in that: include: In the case where there are N photovoltaic arrays in the target area, a first matrix corresponding to the N photovoltaic arrays is determined, wherein the dimension of the first matrix is ​​N*P*M*3, the first matrix records the coordinates of M points corresponding to each photovoltaic array in the i-th array set, the i-th array set includes P photovoltaic arrays adjacent to the i-th photovoltaic array in the N photovoltaic arrays, the coordinates of each point in the first matrix are the coordinates of the photovoltaic array corresponding to the point in the spatial coordinate system, i is an integer greater than or equal to 1 and less than or equal to N, P is an integer greater than or equal to 1 and less than or equal to N-1, and N and M are both integers greater than or equal to 2; Determine N*P*M vectors according to the N*P*M coordinates of the N*P*M points in the first matrix, wherein the j-th vector among the N*P*M vectors is a vector determined according to the coordinates of the j-th point and the altitude angle and azimuth angle of the sun at the current moment, and j is an integer greater than or equal to 1 and less than or equal to N*P*M; determine N*P*M intersections of the N*P*M vectors and the N planes corresponding to the N photovoltaic arrays, and determine the number of target grids among the multiple grids corresponding to the i-th photovoltaic array among the N photovoltaic arrays to obtain the i-th number, wherein the target grid is a grid having the intersection; determine the shadow shielding area of ​​the i-th photovoltaic array according to the i-th number and the area of ​​each grid.

2. The method according to claim 1, characterized in that Determining a first matrix corresponding to the N photovoltaic arrays includes: Determine a second matrix corresponding to the N photovoltaic arrays, wherein the dimension of the second matrix is ​​N*P*M*3, the second matrix records the coordinates of M points corresponding to each photovoltaic array in the i-th array set, and the coordinates of the points in the second matrix are coordinates in a geodetic coordinate system; The coordinates of the points in the second matrix are transformed into a coordinate system to obtain the first matrix.

3. The method according to claim 2, characterized in that Performing coordinate system transformation on the coordinates of the points in the second matrix to obtain the first matrix includes: Determine a third matrix, wherein the third matrix records the three-dimensional vectors of the spatial coordinate system of the photovoltaic array corresponding to each point in the first matrix on the X-axis, the Y-axis and the Z-axis, the dimension of the third matrix is ​​N*P*3*3, the X-axis of the spatial coordinate system of the photovoltaic array is parallel to the line segment where the width of the photovoltaic array is located, the Y-axis of the spatial coordinate system of the photovoltaic array is parallel to the line segment where the length of the photovoltaic array is located, and the Z-axis of the spatial coordinate system of the photovoltaic array is perpendicular to the plane where the photovoltaic array is located; based on the third matrix, the coordinate system of the coordinates of the points in the second matrix is ​​transformed to obtain the first matrix.

4. The method according to claim 3, characterized in that Determine the third matrix, including: According to the central axis of the long side of each photovoltaic array in the N photovoltaic arrays, the vertex coordinates of the two side end columns of each photovoltaic array are obtained to obtain a fourth matrix and a fifth matrix, wherein the fourth matrix photovoltaic has the vertex coordinates of the left end column of each photovoltaic array, and the fifth matrix has the vertex coordinates of the right end column of each photovoltaic array, and the dimensions of the fourth matrix and the fifth matrix are both N*3; Determine a sixth matrix, wherein the sixth matrix includes a deviation angle between each photovoltaic array of the N photovoltaic arrays and a corresponding end pole; Determine a seventh matrix, wherein the seventh matrix has the coordinates of the center point of each of the N photovoltaic arrays, the dimension of the seventh matrix is ​​N*3, and the coordinates in the fourth matrix, the fifth matrix and the seventh matrix are coordinates in a geodetic coordinate system; Determine the three-dimensional vectors of the spatial coordinate system corresponding to each photovoltaic array on the X-axis, the Y-axis and the Z-axis according to the fourth matrix, the fifth matrix, the sixth matrix and the seventh matrix to obtain an eighth matrix, wherein the dimension of the eighth matrix is ​​N*3*3; The third matrix is ​​determined according to the eighth matrix.

5. The method according to claim 2, characterized in that: Determining a second matrix corresponding to the N photovoltaic arrays includes: Determining the maximum length and the maximum width of the N photovoltaic arrays, and constructing a standard rectangle according to the maximum length and the maximum width; Rasterizing the standard rectangle to form M points; Determine M points corresponding to each of the N photovoltaic arrays, and obtain a ninth matrix according to the coordinates of the M points corresponding to each photovoltaic array, wherein the dimension of the ninth matrix is ​​N*M*3, and the coordinates of the M points corresponding to each photovoltaic array are recorded in the ninth matrix. The coordinates of the points in the ninth matrix are coordinates in the geodetic coordinate system. When the size of the photovoltaic array is smaller than the standard rectangle, the coordinates of the points outside the photovoltaic array among the M points are determined to be (0, 0, 0); The second matrix is ​​determined according to the ninth matrix.

6. The method according to claim 1, characterized in that Determining the shadow shielding area of ​​the i-th photovoltaic array according to the i-th number and the area of ​​each grid includes: The shadow shielding area of ​​the i-th photovoltaic array is obtained by multiplying the i-th number by the area of ​​each grid.

7. A device for determining an obstruction area, characterized in that: include: A first determination module is used to determine a first matrix corresponding to the N photovoltaic arrays when there are N photovoltaic arrays in the target area, wherein the dimension of the first matrix is ​​N*P*M*3, the first matrix records the coordinates of M points corresponding to each photovoltaic array in the i-th array set, the i-th array set includes P photovoltaic arrays adjacent to the i-th photovoltaic array in the N photovoltaic arrays, the coordinates of each point in the first matrix are the coordinates of the photovoltaic array corresponding to the point in the space coordinate system, i is an integer greater than or equal to 1 and less than or equal to N, P is an integer greater than or equal to 1 and less than or equal to N-1, and N and M are both integers greater than or equal to 2; A second determination module is used to determine N*P*M vectors according to the N*P*M coordinates of the N*P*M points in the first matrix, wherein the j-th vector among the N*P*M vectors is a vector determined according to the coordinates of the j-th point and the altitude angle and azimuth angle of the sun at the current moment, and j is an integer greater than or equal to 1 and less than or equal to N*P*M; A third determination module is used to determine N*P*M intersections of the N*P*M vectors and the N planes corresponding to the N photovoltaic arrays, and determine the number of target grids in a plurality of grids corresponding to an i-th photovoltaic array in the N photovoltaic arrays to obtain an i-th number, wherein the target grid is a grid having the intersection; The fourth determining module is used to determine the shadow shielding area of ​​the i-th photovoltaic array according to the i-th number and the area of ​​each grid.

8. A computer-readable storage medium, characterized in that: The computer-readable storage medium includes a stored program, wherein the program executes the method according to any one of claims 1 to 6 when executed.

9. An electronic device comprising a memory and a processor, characterized in that: A computer program is stored in the memory, and the processor is configured to execute the method according to any one of claims 1 to 6 through the computer program.

10. A computer program product, comprising a computer program, characterized in that When the computer program is executed by a processor, the method according to any one of claims 1 to 6 is implemented.