CN113820756B - Real-time sunlight detection method and system - Google Patents

Real-time sunlight detection method and system Download PDF

Info

Publication number
CN113820756B
CN113820756B CN202111336167.1A CN202111336167A CN113820756B CN 113820756 B CN113820756 B CN 113820756B CN 202111336167 A CN202111336167 A CN 202111336167A CN 113820756 B CN113820756 B CN 113820756B
Authority
CN
China
Prior art keywords
sunlight
target building
measuring point
point
time
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202111336167.1A
Other languages
Chinese (zh)
Other versions
CN113820756A (en
Inventor
倪钰翔
李星星
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Lufasibang Technology Co ltd
Original Assignee
Nanjing Lufasibang Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing Lufasibang Technology Co ltd filed Critical Nanjing Lufasibang Technology Co ltd
Publication of CN113820756A publication Critical patent/CN113820756A/en
Application granted granted Critical
Publication of CN113820756B publication Critical patent/CN113820756B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01WMETEOROLOGY
    • G01W1/00Meteorology
    • G01W1/12Sunshine duration recorders
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Environmental & Geological Engineering (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Atmospheric Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Ecology (AREA)
  • Environmental Sciences (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)

Abstract

The invention provides a real-time sunlight detection method and a real-time sunlight detection system. The invention firstly calculates the sunshine condition of the target building under each time interval and records the sunshine condition in the form of a bar shadow ray set. And then, according to the sunlight measuring points of the building, respectively calling the rod shadow ray set to the positions of each sunlight measuring point, and according to whether the rod shadow ray set rays are blocked by obstacles, marking the sunlight time. And finally, counting the effective sunlight time of each sunlight measuring point according to a specified calculation mode to obtain the sunlight duration of each sunlight measuring point. The light rays in the rod shadow light ray concentration are matched with the sunlight azimuth of a target building, different sunlight time gaps are marked, real-time illumination conditions of sunlight measuring points under different sunlight time gaps can be obtained directly through whether the light rays are shielded by obstacles, and statistics of effective sunlight time length is accurately achieved flexibly according to the requirement of sunlight time length calculation modes.

Description

Real-time sunlight detection method and system
Technical Field
The invention relates to the technical field of building sunshine detection, in particular to a real-time sunshine detection method and a real-time sunshine detection system.
Background
According to the requirements of urban living area planning and design specifications, parameters such as sunshine standard day, effective sunshine time zone, sunshine time number and the like of a building are required to meet mandatory treaty index requirements. And the sunshine of the building and the site is calculated and obtained according to the standard specification, and whether the sunshine standard of the building is met or not can be measured and calculated.
The existing method for measuring and calculating the sunlight of the building comprises the following steps:
1. sun stick photographs. According to the incident height and the incident direction of solar rays, the shadow length and the shadow direction of a stick are calculated to simultaneously draw the stick shadow in an effective time period, and the shadow length and the shadow direction are used for representing the sunlight influence range of a measuring point, the sunlight duration of the measuring point and the distance and the direction between the measuring point and a surrounding building. However, the sun stick graph identifies a building sun condition with the following disadvantages: firstly, the sun shine stick shadow map has no clear sun shine time length and shielding relation, and is in a sector net shape when used for judging the distance and the direction, so that the relation between the measuring point and the surrounding environment cannot be clearly identified; secondly, the calculation mode of the sunshine duration is single, only all sunshine time can be accumulated, and the requirement of building sunshine standard cannot be met; in addition, in the current CAD plug-in, the measuring points are required to be set manually and specially to display the sunlight stick-figure, and the stick-figure cannot be automatically adjusted along with the adjustment of the three-dimensional building and is required to be set manually again, so that dynamic display cannot be realized.
2. Graph analysis along solar rays. The method comprises the steps of setting sampling precision at a specified elevation of a building contour line on a two-dimensional plane, obtaining a sunlight measuring point to be calculated, calculating sunlight duration of the measuring point, and marking the sunlight duration after calculation at a base point position. But this technique has the following disadvantages when used to view the duration of sunlight on the building contour: firstly, a plurality of unnecessary sunlight measuring points on a north vertical face are calculated, and for a residential building, the necessary sunlight measuring points are only points with specified elevation on the south vertical face of the building, and the algorithm calculates a plurality of unnecessary sunlight points, so that computer resources are seriously wasted; secondly, the sunlight duration and the shielding relation cannot be clearly displayed, and the sunlight shielding relation with surrounding buildings cannot be visually displayed; in addition, the calculation mode of the sunshine duration is single, the sunshine duration cannot be accumulated according to the existing building sunshine specification, and only all sunshine time can be accumulated; furthermore, the sunshine duration is only displayed on the two-dimensional plan, so that the sunshine duration is not easy to view on the three-dimensional building model, and the visual readability is poor; and after the position of the three-dimensional building model is adjusted in the CAD plug-in, the method can not update the sunlight state along with the change of the model position, and needs to manually set and re-consume a large amount of calculation to repeatedly calculate from far.
3. Sunlight multi-point diagram. On a two-dimensional plane, a certain range is selected in a frame mode, the distribution interval and elevation of sunlight measuring points are set, calculation of sunlight duration is conducted on each measuring point, and then calculated sunlight duration numbers and corresponding colors are displayed. But this technique has the following drawbacks when used to view the sun environment between multiple buildings: firstly, calculating a plurality of unnecessary sunlight measuring points positioned on a northbound elevation of a building, wasting computer resources and consuming calculation time; secondly, the calculation mode of the sunshine duration is single, the sunshine duration cannot be accumulated according to the existing building sunshine specification, and only all sunshine time can be accumulated; in addition, the sunlight numerical value has poor visual readability, the sunlight duration can only be displayed in a two-dimensional plan view, and the sunlight duration is not easy to view on a three-dimensional building model, and the display of the sunlight duration is only distinguished by the numerical value and the color, so that the duration of a specific position is not easy to view; meanwhile, the sun shine can not be dynamically displayed, and the sun shine time length parameter can not be automatically updated in direct response to the adjustment of the position of the three-dimensional building model.
4. Solar conic graph. Setting a plurality of moments in a three-dimensional space, taking a datum line and a datum point, displaying the sun position at each moment in the space around the datum point, and connecting the datum point and each sun position by using the datum line to form a conical sector. This technique has the following disadvantages when looking at the light occlusion of the datum point and surrounding building space: firstly, the relationship between sunlight duration and shielding cannot be clearly represented, and the relationship between shielding and spacing between all buildings and measuring points in a conical surface cannot be clearly judged; secondly, the calculation mode of the sunshine duration is single, the sunshine duration cannot be accumulated according to the existing building sunshine specification, and only all sunshine time can be accumulated; meanwhile, the obtained sunlight value is poor in visual readability, sunlight duration at the measuring point is not obvious in display, and light rays which are not shielded temporarily and light rays which are shielded temporarily are easy to interfere visual judgment during display; in addition, dynamic sunlight display cannot be performed, in the plug-in unit of the CAD at present, a measuring point is required to be set specially by people to display the sunlight cone map, and the distance between the detecting point and the shielding building is required to be detected and adjusted one by people.
Therefore, how to provide a reasonable, efficient and accurate scheme for real-time measurement and calculation of the sunlight condition of a building is a technical problem to be solved by those skilled in the art.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides a real-time sunshine detection method and a real-time sunshine detection system, which realize calculation of sunshine shielding by establishing a rod shadow ray set and selecting necessary sunshine measuring points. The invention adopts the following technical scheme.
Firstly, in order to achieve the above objective, a real-time sunlight detection method is provided, which includes the steps of: according to the longitude and latitude of the target building and the height thereof, respectively calculating a rod shadow ray set under the sunlight condition corresponding to the target building in each time interval according to a preset time interval; secondly, respectively determining corresponding sunlight measuring points for each target building, respectively starting a sunlight shielding algorithm thread for each target building, and respectively executing the following third to seventh calculation steps for each sunlight measuring point of the target building through the sunlight shielding algorithm thread; thirdly, comparing the two-dimensional projection range of the rod shadow ray set corresponding to the target building with the longitude and latitude positions of surrounding buildings, if the surrounding buildings are in the two-dimensional projection range of the rod shadow ray set, recording the surrounding buildings as obstacles and jumping to the fourth step, otherwise, jumping to the eighth step; fourth, sorting the barriers according to the distance between the barrier and the target building; fifthly, calling the rod shadow ray set to each sunlight measuring and calculating point of the target building, and respectively marking serial numbers for each ray in the rod shadow ray set according to the sequence of the sun azimuth angle; step six, respectively establishing 3 memories for each sunlight measuring point of the target building, storing the rod shadow ray set in the step five and the serial numbers marked by each ray in the rod shadow ray set by using a first memory, storing the serial numbers of the rays shielded by the barrier by using a second memory, and storing the serial numbers of the rays not shielded by the barrier by using a third memory; seventh, taking each sunlight measuring point as a base point, and sequentially judging each obstacle: acquiring an azimuth angle range between the boundary of the obstacle and the base point, calling all the light rays positioned in the azimuth angle range from a first memory, judging whether the light rays obtained by calling intersect the obstacle, if so, moving the light rays into a second memory, and if not, moving the light rays into a third memory; and eighth, respectively calling the light rays in the third memory corresponding to each sunlight measuring point according to a specified sunlight duration calculation mode, and calculating whether the sunlight measuring point reaches the sunlight duration standard.
Optionally, the real-time insolation detection method according to any one of the above, further comprising: ninth, if the sunlight measuring point reaches the sunlight duration standard, marking the corresponding sunlight duration at the position of the sunlight measuring point in a first display mode; otherwise, marking the corresponding sunshine duration at the position of the sunshine measuring point in a second display mode; the marking fonts, the font sizes and the thicknesses of the sunshine duration are all larger than those of the first display mode in the second display mode, and the marking colors of the sunshine duration are different from those of the first display mode in the second display mode.
Optionally, the real-time insolation detection method according to any one of the above, further comprising: tenth, every time any one of the target building or surrounding building is moved or changed, the second to ninth steps are automatically triggered to be re-executed.
Optionally, the real-time sunlight detection method according to any one of the preceding claims, wherein the rod shadow ray set is obtained by calculating: step 101, according to latitude and longitude (lambda, phi) of a target building, declination angle delta and true solar angle omega under a corresponding time gap, and sinh s Calculating to obtain a solar altitude angle hs corresponding to a target building according to the formula of sin phi sin delta+cos phi cos delta cos omega Calculating to obtain a solar azimuth angle As corresponding to the target building;step 102, calculating a bar-shadow bus terminal coordinate point P (Xp, yp, zp) matched with the maximum height H of a surrounding building, wherein xp=H×sinAs×coths, yp= -1×H×cosas×coths, zp=H; step 103, connecting the coordinate origin O (0, 0) with the terminal coordinate point P of the bar shadow bus to obtain the bar shadow light OP under the sunshine condition corresponding to the target building under the corresponding time gap, wherein the two-dimensional projection of the bar shadow light OP is OP 0 Wherein P is 0 Coordinates (Xp, yp, 0); step 104, executing the steps 101 to 102 respectively according to the predetermined time intervals, and storing the set of rod shadow rays OP under the sunlight condition corresponding to the target building under each time interval as a rod shadow ray set.
Optionally, in any one of the above real-time solar radiation detection methods, in the fifth step, the specific step of calling the rod shadow ray set to each solar radiation measuring point of the target building is: and respectively moving the coordinate origin of each rod shadow ray OP in the rod shadow ray set to each sunlight measuring point of the target building to obtain a rod shadow ray O 'P' matched with the sunlight measuring point, wherein the P 'coordinate is (Xp+Xo, yp+yo, zp+zo) and the sunlight measuring point coordinate is O' (Xo, yo, zo).
Optionally, the real-time sunlight detection method according to any one of the above, wherein the sunlight measuring point of the target building is determined by: step 201, taking the standard floor plan established by using the multi-section lines in anticlockwise sequence as a reference, obtaining a standard floor plan of the target building, extracting the coordinates of each endpoint in the standard floor plan, recording the endpoint at the southwest end as Pt (Xpt, ypt, zpt), recording the endpoint at the westst end as Pt Xmin (Xmin, X0, Z0) in which the end point Pt located at the eastern end is recorded Xmax (Xmax, X1, Z1); step 202, if the number of endpoints at the south end exceeds 1, determining that the deflection angle of the target building is 0, otherwise, determining that the target building deflects; step 203, when it is determined that there is deflection in the target building, extracting an edge line segment L where the start point coordinate value Pts (Xstart, ystart, zstart) and the end point Pt (Xpt, ypt, zpt) located at the south end overlap in the standard floor plan view, and recording the end point of the edge line segment LCoordinates are Pte (Xend, yend, zend), and a vector V (Xend-Xstart, yend-ystart, zend-Zstart) is established according to the direction of the edge line segment L; step 204, calculating the included angle between the vector V and the axial quantity Vx (1, 0) to obtain a deflection value RotAngle corresponding to the target building; step 205, calculating a deflection angle corresponding to the target building according to the deflection value RotAngle: when (when) When the bias Angle is determined to be angle=rotangle-90 °, when +.>Determining that the declination Angle is angle=rotangle; step 206, taking the center of the standard floor type plane graph as a rotation point, deflecting the standard floor type plane graph by an angle of-1 x, and under the state, acquiring each extreme value coordinate point in the rotated standard floor type plane graph as follows: x ' max, X ' min, Y ' max, Y ' min, Z '; step 207, constructing an external rectangle of the rotated standard layer house type plane graph according to each extreme value coordinate point, and marking corresponding virtual measuring points on a southbound edge line of the external rectangle according to the number of measuring points or the interval requirement; and step 208, taking the virtual measuring points as the bottom ends, making a vertical line along a standard house type plane of the target building, marking corresponding coordinates on the vertical line as initial measuring points according to the height requirement of the measuring points or house type floor height, and deflecting each virtual measuring point and each initial measuring point back to an Angle by taking the center of a standard house type plane graph as a rotation center to obtain a sunlight measuring point corresponding to the south of the target building.
Meanwhile, in order to achieve the above object, the present invention also provides a real-time insolation detection system, which includes: the rod shadow ray set generation module is used for respectively calculating rod shadow ray sets under the corresponding sunlight conditions of the target building in each time interval according to the longitude and latitude and the height of the target building and the preset time interval; the sunlight measuring point generating module is used for determining corresponding sunlight measuring points for each target building respectively; the sunlight shielding algorithm module is used for respectively starting a sunlight shielding algorithm thread for each target building, and judging whether the concentrated light of the rod shadow rays is shielded or not for each sunlight measuring point of the target building through the sunlight shielding algorithm thread so as to obtain the light shielding condition of each sunlight measuring point; the sunlight duration calculation module is used for respectively calculating whether the non-shielded light rays corresponding to each sunlight measuring point reach the sunlight duration standard according to the light shielding condition of each sunlight measuring point obtained by the sunlight shielding algorithm module and a specified sunlight duration calculation mode, and marking the corresponding sunlight duration at the position of the sunlight measuring point in a first display mode if the non-shielded light rays reach the sunlight duration standard; otherwise, marking the corresponding sunshine duration at the position of the sunshine measuring point in a second display mode; the marking fonts, the font sizes and the thicknesses of the sunshine duration are all larger than those of the first display mode in the second display mode, and the marking colors of the sunshine duration are different from those of the first display mode in the second display mode.
Optionally, the real-time sunlight detection system according to any one of the preceding claims, wherein the rod shadow ray set generating module specifically calculates the rod shadow ray set under the sunlight condition corresponding to the target building in each time interval according to the following steps: step 101, according to latitude and longitude (lambda, phi) of a target building, declination angle delta and true solar angle omega under a corresponding time gap, and sinh s Calculating to obtain a solar altitude angle hs corresponding to a target building according to the formula of sin phi sin delta+cos phi cos delta cos omegaCalculating to obtain a solar azimuth angle As corresponding to the target building; step 102, calculating a bar-shadow bus terminal coordinate point P (Xp, yp, zp) matched with the maximum height H of a surrounding building, wherein xp=H×sinAs×coths, yp= -1×H×cosas×coths, zp=H; step 103, connecting the coordinate origin O (0, 0) with the terminal coordinate point P of the bar shadow bus to obtain the bar shadow light OP under the sunshine condition corresponding to the target building under the corresponding time gap, wherein the two-dimensional projection of the bar shadow light OP is OP 0 Wherein P is 0 Coordinates (Xp, yp, 0); step 104, respectively executing the following steps according to the preset time intervalThe steps 101 to 102 store the set of rod shadows OP under the sunlight condition corresponding to the target building in each time interval as a rod shadow set.
Optionally, the real-time sunlight detection system according to any one of the preceding claims, wherein the sunlight measuring point generating module specifically determines the sunlight measuring point of the target building according to the following steps: step 201, taking the standard floor plan established by using the multi-section lines in anticlockwise sequence as a reference, obtaining a standard floor plan of the target building, extracting the coordinates of each endpoint in the standard floor plan, recording the endpoint at the southwest end as Pt (Xpt, ypt, zpt), recording the endpoint at the westst end as Pt Xmin (Xmin, X0, Z0) in which the end point Pt located at the eastern end is recorded Xmax (Xmax, X1, Z1); step 202, if the number of endpoints at the south end exceeds 1, determining that the deflection angle of the target building is 0, otherwise, determining that the target building deflects; step 203, when it is determined that there is deflection in the target building, extracting an edge line segment L overlapping a start point coordinate value Pts (Xstart, ystart, zstart) and an end point Pt (Xpt, ypt, zpt) located at the south end in the standard floor plan view, recording an end point coordinate of the edge line segment L as Pte (Xend, yend, zend), and creating a vector V (Xend-Xstart, yend-Zstart) according to the direction of the edge line segment L; step 204, calculating the included angle between the vector V and the axial quantity Vx (1, 0) to obtain a deflection value RotAngle corresponding to the target building; step 205, calculating a deflection angle corresponding to the target building according to the deflection value RotAngle: when (when) When the bias Angle is determined to be angle=rotangle-90 DEG, whenDetermining that the declination Angle is angle=rotangle; step 206, taking the center of the standard floor type plane graph as a rotation point, deflecting the standard floor type plane graph by an angle of-1 x, and under the state, acquiring each extreme value coordinate point in the rotated standard floor type plane graph as follows: x 'max, X'min, Y ' max, Y ' min, Z '; step 207, constructing an external rectangle of the rotated standard layer house type plane graph according to each extreme value coordinate point, and marking corresponding virtual measuring points on a southbound edge line of the external rectangle according to the number of measuring points or the interval requirement; and step 208, taking the virtual measuring points as the bottom ends, making a vertical line along a standard house type plane of the target building, marking corresponding coordinates on the vertical line as initial measuring points according to the height requirement of the measuring points or house type floor height, and deflecting each virtual measuring point and each initial measuring point back to an Angle by taking the center of a standard house type plane graph as a rotation center to obtain a sunlight measuring point corresponding to the south of the target building.
Optionally, the real-time insolation detection system according to any of the above, further comprising: and the automatic updating module automatically triggers the bar shadow ray set generating module, the sunlight measuring point generating module, the sunlight shielding algorithm module, the sunlight duration calculating module and the recalculation to display the new sunlight duration corresponding to each sunlight measuring point when receiving the instruction of moving or changing any target building or surrounding buildings.
Advantageous effects
According to the method, firstly, a rod shadow ray set under the sunlight condition corresponding to a target building in each time interval is calculated according to the longitude and latitude and the height of the target building in the preset time interval. And then respectively calling the rod shadow ray set to the position of each sunlight measuring point according to the sunlight measuring points of the building, and marking the sunlight time of the sunlight measuring points according to whether the rod shadow ray set rays are shielded by the obstacle or not. Finally, according to a specified sunshine duration calculation mode, the sunshine time of each sunshine measuring point of the target building is counted, and the sunshine duration of each sunshine measuring point is correspondingly identified. The light rays in the rod shadow light ray concentration adopted by the method are matched with the sunlight azimuth of a target building, different sunlight time gaps can be marked at the same time, so that real-time illumination conditions of sunlight measuring points under different sunlight time gaps can be obtained directly through whether the light rays are shielded by obstacles, and statistics of effective sunlight time is accurately realized according to different sunlight time calculation modes flexibly.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention.
Drawings
The accompanying drawings are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate the invention and together with the embodiments of the invention, and do not limit the invention. In the drawings:
FIG. 1 is a schematic diagram of a sunlight shielding principle on which the real-time sunlight detection method and system of the present invention is based;
FIG. 2 is a top plan view corresponding to FIG. 1;
FIG. 3 is a side view corresponding to FIG. 1;
FIG. 4 is a schematic diagram of a computing rod shadow ray in accordance with the present invention;
FIG. 5 is a flow chart of the real-time sunlight detection method provided by the invention;
FIG. 6 is a flow chart of a detailed calculation step in the method of the present invention;
FIG. 7 is a plan view of a non-occluded ray set and a solar duration for an unqualified solar measuring point obtained by the present invention
FIG. 8 is a perspective view of a set of non-occluded rays and a length of sunlight for an substandard sunlight measuring point obtained by the present invention;
FIG. 9 is an enlarged partial schematic view of the undershot spot of FIG. 8
FIG. 10 is a side view partial magnified view of a non-occluded ray set and a solar duration of a non-standard solar measurement point obtained by the present invention.
Detailed Description
In order to make the purpose and technical solutions of the embodiments of the present invention more clear, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings of the embodiments of the present invention. It will be apparent that the described embodiments are some, but not all, embodiments of the invention. All other embodiments, which can be made by a person skilled in the art without creative efforts, based on the described embodiments of the present invention fall within the protection scope of the present invention.
It will be understood by those skilled in the art that, unless otherwise defined, all terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. It will be further understood that terms, such as those defined in commonly used dictionaries, should be interpreted as having a meaning that is consistent with their meaning in the context of the prior art and will not be interpreted in an idealized or overly formal sense unless expressly so defined herein.
The meaning of "and/or" in the present invention means that each exists alone or both exist.
The meaning of "up and down" in the present invention means that the direction from the ground to the sun is up or down in the right angle to the building, and is not limited to the specific mechanism of the device of the present invention.
Grasshopper (Grasshopper) is a parameterized platform on Rhino (Rhino) modeling software that is mainly used for various three-dimensional modeling. The method is based on CAD (Computer Aided Design), namely general three-dimensional modeling and drawing software and Grasshopper parameterization platform on the Rhino modeling software, realizes the simulation and calculation of sunlight rays, so as to obtain data of sunlight conditions of different sunlight measuring points of a target building, and further judges and counts the effective sunlight duration of the target building according to the data in a specified calculation mode, so as to obtain the sunlight condition of the building.
The real-time insolation detection system of the present invention may be configured to include:
the rod shadow ray set generation module is used for respectively calculating rod shadow ray sets under the corresponding sunlight conditions of the target building in each time interval according to the longitude and latitude and the height of the target building and the preset time interval;
the sunlight measuring point generating module is used for determining corresponding sunlight measuring points for each target building respectively;
the sunlight shielding algorithm module is used for respectively starting a sunlight shielding algorithm thread for each target building, and judging whether the concentrated light of the rod shadow rays is shielded or not for each sunlight measuring point of the target building through the sunlight shielding algorithm thread so as to obtain the light shielding condition of each sunlight measuring point;
the sunlight duration calculation module is used for respectively calculating whether the non-shielded light rays corresponding to each sunlight measuring point reach the sunlight duration standard according to the light shielding condition of each sunlight measuring point obtained by the sunlight shielding algorithm module and a specified sunlight duration calculation mode, and marking the corresponding sunlight duration at the position of the sunlight measuring point in a first display mode if the non-shielded light rays reach the sunlight duration standard; otherwise, marking the corresponding sunshine duration at the position of the sunshine measuring point in a second display mode; the marking fonts, the font sizes and the thicknesses of the sunshine duration are all larger than those of the first display mode in the second display mode, and the marking colors of the sunshine duration are different from those of the first display mode in the second display mode.
The system can realize statistics and calculation of the sunlight duration corresponding to each sunlight measuring point in a manner and according to the steps shown in fig. 5:
firstly, establishing a rod shadow ray set;
then, determining a household southern sunlight measuring point;
then, a sunlight shielding algorithm thread is executed on each target building in parallel, so that real-time illumination conditions of sunlight time of each sunlight measuring point on each target building are obtained;
finally, the statistics of the effective sunshine duration is accurately realized flexibly according to the requirement of the sunshine duration calculation mode.
In a more specific implementation manner, the rod shadow ray set may be established by the following steps:
1. selecting the maximum height of all buildings to be measured as a building maximum height H; taking 'Year (Year), month (Month), day (Day), time (Hour) and time (Hour) as parameters, and converting minute data into data in hours, geographic longitude (lambda) and geographic latitude (phi)' through a calculation formula sinh of solar altitude angle s =sThe calculation formula of in phi sin delta+cos phi cos delta cos omega and solar azimuth angle comprises the following steps:and after the related parameters are corrected, obtaining a solar altitude hs and a solar azimuth As at the moment, wherein h s Representing the solar altitude (°), A representing the rod shadow azimuth (°), A s Represents the azimuth angle (°), h s Representing a solar altitude angle (°), phi represents a geographical latitude (°), delta represents an declination angle (°), and omega represents a true solar time angle (°;
2. as can be seen from fig. 4, the end coordinates of the bar-shadow ray bus bar corresponding to the maximum height H of the building can be expressed as x=h×sinas×coths, y= -1×h×cosas×coths, z=h, and thus the terminal coordinate point P (Xp, yp, zp) of the bar-shadow bus bar at that time can be calculated, where xp=h×sinas×coths, yp= -1×h×cosas×coths, zp=h, and the point P and the coordinate origin O (0, 0) shown in fig. 4 are connected to obtain the three-dimensional bar-shadow ray OP with the coordinate origin as the sunlight measuring point at that time, and the Z value of the point P is changed to 0, and the two-dimensional projection point P of 0.00 elevation plane in the coordinate system can be obtained 0 (Xp,Yp,0);
3. The method comprises the steps of (1) changing the value of time (Hour) according to the 'time gap precision', repeating the steps (2) and (3) to obtain a series of three-dimensional rod shadow rays and corresponding two-dimensional projection points, wherein the combination of the three-dimensional rod shadow rays and the corresponding two-dimensional projection points is named as a 'rod shadow ray set' (the display effect of the rod shadow ray set is similar to that of ray marks in fig. 7 and 8), connecting all the two-dimensional projection points according to the time sequence, and forming a two-dimensional graph together with a world coordinate origin O (0, 0).
Calculated by a calculation formula sinh of solar altitude s Calculation formula of =sin Φsin δ+cos Φcos δcos Ω and solar azimuth:and P (Xp, yp, zp) and other relevant correction parameters can be known, under the specified geographic longitude and latitude, the solar altitude angle and solar azimuth angle of the region can be calculated from specific momentsAnd can obtain the rod shadow ray with the specified rod height. Therefore, the time can be subdivided into adjacent continuous moments within a period of time, the more the subdivided moments are, the more the rod shadow rays are, the more the incident condition of the sun rays can be accurately simulated, and otherwise, the more inaccurate the incident condition of the sun rays is. The invention names the time difference between each successive moment as "time gap precision", such as: when the time gap precision is 5 minutes (time difference is 5 minutes), the successive moments are 8:05, 8:10, 8:15, 8:20. When the time gap accuracy is 1 minute (time difference is 1 minute), the successive times are 8:01, 8:02, 8:03, 8:04.
However, the time varies every minute and every second, and it is impossible to calculate the rod shadow ray at every moment, so the invention considers that the rod shadow ray at each moment can be regarded as representing a small period of time before and after the moment on the premise of setting the time gap precision, and the value of the small period of time is equal to the time gap precision. Specifically, when the time gap accuracy is 1 minute, the bar shadow ray at the time of 8:01 (8 hours 1 minute 0 seconds) represents a period of 1 minute between 8:00:30 (8 hours 0 minutes 30 seconds) and 8:01:30 (8 hours 1 minute 30 seconds).
Through the above discussion, each rod shadow ray discussed in the present invention has two properties, one is the spatial position of the rod shadow ray at a specific time (the spatial position characteristic of the rod shadow ray reflects on the "1. Occlusion principle" part of d.2.1), and the other is a time period, and the value of the time period is equal to the "time gap precision".
Therefore, on the premise of "time gap precision", a plurality of rod shadows (i.e. "rod shadow ray sets") can be used to represent a long period of time, and the three of time, time and number of rod shadow ray sets are combined, as shown in the following formula:
t=t×n, (formula 13);
wherein: t represents a time period (minutes); t represents time gap accuracy (minutes); n represents the number of rays of the rod shadow ray set.
4. And confirming a sunlight measuring point of the house type south which obviously influences the sunlight duration value.
The three-mode model of the residential building generally adopts a standard floor plane Plan as a bottom surface to establish a model Brep with a height H, and the overall shape is similar to a rectangle. Therefore, the standard layer plane Plan is taken as a target for acquiring the measuring points, the using plane is a default modeling plane Plan.
The method comprises the following steps of:
step 401: establishing a standard floor type plane by using a multi-section line according to a anticlockwise sequence, and detecting the deflection Angle of the floor type plane relative to the (0, -1, 0) direction, wherein anticlockwise is positive and clockwise is negative;
The specific calculation mode is as follows:
4001. establishing a standard layer plane Plan by using a multi-section line according to the anticlockwise sequence;
4002. the coordinates of each endpoint of the standard layer plane Plan are obtained, and the endpoint Pt (X) with the minimum Y value is found out pt ,Y pt ,Z pt ) Point Pt with minimum X value Xmin (X min ,X 0 ,Z 0 ) Point Pt with maximum X value Xmax (X max ,X 1 ,Z 1 );
4003. Judging the number of points with the minimum Y value, if the number is not 1, indicating that the standard layer plane Plan has no deflection relative to the plane XY, wherein the angle is 0 DEG, and if the angle is 1, indicating that the standard layer plane Plan has deflection;
4004. when the number of points with the minimum Y value is 1, the starting point coordinate value Pt is obtained in all line segments of the standard layer plane Plan s (X start ,Y start ,Z start ) Equal to the coordinate value Pt (X pt ,Y pt ,Z pt ) The endpoint coordinates of the line segment L, L are Pt e (X end ,Y end ,Z end ) And creates a vector V (X) in the direction of L end -X start ,Y end -Y tart ,Z end -Z start );
4005. By vector V and world coordinate positive X-axis V X (1, 0) performing included angle analysis to obtain a deflection angle RotAngle;
4006. the deflection direction is adapted to the world coordinate system, taking Rhino/Grasshopper as an example, clockwise is negative, and anticlockwise is positive:
when (when)The declination Angle is angle=rotangle-90 °;
when (when)The declination Angle is angle=rotangle;
step 402: the central point of the house type plane is taken as a rotation point CenterPt, the angle of-1X angle is deflected, the house type plane is rotated to a non-deflection state, namely the deflection angle is 0 DEG, and in the state, the extreme value (X 'of the coordinates of all points of the house type plane is obtained' max ,X’ min ,Y’ max ,Y’ min To construct two vertexes of the house-type circumscribed rectangle in the south direction, namely (X) min ,Y min ,Z)、(X max ,Y min Z), thereby establishing a southbound line, and establishing virtual sunlight measuring points on the line segment according to the number or the interval of the measuring points;
step 403: making a vertical line with a standard house type plane by using a virtual measuring point, wherein the vertical point is the plane
The original measuring points of (plane. WorldXY) are adjusted to adjust Z values to construct measuring points at different elevations, and the measuring points are deflected by Angle by taking center Pt as a center point to return to an initial state, so that a household type southern sunlight measuring point is obtained "
5. Comparing the two-dimensional projection range of the rod shadow ray set corresponding to the target building with the longitude and latitude positions of surrounding buildings, if the surrounding buildings are in the two-dimensional projection range of the rod shadow ray set, recording the surrounding buildings as obstacles and jumping to the step 6, otherwise, jumping to the step 10;
6. sequencing the obstacles according to the distance between the obstacle and the target building;
7. invoking the rod shadow ray set to each sunlight measuring and calculating point of the target building, and respectively marking serial numbers for each ray in the rod shadow ray set according to the sequence of the sun azimuth angle;
8. respectively establishing 3 memories for each sunlight measuring point of a target building, storing a rod shadow ray set in the fifth step and a sequence number marked by each ray in the rod shadow ray set by using a first memory, storing a sequence number of a ray shielded by an obstacle by using a second memory, and storing a sequence number of a ray not shielded by the obstacle by using a third memory;
9. Taking each sunlight measuring point as a base point, sequentially judging each obstacle: acquiring an azimuth angle range between the boundary of the obstacle and the base point, calling all the light rays positioned in the azimuth angle range from a first memory, judging whether the light rays obtained by calling intersect the obstacle, if so, moving the light rays into a second memory, and if not, moving the light rays into a third memory;
10. and respectively calling the light rays in the third memory corresponding to each sunlight measuring point according to the requirement of a specified sunlight duration calculation mode, and calculating whether the sunlight measuring point reaches the sunlight duration standard.
The principle of the sunlight shielding algorithm adopted by the scheme is described below.
Reference is made to fig. 1, 2 and 3. The condition that no light is incident on the sunlight measuring point (light is blocked) needs to meet the following two formulas:
the minus is not less than As and not more than maxAs; (1)
Wherein:
minus, the minimum angle formed by the buildings in the south of the measuring point and the measuring point relative to the forward south, namely
Angle COS, in fig. 2, is the angle between OC and the forward and the south;
maxAs: the maximum included angle formed by the buildings on the south of the measuring point and the measuring point relative to the forward and the south, namely
Angle DOS, in FIG. 2, is the angle between OD and the forward and the south;
as: the earth is at all solar azimuth angles, in figure 1 +.bos;
measuring and calculating the height of a south building, wherein the height of a No. 2 building is shown in fig. 1;
H p2 measuring and calculating the relative elevation from the bottom of the building on the south to the 0.00 elevation, wherein the elevation is the elevation of the bottom of the building 2# in the figure 1;
H p1 measuring the relative elevation from the bottom of the building where the point is located to the 0.00 elevation, wherein the elevation of the bottom of the 1# building is shown in fig. 1;
h o measuring and calculating the height from the point to the bottom of the building, wherein in fig. 3, the height from the point O to the plane of the bottom of the 1# building is shown;
L As : the shortest distance from the measuring point to the south building is measured in the vertical plane of the solar azimuth As at the moment of the ground, and the length of the OB is shown in fig. 3;
hs: the solar altitude at this point in time is +.LOB in FIG. 3.
Each letter meaning in fig. 1 to 3 is explained:
two buildings with No. 1 and No. 2, and No. 1 and No. 2
O: solar measurement points in the south of the 1# building;
a: solar ray OL passes through the intersection point on the side closer to the O point after the 2# building;
b: the projection point of the point A on the plane where the point O is located;
c: the line segment connected with the building boundary No. 2 and the O point forms a minimum included angle with the forward direction;
d: the line segment connected with the building boundary No. 2 and the O point forms the point of the maximum included angle with the forward direction and the south direction;
H: height of the 2# building;
l: virtual sun position at a specific moment;
s: the direction of the plane where the O point is located
East/West/South/North: east, west, south, north orientation line of + -0.00 standard plane;
and (3) OL: solar rays at that time of the earth;
the following specifically describes a calculation process for realizing the sunlight shielding situation by applying the above principle in the scenes shown in fig. 7 and 8.
Referring to fig. 6, it is shown:
step one: the method comprises the steps of inputting date information, longitude and latitude information, sunlight standard duration, time gap precision and residential building height information, establishing a local solar azimuth angle and solar altitude data set under the condition of the precision, and calculating according to the azimuth angle data set, the altitude data set and building parameters to obtain a three-dimensional rod shadow ray set of the building and two-dimensional projection of the three-dimensional rod shadow ray set, so as to establish a two-dimensional projection boundary line.
Step two: parallel threads are started for each building. Each thread executes the following steps three to eight respectively. The threads of each building are mutually independent, and the threads can correspondingly calculate the light which is not shielded for each sunlight measuring point of each building according to the rod shadow ray set obtained in the step.
Step three: setting elevation, clearance or number of measuring points for each building, and generating necessary sunlight measuring points on the south of the building.
Step four: detecting whether a building has an obstacle or not by using a two-dimensional projection boundary line of a bar shadow ray set, outputting all measuring points of the building and sunlight duration data thereof if the building has no obstacle, setting a duration value as a standard sunlight duration required locally, and jumping to a step nine; if there is an obstacle, continuing the subsequent steps.
Step five: reordering the obstacles of the building from the near to the far according to the distance;
step six: and (3) calling the rod shadow ray set with the corresponding height in the first step to the measuring point of each building, marking serial numbers for each ray of the point, and enabling serial number data to correspond to the serial numbers of the solar azimuth angle data set in the first step one by one.
Step seven: establishing three data storages for each measuring point, wherein the first storage is used for storing the bar shadow ray set and the serial numbers in the step six, the second storage is used for storing the serial numbers of all the shielded rays, and the third storage is used for storing the rays which are not shielded and the serial numbers thereof;
step eight: and D, taking the measuring point as a base point, acquiring the azimuth angle range between the boundary of the first obstacle and the base point in the step five, reading all the rod shadow rays and serial numbers in the range in a first storage, judging by intersecting the rays with the obstacle, if intersecting, indicating that the rays are blocked, storing the serial numbers of the blocked rays in a second storage, if not intersecting, indicating that the rays are not blocked, and storing the rays and the serial numbers in a third storage. And (3) carrying out the same solar azimuth angle range acquisition flow on the second and subsequent obstacles, acquiring the rod shadow rays and the serial numbers in each azimuth angle range in a first storage, judging whether the serial numbers of the rays are in a second storage, if so, not carrying out intersection judgment any more, and if not, continuing to carry out intersection judgment, and storing the judgment results in the second storage or a third storage respectively. And repeating the flow of the step until all the obstacles are judged, and finally obtaining the data of the third storage of the measuring point, namely, the light rays which are not blocked and the serial numbers.
Step nine: and processing the sunlight data output by each measuring point according to a local required sunlight duration accumulation mode to obtain the sunlight duration meeting the requirements. Examples: the local required sunlight duration accumulating mode is as follows: the sun exposure standard of 2 hours on the cold day is accumulated for not more than 2 sections, and each section is not less than 30 minutes. In the actual sunshine process, assuming that the time period of being irradiated by solar rays at a certain measuring point is 9:00-9:25, 10:20-11:15, 14:40-15:25, the continuous time is 25 minutes, 55 minutes and 45 minutes respectively, and 3 actual total time periods are accumulated for 2.08 hours (25+55+45=125 minutes). However, according to the foregoing accumulating mode, the time period of greater 2 periods and more than or equal to 30 minutes should be calculated, and the specified time periods are 10:20-11:15 and 14:40-15:25, and the total time is 1.67 hours (55+45=100 minutes) according to the specification, and at this time, the 2-hour sunlight standard required by the specification is not satisfied. The specifications are more stringent than the actual situation.
The specific calculation mode of the duration T is as follows:
901. acquiring a third storage after the step seven is finished, wherein the third storage stores the light rays which are not shielded and serial numbers thereof;
902. grouping the third storage according to the standard that the difference value between two adjacent serial numbers is 1 to obtain K 0 Group data;
903. counting the number N of light rays of each group k Data N according to the calculation method of t=t×n k Time length calculation T k If the duration T k Discarding if the time is less than 30 minutes, otherwise, reserving; wherein t represents time gap precision (minutes), and N represents the number of rays of the rod shadow ray set;
904. according to T k Is to take its top 2 sets of data T 1 、T 2 Less than 2 groups are 1 group T 1 Accumulate T 1 、T 2 Or T 2 Obtaining a final duration.
Step ten: judging whether the sunshine duration of each measuring point in the step nine accords with the sunshine standard required by the local standard, if so, displaying the sunshine duration at the measuring point only, wherein the display attribute is a black numerical value with fixed size; if the display is not consistent with the window size, displaying the unoccluded ray set at the measuring point, wherein the display attribute of the sunshine numerical value is red, the font size is scaled with the window size in the same proportion, and the display is shown as the effect of fig. 9 or fig. 10.
Step eleven: on the CAD platform, a program can be restarted by setting the program to acquire a mouse or keyboard signal, and the steps one to ten are automatically repeated; on the Grasshopper platform, the procedure is restarted by detecting changes in any building parameters or movement of any building model, and steps one through ten are automatically repeated.
In summary, the invention first calculates the sunlight condition of the target building at each time interval and records the sunlight condition in the form of a bar shadow ray set. And then, according to the sunlight measuring points of the building, respectively calling the rod shadow ray set to the positions of each sunlight measuring point, and according to whether the rod shadow ray set rays are blocked by obstacles, marking the sunlight time. And finally, counting the effective sunlight time of each sunlight measuring point according to a specified calculation mode to obtain the sunlight duration of each sunlight measuring point. The light rays in the rod shadow light ray concentration are matched with the sunlight azimuth of a target building, different sunlight time gaps are marked, real-time illumination conditions of sunlight measuring points under different sunlight time gaps can be obtained directly through whether the light rays are shielded by obstacles, and statistics of effective sunlight time length is accurately achieved flexibly according to the requirement of sunlight time length calculation modes.
The application has the advantages that:
1. only calculate the sunlight shielding condition of the necessary sunlight measuring point, save the computing resource
The invention establishes a set of algorithm which can generate a specified number or interval of sunlight measuring points at a specified elevation of the south of the residence, and the algorithm does not calculate unnecessary sunlight points, thereby greatly saving calculation resources; and the multi-thread coding is adopted, so that CPU resources of a computer are fully utilized, and the operation speed is increased.
2. Dynamic sun display
The algorithm program can synchronously update and calculate with the adjustment of the building model (such as plane change, altitude change, position change and the like) and synchronously display the building model by adding a specific instruction to the CAD/Rhino (Grasshopper) three-dimensional modeling platform. The repeated operation process and time of operators can be simplified, and the operation is smoother and more convenient.
3. Good sunlight numerical vision readability
The invention establishes a concise sunshine duration display mode: and the sunlight display is built only at necessary measuring points, so that the visual interference of unnecessary points is eliminated. Each measuring point calculates and displays the sunlight value of the point according to a sunlight duration accumulation mode, wherein the value has two attributes of color and size; when the sunshine numerical value does not reach the standard sunshine duration, the numerical value is displayed in red, the font size is in fixed proportion to the size of the operation window, and the font size scales with the window size in the same proportion; when the window reaches the standard, the numerical value is displayed in black, and the font size is a fixed value (the value can be manually adjusted) and does not change along with the change of the window size. The visual influence of unnecessary points can be eliminated, the interference of standard reaching points is eliminated, the readability of sunlight data is enhanced, and the judgment error caused by artificial factors is avoided.
4. Multiple sunshine duration calculation modes
A plurality of sunlight duration accumulating modes are established: mode 1 is to accumulate all sunlight time in the effective period, mode 2 is the requirement of sunlight standard of most local residential buildings at present, namely, continuously accumulating sunlight time in a segmented mode, wherein if sunlight standard is 2 hours on a cold day, the continuous sunlight time period is not more than 2 segments, and each segment is not less than 30 minutes. Wherein "day of cold", "2 hours", "2 segments", "30 minutes" are variables, and different regions may be different. The invention can conveniently and rapidly meet the sunlight standards and various calculation rules in various places, and greatly reduces sunlight re-detection work caused by different algorithms.
5. Clear sunshine duration and shielding relation. According to the invention, the light shielding condition and the sunshine duration value are intuitively related, so that the relationship between the light shielding condition and the sunshine duration value can be clearly and conveniently displayed, and the efficiency of adjusting the building layout by an operator is improved.
5.1, in principle of a sun stick shadow map and a sun cone map, a constant value length bus in the sun cone map is improved to be an indefinite value length line generated based on measuring and calculating the maximum height in a building, and then a fan mesh shape is improved to be a fan ray set consisting of indefinite value length lines at each moment based on the sun stick shadow map, and the ray set is named as a 'stick shadow ray set', wherein the precision of the ray set is determined by the quantity of rays, and the more the quantity of rays is, the higher the precision is;
5.2, applying the rod shadow ray set to each sunlight measuring point, obtaining the accurate shielding relation between the ray set of the point and surrounding buildings, and simultaneously calculating according to a sunlight duration accumulating mode to obtain the sunlight duration of the measuring point, and displaying the sunlight value according to the display mode established by the invention when the duration does not reach the standard;
and 5.3, displaying the non-shielded bar shadow ray set at the measuring point only when the sunlight duration of the measuring point does not reach the standard, and displaying the non-shielded bar shadow ray set when the sunlight duration of the measuring point reaches the standard, thereby eliminating the visual interference of the measuring point reaching the standard and clearly displaying the shielding relation.
The invention combines the advantages of both the sun stick shadow map and the sun cone map. a. As with the traditional sunlight cone diagram, the space position and angle of the sunlight incident to the measuring point at each moment can be still displayed, the difference is that each line with different length is essentially displayed at different moments, the area range of the sunlight duration of the measuring point can be influenced, namely if the sunlight cone diagram has an intersection with a certain line, the sunlight is blocked, otherwise, the sunlight is not blocked if the sunlight cone diagram has no intersection; b. the lines at each moment form a continuous surface in space, the overall form is similar to that of the traditional sun-stick photo, and the distinguishing form is formed by lines instead of curved surfaces, so that the non-occluded light rays can be calculated and displayed clearly, and the traditional sun-stick photo cannot display the occlusion relation of each specific light ray
The foregoing is a description of embodiments of the invention, which are specific and detailed, but are not to be construed as limiting the scope of the invention. It should be noted that it will be apparent to those skilled in the art that several variations and modifications can be made without departing from the spirit of the invention, which are all within the scope of the invention.

Claims (10)

1. The real-time sunlight detection method is characterized by comprising the following steps of:
according to the longitude and latitude of the target building and the height thereof, respectively calculating a rod shadow ray set under the sunlight condition corresponding to the target building in each time interval according to a preset time interval;
secondly, respectively determining corresponding sunlight measuring points for each target building, respectively starting a sunlight shielding algorithm thread for each target building, and respectively executing the following third to seventh calculation steps for each sunlight measuring point of the target building through the sunlight shielding algorithm thread;
thirdly, comparing the two-dimensional projection range of the rod shadow ray set corresponding to the target building with the longitude and latitude positions of surrounding buildings, if the surrounding buildings are in the two-dimensional projection range of the rod shadow ray set, recording the surrounding buildings as obstacles and jumping to the fourth step, otherwise, jumping to the eighth step;
Fourth, sorting the barriers according to the distance between the barrier and the target building;
fifthly, calling the rod shadow ray set to each sunlight measuring and calculating point of the target building, and respectively marking serial numbers for each ray in the rod shadow ray set according to the sequence of the sun azimuth angle;
step six, respectively establishing 3 memories for each sunlight measuring point of the target building, storing the rod shadow ray set in the step five and the serial numbers marked by each ray in the rod shadow ray set by using a first memory, storing the serial numbers of the rays shielded by the barrier by using a second memory, and storing the serial numbers of the rays not shielded by the barrier by using a third memory;
seventh, taking each sunlight measuring point as a base point, and sequentially judging each obstacle: acquiring an azimuth angle range between the boundary of the obstacle and the base point, calling all the light rays positioned in the azimuth angle range from a first memory, judging whether the light rays obtained by calling intersect the obstacle, if so, moving the light rays into a second memory, and if not, moving the light rays into a third memory;
and eighth, respectively calling the light rays in the third memory corresponding to each sunlight measuring point according to a specified sunlight duration calculation mode, and calculating whether the sunlight measuring point reaches the sunlight duration standard.
2. The real-time insolation detection method according to claim 1, further comprising: ninth, if the sunlight measuring point reaches the sunlight duration standard, marking the corresponding sunlight duration at the position of the sunlight measuring point in a first display mode; otherwise, marking the corresponding sunshine duration at the position of the sunshine measuring point in a second display mode;
the marking fonts, the font sizes and the thicknesses of the sunshine duration are all larger than those of the first display mode in the second display mode, and the marking colors of the sunshine duration are different from those of the first display mode in the second display mode.
3. The real-time insolation detection method according to claim 2, further comprising: tenth, every time any one of the target building or surrounding building is moved or changed, the second to ninth steps are automatically triggered to be re-executed.
4. The real-time insolation detection method according to claim 3, wherein the bar shadow ray set is calculated by:
step 101, according to latitude and longitude (lambda, phi) of a target building, declination angle delta and true solar angle omega under a corresponding time gap, and sinh s Calculating to obtain a solar altitude angle hs corresponding to a target building according to the formula of sin phi sin delta+cos phi cos delta cos omega Calculating to obtain a solar azimuth angle As corresponding to the target building;
step 102, calculating a bar shadow bus bar terminal coordinate point P (Xp,
yp, zp), wherein xp=h×sinas×coths, yp= -1×h×cosas×coths, zp=h;
step 103, connecting the coordinate origin O (0, 0) with the terminal coordinate point P of the bar shadow bus to obtain the bar shadow light OP under the sunshine condition corresponding to the target building under the corresponding time gap, wherein the two-dimensional projection of the bar shadow light OP is OP 0 Wherein P is 0 Coordinates (Xp, yp, 0);
step 104, executing the steps 101 to 102 respectively according to the predetermined time intervals, and storing the set of rod shadow rays OP under the sunlight condition corresponding to the target building under each time interval as a rod shadow ray set.
5. The method for real-time sunlight detection method according to claim 4, wherein the step of calling the rod shadow ray set to each sunlight measuring point of the target building comprises the following specific steps: and respectively moving the coordinate origin of each rod shadow ray OP in the rod shadow ray set to each sunlight measuring point of the target building to obtain a rod shadow ray O 'P' matched with the sunlight measuring point, wherein the P 'coordinate is (Xp+Xo, yp+yo, zp+zo) and the sunlight measuring point coordinate is O' (Xo, yo, zo).
6. The real-time insolation detection method according to claim 3, wherein the insolation measuring point of the target building is determined by:
step 201, taking the standard floor plan established by using the multi-section lines in anticlockwise sequence as a reference, obtaining a standard floor plan of the target building, extracting the coordinates of each endpoint in the standard floor plan, recording the endpoint at the southwest end as Pt (Xpt, ypt, zpt), recording the endpoint at the westst end as Pt Xmin (Xmin, X0, Z0) in which the end point Pt located at the eastern end is recorded Xmax (Xmax, X1, Z1); step 202, if the number of endpoints at the south end exceeds 1, determining that the deflection angle of the target building is 0, otherwise, determining that the target building deflects;
step 203, when it is determined that there is deflection in the target building, extracting an edge line segment L where a start point coordinate value Pts (Xstart, ystart, zstart) in the standard floor plan view overlaps with an end point Pt (Xpt, ypt, zpt) located at the south end, recording the end point coordinate of the edge line segment L as Pte (Xend, yend, zend), creating a vector V (Xend-Xstart, yend-yttart) according to the direction of the edge line segment L,
Zend-Zstart);
step 204, calculating the included angle between the vector V and the axial quantity Vx (1, 0) to obtain a deflection value RotAngle corresponding to the target building;
Step 205, calculating a deflection angle corresponding to the target building according to the deflection value RotAngle:
when (when)When the deflection Angle is determined to be angle=rotangle-90 °, when +.>Determining that the deflection Angle is angle=rotangle;
step 206, taking the center of the standard floor type plane graph as a rotation point, deflecting the standard floor type plane graph by an angle of-1 x, and under the state, acquiring each extreme value coordinate point in the rotated standard floor type plane graph as follows: x ' max, X ' min, Y ' max, Y ' min, Z ';
step 207, constructing an external rectangle of the rotated standard layer house type plane graph according to each extreme value coordinate point, and marking corresponding virtual measuring points on a southbound edge line of the external rectangle according to the number of measuring points or the interval requirement;
and step 208, taking the virtual measuring points as the bottom ends, making a vertical line along a standard house type plane of the target building, marking corresponding coordinates on the vertical line as initial measuring points according to the height requirement of the measuring points or house type floor height, and deflecting each virtual measuring point and each initial measuring point back to an Angle by taking the center of a standard house type plane graph as a rotation center to obtain a sunlight measuring point corresponding to the south of the target building.
7. A real-time insolation detection system, characterized in that the real-time insolation detection method according to any one of claims 1-6 is used, the real-time insolation detection system comprising:
the rod shadow ray set generation module is used for respectively calculating rod shadow ray sets under the corresponding sunlight conditions of the target building in each time interval according to the longitude and latitude and the height of the target building and the preset time interval;
the sunlight measuring point generating module is used for determining corresponding sunlight measuring points for each target building respectively;
the sunlight shielding algorithm module is used for respectively starting a sunlight shielding algorithm thread for each target building, and judging whether the concentrated light of the rod shadow rays is shielded or not for each sunlight measuring point of the target building through the sunlight shielding algorithm thread so as to obtain the light shielding condition of each sunlight measuring point;
the sunlight duration calculation module is used for respectively calculating whether the non-shielded light rays corresponding to each sunlight measuring point reach the sunlight duration standard according to the light shielding condition of each sunlight measuring point obtained by the sunlight shielding algorithm module and a specified sunlight duration calculation mode, and marking the corresponding sunlight duration at the position of the sunlight measuring point in a first display mode if the non-shielded light rays reach the sunlight duration standard; otherwise, marking the corresponding sunshine duration at the position of the sunshine measuring point in a second display mode; the marking fonts, the font sizes and the thicknesses of the sunshine duration are all larger than those of the first display mode in the second display mode, and the marking colors of the sunshine duration are different from those of the first display mode in the second display mode.
8. The real-time insolation detection system according to claim 7, wherein the bar shadow ray set generating module calculates the bar shadow ray set under the insolation condition corresponding to the target building in each time interval according to the following steps:
step 101, according to latitude and longitude (lambda, phi) of a target building, declination angle delta and true solar angle omega under a corresponding time gap, and sinh s Calculating to obtain a solar altitude angle hs corresponding to a target building according to the formula of sin phi sin delta+cos phi cos delta cos omegaCalculating to obtain a solar azimuth angle As corresponding to the target building;
step 102, calculating a bar shadow bus bar terminal coordinate point P (Xp,
yp, zp), wherein xp=h×sinas×coths, yp= -1×h×cosas×coths, zp=h;
step 103, connecting the coordinate origin O (0, 0) with the terminal coordinate point P of the bar shadow bus to obtain the bar shadow light OP under the sunshine condition corresponding to the target building under the corresponding time gap, wherein the two-dimensional projection of the bar shadow light OP is OP 0 Wherein P is 0 Coordinates (Xp, yp, 0);
step 104, executing the steps 101 to 102 respectively according to the predetermined time intervals, and storing the set of rod shadow rays OP under the sunlight condition corresponding to the target building under each time interval as a rod shadow ray set.
9. The real-time insolation detection system according to claim 8, wherein the insolation measuring point generating module determines the insolation measuring point of the target building according to the following steps:
step 201, taking the standard floor plan established by using the multi-section lines in anticlockwise sequence as a reference, obtaining a standard floor plan of the target building, extracting the coordinates of each endpoint in the standard floor plan, recording the endpoint at the southwest end as Pt (Xpt, ypt, zpt), recording the endpoint at the westst end as Pt Xmin (Xmin, X0, Z0) in which the end point Pt located at the eastern end is recorded Xmax (Xmax, X1, Z1); step 202, if the number of endpoints at the south end exceeds 1, determining that the deflection angle of the target building is 0, otherwise, determining that the target building deflects;
step 203, when it is determined that there is deflection in the target building, extracting an edge line segment L where a start point coordinate value Pts (Xstart, ystart, zstart) in the standard floor plan view overlaps with an end point Pt (Xpt, ypt, zpt) located at the south end, recording the end point coordinate of the edge line segment L as Pte (Xend, yend, zend), creating a vector V (Xend-Xstart, yend-yttart) according to the direction of the edge line segment L,
Zend-Zstart);
Step 204, calculating the included angle between the vector V and the axial quantity Vx (1, 0) to obtain a deflection value RotAngle corresponding to the target building;
step 205, calculating a deflection angle corresponding to the target building according to the deflection value RotAngle:
when (when)When the deflection Angle is determined to be angle=rotangle-90 °, when +.>Determining that the deflection Angle is angle=rotangle;
step 206, taking the center of the standard floor type plane graph as a rotation point, deflecting the standard floor type plane graph by an angle of-1 x, and under the state, acquiring each extreme value coordinate point in the rotated standard floor type plane graph as follows: x ' max, X ' min, Y ' max, Y ' min, Z ';
step 207, constructing an external rectangle of the rotated standard layer house type plane graph according to each extreme value coordinate point, and marking corresponding virtual measuring points on a southbound edge line of the external rectangle according to the number of measuring points or the interval requirement;
and step 208, taking the virtual measuring points as the bottom ends, making a vertical line along a standard house type plane of the target building, marking corresponding coordinates on the vertical line as initial measuring points according to the height requirement of the measuring points or house type floor height, and deflecting each virtual measuring point and each initial measuring point back to an Angle by taking the center of a standard house type plane graph as a rotation center to obtain a sunlight measuring point corresponding to the south of the target building.
10. The real-time insolation detection system according to claim 9, further comprising: and the automatic updating module automatically triggers the bar shadow ray set generating module, the sunlight measuring point generating module, the sunlight shielding algorithm module, the sunlight duration calculating module and the recalculation to display the new sunlight duration corresponding to each sunlight measuring point when receiving the instruction of moving or changing any target building or surrounding buildings.
CN202111336167.1A 2021-04-18 2021-11-12 Real-time sunlight detection method and system Active CN113820756B (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202110415316 2021-04-18
CN2021104153167 2021-04-18

Publications (2)

Publication Number Publication Date
CN113820756A CN113820756A (en) 2021-12-21
CN113820756B true CN113820756B (en) 2024-01-30

Family

ID=78919605

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111336167.1A Active CN113820756B (en) 2021-04-18 2021-11-12 Real-time sunlight detection method and system

Country Status (1)

Country Link
CN (1) CN113820756B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114693513B (en) * 2022-03-28 2024-03-26 西南石油大学 House sunshine analysis method based on visual image

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007052029A (en) * 2006-09-20 2007-03-01 Asahi Kasei Homes Kk Method of sunshine simulation for housing
CN108415103A (en) * 2018-01-15 2018-08-17 马步真 A kind of foundation of building sunshine tester and its application method
CN111505743A (en) * 2020-03-31 2020-08-07 上海天华崧易建筑设计有限公司 Building sunshine duration measuring instrument
CN211717431U (en) * 2020-05-09 2020-10-20 上海天华崧易建筑设计有限公司 Slope ruler for sunlight elevation height
CN211905746U (en) * 2020-04-14 2020-11-10 上海天华崧易建筑设计有限公司 Building sunlight sweeping angle control ruler

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007052029A (en) * 2006-09-20 2007-03-01 Asahi Kasei Homes Kk Method of sunshine simulation for housing
CN108415103A (en) * 2018-01-15 2018-08-17 马步真 A kind of foundation of building sunshine tester and its application method
CN111505743A (en) * 2020-03-31 2020-08-07 上海天华崧易建筑设计有限公司 Building sunshine duration measuring instrument
CN211905746U (en) * 2020-04-14 2020-11-10 上海天华崧易建筑设计有限公司 Building sunlight sweeping angle control ruler
CN211717431U (en) * 2020-05-09 2020-10-20 上海天华崧易建筑设计有限公司 Slope ruler for sunlight elevation height

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于棒影法的居室有效日照时数分析——以清远市区为例;彭窈;清远职业技术学院学报;第12卷(第5期);全文 *

Also Published As

Publication number Publication date
CN113820756A (en) 2021-12-21

Similar Documents

Publication Publication Date Title
CN110111414B (en) Orthographic image generation method based on three-dimensional laser point cloud
CN113034689B (en) Laser point cloud-based terrain three-dimensional model, terrain map construction method and system, and storage medium
CN108010092A (en) A kind of city high density area Solar use potential evaluation method based on low altitude photogrammetry
US10115165B2 (en) Management of tax information based on topographical information
CN110335355A (en) A kind of high automatic calculating method of large-scale shallow water lake
CN110070038B (en) Method for measuring and calculating urban green land rate based on urban greening coverage area
CN115345881B (en) Pavement disease detection method based on computer vision
CN113820756B (en) Real-time sunlight detection method and system
CN111723464A (en) Typhoon elliptic wind field parametric simulation method based on remote sensing image characteristics
CN106780586B (en) A kind of solar energy potential evaluation method based on ground laser point cloud
Dai et al. Road extraction from high-resolution satellite images based on multiple descriptors
CN115063551B (en) Method and device for generating slice orthoimage based on oblique photography three-dimensional model
CN114117702A (en) Point cloud-based automatic reverse modeling method for power transmission line
Stasinopoulos A survey of solar envelope properties using solid modelling
CN111982077B (en) Electronic map drawing method and system and electronic equipment
CN114937137B (en) BIM and GIS-based building environment intelligent analysis method
CN115760885B (en) High-closure-degree wetland forest parameter extraction method based on consumer-level unmanned aerial vehicle image
CN110135049B (en) Distributed photovoltaic power station automatic configuration design method based on three-dimensional point cloud
CN116778098A (en) Plant disease and pest data set construction method based on deep learning
CN114463494B (en) Automatic topographic feature line extraction method
CN115984490A (en) Modeling analysis method and system for wind field characteristics, unmanned aerial vehicle equipment and storage medium
Xu et al. Methods for the construction of DEMs of artificial slopes considering morphological features and semantic information
CN113761628B (en) Building street facade transformation parameterized design generation method
CN114912171A (en) Building shielding-based sunlight radiation calculation method
CN114581775A (en) Scene matching feature data processing method and system based on city three-dimensional space form

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant