CN103065055B - Light distribation computing method in heliogreenhouse - Google Patents
Light distribation computing method in heliogreenhouse Download PDFInfo
- Publication number
- CN103065055B CN103065055B CN201310015165.1A CN201310015165A CN103065055B CN 103065055 B CN103065055 B CN 103065055B CN 201310015165 A CN201310015165 A CN 201310015165A CN 103065055 B CN103065055 B CN 103065055B
- Authority
- CN
- China
- Prior art keywords
- heliogreenhouse
- bin
- light
- direct
- distribation
- 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
Links
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
The present invention relates to the heliogreenhouse technical field in industrialized agriculture, be specifically related to Light distribation computing method in a kind of heliogreenhouse.In this heliogreenhouse, Light distribation computing method comprise step: S1. obtains the distributed intelligence of heliogreenhouse exterior light; S2. construct initial luminous environment in heliogreenhouse in conjunction with described exterior light distributed intelligence and calculate the initial beam intensity of each bin; S3. in conjunction with described initial beam intensity, heliogreenhouse inscattering Light distribation is calculated according to radiosity method; S4. by direct light distribution in heliogreenhouse and scattering light distribution superposition, light distribution information in heliogreenhouse is obtained.The present invention accurately can calculate the light distribution information in heliogreenhouse, for the research of sunlight Interior Illumination Environment provides Data support, for the development of industrialized agriculture provides power-assisted.
Description
Technical field
The present invention relates to the heliogreenhouse technical field in industrialized agriculture, be specifically related to Light distribation computing method in a kind of heliogreenhouse.
Background technology
When the contradiction of " population, resource, environment and food " is day by day serious, Facility Agriculture Development is in the ascendant, has become one of the most great-hearted new industry in the world today.Greenhouse production is as the important component part of modern installations agricultural, and have high in technological content, added value of product is high, land productivity and labour productivity high, is important channel that is populous, the sustainable development of shortage of resources National agricultural.At present, focus mostly on in Structural Performance of Greenhouse, building materials, daylighting, warming, heat supply, ventilation, balance fertilizing, water-saving irrigation, prevention and control of plant diseases, pest control etc. about the research of heliogreenhouse, the research in sunlight Interior Illumination Environment then rarely has report.
Although sunlight Interior Illumination Environment has material impact to plant growth, Light distribation is carried out for heliogreenhouse and calculates that to calculate with insulation be the important evaluation criteria understanding the whether applicable actual agricultural production of this heliogreenhouse; But sunlight Interior Illumination Environment is compared with outdoor, more complex, not only light conditions outer with heliogreenhouse is relevant for it, is also subject to the impact of the factors such as the optical characteristics of the shape of heliogreenhouse, orientation, cladding material and structured material; Therefore, heliogreenhouse luminous environment is the difficult point of research always.The research of current sunlight Interior Illumination Environment, main with day solar radiation semi-invariant for index, investigate the impact on direct light in heliogreenhouse and scattered light transmissivity and distribution of light etc. such as heliogreenhouse orientation, structured material, and on this basis layouting of optical sensor is studied, make measured value more accurately with representative.
In sum, in actual, the heliogreenhouse that possesses high practicability of a kind of combination, Light distribation computing method urgently provide.
Summary of the invention
(1) technical matters that will solve
The object of the present invention is to provide Light distribation computing method in a kind of heliogreenhouse, for accurately calculating the light distribution information in heliogreenhouse, for the research of sunlight Interior Illumination Environment provides Data support.
(2) technical scheme
Technical solution of the present invention is as follows:
Light distribation computing method in a kind of heliogreenhouse, comprise step:
S1. the distributed intelligence of heliogreenhouse exterior light is obtained;
S2. construct initial luminous environment in heliogreenhouse in conjunction with described exterior light distributed intelligence and calculate the initial beam intensity of each bin;
S3. in conjunction with described initial beam intensity, heliogreenhouse inscattering Light distribation is calculated according to radiosity method;
S4. by direct light distribution in heliogreenhouse and scattering light distribution superposition, light distribution information in heliogreenhouse is obtained.
Preferably, described step S1 comprises:
S101. altitude of the sun and position angle is calculated according to current time and heliogreenhouse geographic location information;
S102. the outside direct projection intensity of illumination of heliogreenhouse and light scattering intensity is calculated.
Preferably, described step S102 comprises:
The outside total light intensity of note heliogreenhouse is I, and the number percent of direct light is q;
Then the outside direct light of heliogreenhouse is according to intensity I
direct=I*q,
Heliogreenhouse externals stray intensity of illumination I
diffuse=I* (1-q).
Preferably, described step S2 comprises:
S201. heliogreenhouse all bins projection plane on direct light direction except front domatic bin is divided into multiple subregion, every sub regions at least comprises the view field of a bin;
S202. the bin that this subregion of chosen distance is nearest from bin corresponding to every sub regions is light source bin and calculates the initial beam intensity of light source bin, and other bin initial strengths that this subregion is corresponding are set to zero;
S203. the initial beam intensity of domatic bin before calculating.
Preferably, by Z-Buffer projecting method, heliogreenhouse all bins except front domatic bin are projected on direct light direction in described step S201.
Preferably, in described step S202, the initial beam intensity calculating of light source bin comprises:
Judge whether this light source bin and direct light angular separation cosine value cos θ are greater than 0:
Be, then this light source bin initial beam intensity I
d=I
direct* cos θ * τ;
No, then this light source bin initial beam intensity I
d=0;
Wherein, τ is the direct light transmittance of front domatic bin.
Preferably, domatic bin initial beam intensity I before in described step S203
f=I
direct* (-cos θ) * τ+I
diffuse* ω; Wherein, ω is the scattered light transmittance of front domatic bin.
Preferably, described step S3 comprises:
S301. in conjunction with described initial beam intensity, radiometric algorithms's model is built;
S302. radiometric algorithms's model described in iterative, obtains the scattering light distribution of diverse location in heliogreenhouse.
Preferably, in described step S301, radiometric algorithms's model is:
Wherein, B
ifor the radiancy of bin i, N is bin sum, E
ifor the initial illumination of bin i, r
ifor the reflectivity of bin i, F
ijfor the form factor between bin i and bin j.
Preferably, initial illumination
(3) beneficial effect
Light distribation computing method in a kind of heliogreenhouse of the present invention, construct initial luminous environment in heliogreenhouse in conjunction with the distributed intelligence of heliogreenhouse exterior light and calculate the initial beam intensity of each bin, heliogreenhouse inscattering Light distribation is calculated according to radiosity method, by direct light distribution in heliogreenhouse and scattering light distribution superposition, light distribution information in heliogreenhouse can be obtained.The present invention accurately can calculate the light distribution information in heliogreenhouse, for the research of sunlight Interior Illumination Environment provides Data support, for the development of industrialized agriculture provides power-assisted.
Accompanying drawing explanation
Fig. 1 is the process flow diagram of Light distribation computing method in a kind of heliogreenhouse of the present invention.
Embodiment
Below in conjunction with drawings and Examples, the embodiment of invention is described further.Following examples only for illustration of the present invention, but are not used for limiting the scope of the invention.
Light distribation computing method in process flow diagram a kind of heliogreenhouse as shown in Figure 1, mainly comprise step:
S1. the distributed intelligence of heliogreenhouse exterior light is obtained, for the structure of luminous environment initial in heliogreenhouse provides data; This step mainly comprises:
S101. altitude of the sun and solar azimuth is calculated according to current time (date Hour Minute Second) and heliogreenhouse geographic location (longitude and latitude) information;
S102. the outside direct projection intensity of illumination of heliogreenhouse and light scattering intensity is calculated; Be specially: the outside total light intensity of note heliogreenhouse is I, and the number percent wherein shared by direct light is q; Then the outside direct light of heliogreenhouse is according to intensity I
direct=I*q, heliogreenhouse externals stray intensity of illumination I
diffuse=I* (1-q).
S2. initial luminous environment in the exterior light distributed intelligence structure heliogreenhouse obtained in integrating step S1, and by calculating direct light distribution, determine the initial beam intensity of each bin in radiancy computation process further; This step mainly comprises:
S201. heliogreenhouse comprises multiple bin, and the reflection and transmission attribute of each bin can be obtained by actual measurement; Wherein, front domatic bin is light transmissive material, and it can be plastic foil, glass etc.; Each bin all has normal direction, and the enclosed environment radiancy Light distribation conveniently carried out in next step calculates, and the normal direction of all bins is all pointed in heliogreenhouse; Projected on direct light direction by heliogreenhouse all bins except front domatic bin by Z-Buffer projecting method (a kind of three-dimensional picture Hidden surface removal algorithm), view plane normal points to heliogreenhouse; Because front domatic bin is light transmissive material, a catching portion spectrophotometric intensity, does not therefore participate in projection and calculates; Above-mentioned projection plane is divided into multiple subregion, and every sub regions at least comprises the view field of a bin;
S202. the bin that this subregion of chosen distance is nearest from bin corresponding to every sub regions is light source bin, and calculate the initial beam intensity of light source bin, other bin initial strengths that this subregion is corresponding are set to zero; Wherein, the initial beam intensity of light source bin calculates and is specially:
Judge whether this light source bin and direct light angular separation cosine value cos θ are greater than 0:
Be, then this light source bin initial beam intensity I
d=I
direct* cos θ * τ;
No, then this light source bin initial beam intensity I
d=0;
Wherein, τ is the direct light transmittance of front domatic bin.
S203. the initial beam intensity of domatic bin before calculating; Before domatic bin initial beam intensity be using scattered light and direct light separately through part as initial beam intensity, be specially: front domatic bin initial beam intensity I
f=I
direct* (-cos θ) * τ+I
diffuse* ω; Wherein, ω is the scattered light transmittance of front domatic bin, and because front domatic panel method is to inwardly, therefore included angle cosine adds negative sign.
So far, in heliogreenhouse, initial luminous environment structure is complete, and the interior light source bin illuminated by direct light of heliogreenhouse is I to the light intensity value that heliogreenhouse is contributed
direct* co s θ * τ * r
i, wherein r
iit is the reflectivity of i-th bin; Front domatic bin is I to the light intensity value that heliogreenhouse is contributed
f.
The initial beam intensity of each bin S3. calculated in integrating step S2, calculates heliogreenhouse inscattering Light distribation according to radiosity method; Radiosity method is after calculating the initial luminous environment of heliogreenhouse, and the thought based on energy equilibrium adopts iterative process to make sunlight Interior Illumination Environment reach the state of a natural equilibrium; This step mainly comprises:
S301. in conjunction with described initial beam intensity, build radiometric algorithms RGM(Radiosity-Graphics Model, radiation is painted) model; Specific as follows:
Wherein, B
ifor the radiancy of bin i, N is bin sum, E
ifor the initial illumination of bin i, r
ifor the reflectivity of bin i, F
ijfor the form factor between bin i and bin j; Wherein, initial illumination
the initial beam intensity of each bin namely in above-mentioned steps S2; Form factor F between bin i and bin j
ij, be and leave the radiancy ratio that a jth bin arrives i-th bin, can calculate by half cube method;
S302. radiometric algorithms's model described in iterative, after solving the threshold value reaching certain, iteration completes; The B finally obtained
ivalue is the scattering light distribution of diverse location in heliogreenhouse.
S4. by direct light distribution in heliogreenhouse and scattering light distribution superposition, light distribution information in heliogreenhouse is obtained; To calculate the light intensity of assigned address in heliogreenhouse, the bin being zero at assigned address structure reflectivity participates in calculating.
In a kind of heliogreenhouse of the present invention, Light distribation computing method are for domestic modal heliogreenhouse, under the prerequisite of given extraneous illumination condition, utilize the heliogreenhouse light transmissive material after direct light calculates as the light source of enclosed environment, the radiancy of carrying out enclosed environment calculates, and then obtain the luminous environment distribution of heliogreenhouse inside, the luminous environment distribution of this heliogreenhouse also can be provided when designing heliogreenhouse; In the heliogreenhouse that method of the present invention calculates, Light distribation is with to measure the Light distribation obtained in real daylight greenhouse substantially identical.
Above embodiment is only for illustration of the present invention; and be not limitation of the present invention; the those of ordinary skill of relevant technical field; without departing from the spirit and scope of the present invention; can also make a variety of changes and modification, therefore all equivalent technical schemes also belong to protection category of the present invention.
Claims (9)
1. Light distribation computing method in heliogreenhouse, is characterized in that, comprise step:
S1. the distributed intelligence of heliogreenhouse exterior light is obtained;
S2. construct initial luminous environment in heliogreenhouse in conjunction with described exterior light distributed intelligence and calculate the initial beam intensity of each bin;
S3. in conjunction with described initial beam intensity, heliogreenhouse inscattering Light distribation is calculated according to radiosity method;
S4. by direct light distribution in heliogreenhouse and scattering light distribution superposition, light distribution information in heliogreenhouse is obtained;
Described step S2 comprises:
S201. heliogreenhouse all bins projection plane on direct light direction except front domatic bin is divided into multiple subregion, every sub regions at least comprises the view field of a bin;
S202. the bin that this subregion of chosen distance is nearest from bin corresponding to every sub regions is light source bin and calculates the initial beam intensity of light source bin, and other bin initial strengths that this subregion is corresponding are set to zero;
S203. the initial beam intensity of domatic bin before calculating.
2. Light distribation computing method in heliogreenhouse according to claim 1, is characterized in that, described step S1 comprises:
S101. altitude of the sun and position angle is calculated according to current time and heliogreenhouse geographic location information;
S102. the outside direct projection intensity of illumination of heliogreenhouse and light scattering intensity is calculated.
3. Light distribation computing method in heliogreenhouse according to claim 2, is characterized in that, described step S102 comprises:
The outside total light intensity of note heliogreenhouse is I, and the number percent of direct light is q;
Then the outside direct light of heliogreenhouse is according to intensity I
direct=I*q,
Heliogreenhouse externals stray intensity of illumination I
diffuse=I* (1-q).
4. Light distribation computing method in heliogreenhouse according to claim 3, is characterized in that, are projected by heliogreenhouse all bins except front domatic bin in described step S201 by Z-Buffer projecting method on direct light direction.
5. Light distribation computing method in heliogreenhouse according to claim 4, is characterized in that, in described step S202, the initial beam intensity of light source bin calculates and comprises:
Judge whether the cosine value cos θ of this light source bin and direct light angular separation θ is greater than 0:
Be, then this light source bin initial beam intensity I
d=I
direct* cos θ * τ;
No, then this light source bin initial beam intensity I
d=0;
Wherein, τ is the direct light transmittance of front domatic bin.
6. Light distribation computing method in heliogreenhouse according to claim 5, is characterized in that, domatic bin initial beam intensity I before in described step S203
f=I
direct* (-cos θ) * τ+I
diffuse* ω; Wherein, ω is the scattered light transmittance of front domatic bin.
7. Light distribation computing method in heliogreenhouse according to claim 6, is characterized in that, described step S3 comprises:
S301. in conjunction with described initial beam intensity, radiometric algorithms's model is built;
S302. radiometric algorithms's model described in iterative, obtains the scattering light distribution of diverse location in heliogreenhouse.
8. Light distribation computing method in heliogreenhouse according to claim 7, is characterized in that, in described step S301, radiometric algorithms's model is:
Wherein, B
ifor the radiancy of bin i, N is bin sum, E
ifor the initial illumination of bin i, r
ifor the reflectivity of bin i, F
ijfor the form factor between bin i and bin j.
9. Light distribation computing method in heliogreenhouse according to claim 8, is characterized in that, described initial illumination
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310015165.1A CN103065055B (en) | 2013-01-15 | 2013-01-15 | Light distribation computing method in heliogreenhouse |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310015165.1A CN103065055B (en) | 2013-01-15 | 2013-01-15 | Light distribation computing method in heliogreenhouse |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103065055A CN103065055A (en) | 2013-04-24 |
| CN103065055B true CN103065055B (en) | 2015-08-12 |
Family
ID=48107683
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310015165.1A Active CN103065055B (en) | 2013-01-15 | 2013-01-15 | Light distribation computing method in heliogreenhouse |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103065055B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114637351B (en) * | 2022-03-14 | 2023-03-28 | 江西省农业科学院园艺研究所 | Greenhouse environment regulation and control method and system for facility crops |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5313568A (en) * | 1990-05-31 | 1994-05-17 | Hewlett-Packard Company | Three dimensional computer graphics employing ray tracing to compute form factors in radiosity |
| CN101900602A (en) * | 2010-06-24 | 2010-12-01 | 北京农业信息技术研究中心 | Calculation method for distribution of scattered light on crop canopy |
| CN102680084A (en) * | 2012-04-25 | 2012-09-19 | 北京农业信息技术研究中心 | Light distribution measuring device and method for crop canopy |
-
2013
- 2013-01-15 CN CN201310015165.1A patent/CN103065055B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5313568A (en) * | 1990-05-31 | 1994-05-17 | Hewlett-Packard Company | Three dimensional computer graphics employing ray tracing to compute form factors in radiosity |
| CN101900602A (en) * | 2010-06-24 | 2010-12-01 | 北京农业信息技术研究中心 | Calculation method for distribution of scattered light on crop canopy |
| CN102680084A (en) * | 2012-04-25 | 2012-09-19 | 北京农业信息技术研究中心 | Light distribution measuring device and method for crop canopy |
Non-Patent Citations (3)
| Title |
|---|
| Qinhuo Liu,Huaguo Huang,Wenhan Qin等.《An Extended 3-D Radiosity-Graphics Combined Model for Studying Thermal-Emission Directionality of Crop Canopy》.《IEEE Transaction on GeoScience and Remote Sensing》.2007,第45卷(第9期), * |
| 温维亮,孟军,郭新宇等.《基于辐射照度的作物冠层光分布计算系统设计》.《农业机械学报》.2009,第40卷(第S1期), * |
| 王锡平,郭焱,李保国.《植物冠层内光合有效辐射空间分布模拟模型》.《中国数字农业与农村信息化学术研究研讨会论文集》.2005, * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103065055A (en) | 2013-04-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wang et al. | Temperature error correction based on BP neural network in meteorological wireless sensor network | |
| Zhu et al. | The effect of urban morphology on the solar capacity of three-dimensional cities | |
| Escobar et al. | Solar energy resource assessment in Chile: Satellite estimation and ground station measurements | |
| CN103955749B (en) | Direct solar radiation value Forecasting Methodology and system | |
| CN105373997A (en) | Real-time weather data-based evacuation method during occurrence of dangerous chemical leakage accident | |
| Anurag et al. | General design procedures for airport-based solar photovoltaic systems | |
| CN101915914A (en) | Lookup table based pixel-by-pixel atmospheric correction method of remote sensing images | |
| CN102435586A (en) | Method and system for generating earth surface albedo product | |
| Hofierka et al. | The spatial distribution of photovoltaic power plants in relation to solar resource potential: the case of the Czech Republic and Slovakia | |
| Zhao et al. | Coupling meteorological stations data and satellite data for prediction of global solar radiation with machine learning models | |
| Bernabé et al. | Radiative and heat storage properties of the urban fabric derived from analysis of surface forms | |
| CN103065055B (en) | Light distribation computing method in heliogreenhouse | |
| Arias-Rosales et al. | Shadow modeling in urban environments for solar harvesting devices with freely defined positions and orientations | |
| CN103823976B (en) | Heliogreenhouse light thermal environment computational methods | |
| Govehovitch et al. | Numerical validation of the radiative model for the solar cadaster developed for greater Geneva | |
| Manni et al. | Solar energy digitalization at high latitudes: A model chain combining solar irradiation models, a LiDAR scanner, and high-detail 3D building model | |
| Yoon et al. | Window material daylighting performance assessment algorithm: Comparing radiosity and split-flux methods | |
| Danks et al. | Urban scale simulations of solar reflections in the built environment: methodology & validation | |
| Salgado-Conrado et al. | Review of heliodon developments and computational tools for building shadow analysis | |
| JP5724102B2 (en) | Indoor illumination simulation method and system | |
| Liu et al. | Determination of the optimal orientation of Chinese solar greenhouses using 3D light environment simulations | |
| CN102331995A (en) | Point source-based information acquisition method and system for three-dimensional geo-information model | |
| Park et al. | Simulation analysis of daylight characteristics and cooling load based on performance test of covering materials used in smart farms | |
| CN104933302A (en) | Method and system for realizing aviation remote sensing simulation of unmanned aerial vehicle load | |
| Freitas | Photovoltaic potential in building façades |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |