Spherical coordinate-based sunshine duration calculation method and system for building regeneration

Abstract

The invention discloses a spherical coordinate-based sunshine duration calculation method and system for building regeneration, and the method comprises the steps: taking a movement track of the sun in one day as a bottom surface, taking a measurement point as a vertex, obtaining a sunshine conical surface under a spherical coordinate system, and obtaining sunshine conical surface parameters, including a cone vector, a cone angle and a sunshine analysis time interval of the conical vertex; performing triangular meshing on the building model, and converting vertexes of a triangular mesh into spherical coordinates which are in the same spherical coordinate system with the sunlight conical surface; finding out a triangular mesh colliding with the sunlight conical surface; calculating the longitude of the intersection point of the triangular mesh and the sunlight conical surface to obtain a shielding interval of the triangular mesh; carrying out interval combination on the shielding intervals of all the triangular meshes to obtain a shielding interval of the building model; and subtracting the shielding interval of the building model from the sunshine analysis time interval to obtain the sunshine interval of the measuring point, and calculating the sunshine duration. The method solves the problems that the calculation precision and the calculation speed are difficult to consider at the same time, and the calculation difficulty exists for complex buildings.

Description

Technical Field

[0001] The present invention relates to the technical field of computer-aided design, and in particular to a method and system for calculating sunshine duration for building regeneration based on spherical coordinates. Background Art

[0002] A large number of CAD architectural drawing software, intelligent planning software, etc. need to calculate the sunshine duration of buildings. For example, the design model of the building after demolition and modification needs to calculate the sunshine duration, so as to design the distance and shape between multiple buildings more appropriately, obtain better lighting conditions, and avoid excessive light blocking that affects the living experience of residents. Most of the existing sunshine calculation software uses the method of sampling and subdividing time to calculate the sunshine duration. The calculation accuracy depends on the sampling density. The problem caused by this is that the sunshine accuracy and calculation speed are highly related. According to GB / T50947-2014 "Standard for Building Sunlight Calculation Parameters", the sunshine sampling interval must be less than 1min, and the overall error must be less than 3min; however, the accuracy of the sunshine calculation software based on sampling is limited, and if the calculation accuracy is improved, the cost of computing resources is too high, which will greatly increase the calculation time, and it is difficult to balance the calculation accuracy and speed. The existing sunshine calculation software based on the occlusion interval can effectively speed up the calculation speed, but the existing algorithm cannot perform sunshine calculation for buildings with complex geometric shapes. Summary of the invention

[0003] 1. Technical issues to be resolved Based on the above problems, the present invention provides a sunshine duration calculation method and system for building regeneration based on spherical coordinates, which solves the problem that the existing technology relies on sunshine sampling density, it is difficult to balance calculation accuracy and calculation speed, and there is difficulty in calculation for complex buildings.

[0004] (II) Technical solution Based on the above technical problems, the present invention provides a method for calculating sunshine duration of a building model based on spherical coordinate transformation, comprising: S1. Taking the moving track of the sun in one day as the base and the measuring point as the vertex, a sunshine cone in the spherical coordinate system is obtained, and sunshine cone parameters are obtained according to the sunshine analysis date and time period. The sunshine cone parameters include the cone direction vector, cone angle, and sunshine analysis time interval of the cone vertex; S2, triangulating the building model, and converting the spatial rectangular coordinates of the vertices of the triangulated mesh into spherical coordinates in the same spherical coordinate system as the sunlight cone surface according to the cone direction vector of the cone vertex; S3, finding the triangular mesh that collides with the sunlight cone according to the cone angle; S4, calculating the longitudes of two intersection points between the triangular mesh colliding with the sunlight cone and the sunlight cone, and the range of the longitudes of the two intersection points is the occlusion interval of the triangular mesh; S5, performing interval merging Boolean operation on the occlusion intervals corresponding to all the triangular meshes that collide with the sunlight cone to obtain the occlusion intervals of the building model; S6. Subtract the shielding interval of the building model from the sunshine analysis period to obtain the sunshine interval of the measuring point, and calculate the sunshine duration according to the sunshine interval of the measuring point.

[0005] Furthermore, S3 includes: determining whether the azimuth coordinates of the three vertices of the triangular mesh are all greater than the cone angle or all less than the cone angle; if so, there is no collision between the triangular mesh and the sunlight cone; otherwise, there is a collision.

[0006] Furthermore, S4 includes: if the edges intersecting with the sunlight cone are AB and AC, the spherical coordinates of vertex A are (r1, θ1, φ1), the spherical coordinates of vertex B are (r2, θ2, φ2), the spherical coordinates of vertex C are (r3, θ3, φ3), and α is the cone angle, then the longitude of the intersection point D of edge AB and the sunlight cone is d= (α-θ2)*(φ1-φ2) / (θ1-θ2)+φ2, and the longitude of the intersection point E of edge AC and the sunlight cone is e= (α-θ3)*(φ1-φ3) / (θ1-θ3)+φ3. Assuming e>d, the occlusion interval of the triangular mesh ABC is (d, e).

[0007] Furthermore, in S5, the interval merging Boolean operation is to obtain a union.

[0008] Further, in S6, calculating the sunshine duration according to the sunshine interval of the measuring point includes: obtaining the sunshine duration according to the difference between the end values ​​of the sunshine interval of the measuring point.

[0009] The present invention also discloses a sunshine duration calculation system for building regeneration based on spherical coordinates, comprising: at least one processor; and at least one memory in communication with the processor, wherein: The memory stores program instructions that can be executed by the processor, and the processor can execute the method by calling the program instructions.

[0010] The present invention also discloses a non-transitory computer-readable storage medium, wherein the non-transitory computer-readable storage medium stores computer instructions, wherein the computer instructions enable the computer to execute the method.

[0011] (III) Beneficial effects The above technical solution of the present invention has the following advantages: (1) The present invention converts the sunshine trajectory into a sunshine cone in a spherical coordinate system, obtains the occlusion interval of the triangular mesh by calculating the longitude of the intersection of the sunshine cone and each colliding triangular mesh of the building model, and then obtains the occlusion interval of the building model through interval merging Boolean operation, that is, the sunshine interval and sunshine duration of the measuring point. It does not need to rely on the time sampling density, and can still obtain millisecond-level sunshine duration accuracy. Therefore, the present invention reduces the difficulty of sunshine duration calculation while greatly improving the calculation accuracy and speed; (2) In the sunshine duration calculation method of the present invention, the method of converting the building model into a triangular mesh is extremely mature and applicable to various building models including complex buildings. The sunshine duration calculation method can be used to calculate the sunshine duration of building models on all geometric building software, which reduces the difficulty of calculation and has excellent versatility. BRIEF DESCRIPTION OF THE DRAWINGS

[0012] The features and advantages of the present invention will be more clearly understood by referring to the accompanying drawings, which are schematic and should not be construed as limiting the present invention in any way. In the accompanying drawings: Figure 1 A schematic flow chart of a method for calculating sunshine duration for building regeneration based on spherical coordinates according to an embodiment of the present invention; Figure 2 is a schematic diagram of a sunlight cone according to an embodiment of the present invention; Figure 3 A schematic diagram of a collision between a building model and a sunlight cone according to an embodiment of the present invention; Figure 4 It is a schematic diagram of the principle of the occlusion interval of the triangular mesh colliding with the sunlight cone according to an embodiment of the present invention. DETAILED DESCRIPTION

[0013] The specific implementation of the present invention is further described in detail below in conjunction with the accompanying drawings and examples. The following examples are used to illustrate the present invention, but are not intended to limit the scope of the present invention.

[0014] The embodiment of the present invention is a method for calculating sunshine duration for building regeneration based on spherical coordinates. Figure 1 As shown, the following steps are included: S1. Taking the moving track of the sun in one day as the base and the measuring point as the vertex, a sunshine cone in the spherical coordinate system is obtained, and sunshine cone parameters are obtained according to the sunshine analysis date and time period. The sunshine cone parameters include the cone direction vector, cone angle, and sunshine analysis time interval of the cone vertex; According to the laws of solar motion and the specifications of sunshine calculation parameters, the sun's movement trajectory in a day can be approximately considered as a circle on the celestial sphere with the measuring point as the center of the sphere. The trajectory and the center of the circle can form a cone, that is, a sunshine cone with the sun's movement trajectory in a day as the base and the measuring point as the vertex.

[0015] like Figure 2 As shown in the figure, the sunshine cone parameters can be obtained according to the sunshine analysis date and analysis period. The sunshine cone parameters include the cone direction vector of the cone vertex: =( v x , v y , v z ), the cone angle is α, and the sunshine analysis period interval is (T1, T2), where T1 and T2 are the start time and end time of the sunshine analysis period respectively, T1 and T2 are both true solar time, the cone angle α also represents the latitude of the sunshine cone, and the cone direction vector is based on the spatial rectangular coordinate system, pointing from the cone vertex to the center of the circle.

[0016] S2, triangulating the building model, and converting the spatial rectangular coordinates of the vertices of the triangulated mesh into spherical coordinates in the same spherical coordinate system as the sunlight cone surface according to the cone direction vector of the cone vertex; According to the cone vertex cone vector =( v x , v y , v z ) is used to perform spherical coordinate conversion, thereby converting the spatial rectangular coordinates of the vertices of the triangular mesh into spherical coordinates based on the spherical coordinate center of the cone vertex, so that the sunlight cone surface and the building model are both in the same spherical coordinates.

[0017] S3. Find the triangular mesh that collides with the sunlight cone surface according to the cone angle: determine whether the azimuth coordinates of the three vertices of the triangular mesh are all greater than the cone angle or all less than the cone angle. If so, there is no collision between the triangular mesh and the sunlight cone surface, otherwise, there is a collision.

[0018] like Figure 3 As shown, even a complex building model is composed of multiple triangular meshes, and the occlusion interval problem of the complex building model can be converted into the occlusion interval problem of multiple triangular meshes.

[0019] When there is no collision between the triangular mesh and the sunlight cone, if the triangular mesh is above the sunlight cone, the azimuth angles of the three vertices of the triangular mesh are all greater than the cone angle of the sunlight cone; if the triangular mesh is below the sunlight cone, the azimuth angles of the three vertices of the triangular mesh are all less than the cone angle of the sunlight cone; when there is a collision between the triangular mesh and the sunlight cone, the three vertices of the triangular mesh must not be on the same side of the sunlight cone, then the azimuth angle of one or two vertices must be greater than the cone angle of the sunlight cone, that is, the azimuth coordinates of the three vertices of the triangular mesh do not satisfy the requirement of being greater than the cone angle or less than the cone angle.

[0020] S4, calculating the longitudes of two intersection points between the triangular mesh colliding with the sunlight cone and the sunlight cone, and the range of the longitudes of the two intersection points is the occlusion interval of the triangular mesh; When a triangle mesh collides with a sunlight cone, the three vertices of the triangle mesh must not be on the same side of the sunlight cone. Then the triangle mesh with collision must have two edges intersecting with the sunlight cone. The interval obtained by the longitude of the two intersection points of the two edges of the triangle mesh with collision and the sunlight cone is the collision interval of the triangle mesh and the sunlight cone, that is, the occlusion interval of the triangle mesh. The method for solving the longitude of the two intersection points of the two edges of the triangle mesh with collision and the sunlight cone is as follows: like Figure 4 As shown, if the edges intersecting with the sunshine cone are AB and AC, the spherical coordinates of vertex A are (r1, θ1, φ1), the spherical coordinates of vertex B are (r2, θ2, φ2), and the spherical coordinates of vertex C are (r3, θ3, φ3), θ1, θ2, and θ3 are the latitudes of the corresponding vertices, φ1, φ2, and φ3 are the longitudes of the corresponding vertices, r1, r2, and r3 are the altitudes of the corresponding vertices, but they will not be used when solving the longitude. The intersection points D and E of the edges AB and AC intersecting with the sunshine cone are on the sunshine cone, and the latitudes of the intersection points D and E are consistent with the latitude of the sunshine cone, both of which are cone angles α. Then the longitude of the intersection point D can be calculated as d= (α-θ2)*(φ1-φ2) / (θ1-θ2)+φ2, and the longitude of the intersection point E is e= (α-θ3)*(φ1-φ3) / (θ1-θ3)+φ3, assuming e>d, the occlusion interval of the triangle mesh ABC is (d, e).

[0021] S5, performing interval merging Boolean operation on the occlusion intervals corresponding to all the triangular meshes that collide with the sunlight cone to obtain the occlusion intervals of the building model; From S4, the occlusion interval list of all triangular meshes can be obtained. All triangular meshes constitute the building model. Therefore, the occlusion intervals of all triangular meshes are subjected to interval merging Boolean operation, that is, the union is calculated to obtain the occlusion intervals of the building model.

[0022] S6, subtracting the shielding interval of the building model from the sunshine analysis period to obtain the sunshine interval of the measuring point, and calculating the sunshine duration according to the sunshine interval of the measuring point; The sunshine analysis period interval (T1, T2) is subtracted from the shading interval of the building model to obtain the sunshine interval of the measuring point, and the sunshine duration is obtained according to the difference between the end values ​​of the sunshine interval of the measuring point, thereby obtaining the accurate sunshine interval and sunshine duration of the measuring point.

[0023] Finally, it should be noted that the above-mentioned calculation method can be converted into software program instructions, which can be implemented by using a control system including a processor and a memory, or by computer instructions stored in a non-transitory computer-readable storage medium. The above-mentioned integrated unit implemented in the form of a software functional unit can be stored in a computer-readable storage medium. The above-mentioned software functional unit is stored in a storage medium, including a number of instructions to enable a computer device (which can be a personal computer, a server, or a network device, etc.) or a processor to perform some steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (Read-Only Memory, ROM), random access memory (Random Access Memory, RAM), disk or optical disk and other media that can store program codes.

[0024] In summary, the above-mentioned sunshine duration calculation method and system for building regeneration based on spherical coordinates have the following beneficial effects: (1) The present invention converts the sunshine trajectory into a sunshine cone in a spherical coordinate system, obtains the occlusion interval of the triangular mesh by calculating the longitude of the intersection of the sunshine cone and each colliding triangular mesh of the building model, and then obtains the occlusion interval of the building model through interval merging Boolean operation, that is, the sunshine interval and sunshine duration of the measuring point. It does not need to rely on the time sampling density, and can still obtain millisecond-level sunshine duration accuracy. Therefore, the present invention reduces the difficulty of sunshine duration calculation while greatly improving the calculation accuracy and speed; (2) In the sunshine duration calculation method of the present invention, the method of converting the building model into a triangular mesh is extremely mature and applicable to various building models including complex buildings. The sunshine duration calculation method can be used to calculate the sunshine duration of building models on all geometric building software, which reduces the difficulty of calculation and has excellent versatility.

[0025] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the embodiments of the present invention are described in conjunction with the accompanying drawings, those skilled in the art may make various modifications and variations without departing from the spirit and scope of the present invention, and such modifications and variations shall fall within the scope defined by the appended claims.

Claims

1. A method for calculating sunshine duration for building regeneration based on spherical coordinates, characterized in that: include: S1. Taking the moving track of the sun in one day as the base and the measuring point as the vertex, a sunshine cone in the spherical coordinate system is obtained, and sunshine cone parameters are obtained according to the sunshine analysis date and time period. The sunshine cone parameters include the cone direction vector, cone angle, and sunshine analysis time interval of the cone vertex; S2, triangulating the building model, and converting the spatial rectangular coordinates of the vertices of the triangulated mesh into spherical coordinates in the same spherical coordinate system as the sunlight cone surface according to the cone direction vector of the cone vertex; S3, finding the triangular mesh that collides with the sunlight cone according to the cone angle; S4, calculating the longitudes of two intersection points between the triangular mesh colliding with the sunlight cone and the sunlight cone, and the range of the longitudes of the two intersection points is the occlusion interval of the triangular mesh; S5, performing interval merging Boolean operation on the occlusion intervals corresponding to all the triangular meshes that collide with the sunlight cone to obtain the occlusion intervals of the building model; S6. Subtract the shielding interval of the building model from the sunshine analysis period to obtain the sunshine interval of the measuring point, and calculate the sunshine duration according to the sunshine interval of the measuring point.

2. The sunshine duration calculation method for building regeneration based on spherical coordinates according to claim 1 is characterized in that: The S3 includes: determining whether the azimuth coordinates of the three vertices of the triangular mesh are all greater than the cone angle or all less than the cone angle; if so, there is no collision between the triangular mesh and the sunlight cone surface; otherwise, there is a collision.

3. The sunshine duration calculation method for building regeneration based on spherical coordinates according to claim 1 is characterized in that: The S4 includes: if the edges intersecting with the sunlight cone are AB and AC, the spherical coordinates of vertex A are (r1, θ1, φ1), the spherical coordinates of vertex B are (r2, θ2, φ2), the spherical coordinates of vertex C are (r3, θ3, φ3), α is the cone angle, then the longitude of the intersection point D of the edge AB and the sunlight cone is d= (α-θ2)*(φ1-φ2) / (θ1-θ2)+φ2, the longitude of the intersection point E of the edge AC and the sunlight cone is e=(α-θ3)*(φ1-φ3) / (θ1-θ3)+φ3, assuming e>d, the occlusion interval of the triangular mesh ABC is (d, e).

4. The sunshine duration calculation method for building regeneration based on spherical coordinates according to claim 1 is characterized in that: In S5, the interval merging Boolean operation is to obtain a union.

5. The sunshine duration calculation method for building regeneration based on spherical coordinates according to claim 1 is characterized in that: In S6, the calculating of the sunshine duration according to the sunshine interval of the measuring point includes: obtaining the sunshine duration according to the difference between the end values ​​of the sunshine interval of the measuring point.

6. A sunshine duration calculation system for building regeneration based on spherical coordinates, characterized in that: include: at least one processor; and at least one memory in communication with the processor, wherein: The memory stores program instructions executable by the processor, and the processor can execute the method according to any one of claims 1 to 5 by calling the program instructions.

7. A non-transitory computer-readable storage medium, characterized in that: The non-transitory computer-readable storage medium stores computer instructions, which cause the computer to execute the method according to any one of claims 1 to 5.

Images