CN113469537A - Urban greenbelt sporopollen sensitization risk value calculation method based on weather and population density - Google Patents
Urban greenbelt sporopollen sensitization risk value calculation method based on weather and population density Download PDFInfo
- Publication number
- CN113469537A CN113469537A CN202110767674.4A CN202110767674A CN113469537A CN 113469537 A CN113469537 A CN 113469537A CN 202110767674 A CN202110767674 A CN 202110767674A CN 113469537 A CN113469537 A CN 113469537A
- Authority
- CN
- China
- Prior art keywords
- urban
- sensitization
- value
- population density
- spore powder
- 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.)
- Granted
Links
- 206010070834 Sensitisation Diseases 0.000 title claims abstract description 145
- 230000008313 sensitization Effects 0.000 title claims abstract description 145
- 238000004364 calculation method Methods 0.000 title abstract description 6
- 238000000034 method Methods 0.000 claims abstract description 30
- 238000012545 processing Methods 0.000 claims abstract description 6
- 239000000843 powder Substances 0.000 claims description 84
- 241000196324 Embryophyta Species 0.000 claims description 55
- 230000010152 pollination Effects 0.000 claims description 27
- 230000017260 vegetative to reproductive phase transition of meristem Effects 0.000 claims description 21
- 239000002609 medium Substances 0.000 claims description 12
- 239000002245 particle Substances 0.000 claims description 12
- 241000894007 species Species 0.000 claims description 12
- 238000005286 illumination Methods 0.000 claims description 7
- 239000012194 insect media Substances 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000001932 seasonal effect Effects 0.000 claims description 6
- 238000010586 diagram Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 238000010998 test method Methods 0.000 claims description 3
- 206010020751 Hypersensitivity Diseases 0.000 abstract description 7
- 208000026935 allergic disease Diseases 0.000 abstract description 7
- 230000007815 allergy Effects 0.000 abstract description 7
- 208000035285 Allergic Seasonal Rhinitis Diseases 0.000 abstract 3
- 206010048908 Seasonal allergy Diseases 0.000 abstract 3
- 201000004338 pollen allergy Diseases 0.000 abstract 3
- 230000009286 beneficial effect Effects 0.000 description 9
- 230000036541 health Effects 0.000 description 6
- 241000282414 Homo sapiens Species 0.000 description 4
- 238000010276 construction Methods 0.000 description 2
- 230000033001 locomotion Effects 0.000 description 2
- 230000008520 organization Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 206010052613 Allergic bronchitis Diseases 0.000 description 1
- 206010012434 Dermatitis allergic Diseases 0.000 description 1
- 206010061218 Inflammation Diseases 0.000 description 1
- 206010039085 Rhinitis allergic Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000004887 air purification Methods 0.000 description 1
- 201000009961 allergic asthma Diseases 0.000 description 1
- 230000000172 allergic effect Effects 0.000 description 1
- 201000010105 allergic rhinitis Diseases 0.000 description 1
- 208000006673 asthma Diseases 0.000 description 1
- 208000010668 atopic eczema Diseases 0.000 description 1
- 210000000795 conjunctiva Anatomy 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000004054 inflammatory process Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 231100000915 pathological change Toxicity 0.000 description 1
- 230000036285 pathological change Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000005180 public health Effects 0.000 description 1
- 210000002345 respiratory system Anatomy 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A30/00—Adapting or protecting infrastructure or their operation
- Y02A30/60—Planning or developing urban green infrastructure
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- Data Mining & Analysis (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Databases & Information Systems (AREA)
- Operations Research (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Game Theory and Decision Science (AREA)
- Computational Mathematics (AREA)
- Quality & Reliability (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Bioinformatics & Computational Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Algebra (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Software Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Remote Sensing (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Cultivation Of Plants (AREA)
Abstract
本发明公开了一种基于气象和人口密度的城市绿地孢粉致敏风险值计算方法，包括以下步骤：S1、根据城市绿地的植物特征，计算每种植物的致敏力；S2、根据每种植物的致敏力，计算城市各类绿地的孢粉致敏风险值；S3、将城市各类绿地的孢粉致敏风险值的分布图进行网格化处理，绘制城市绿地孢粉致敏风险加权值地图；S4、根据城市街道的人口密度，基于人口密度分布图，得到城市人口密度加权值；S5、采用移动式气象站对城市的气象数据进行采集，得到日气象均值；S6、根据城市人口密度加权值和日气象均值，计算基于气象和人口密度的城市绿地孢粉致敏风险值；本发明解决了现有技术中缺乏城市绿地孢粉致敏风险值计算方法的问题。
The invention discloses a method for calculating the risk value of pollen allergy in urban green space based on meteorology and population density. The sensitization power of plants, calculate the sporopollen allergy risk value of various urban green spaces; S3. Perform grid processing on the distribution map of sporopollen allergy risk values of various urban green spaces to draw the sporopollen allergy risk of urban green space Weighted value map; S4. According to the population density of the city streets, based on the population density distribution map, the weighted value of the urban population density is obtained; S5. Use a mobile weather station to collect the meteorological data of the city to obtain the daily meteorological average value; S6. According to the city The population density weighted value and the daily meteorological mean value are used to calculate the urban green space pollen allergy risk value based on meteorology and population density; the invention solves the problem that the prior art lacks the urban green space pollen allergy risk value calculation method.
Description
Technical Field
The invention relates to a method for calculating a spore powder sensitization risk value, in particular to a method for calculating an urban green space spore powder sensitization risk value based on weather and population density.
Background
The health report of the world health organization 2016 indicates that the urban ecological environment is an important factor in determining public health. The urban green land is the region with the most concentrated urban species diversity and is the green carrier which is most directly acted on the human beings and forms the healthy space for the human beings. The plants play a vital role in maintaining the ecological sustainability of urban development and building a livable urban space. Currently, in the current urban green land construction process, most attention points are focused on positive benefits brought to cities by garden plants, such as ornamental value, city air purification, microclimate adjustment, city biological diversity maintenance and the like. The risks brought by vegetation selection to human health and even urban ecosystems in the process of urbanization, such as species invasion, capital construction damage, sporopollen (spores and pollen of plants) sensitization and the like, are ignored.
The spore powder allergy is mainly manifested by catarrhal inflammation of respiratory tract and conjunctiva, such as allergic rhinitis, allergic rash, bronchitis, asthma and the like, can be accompanied by pathological changes of skin and other organs, and seriously threatens human health. The data of the world health organization predicts that 35% of the world population will be allergic to sporopollen in the next 20 years. The incidence rate of sporopollen allergy in China is increased year by year, and the sporopollen allergy is a famous and genuine epidemic disease. However, the urban greenbelt sporopollen sensitization risk assessment and prediction have not been carried out domestically. Therefore, the evaluation of the spore allergy risk of urban green lands is imminent.
Disclosure of Invention
Aiming at the defects in the prior art, the method for calculating the sensitization risk value of the urban greenbelt spore powder based on weather and population density solves the problem that the method for calculating the sensitization risk value of the urban greenbelt spore powder is lacked in the prior art.
In order to achieve the purpose of the invention, the invention adopts the technical scheme that: a method for calculating an urban greenbelt sporopollen sensitization risk value based on weather and population density comprises the following steps:
s1, calculating the allergenicity of each plant according to the plant characteristics of the urban green land;
s2, calculating the spore powder sensitization risk value of various urban greenbelts according to the sensitization force of each plant;
s3, gridding the distribution map of the spore powder sensitization risk values of various greenbelts in the city, and drawing a map of the spore powder sensitization risk weighted value of the greenbelts in the city;
s4, obtaining an urban population density weighted value based on the central urban residential land population density distribution map according to the population density of urban streets;
s5, collecting urban meteorological data by a mobile meteorological station according to the urban greenbelt spore powder sensitization risk weighted value map to obtain a solar meteorological mean value;
and S6, calculating the urban greenbelt sporopollen sensitization risk value based on weather and population density according to the urban population density weighted value and the daily weather mean value.
Further, the plant characteristics in step S1 include: spore powder sensitization capacity, pollination type and flowering time.
Further, step S1 includes the following substeps:
s11, dividing the spore powder sensitization capability of the plants into 4 grades according to the spore powder particle diameter of the urban green land plants;
s12, assigning values of 0-3 to the sensitization capacity of 4 grades of sporopollen in sequence to obtain sporopollen with the sporopollen particle diameter larger than 200 microns, wherein the sensitization capacity is assigned to be 0, the sensitization capacity is assigned to be 1, the sensitization capacity is assigned to be 2, the sporopollen particle diameter is smaller than 20 microns or is 50-100 microns, and the sensitization capacity is assigned to be 3, wherein the sporopollen particle diameter is 20-50 microns;
s13, assigning values to the pollination types by adopting values of 1-3 according to the pollination types of the urban green land plants to obtain the assignment value of the plants of the insect media pollination type as 1, the assignment values of the plants of the insect media pollination and wind media pollination types as 2 and the assignment values of the plants of the wind media pollination types as 3;
s14, measuring the flowering time required by 5% of flowers of urban green land plants to 95% of flowers to wither, assigning the flowering time of 1-3 weeks to 1, assigning the flowering time of 4-6 weeks to 2, and assigning the flowering time of 7 weeks or more to 3;
and S15, calculating the sensitization force of each plant according to the assignment of the sensitization capability of the sporopollen, the assignment of the pollination type and the assignment of the flowering time.
Further, the formula of the sensitization power of each plant in step S15 is:
VPA＝tp·dpp·ap
wherein VPA is the sensitization power of each plant, tp is the assignment of pollination type, dpp is the assignment of flowering time, and ap is the assignment of spore powder sensitization power.
The beneficial effects of the above further scheme are: the method is beneficial to understanding the allergenicity of various plants in the urban green land, and provides a basis for the safety and scientificity of urban green land species selection.
Further, the formula for calculating the spore powder sensitization risk value of various greenbelts in the city in step S2 is as follows:
wherein, IUGZAThe spore powder sensitization risk value of various urban greenbelts, VPA is the sensitization power of each plant, SiIs the floor area of the i-th species, HiHeight of species of type i, STIs the area of a certain urban green land, I is the total species class, and max VPA is the maximum value of VPA in all plants of a certain urban green land.
The beneficial effects of the above further scheme are: the sensitization force of each plant in the green land is comprehensively considered, and superposition calculation is carried out on the basis, so that the calculation of the sensitization risk of the green land is more comprehensive and scientific.
Further, step S3 includes the following substeps:
s31, adopting ArcGIS software to perform gridding processing on the distribution diagram of the spore powder sensitization risk values of various greenbelts in the city to obtain the average spore powder sensitization risk value of unit grids in different seasons;
s32, dividing unit grids into low, medium and high sensitization risk grades according to the average spore powder sensitization risk value;
s33, assigning a value of 1 to the unit grid corresponding to the low sensitization risk level, assigning a value of 2 to the unit grid corresponding to the medium sensitization risk level, and assigning a value of 3 to the unit grid corresponding to the high sensitization risk level;
s34, calculating the seasonal spore powder sensitization risk weighted value of the unit grid by adopting a Sudoku weighting method according to the assignment of the unit grid;
and S35, drawing an urban greenbelt spore powder sensitization risk weighted value map based on the seasonal spore powder sensitization risk weighted value of the unit grid.
The beneficial effects of the above further scheme are: in cities, sporopollen can move for a short distance under the action of wind, so that the sporopollen sensitization risk between adjacent areas easily generates an interactive effect, and the interaction effect of peripheral grids on a core area can be comprehensively considered by adopting a Sudoku weighting method, so that the drawing of an urban greenbelt sporopollen sensitization risk map is more scientific.
Further, step S4 includes the following substeps:
s41, dividing the population density distribution map of the residential land in the central urban area into 5 levels of extremely high density, medium density, low density and extremely low density according to the population density of the urban streets;
s42, assigning the extremely-high density area to be 5, assigning the high density area to be 4, assigning the medium density area to be 3, assigning the low density area to be 2, and assigning the extremely-low density area to be 1 to obtain the population density distribution map of the residential land of the central urban area after grading assignment;
s43, meshing the population density distribution map of the residential land of the central urban area after grading assignment by adopting ArcGIS software to obtain the average population density value of the unit grid;
and S44, obtaining the weighted value of the urban population density according to the average population density value of the unit grid.
The beneficial effects of the above further scheme are: the urban greenbelt spore powder sensitization risk value is calculated to enable the plants to better serve the users, and if the population density of the area where the greenbelt is located is small, health loss caused by greenbelt spore powder sensitization risk is relatively small, and vice versa. Therefore, the urban population density weighted value is calculated, and the problem that different areas of the city can be exposed in a more targeted manner by calculating the urban greenbelt sensitization risk value can be solved.
Further, step S5 includes the following substeps:
s51, setting an urban mobile weather station according to the urban greenbelt sporopollen sensitization risk weighted value map;
s52, collecting air temperature data, relative humidity data, wind speed data, rainfall data and illumination intensity data through the urban mobile weather station;
and S53, calculating the daily average value of the air temperature data, the relative humidity data, the wind speed data, the rainfall data and the illumination intensity data to obtain the daily weather average value of the 5 kinds of weather data.
The beneficial effects of the above further scheme are: air temperature data, relative humidity data, wind speed data, rainfall data and illumination intensity data in urban greenbelts influence the generation and movement of the greenbelt sporopollen, and the scheme is beneficial to more scientifically estimating the generation amount and the movement track of the greenbelt sporopollen.
Further, step S6 includes the following substeps:
s61, standardizing the urban population density weighted value, the spore powder sensitization risk value of various urban greenbelts and the solar meteorological mean value of 5 meteorological data to obtain 7 standardized data;
s62, calculating a component matrix of 7 kinds of standardized data by a KMO method and a Bartlett sphericity test method;
s63, calculating the eigenvalue of the component matrix to obtain the eigenvalue of 7 kinds of standardized data;
s64, calculating the square of the characteristic value of each standardized data to obtain the weight value of each standardized data;
s65, carrying out normalization processing on the weight value of each kind of standardized data to enable the sum of the weight values of the 7 kinds of standardized data to be 1, and obtaining the normalized weight values of the 7 kinds of standardized data;
and S66, calculating the sensitization risk value of the urban greenbelt spore powder based on weather and population density according to the normalized weight values of the 7 kinds of standardized data.
Further, the formula for calculating the sensitization risk value of the urban greenswamp spore powder based on weather and population density in step S66 is as follows:
RUGZA＝a1·IUGZA+a2·x1+a3·x2+a4·x3+a5·x4+a6·x5+a7·x6
wherein R isUGZAIs an urban greenswamp spore powder sensitization risk value based on weather and population density, a1、a2、 a3、a4、a5、a6And a7Normalized weight values, x, for 7 normalized data1、x2、x3、x4、x5、 x6The corresponding urban population density weighted value and the mean value of the solar weather of 5 kinds of weather data, IUGZAThe spore powder sensitization risk value of various greenbelts in cities.
The beneficial effects of the above further scheme are: the influence of main microclimate factors of urban greenbelts on the sensitization risk value of the greenbelts is fully considered, so that the calculation of the risk value is more scientific and comprehensive.
In conclusion, the beneficial effects of the invention are as follows: the method calculates the sensitization power of each plant, further calculates to obtain the spore powder sensitization risk values of various greenbelts in the city, draws a map of the spore powder sensitization risk weighted values of the greenbelts in the city, and then calculates to obtain the spore powder sensitization risk values of the greenbelts in the city based on weather and population density by combining the urban population density weighted value and the weather data of the city.
Drawings
FIG. 1 is a flow chart of a method for calculating the sensitization risk value of urban greenswamp spore powder based on weather and population density.
Detailed Description
The following description of the embodiments of the present invention is provided to facilitate the understanding of the present invention by those skilled in the art, but it should be understood that the present invention is not limited to the scope of the embodiments, and it will be apparent to those skilled in the art that various changes may be made without departing from the spirit and scope of the invention as defined and defined in the appended claims, and all matters produced by the invention using the inventive concept are protected.
As shown in fig. 1, a method for calculating an urban greenswamp spore powder sensitization risk value based on weather and population density comprises the following steps:
s1, calculating the allergenicity of each plant according to the plant characteristics of the urban green land;
the plant characteristics in step S1 include: spore powder sensitization capacity, pollination type and flowering time.
Step S1 includes the following substeps:
s11, dividing the spore powder sensitization capability of the plants into 4 grades according to the spore powder particle diameter of the urban green land plants;
s12, assigning values of 0-3 to the sensitization capacity of 4 grades of sporopollen in sequence to obtain sporopollen with the sporopollen particle diameter larger than 200 microns, wherein the sensitization capacity is assigned to be 0, the sensitization capacity is assigned to be 1, the sensitization capacity is assigned to be 2, the sporopollen particle diameter is smaller than 20 microns or is 50-100 microns, and the sensitization capacity is assigned to be 3, wherein the sporopollen particle diameter is 20-50 microns;
s13, assigning values to the pollination types by adopting values of 1-3 according to the pollination types of the urban green land plants to obtain the assignment value of the plants of the insect media pollination type as 1, the assignment values of the plants of the insect media pollination and wind media pollination types as 2 and the assignment values of the plants of the wind media pollination types as 3;
s14, measuring the flowering time required by 5% of flowers of urban green land plants to 95% of flowers to wither, assigning the flowering time of 1-3 weeks to 1, assigning the flowering time of 4-6 weeks to 2, and assigning the flowering time of 7 weeks or more to 3;
and S15, calculating the sensitization force of each plant according to the assignment of the sensitization capability of the sporopollen, the assignment of the pollination type and the assignment of the flowering time.
The formula of the sensitization power of each plant in the step S15 is as follows:
VPA＝tp·dpp·ap
wherein VPA is the sensitization power of each plant, tp is the assignment of pollination type, dpp is the assignment of flowering time, and ap is the assignment of spore powder sensitization power.
S2, calculating the spore powder sensitization risk value of various urban greenbelts according to the sensitization force of each plant;
in step S2, the formula for calculating the spore powder sensitization risk value of various greenbelts in the city is:
wherein, IUGZAThe spore powder sensitization risk value of various urban greenbelts, VPA is the sensitization power of each plant, SiIs the floor area of the i-th species, HiHeight of species of type i, STIs the area of a certain urban green land, I is the total species class, and max VPA is the maximum value of VPA in all plants of a certain urban green land.
S3, gridding the distribution map of the spore powder sensitization risk values of various greenbelts in the city, and drawing a map of the spore powder sensitization risk weighted value of the greenbelts in the city;
step S3 includes the following substeps:
s31, adopting ArcGIS software to perform gridding processing on the distribution diagram of the spore powder sensitization risk values of various greenbelts in the city to obtain the average spore powder sensitization risk value of unit grids in different seasons;
in the present embodiment, the size of the unit cell in step S31 is 500m × 500 m.
S32, dividing unit grids into low, medium and high sensitization risk grades according to the average spore powder sensitization risk value;
s33, assigning a value of 1 to the unit grid corresponding to the low sensitization risk level, assigning a value of 2 to the unit grid corresponding to the medium sensitization risk level, and assigning a value of 3 to the unit grid corresponding to the high sensitization risk level;
s34, calculating the seasonal spore powder sensitization risk weighted value of the unit grid by adopting a Sudoku weighting method according to the assignment of the unit grid;
and S35, drawing an urban greenbelt spore powder sensitization risk weighted value map based on the seasonal spore powder sensitization risk weighted value of the unit grid.
S4, obtaining an urban population density weighted value based on the central urban residential land population density distribution map according to the population density of urban streets;
step S4 includes the following substeps:
s41, dividing the population density distribution map of the residential land in the central urban area into 5 levels of extremely high density, medium density, low density and extremely low density according to the population density of the urban streets;
s42, assigning the extremely-high density area to be 5, assigning the high density area to be 4, assigning the medium density area to be 3, assigning the low density area to be 2, and assigning the extremely-low density area to be 1 to obtain the population density distribution map of the residential land of the central urban area after grading assignment;
s43, meshing the population density distribution map of the residential land of the central urban area after grading assignment by adopting ArcGIS software to obtain the average population density value of the unit grid;
in the present embodiment, the size of the unit cell in step S43 is 500m × 500 m.
And S44, obtaining the weighted value of the urban population density according to the average population density value of the unit grid.
S5, collecting urban meteorological data by a mobile meteorological station according to the urban greenbelt spore powder sensitization risk weighted value map to obtain a solar meteorological mean value;
step S5 includes the following substeps:
s51, setting an urban mobile weather station according to the urban greenbelt sporopollen sensitization risk weighted value map;
s52, collecting air temperature data, relative humidity data, wind speed data, rainfall data and illumination intensity data through the urban mobile weather station;
in this embodiment, the urban mobile weather station in step S52 may collect weather data in 24 hours per 2 days per month.
And S53, calculating the daily average value of the air temperature data, the relative humidity data, the wind speed data, the rainfall data and the illumination intensity data to obtain the daily weather average value of the 5 kinds of weather data.
And S6, calculating the urban greenbelt sporopollen sensitization risk value based on weather and population density according to the urban population density weighted value and the daily weather mean value.
Step S6 includes the following substeps:
s61, standardizing the urban population density weighted value, the spore powder sensitization risk value of various urban greenbelts and the solar meteorological mean value of 5 meteorological data to obtain 7 standardized data;
s62, calculating a component matrix of 7 kinds of standardized data by a KMO method and a Bartlett sphericity test method;
s63, calculating the eigenvalue of the component matrix to obtain the eigenvalue of 7 kinds of standardized data;
s64, calculating the square of the characteristic value of each standardized data to obtain the weight value of each standardized data;
s65, carrying out normalization processing on the weight value of each kind of standardized data to enable the sum of the weight values of the 7 kinds of standardized data to be 1, and obtaining the normalized weight values of the 7 kinds of standardized data;
and S66, calculating the sensitization risk value of the urban greenbelt spore powder based on weather and population density according to the normalized weight values of the 7 kinds of standardized data.
In step S66, the formula for calculating the sensitization risk value of the urban greenbelt spore powder based on weather and population density is as follows:
RUGZA＝a1·IUGZA+a2·x1+a3·x2+a4·x3+a5·x4+a6·x5+a7·x6
wherein R isUGZAIs an urban greenswamp spore powder sensitization risk value based on weather and population density, a1、a2、 a3、a4、a5、a6And a7Normalized weight values, x, for 7 normalized data1、x2、x3、x4、x5、 x6The corresponding urban population density weighted value and the mean value of the solar weather of 5 kinds of weather data, IUGZAThe spore powder sensitization risk value of various greenbelts in cities.
Claims (10)
1. A method for calculating an urban greenbelt sporopollen sensitization risk value based on weather and population density is characterized by comprising the following steps:
s1, calculating the allergenicity of each plant according to the plant characteristics of the urban green land;
s2, calculating the spore powder sensitization risk value of various urban greenbelts according to the sensitization force of each plant;
s3, gridding the distribution map of the spore powder sensitization risk values of various greenbelts in the city, and drawing a map of the spore powder sensitization risk weighted value of the greenbelts in the city;
s4, obtaining an urban population density weighted value based on the central urban residential land population density distribution map according to the population density of urban streets;
s5, collecting urban meteorological data by a mobile meteorological station according to the urban greenbelt spore powder sensitization risk weighted value map to obtain a solar meteorological mean value;
and S6, calculating the urban greenbelt sporopollen sensitization risk value based on weather and population density according to the urban population density weighted value and the daily weather mean value.
2. The method for calculating the sensitization risk value of the urban greenswamp spore powder according to the weather and the population density as claimed in claim 1, wherein the plant characteristics in the step S1 comprise: spore powder sensitization capacity, pollination type and flowering time.
3. The method for calculating the sensitization risk value of the urban greenswamp spore powder according to the weather and the population density as claimed in claim 2, wherein the step S1 comprises the following substeps:
s11, dividing the spore powder sensitization capability of the plants into 4 grades according to the spore powder particle diameter of the urban green land plants;
s12, assigning values of 0-3 to the sensitization capacity of 4 grades of sporopollen in sequence to obtain the sporopollen with the sporopollen particle diameter larger than 200 microns, wherein the sensitization capacity is assigned to be 0, the sensitization capacity is assigned to be 1, the sensitization capacity is assigned to be 2, the sporopollen particle diameter is smaller than 20 microns or is 50-100 microns, and the sensitization capacity is assigned to be 3, wherein the sporopollen particle diameter is 20-50 microns;
s13, assigning values to the pollination types by adopting values of 1-3 according to the pollination types of the urban green land plants to obtain the assignment value of the plants of the insect media pollination type as 1, the assignment values of the plants of the insect media pollination and wind media pollination types as 2 and the assignment values of the plants of the wind media pollination types as 3;
s14, measuring the flowering time required by 5% of flowers of urban green land plants to 95% of flowers to wither, assigning the flowering time of 1-3 weeks to 1, assigning the flowering time of 4-6 weeks to 2, and assigning the flowering time of 7 weeks or more to 3;
and S15, calculating the sensitization force of each plant according to the assignment of the sensitization capability of the sporopollen, the assignment of the pollination type and the assignment of the flowering time.
4. The method for calculating the sensitization risk value of the urban greenswamp spore powder according to the weather and the population density as claimed in claim 3, wherein the formula of the sensitization power of each plant in the step S15 is as follows:
VPA＝tp·dpp·ap
wherein VPA is the sensitization power of each plant, tp is the assignment of pollination type, dpp is the assignment of flowering time, and ap is the assignment of spore powder sensitization power.
5. The method for calculating the sensitization risk value of the spore powder of the urban greenbelt based on weather and population density according to claim 1, wherein the formula for calculating the sensitization risk value of the spore powder of each type of greenbelt of the urban in the step S2 is as follows:
wherein, IUGZAThe spore powder sensitization risk value of various urban greenbelts, VPA is the sensitization power of each plant, SiIs the floor area of the i-th species, HiHeight of species of type i, STIs the area of a certain urban green land, I is the total species class, and max VPA is the maximum value of VPA in all plants of a certain urban green land.
6. The method for calculating the sensitization risk value of the urban greenswamp spore powder according to the weather and the population density as claimed in claim 1, wherein the step S3 comprises the following substeps:
s31, gridding the distribution diagram of the spore powder sensitization risk values of various greenbelts in the city by adopting ArcGIS software to obtain the average spore powder sensitization risk value of unit grids in different seasons;
s32, dividing unit grids into low, medium and high sensitization risk grades according to the average spore powder sensitization risk value;
s33, assigning a value of 1 to the unit grid corresponding to the low sensitization risk level, assigning a value of 2 to the unit grid corresponding to the medium sensitization risk level, and assigning a value of 3 to the unit grid corresponding to the high sensitization risk level;
s34, calculating the seasonal spore powder sensitization risk weighted value of the unit grid by adopting a Sudoku weighting method according to the assignment of the unit grid;
and S35, drawing an urban greenbelt spore powder sensitization risk weighted value map based on the seasonal spore powder sensitization risk weighted value of the unit grid.
7. The method for calculating the sensitization risk value of the urban greenswamp spore powder according to the weather and the population density as claimed in claim 1, wherein the step S4 comprises the following substeps:
s41, dividing the population density distribution map of the residential land in the central urban area into 5 levels of extremely high density, medium density, low density and extremely low density according to the population density of the urban streets;
s42, assigning the extremely-high density area to be 5, assigning the high density area to be 4, assigning the medium density area to be 3, assigning the low density area to be 2, and assigning the extremely-low density area to be 1 to obtain the population density distribution map of the residential land of the central urban area after grading assignment;
s43, meshing the population density distribution map of the residential land of the central urban area after grading assignment by adopting ArcGIS software to obtain the average population density value of the unit grid;
and S44, obtaining the weighted value of the urban population density according to the average population density value of the unit grid.
8. The method for calculating the sensitization risk value of the urban greenswamp spore powder according to the weather and the population density as claimed in claim 1, wherein the step S5 comprises the following substeps:
s51, setting an urban mobile weather station according to the urban greenbelt sporopollen sensitization risk weighted value map;
s52, collecting air temperature data, relative humidity data, wind speed data, rainfall data and illumination intensity data through the urban mobile weather station;
and S53, calculating the daily average value of the air temperature data, the relative humidity data, the wind speed data, the rainfall data and the illumination intensity data to obtain the daily weather average value of the 5 kinds of weather data.
9. The method for calculating the sensitization risk value of the urban greenswamp spore powder according to the weather and the population density as claimed in claim 8, wherein the step S6 comprises the following substeps:
s61, standardizing the urban population density weighted value, the spore powder sensitization risk value of various urban greenbelts and the solar meteorological mean value of 5 meteorological data to obtain 7 standardized data;
s62, calculating a component matrix of 7 kinds of standardized data by a KMO method and a Bartlett sphericity test method;
s63, calculating the eigenvalue of the component matrix to obtain the eigenvalue of 7 kinds of standardized data;
s64, calculating the square of the characteristic value of each standardized data to obtain the weight value of each standardized data;
s65, carrying out normalization processing on the weight value of each kind of standardized data to enable the sum of the weight values of the 7 kinds of standardized data to be 1, and obtaining the normalized weight values of the 7 kinds of standardized data;
and S66, calculating the urban greenbelt spore powder sensitization risk value based on weather and population density according to the normalized weight values of the 7 kinds of standardized data.
10. The method for calculating the sensitization risk value of urban greenswamp spore powder based on weather and population density according to claim 9, wherein the formula for calculating the sensitization risk value of urban greenswamp spore powder based on weather and population density in step S66 is as follows:
RUGZA＝a1·IUGZA+a2·x1+a3·x2+a4·x3+a5·x4+a6·x5+a7·x6
wherein R isUGZAIs an urban greenswamp spore powder sensitization risk value based on weather and population density, a1、a2、a3、a4、a5、a6And a7Normalized weight values, x, for 7 normalized data1、x2、x3、x4、x5、x6The corresponding urban population density weighted value and the mean value of the solar weather of 5 kinds of weather data, IUGZAThe method is the spore powder sensitization risk value of various urban greenbelts.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110767674.4A CN113469537B (en) | 2021-07-07 | 2021-07-07 | Urban greenbelt sporopollen sensitization risk value calculation method based on weather and population density |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110767674.4A CN113469537B (en) | 2021-07-07 | 2021-07-07 | Urban greenbelt sporopollen sensitization risk value calculation method based on weather and population density |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113469537A true CN113469537A (en) | 2021-10-01 |
| CN113469537B CN113469537B (en) | 2023-04-07 |
Family
ID=77879218
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110767674.4A Active CN113469537B (en) | 2021-07-07 | 2021-07-07 | Urban greenbelt sporopollen sensitization risk value calculation method based on weather and population density |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113469537B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114707585A (en) * | 2022-03-22 | 2022-07-05 | 杭州睿胜软件有限公司 | Allergic species broadcasting method and system and readable storage medium |
| CN117253225A (en) * | 2023-09-05 | 2023-12-19 | 西南交通大学 | Sensitization risk judging method based on garden pollen morphology observation |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120066005A1 (en) * | 2010-09-10 | 2012-03-15 | State Farm Mutual Automobile Insurance Company | Systems and methods for grid-based insurance rating |
| CN106372807A (en) * | 2016-09-14 | 2017-02-01 | 国家测绘地理信息局第六地形测量队 | County region based resource environment carrying capacity evaluation method |
| WO2017032210A1 (en) * | 2015-08-24 | 2017-03-02 | 国家电网公司 | Cluster analysis based power transmission line mountain fire risky area division method |
| CN110738404A (en) * | 2019-09-29 | 2020-01-31 | 天津大学 | Improved ecological system ecological service value evaluation method |
-
2021
- 2021-07-07 CN CN202110767674.4A patent/CN113469537B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120066005A1 (en) * | 2010-09-10 | 2012-03-15 | State Farm Mutual Automobile Insurance Company | Systems and methods for grid-based insurance rating |
| WO2017032210A1 (en) * | 2015-08-24 | 2017-03-02 | 国家电网公司 | Cluster analysis based power transmission line mountain fire risky area division method |
| CN106372807A (en) * | 2016-09-14 | 2017-02-01 | 国家测绘地理信息局第六地形测量队 | County region based resource environment carrying capacity evaluation method |
| CN110738404A (en) * | 2019-09-29 | 2020-01-31 | 天津大学 | Improved ecological system ecological service value evaluation method |
Non-Patent Citations (1)
| Title |
|---|
| PALOMA CARI˜NANOS等: "Estimating the allergenic potential of urban green spaces: A case-study in Granada, Spain", 《LANDSCAPE AND URBAN PLANNING》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114707585A (en) * | 2022-03-22 | 2022-07-05 | 杭州睿胜软件有限公司 | Allergic species broadcasting method and system and readable storage medium |
| CN117253225A (en) * | 2023-09-05 | 2023-12-19 | 西南交通大学 | Sensitization risk judging method based on garden pollen morphology observation |
| CN117253225B (en) * | 2023-09-05 | 2024-08-27 | 西南交通大学 | A method for judging allergenic risk based on morphological observation of garden pollen |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113469537B (en) | 2023-04-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Han et al. | Prediction and evaluation of spatial distributions of ozone and urban heat island using a machine learning modified land use regression method | |
| Yang et al. | Quantifying air pollution removal by green roofs in Chicago | |
| Peper et al. | Equations for predicting diameter, height, crown width, and leaf area of San Joaquin Valley street trees | |
| Yilmaz et al. | Assessing the effects of different urban landscapes and built environment patterns on thermal comfort and air pollution in Erzurum city, Turkey | |
| Zhong et al. | Pedestrian-level gust wind flow and comfort around a building array–influencing assessment on the pocket park | |
| CN113469537B (en) | Urban greenbelt sporopollen sensitization risk value calculation method based on weather and population density | |
| CN115935855B (en) | Urban greening method and device based on optimized tree pollen concentration index | |
| NL2032402A (en) | Method for evaluating implementation effect of regional ecological protection and restoration project | |
| CN113740221B (en) | Air pollution source area determination method based on fine particulate matter and ozone | |
| KR20140021179A (en) | Solar radiation prediction method and apparatus | |
| CN114578457A (en) | Atmospheric pollutant concentration space-time prediction method based on evolution ensemble learning | |
| CN106845080A (en) | Scene Tourist meteorological disaster intelligent Forecasting based on difference amendment | |
| CN114757387A (en) | Forest fire danger state monitoring and analyzing system and method based on artificial intelligence | |
| CN115631077A (en) | Rainstorm disaster risk assessment method based on risk degree | |
| Dunichkin et al. | Studies and evaluation of bioclimatic comfort of residential areas for improving the quality of environment | |
| Perini et al. | Vertical greening systems for pollutants reduction | |
| CN116883215A (en) | A method for ecological and hydrological control of lakes and wetlands that considers the dynamic needs of migratory bird habitats | |
| CN118036891A (en) | Urban park greening construction planning evaluation system based on Internet of things | |
| CN106940830B (en) | Comprehensive evaluation technology for influences and risks of future climate change on biodiversity | |
| Yin et al. | Structure of an urban green space indirectly affects the distribution of airborne particulate matter: a study based on structural equation modelling | |
| Percy et al. | New exposure-based metric approach for evaluating O3 risk to North American aspen forests | |
| Cao et al. | Quantifying the vertical distribution pattern of land-use conversion in the loess hilly region of northern Shaanxi Province 1995–2015 | |
| CN118036889B (en) | Multi-time-space scale evolution feature analysis method and system for different types of drought | |
| Chen et al. | Spatial-temporal pattern of ecosystem services and sustainable development in representative mountainous cities: A case study of Chengdu-Chongqing Urban Agglomeration | |
| CN113642842A (en) | Department report based ecological system production value business accounting method |
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 |