CN113469537B - 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
- CN113469537B CN113469537B CN202110767674.4A CN202110767674A CN113469537B CN 113469537 B CN113469537 B CN 113469537B CN 202110767674 A CN202110767674 A CN 202110767674A CN 113469537 B CN113469537 B CN 113469537B
- Authority
- CN
- China
- Prior art keywords
- urban
- sensitization
- value
- sporopollen
- population density
- 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
- 206010070834 Sensitisation Diseases 0.000 title claims abstract description 138
- 230000008313 sensitization Effects 0.000 title claims abstract description 138
- 238000004364 calculation method Methods 0.000 title description 7
- 239000000843 powder Substances 0.000 claims abstract description 67
- 238000000034 method Methods 0.000 claims abstract description 22
- 241000196324 Embryophyta Species 0.000 claims description 65
- 230000010152 pollination Effects 0.000 claims description 24
- 230000017260 vegetative to reproductive phase transition of meristem Effects 0.000 claims description 21
- 239000002245 particle Substances 0.000 claims description 13
- 241000894007 species Species 0.000 claims description 12
- 238000005286 illumination Methods 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000001932 seasonal effect Effects 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000010998 test method Methods 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 description 9
- 230000036541 health Effects 0.000 description 6
- 206010020751 Hypersensitivity Diseases 0.000 description 5
- 208000026935 allergic disease Diseases 0.000 description 4
- 230000007815 allergy Effects 0.000 description 4
- 238000010276 construction Methods 0.000 description 2
- 230000033001 locomotion Effects 0.000 description 2
- 238000004519 manufacturing process Methods 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
- 201000010105 allergic rhinitis Diseases 0.000 description 1
- 208000006673 asthma 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)
- Data Mining & Analysis (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Databases & Information Systems (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- General Business, Economics & Management (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Game Theory and Decision Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Quality & Reliability (AREA)
- Probability & Statistics with Applications (AREA)
- Evolutionary Biology (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Algebra (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Software Systems (AREA)
- Primary Health Care (AREA)
- Remote Sensing (AREA)
- Cultivation Of Plants (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method for calculating an urban greenbelt sporopollen sensitization risk value based on weather and population density, which comprises the following steps of: s1, calculating the allergenicity of each plant according to the plant characteristics of urban green lands; s2, calculating the spore powder sensitization risk values of various urban greenbelts according to the sensitization force of each plant; s3, gridding distribution graphs of spore powder sensitization risk values of various greenbelts in the city, and drawing a map of spore powder sensitization risk weighted values of the greenbelts in the city; s4, according to the population density of the urban street, based on the population density distribution map, obtaining an urban population density weighted value; s5, acquiring meteorological data of the city by adopting a mobile meteorological station to obtain a solar meteorological mean value; 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; the invention solves the problem that the method for calculating the sensitization risk value of the urban greenbelt sporopollen in the prior art is lacked.
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 by nature and human and forms a human living healthy space. 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 is focused on the positive benefits of garden plants to cities, such as ornamental value, city air purification, microclimate regulation, 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 develop allergic reactions 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 genuine epidemic disease. However, the urban greenbelt sporopollen sensitization risk assessment and prediction have not been carried out in China. 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 sensitization of each plant according to the plant characteristics of urban green land;
s2, calculating the pollen sensitization risk value of various greenbelts in the city according to the sensitization force of each plant;
s3, gridding distribution graphs of spore powder sensitization risk values of various greenbelts in the city, and drawing a map of spore powder sensitization risk weighted values of the greenbelts in the city;
s4, according to the population density of the urban street, based on the population density distribution map of the residential land of the central urban area, obtaining an urban population density weighted value;
s5, acquiring weather data of the city by adopting a mobile weather station according to the map of the sensitization risk weighted value of the spore powder of the green land of the city to obtain a solar weather 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 comprises the following substeps:
s11, dividing the spore powder sensitization capability of the plants into 4 grades according to the diameter of the spore powder particles 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 particle diameter of the sporopollen of the plant larger than 200 microns, assigning the sensitization capacity of the sporopollen with the particle diameter of the sporopollen of the plant of 100-200 microns as 0, assigning the sensitization capacity of the sporopollen with the particle diameter of the sporopollen of the plant smaller than 20 microns or between 50-100 microns as 2, and assigning the sensitization capacity of the sporopollen with the particle diameter of the sporopollen of the plant between 20-50 microns as 3;
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-borne pollination type as 1, the assignment value of the plants of the insect-borne and wind-borne pollination types as 2 and the assignment value of the plants of the wind-borne pollination type 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: is beneficial to understanding the allergenicity of various plants in urban green land and provides a basis for the safety and scientificity of urban green land species selection.
Further, the formula for calculating the sporopollen sensitization risk value of each type of green land in the city in the step S2 is as follows:
wherein, I UGZA The spore powder sensitization risk value of various urban greenbelts, VPA is the sensitization power of each plant, S i Is the floor area of the i-th species, H i Height of species of type i, S T Is 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 performed on the basis, so that the calculation of the sensitization risk of the green land is more comprehensive and scientific.
Further, step S3 comprises the following sub-steps:
s31, gridding distribution graphs of sporopollen sensitization risk values of various greenbelts in cities by adopting ArcGIS software to obtain average sporopollen sensitization risk values of unit grids in different seasons;
s32, dividing unit grids into low, medium and high sensitization risk levels according to the average sporopollen 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 in 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 of 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 region to be 5, assigning the high density region to be 4, assigning the medium density region to be 3, assigning the low density region to be 2, and assigning the extremely low density region to be 1 to obtain a 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 subjected to grading assignment by adopting ArcGIS software to obtain the average population density value of the unit grid;
and S44, obtaining a city population density weighted value 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 users, if population density of a region where greenbelts are 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 an 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 an urban mobile weather station;
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 greenbelt all influence the production and the movement of greenbelt sporopollen, and the scheme is beneficial to estimating the production quantity and the movement track of the greenbelt sporopollen more scientifically.
Further, step S6 includes the following substeps:
s61, standardizing the weighted value of the population density of the city, the spore powder sensitization risk values of various greenbelts of the city and the daily 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 type of standardized data to enable the sum of the weight values of the 7 types of standardized data to be 1, and obtaining the normalized weight values of the 7 types 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.
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:
R UGZA ＝a 1 ·I UGZA +a 2 ·x 1 +a 3 ·x 2 +a 4 ·x 3 +a 5 ·x 4 +a 6 ·x 5 +a 7 ·x 6
wherein R is UGZA Is an urban greenswamp spore powder sensitization risk value based on weather and population density, a 1 、a 2 、a 3 、a 4 、a 5 、a 6 And a 7 Normalized weight values, x, for 7 normalized data 1 、x 2 、x 3 、x 4 、x 5 、x 6 The corresponding urban population density weighted value and the mean value of the solar weather of 5 kinds of weather data, I UGZA For various green areas of citiesPollen sensitization risk value.
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: according to the method, the allergenicity of each plant is calculated, so that the spore powder allergenicity risk values of various greenbelts of a city are obtained through calculation, a city greenbelt spore powder allergenicity risk weighted value map is drawn, and then the city greenbelt spore powder allergenicity risk values based on weather and population density are obtained through calculation by combining the city population density weighted value and the weather data of the city.
Drawings
FIG. 1 is a flow chart of a method for calculating an urban greenswamp spore powder sensitization risk value 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 urban green lands;
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 diameter of the spore powder particles of the urban green land plants;
s12, assigning values to the sensitization capacity of the sporopollen of 4 grades by adopting values of 0-3 in sequence to obtain a sporopollen with the sporopollen particle diameter larger than 200 microns of the plant, wherein the sensitization capacity is assigned to 0, the sensitization capacity is assigned to 1, the sensitization capacity is assigned to 2, and the sensitization capacity is assigned to 3, wherein the sporopollen particle diameter of the plant is smaller than 20 microns or is between 50 microns and 100 microns, and the sensitization capacity is assigned to 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-borne pollination type as 1, the assignment value of the plants of the insect-borne and wind-borne pollination types as 2 and the assignment value of the plants of the wind-borne pollination type 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 values of various urban greenbelts according to the sensitization force of each plant;
in the step S2, the formula for calculating the spore powder sensitization risk value of various greenbelts in the city is as follows:
wherein, I UGZA The spore powder sensitization risk value of various urban greenbelts, VPA is the sensitization power of each plant, S i Is the floor area of the i-th species, H i Height of species of type i, S T Is 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 distribution graphs of spore powder sensitization risk values of various greenbelts in the city, and drawing a map of spore powder sensitization risk weighted values of the greenbelts in the city;
step S3 comprises the following substeps:
s31, gridding distribution graphs of spore powder sensitization risk values of various greenbelts in cities by adopting ArcGIS software to obtain average spore powder sensitization risk values of unit grids in different seasons;
in the present embodiment, the size of the unit cell in step S31 is 500m × 500m.
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 a unit grid corresponding to a low sensitization risk level, assigning a value of 2 to a unit grid corresponding to a medium sensitization risk level, and assigning a value of 3 to a unit grid corresponding to a 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, according to the population density of the urban street, based on the population density distribution map of the residential land of the central urban area, obtaining an urban population density weighted value;
step S4 includes the following substeps:
s41, dividing the population density distribution map of the residential land of 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 a 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 subjected to 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 × 500m.
And S44, obtaining an urban population density weighted value according to the average population density value of the unit grid.
S5, acquiring meteorological data of the city by adopting 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 an 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 an urban mobile weather station;
in this embodiment, the urban mobile weather station in step S52 may collect the weather data in 24 hours every 2 days every month.
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 weighted value of the population density of the city, the spore powder sensitization risk values of various greenbelts of the city and the daily 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 characteristic values of the component matrix to obtain the characteristic values 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 type of standardized data to enable the sum of the weight values of the 7 types of standardized data to be 1, and obtaining the normalized weight values of the 7 types 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.
In the step S66, the formula for calculating the sensitization risk value of the urban greenbelt sporopollen based on weather and population density is as follows:
R UGZA ＝a 1 ·I UGZA +a 2 ·x 1 +a 3 ·x 2 +a 4 ·x 3 +a 5 ·x 4 +a 6 ·x 5 +a 7 ·x 6
wherein R is UGZA Is an urban greenswamp spore powder sensitization risk value based on weather and population density, a 1 、a 2 、a 3 、a 4 、a 5 、a 6 And a 7 Normalized weight values, x, for 7 normalized data 1 、x 2 、x 3 、x 4 、x 5 、x 6 The corresponding urban population density weighted value and the mean value of the solar weather of 5 kinds of weather data, I UGZA The method is the spore powder sensitization risk value of various urban greenbelts.
Claims (1)
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 urban green lands; the plant characteristics in the step S1 comprise: pollen sensitization ability, pollination type and flowering time;
the step S1 comprises the following sub-steps:
s11, dividing the spore powder sensitization capability of the plants into 4 grades according to the diameter of the spore powder particles 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 particle diameter of the sporopollen of the plant larger than 200 microns, assigning the sensitization capacity of the sporopollen with the particle diameter of the sporopollen of the plant of 100-200 microns as 0, assigning the sensitization capacity of the sporopollen with the particle diameter of the sporopollen of the plant smaller than 20 microns or between 50-100 microns as 2, and assigning the sensitization capacity of the sporopollen with the particle diameter of the sporopollen of the plant between 20-50 microns as 3;
s13, assigning 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 insect-borne pollination type plant as 1, the assignment values of the insect-borne and wind-borne pollination types plant as 2 and the assignment value of the wind-borne pollination type plant 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;
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·ao
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 capacity;
s2, calculating the spore powder sensitization risk values of various urban greenbelts according to the sensitization force of each plant;
the formula for calculating the spore powder sensitization risk value of various greenbelts in the city in the step S2 is as follows:
wherein, I UGZA Is the spore powder sensitization risk value of various urban greenbelts, VPA is the sensitization power of each plant,S i is the floor area of the i-th species, H i Height of species of type i, S T The area of a certain urban green land, I is the total species category, and max VPA is the maximum value of VPA in all plants of a certain urban green land;
s3, gridding distribution graphs of spore powder sensitization risk values of various greenbelts in the city, and drawing a map of spore powder sensitization risk weighted values of the greenbelts in the city;
the step S3 comprises the following sub-steps:
s31, gridding distribution graphs of sporopollen sensitization risk values of various greenbelts in cities by adopting ArcGIS software to obtain average sporopollen sensitization risk values 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 a unit grid corresponding to a low sensitization risk level, assigning a value of 2 to a unit grid corresponding to a medium sensitization risk level, and assigning a value of 3 to a unit grid corresponding to a 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;
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, according to the population density of the urban street, based on the population density distribution map of the residential land of the central urban area, obtaining an urban population density weighted value;
the step S4 comprises the following sub-steps:
s41, dividing the population density distribution map of the residential land of 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 a 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 in the central urban area after grading assignment by adopting ArcGIS software to obtain an average population density value of a unit grid;
s44, obtaining a weighted value of the urban population density according to the average population density value of the unit grid;
s5, acquiring meteorological data of the city by adopting a mobile meteorological station according to the urban greenbelt spore powder sensitization risk weighted value map to obtain a solar meteorological mean value;
the step S5 comprises the following sub-steps:
s51, setting an urban mobile weather station according to an 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 an urban mobile weather station;
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;
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;
the step S6 comprises the following sub-steps:
s61, standardizing the weighted value of the population density of the city, the spore powder sensitization risk values of various greenbelts of the city and the daily 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 characteristic values of the component matrix to obtain the characteristic values 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 type of standardized data to enable the sum of the weight values of the 7 types of standardized data to be 1, and obtaining the normalized weight values of the 7 types of standardized data;
s66, calculating an 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;
in the step S66, the formula for calculating the urban greenbelt sporopollen sensitization risk value based on weather and population density is as follows:
R UGZA ＝a 1 ·I UGZA +a 2 ·x 1 +a 3 ·x 2 +a 4 ·x 3 +a 5 ·x 4 +a 6 ·x 5 +a 7 ·x 6
wherein R is UGZA Is an urban greenswamp spore powder sensitization risk value based on weather and population density, a 1 、a 2 、a 3 、a 4 、a 5 、a 6 And a 7 Normalized weight values, x, for 7 normalized data 1 、x 2 、x 3 、x 4 、x 5 、x 6 The corresponding urban population density weighted value and the mean value of the solar weather of 5 kinds of weather data, I UGZA The 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 CN113469537A (en) | 2021-10-01 |
| CN113469537B true 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) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114707585B (en) * | 2022-03-22 | 2024-08-02 | 杭州睿胜软件有限公司 | Allergy species broadcasting method, system and readable storage medium |
| CN117253225B (en) * | 2023-09-05 | 2024-08-27 | 西南交通大学 | Sensitization risk judging method based on garden pollen morphology observation |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8504393B2 (en) * | 2010-09-10 | 2013-08-06 | State Farm Mutual Automobile Insurance Company | Systems and methods for grid-based insurance rating |
| CN105184668A (en) * | 2015-08-24 | 2015-12-23 | 国家电网公司 | Forest fire risk area dividing method for power transmission line based on cluster analysis |
-
2021
- 2021-07-07 CN CN202110767674.4A patent/CN113469537B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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ñanos等.Estimating the allergenic potential of urban green spaces: A case-study in Granada, Spain.《Landscape and Urban Planning》.2014,第134-144页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113469537A (en) | 2021-10-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Diener et al. | How can vegetation protect us from air pollution? A critical review on green spaces' mitigation abilities for air-borne particles from a public health perspective-with implications for urban planning | |
| CN113469537B (en) | Urban greenbelt sporopollen sensitization risk value calculation method based on weather and population density | |
| Yang et al. | Quantifying air pollution removal by green roofs in Chicago | |
| Jamei et al. | Urban development and pedestrian thermal comfort in Melbourne | |
| Han et al. | Prediction and evaluation of spatial distributions of ozone and urban heat island using a machine learning modified land use regression method | |
| WO2024098444A1 (en) | Rural ecosystem carbon stock prediction method | |
| Zhang et al. | Continuous monitoring, compositions analysis and the implication of regional transport for submicron and fine aerosols in Beijing, China | |
| Yilmaz et al. | Assessing the effects of different urban landscapes and built environment patterns on thermal comfort and air pollution in Erzurum city, Turkey | |
| Tan et al. | Factors influencing street tree health in constrained planting spaces: Evidence from Kyoto City, Japan | |
| CN115935855B (en) | Urban greening method and device based on optimized tree pollen concentration index | |
| Ksiazek-Mikenas et al. | Traits for stress-tolerance are associated with long-term plant survival on green roofs | |
| CN114565223A (en) | Method for evaluating implementation effect of regional ecological protection restoration project | |
| Chauhan et al. | Assessing inherent vulnerability of farming communities across different biogeographical zones in Himachal Pradesh, India | |
| Dunichkin et al. | Studies and evaluation of bioclimatic comfort of residential areas for improving the quality of environment | |
| CN115631077A (en) | Rainstorm disaster risk assessment method based on risk degree | |
| Zhou et al. | Inter-annual variation in grassland net ecosystem productivity and its coupling relation to climatic factors in China | |
| Yang et al. | General air temperature and humidity features of local climate zones: a multi-city observational study in eastern China | |
| Cao et al. | Quantifying the vertical distribution pattern of land-use conversion in the loess hilly region of northern Shaanxi Province 1995–2015 | |
| Ersoy Tonyaloğlu et al. | Impact of land cover change on urban tree cover and potential regulating ecosystem services: The case of Aydın/Turkey | |
| Lovynska et al. | Study on accumulation of heavy metals by green plantations in the conditions of industrial cities | |
| Chen et al. | Spatial-temporal pattern of ecosystem services and sustainable development in representative mountainous cities: A case study of Chengdu-Chongqing Urban Agglomeration | |
| Wang et al. | Spatio-temporal Evaluation of Multi-scale Cultivated Land System Resilience in Black Soil Region from 2000 to 2019: A Case Study of Liaoning Province, Northeast China | |
| Song et al. | Coupled water-habitat-carbon nexus and driving mechanisms in the Tarim River Basin: A multi-scenario simulation perspective | |
| Madushanka et al. | Strategies to promote greenery in urban boundary wall facades: A case study in residential areas of Colombo district, Sri Lanka | |
| CHEN et al. | Effects of spatial forms of green infrastructure in block scale on PM10 and PM2. 5 removal—a case study of the main city of Wuhan |
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 |