CN112418596A - Biodiversity protection planning method - Google Patents
Biodiversity protection planning method Download PDFInfo
- Publication number
- CN112418596A CN112418596A CN202011084995.6A CN202011084995A CN112418596A CN 112418596 A CN112418596 A CN 112418596A CN 202011084995 A CN202011084995 A CN 202011084995A CN 112418596 A CN112418596 A CN 112418596A
- Authority
- CN
- China
- Prior art keywords
- planning
- granularity
- phylogenetic
- species
- range
- 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
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000009826 distribution Methods 0.000 claims abstract description 17
- 235000019580 granularity Nutrition 0.000 claims abstract 40
- 241000894007 species Species 0.000 claims description 82
- 241000422846 Sequoiadendron giganteum Species 0.000 claims description 19
- 238000013507 mapping Methods 0.000 claims description 11
- 230000011218 segmentation Effects 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 7
- 238000009877 rendering Methods 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 239000002245 particle Substances 0.000 claims 1
- 241000196324 Embryophyta Species 0.000 description 7
- 238000010276 construction Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 239000003086 colorant Substances 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000018109 developmental process Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000004040 coloring Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 238000007639 printing Methods 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 240000007087 Apium graveolens Species 0.000 description 1
- 235000015849 Apium graveolens Dulce Group Nutrition 0.000 description 1
- 235000010591 Appio Nutrition 0.000 description 1
- 240000003243 Thuja occidentalis Species 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 241000902900 cellular organisms Species 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 238000013467 fragmentation Methods 0.000 description 1
- 238000006062 fragmentation reaction Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 125000001997 phenyl group Chemical group [H]C1=C([H])C([H])=C(*)C([H])=C1[H] 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 230000001228 trophic effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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/0637—Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
- G06Q10/06375—Prediction of business process outcome or impact based on a proposed change
-
- 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
- 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
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Tourism & Hospitality (AREA)
- Economics (AREA)
- Educational Administration (AREA)
- Strategic Management (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Databases & Information Systems (AREA)
- General Business, Economics & Management (AREA)
- Development Economics (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Remote Sensing (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method for biodiversity protection planning, which comprises the following steps: 1) determining a planning range and a planning granularity according to a set biodiversity protection planning condition; each planning range comprises one or more planning granularities, a planning granularity list is established according to each obtained planning granularity, and the species to be protected in the corresponding planning granularity are determined according to the range of each planning granularity; 2) constructing a phylogenetic tree according to the species in each planning granularity to obtain a plurality of phylogenetic trees; 3) calculating the phylogenetic diversity index of each phylogenetic tree, and generating a diversity index key value pair in a planning range; 4) and calculating the species diversity distribution of each planning range according to the diversity index key values, and determining a priority protection area according to the species diversity distribution. The method can quickly calculate the phylogenetic diversity index distribution of small-granularity units and determine a reasonable planning result.
Description
Technical Field
The invention belongs to the technical field of applied ecological information, and relates to a large-scale biodiversity protection planning method based on phylogenetic diversity indexes, which is mainly applied to the fields of dynamic planning and ecological protection of national parks, protected areas and protected areas.
Background
Biodiversity is the sum of the characteristics and functions of all the species and their variability in a region, the genetic resources carried by each species, and the population (population), community (community), community (biome), phylum (biota), ecosystem (ecosystem) formed between species, etc., which are the basis for supporting the earth's life system. The natural conservation place refers to a land area or a sea area (Pengzongjing and the like, 2018) which implements long-term protection on important natural ecosystem, natural trails, natural landscapes and natural resources, ecological functions and cultural values borne by the natural trails.
Since the first natural protection area is established in 1956, natural protection lands in China are rapidly constructed and developed, 12000 places of various natural protection lands at all levels occupy about 18% of the land area, and play an important role in protecting biological diversity, guaranteeing ecological safety and the like in China. However, there are problems that the overall development strategy and planning of the natural protection area is lacking, the space layout of the natural protection area is unreasonable, the fragmentation and islanding phenomena of the natural protection area are serious, the protection effect is low, and the like, and the overall effect of the natural protection area is influenced (Xu, Weihua, et al.2019). In order to solve the problems, the nation steadily advances the reformation work of natural conservation ground under the background of ecological civilization construction, perfects the construction of a conservation ground system, proposes to establish a natural conservation ground system taking national parks as a main body, takes the establishment of national parks as a trigger, and solves the contradiction of natural conservation ground in China.
The definition of protection regions and the formulation of protection strategies are often based on the idea that all species have equal protection value (Vane-Wright et al, 1991; Posadas et al, 2001), however, the protection value is not the same among species due to historical differences (Faith, 1992, 2015; Laity et al,2015), and it is not comprehensive and accurate to evaluate and formulate protection strategies based solely on species diversity (Knapp et al, 2008; Daru & le Roux, 2016; wang et al, 2016). As phylogenetic diversity is not influenced by the change of classification status, the functions of the ecosystem (Cadotte et al, 2008; Flynn et al, 2011; Srivastava et al, 2012;) can be relatively well reflected, and the loss rate (Sechrest et al,2012) which is faster than the species diversity is added, so that the research and the evaluation of the phylogenetic diversity are more urgent. Precisely because of the fact that communities with high phylogenetic diversity are more stable, show higher productivity, and possess more species of different trophic levels (Cadotte et al,2012, 2013; Dinnange et al,2012), when establishing a preferential protection zone, attention is first paid to the zone with high phylogenetic diversity to protect as many evolutionary units as possible, thereby ensuring that the evolution history of the species within the region is effectively preserved (Sechrest et al, 2002; Jetz et al, 2012; shptot et al, 2015). In addition, phylogenetic diversity of the evolutionary dimension can clearly reveal branches with special evolutionary histories or endangered (Isaac et al, 2007; Rosauer et al,2009), which provides scientific basis for establishing a priority protection region or optimizing a protection strategy. Biodiversity protection is to protect not only the region rich in existing species, but also a region unique and rich in evolutionary history [ Chen end, Li & Rre, 2018 ].
In recent years, with the rapid development of DNA barcode data, new generation sequencing technology and computational tools, studies on the change, origin and diversity protection of plant region compositions based on the evolutionary relationship among species and depending on phylogenetic diversity and related parameters are developed, and studies on the distribution patterns and formation reasons of biodiversity in different regions by integrating environmental factors or climate factors based on the spatial distribution of species are also gradually carried out (e.g., Pio et al, 2014; Zhang et al, 2015). However, due to the limitation of research foundation, the construction level of species distribution models in different regions and the evaluation book of environmental factors are inconsistent, some earlier-developed researches and technical methods cannot be well integrated into the research of phylogenetic diversity, the sharing and integration of data still have great obstacles in China, and the distribution mode of researching biological diversity on a large scale is limited (xixixixiu celery & li, 2017).
Therefore, in practical applications, particularly in implementation of national parks and protected areas implemented in large-scale ranges and small-granularity units, how to quickly calculate a reliable phylogenetic tree, how to quickly calculate the phylogenetic diversity index distribution of small-granularity units, and how to quickly present planning results are key problems to be solved. The invention provides a method for carrying out biological diversity protection planning based on phylogenetic diversity index, and implements related schemes based on Chinese plant species distribution and Chinese land plant phylogenetic big trees.
Disclosure of Invention
Aiming at the technical problems in the prior art, the invention aims to provide a method for carrying out biological diversity protection planning based on a phylogenetic diversity index, which comprises the following specific processes:
step A: determining the biodiversity protection planning condition: carrying out condition setting on the planned range, the planned granularity and the like, determining the planned granularity individual and the corresponding species in the range, and then executing the step B;
(1) setting biodiversity protection conditions: setting a planned range, a planned granularity and the like;
the specification scope may be an administrative region type, such as province, city, or a combination thereof; it may also be manually selected, such as for a particular water area, mountain area, or other area established for planning.
The planning granularity can be administrative region types, such as county-level administrative districts and the like; it may also be a grid, such as a 10KM by 10KM grid or other grid granularity;
(2) screening and planning granularity individuals: calculating the inclusion range of all planning granularity individuals in the biodiversity planning enclosure to obtain a list (G1, G2, … Gi) containing each planning granularity individual, wherein i is the number of the planning granularity individual; gi is the area range of the ith planning granularity;
when the inclusion range of all planning granularity individuals in the biodiversity planning enclosure is calculated, the planning boundary processing method comprises the following steps: if the individual is based on the check area, the granularity individual at the planning boundary can be the planning strength individual if the individual is in the check area with the area more than or equal to 50 percent. And if the method is based on the administrative region, determining planning granularity individuals according to the division data of the latest administrative region of the country.
(3) Calculating the species in the planning granularity individual range: calculating all species names under the range of each planning granularity individual;
when each planning granularity individual contains all the names of the wild species in the range, the species distribution is based on the distribution data of the wild species collection area, and the relevant projection transformation can be carried out according to the selected planning granularity mode (such as a grid mode or an administrative district mode); the format of the output wild species name is 'generic name + species name'.
And B: and (3) construction of a phylogenetic tree: and constructing a phylogenetic tree for all species of each planning granularity individual, and constructing the phylogenetic tree by taking the planning granularity as a unit. Then executing step C;
(1) standardization of wild species directory: carrying out object name record standardization on all species names under the inclusion range of each planning granularity individual, and judging whether the species names are wild species or not; if the wild species is the local wild species, the wild species is reserved, and if the wild species is the non-local wild species, the wild species is deleted;
(2) matching node elements: inquiring a basic phylogenetic big tree database for a standardized wild species directory to obtain a node element list and a mapping relation between the species directory list and the node element list by matching (according to the steps from a low-order element to a high-order element, for example, matching is carried out by firstly 'planting' and then 'belonging' and then 'subject');
(3) acquiring a node element relation chain: retrieving the relationship chain of each node element from the basic phylogenetic big tree database to obtain the relationship chain of each node element;
(4) calculating node element grouping relation: pairwise comparison and matching are carried out on all the obtained relationship chains from leaf nodes to form a grouping relationship of the node elements;
(5) calculating the branch length of the node elements: calculating each node element and the corresponding new branch length of the group thereof in the phylogenetic tree according to the node weight and the species grouping relation of each node element in the node element list in the phylogenetic base big tree, and generating the phylogenetic tree corresponding to the node element list;
(6) generating a phylogenetic tree of planning granularity individuals: and generating a phylogenetic tree of planning granularity individuals according to the phylogenetic tree corresponding to the node element list and the mapping relation between the species list and the node element list.
And (5) repeating the steps (1) to (6) to obtain the phylogenetic tree with branch length of each planning granularity individual.
The preparation process of the basic phylogenetic big tree database comprises the following steps: firstly, segmenting and reconstructing a file of a system evolution basic big tree, recording segmentation times and segmentation weight, converting the segmentation times and the segmentation weight into a binary tree data structure, and recording a relationship chain of each node element and storing the relationship chain into a database.
And C: and (3) constructing a phylogenetic diversity index:
(1) calculating the individual diversity index of the planning granularity: firstly, respectively aiming at a phylogenetic tree of each species in planning granularity individual, calculating a phylogenetic diversity index: the phylogenetic tree for each planning granularity individual expresses the length of the species evolution history by the length of the branch, and the phylogenetic diversity of one planning granularity is the sum PDi of the branch lengths of all the species in the region on the phylogenetic tree.
(2) Generating a diversity index key value pair in a planning range: then summarizing to obtain a set of parameters including each planning granularity individual and a corresponding phylogenetic diversity index key value pair { Gi: PDi };
step D: biodiversity planning online mapping: and visually presenting the phylogenetic diversity index of each planning granularity in a geographic information system:
(1) establishing a planning granularity individual planning index: firstly, the planning index is standardized to obtain the planning index PLi of each planning granularity unit area. B, calculating to obtain the type of the planning granularity set according to the step A, and judging whether the area is a uniform area, if the area of the planning granularity is the equal area A, and standardizing each planning granularity individual and a key value pair of a planning index thereof, wherein the key value pair is { Gi: PLi } - { Gi: PDi/A }; if the areas of the planning granularity are not equal, the area Ai of each planning granularity individual is used for obtaining a normalized key value pair { Gi: PLi } - { Gi (PDi/Ai) } of each planning granularity individual and a planning index thereof;
(2) normalization process planning index: normalizing all the values of the PLi to obtain a { Gi: PLNi } key value pair, and finally obtaining the value of PLNi within the range of 0-1;
(3) and (3) establishing a rendering color: and calculating the color of each planning granularity individual in the gradient value according to a color gradient algorithm, and finally obtaining the mapping relation between each planning granularity individual and the corresponding RGB color to form a { Gi: RGbi } key value pair set.
(4) Constructing a planning granularity individual boundary service: b, constructing a vector map service of each individual planning granularity, and reading the individual boundary of each planning granularity according to the planning range and the planning granularity in the step A to form the map service;
(5) rendering into a graph: and according to the boundary range, performing color filling on each planning granularity individual according to { Gi: RGBi } coloring data, and rendering into a picture.
Step E: establishing a priority protection area: and according to different planning targets and coloring gradients, carrying out planning and establishment of the biodiversity protection area and generation of a corresponding data product. On one hand, a protection area can be determined and mapped according to the color distribution and the continuity degree of the coloring gradient through the visual map obtained in the step D; on the other hand, the proportion of the protection area in the planned area can be set, and a corresponding protection area is selected to form a diagram according to the value of PLi.
(1) Establishing a planning proportion: when the planning target is the protection area proportion S in the planning region, firstly, carrying out descending order arrangement according to the PLi to obtain a list of planning granularity arranged according to the PLi descending order;
(2) calculation of the priority protection list: then, the corresponding planning granularity is selected from the above to be added to obtain the total area Ab, so that the ratio of the total area Ab to the total planning area Az is greater than or equal to S, namely
Ak is the area of any one of k priority protection planning granularity individuals selected to meet the conditions, and Ai is the area of any one planning granularity individual in the planning range to obtain a priority protection planning granularity individual list;
(3) priority protection granular individual boundary services: reading boundary services of k priority protection planning granularity individuals, and performing highlighted color display on the boundaries;
(4) and (3) generation of a priority protection map product: and finally, generating a map data product of the established priority protection planning area, and printing the map data product into a png format data file for storage.
The invention also provides a server characterized by comprising a memory and a processor, the memory storing a computer program configured to be executed by the processor, the computer program comprising instructions for carrying out the steps of the above method.
The invention also provides a computer-readable storage medium, on which a computer program is stored, characterized in that the computer program realizes the steps of the above-mentioned method when being executed by a processor.
Compared with the prior art, the invention has the following positive effects:
the invention provides a method for carrying out biological diversity protection planning based on a phylogenetic diversity index, which is characterized in that a wild species phylogenetic tree is applied, a phylogenetic tree with planning granularity as a unit is quickly constructed, the phylogenetic diversity index is calculated, and online drawing of biological diversity planning is further carried out.
Drawings
FIG. 1 is a workflow diagram of a method for biodiversity protection planning based on phylogenetic diversity index;
FIG. 2 is a process flow diagram of a method of performing a biodiversity protection plan in an embodiment.
Detailed Description
The invention is further described by the following specific embodiments in conjunction with the accompanying drawings.
Referring to fig. 1, the specific process of the method for planning protection against biodiversity based on phylogenetic diversity index described in this example is:
step A: determining the biodiversity protection planning condition: carrying out condition setting on the planned range, the planned granularity and the like, determining the planned granularity individual and the corresponding species in the range, and then executing the step B;
(1) setting biodiversity protection conditions: the present embodiment is exemplified by a planning setting in which the planning range is yunnan province and the planning granularity is a county-level administrative district.
(2) Screening and planning granularity individuals: the inclusion ranges of all planning-granularity individuals in the biodiversity planning enclosure are calculated, and a list { G1, G2, … Gi }, i being the number of the planning-granularity individual, containing each planning-granularity individual is obtained, and the number is sorted by the first letter of the names of counties in Yunnan province, and the data is shown in table 1.
(3) Calculating the species in the planning granularity individual range: the names of all species in the range included in each county are calculated, and the output is output in the format of "generic name + species name" in table 2.
When calculating all the names of the wild species in the range contained by each planned granularity individual, the species distribution of the embodiment is based on the distribution data of the wild species collection place, the wild species distribution data used in the embodiment is sorted by Chinese plant journal, and the administrative region division is based on a national administrative region information query platform (http:// xzqh.mca.gov.cn/map); and carrying out regional projection transformation aiming at the historical distribution regional name related to Chinese plant annals and the existing county-level administrative districts.
And B: and (3) construction of a phylogenetic tree: and constructing a phylogenetic tree for all species of each planning granularity individual, and constructing the phylogenetic tree by taking the planning granularity as a unit. Then executing step C; in this example, phylogenetic trees of wild species were constructed in units of county-level row and district in Yunnan province.
(1) Standardization of wild species directory: standardizing the name of all species in the range of each county to judge whether the species are wild species (the judging condition of the embodiment is that the species are not included in the name of the species in Chinese plant journal, namely judge not to be the wild species [ the name of the plant species which can be cultivated or is foreign ]);
(2) matching node elements: the method comprises the steps of inquiring a basic Phylogenetic big tree database for a standardized wild species directory to obtain a node element list and a mapping relation between a species directory list and the node element list (matching is carried out according to the steps from a low-order element to a high-order element (such as 'kind of' then 'belongs to' and then 'family'), wherein the basic evolutionary big tree adopted by the embodiment is full coverage of a subordinate unit (Hu H, Liu B, Liang Y S, et al. an updated Chinese vascular plant tree of life: a marine diversity family visited [ J ]. Journal of Systematics and Evolution,2020.), and the obtained mapping relation is 'complete species matching' or 'typical genus matching');
(3) acquiring a node element relation chain: retrieving the relationship chain of each node element from the basic phylogenetic big tree database to obtain the relationship chain of each node element;
(4) calculating node element grouping relation: comparing and matching all the obtained relationship chains from leaf nodes to form a grouping relationship of the node elements;
(5) calculating the branch length of the node elements: calculating the new branch length of each node element and the corresponding group thereof in the phylogenetic development tree according to the node weight and the species grouping relation of each corresponding node element in the phylogenetic basic big tree;
(6) generating a phylogenetic tree: the phylogenetic trees with branch length obtained by the cyclic recursion steps (3) - (5) for each planning granularity individual are shown in table 3.
And (5) repeating the steps (1) to (6) to obtain the phylogenetic tree with branch length of each planning granularity individual.
The preparation process of the basic phylogenetic big tree database comprises the following steps: firstly, carrying out segmentation reconstruction on a file of a system evolution basic large tree, recording segmentation times and segmentation weights, converting the file into a binary tree data structure, and recording a relationship chain of each element to be stored in a database, wherein the specific method refers to SoTree (Meng Z, Li JH, Yang T, Lin L, et al. SoTree: An automated phenyl assembly sugar for each of the entities from big tree of the institute of technology on City, Chengdu 2015,792 797.);
and C: and (3) constructing a phylogenetic diversity index:
(1) calculating the individual diversity index of the planning granularity: firstly, respectively aiming at a phylogenetic tree of each species in planning granularity individual, calculating a phylogenetic diversity index: the phylogenetic tree for each planning granularity individual expresses the length of the species evolution history by the length of the branch, and the phylogenetic diversity of one planning granularity is the sum PDi of the branch lengths of all the species in the region on the phylogenetic tree.
(2) Generating a diversity index key value pair in a planning range: then summarizing to obtain a set of parameters including each planning granularity individual and a corresponding phylogenetic diversity index key value pair { Gi: PDi };
step D: biodiversity planning online mapping: and visually presenting the individual biodiversity index of each planning granularity in a geographic information system:
(1) establishing a planning granularity individual planning index: firstly, carrying out standardized calculation on a planning index PLi, acquiring the type of the planning granularity set in the step A, judging whether the planning granularity is a uniform area, if so, carrying out planning in a grid mode, wherein each standardized individual planning granularity and a key value pair of the planning index thereof are { Gi: PLi } - { Gi: PDi }; if the area is non-uniform, assuming that the area of each individual planning granularity is Ai, obtaining a key value pair { Gi: PLi } - { Gi (PDi/Ai) } of each individual planning granularity and the planning index thereof, where the data sample in this embodiment is shown in table 4, for example;
(2) normalization process planning index: normalizing all the values of the PLi to obtain a { Gi: PLNi } key value pair, wherein the value of PLNi is in the range of 0-1; normalized by PLNi ═ PLNi-min (pln))/(max (pln) -min (pln));
(3) and (3) establishing a rendering color: and calculating the corresponding color of each planning granularity individual in the color gradient according to a color gradient algorithm, and finally obtaining the RGB color mapping relation of each planning granularity individual to form a { Gi: RGBi } key value pair set.
Firstly, in the range of 0-1, taking several values according to gradient to form a set C ═ { C1, C2, C3, … and Cn } to respectively color, so as to gradually change in the colors;
then, determining the values of x and y according to the value of PLNi (x and y are from the set C and PLNi is between x and y), namely determining a color value interval [ x, y ];
and finally, calculating the color value corresponding to each Gi: the color corresponding to the minimum value of the interval [ x, y ] is represented by Colorx ═ rgb (rx, gx, bx), and the color corresponding to the maximum value y of the interval is represented by Colory ═ rgb (ry, gy, by), then the color corresponding to the PLNi in the interval is RGBi ═ rgb { rx + factor (ry-rx), gx + factor (gy-gx), bx + factor (by-bx) }, wherein the factor represents the relative position of the PLNi in the interval [ x, y ], and can be obtained by (PLNi-x)/(y-x).
In this embodiment, several values are taken in a gradient range of 0 to 1, such as C ═ 0,0.25,0.5,0.75,1, and the colors are given: rgb {44,123,182}, rgb {171,217,233}, rgb {255,255,191}, rgb {253,174,97}, rgb {215,25,28} to fade in these colors;
if the PLN value in bingchuan is 0.6174, then [ x, y ] ═ 0.5,0.75], then Colorx ═ rgb {255,255,191 }; color rgb {253,174,97 }; factor (0.6174-0.5)/(0.75-0.5) ═ 0.4696; then, the corresponding color value is RGB bingchuan county (RGB { round _ half _ up (255+0.4696 (253) 255)), round _ half _ up (255+0.4696 (174) 255)), round _ half _ up (191+0.4696 (97-191)) } (RGB (254,217,147) according to the formula. Where round _ half _ up is rounded up.
(4) Constructing a planning granularity individual boundary service: constructing vector map services of planning granularity individuals, reading each planning granularity individual boundary and a corresponding coloring value { Gi: RGbi } to form the map services according to the planning region and the planning granularity in the step A, and constructing the vector map county-level boundary services by applying mapnik V3(https:// mapnik.org);
(5) rendering into a graph: layers built by adding map services by Openlayers V6(https:// openlayers.org) can display gradient colorized planning layers at a browser end, and can be magnified or zoomed out for viewing at different resolution scales.
Step E: establishing a priority protection area: and according to different planning targets and coloring gradients, carrying out planning and establishment of the biodiversity protection area and generation of a corresponding data product. On one hand, a protection area can be determined and mapped according to the color distribution and the continuity degree of the coloring gradient through the visual map obtained in the step D; on the other hand, the proportion of the protection area in the planned area can be set, and a corresponding protection area is selected to form a diagram according to the value of PLi. The present embodiment generally describes the yunnan provincial protection area ratio as 18%, and includes the following steps:
(1) establishing a planning proportion: when the planning target is that the proportion S of the protection area in the planning region is 18%, firstly performing descending arrangement according to the PLi to obtain a list of planning granularity arranged according to the PLi in the descending arrangement;
(2) calculation of the priority protection list: then, the corresponding planning granularity is selected from the above to be added to obtain the total area Ab, so that the ratio of the total area Ab to the total planning area Az is more than or equal to 18 percent, namely
Ak is the area of any one of k priority protection planning granularity individuals selected to meet the conditions, and Ai is the area of any one planning granularity individual in the planning range to obtain a priority protection planning granularity individual list;
(3) priority protection granular individual boundary services: reading k priority protection planning granularity individuals, namely boundary services of priority protection counties, and performing highlighted color display on the boundaries;
(4) and (3) generation of a priority protection map product: and finally, generating a map data product of the established priority protection planning area, and printing the map data product into a png format data file for storage.
Table 1: output data example for setting planning range and planning granularity in embodiment
Table 2: within-individual species output data examples per planning granularity in an embodiment
Table 3: examples of phylogenetic tree data with branch lengths for each planning granularity individual
Table 4: in embodiments, PD output data examples within an individual per planning granularity
The above embodiments are only intended to illustrate the technical solution of the present invention and not to limit the same, and a person skilled in the art can modify the technical solution of the present invention or substitute the same without departing from the spirit and scope of the present invention, and the scope of the present invention should be determined by the claims.
Claims (8)
1. A method of biodiversity protection planning, comprising the steps of:
1) determining a planning range and a planning granularity according to a set biodiversity protection planning condition; each planning range comprises one or more planning granularities, a planning granularity list { G1, G2, … Gi } is established according to each obtained planning granularity, and the species to be protected in the corresponding planning granularity is determined according to the range of each planning granularity; wherein Gi is the region range of the ith planning granularity;
2) constructing a phylogenetic tree according to the species in each planning granularity to obtain i phylogenetic trees;
3) calculating the phylogenetic diversity index of each phylogenetic tree, and generating a diversity index key value pair in a planning range; wherein, the diversity index key value pair corresponding to the ith planning granularity Gi is { Gi: PDi }, and PDi is a phylogenetic diversity index calculated according to the phylogenetic tree of the planning granularity Gi;
4) and calculating the species diversity distribution of each planning range according to the diversity index key values, and determining a priority protection area according to the species diversity distribution.
2. The method of claim 1, wherein the phylogenetic tree is constructed by:
21) carrying out object name record standardization on the species names in the planned particle size range, and judging whether the species names are wild species or not; if the species is a local wild species, the species is reserved, otherwise, the species is deleted; obtaining a list of species lists; carrying out segmentation reconstruction on the file of the phylogenetic basic big tree, recording the segmentation times and the segmentation weight, converting the phylogenetic basic big tree into a binary tree data structure, and storing a relationship chain of each node element in the phylogenetic basic big tree into a basic phylogenetic big tree database;
22) querying the basic phylogenetic big tree database according to the standardized wild species name in the species directory list, matching to obtain a node element list, and recording the mapping relation between the species directory list and the node element list;
23) searching and acquiring a relationship chain of each node element in the node element list from the basic phylogenetic big tree database;
24) comparing the obtained relation chains pairwise, and if the two relation chains are matched, dividing the node elements corresponding to the two matched relation chains into a group to obtain the grouping relation of the node elements;
25) calculating each node element and the corresponding new branch length of the group thereof in the phylogenetic tree according to the node weight and the species grouping relation of each node element in the node element list in the phylogenetic base big tree, and generating the phylogenetic tree corresponding to the node element list;
26) and generating a phylogenetic tree corresponding to the planning granularity individuals according to the phylogenetic tree corresponding to the node element list and the mapping relation.
3. The method of claim 2, wherein said phylogenetic diversity index for each of said phylogenetic trees is calculated by: and expressing the length of species evolution history by the length of the branch in the phylogenetic tree, and calculating the sum PDi of the lengths of the branches on the phylogenetic tree as the phylogenetic diversity index of the phylogenetic tree.
4. The method of claim 2, wherein when pairwise aligning the resulting relationship chains, matching node elements in the two relationship chains occurs in order from a lower order element to a higher order element.
5. The method of claim 1, wherein the phylogenetic diversity index for each planning granularity is visualized in a geographic information system according to diversity index key assignments by:
41) calculating the planning index of each planning granularity unit area to obtain a standardized diversity index key value pair of each planning granularity; the normalized diversity index key value pair corresponding to the planning granularity Gi is { Gi: PLi }; pli is a planning index of a unit area of planning granularity Gi;
42) normalizing the normalized diversity index key value pair of each planning granularity;
43) calculating the color corresponding to each planning granularity in the gradient value according to a color gradient algorithm to obtain the mapping relation between each planning granularity and the RGB color and form a key value pair set; wherein, the key value pair { Gi: RGBi } corresponding to the planning granularity Gi in the key value pair set is RGBi which is the RGB color corresponding to the normalization value PlNi corresponding to the Pli;
44) reading each individual boundary of the planning granularity according to the planning range and the planning granularity to form vector map service corresponding to the planning granularity;
45) and according to the boundary range of each planning granularity, carrying out color filling on the individual corresponding to each planning granularity in the set according to the color-giving data corresponding to each planning granularity by the key value and rendering the individual into a picture.
6. The method as claimed in claim 1, wherein in step 4), the method for determining the preferential protection area comprises:
1) performing descending arrangement on the planning granularity according to the planning index of the unit area of the planning granularity to obtain a planning granularity list of the descending arrangement;
2) selecting k priority protection plan granularities according to the sequence and summing the areas to obtain a total area Ab, so that the ratio of the total area Ab to the total area Az of the plan range is greater than or equal to S, namely
Ak is the area of any one of k priority protection planning granularity individuals selected to meet the conditions, and Ai is the area of any one planning granularity individual in the planning range; s is a set protection area proportion in a planning area;
3) and reading the boundary services of the k priority protection planning granularity individuals, performing highlighted color display on the boundary, and establishing a priority protection area.
7. A server, comprising a memory and a processor, the memory storing a computer program configured to be executed by the processor, the computer program comprising instructions for carrying out the steps of the method according to any one of claims 1 to 6.
8. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method of any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011084995.6A CN112418596B (en) | 2020-10-12 | 2020-10-12 | Biodiversity protection planning method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011084995.6A CN112418596B (en) | 2020-10-12 | 2020-10-12 | Biodiversity protection planning method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112418596A true CN112418596A (en) | 2021-02-26 |
| CN112418596B CN112418596B (en) | 2022-10-11 |
Family
ID=74854404
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011084995.6A Active CN112418596B (en) | 2020-10-12 | 2020-10-12 | Biodiversity protection planning method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112418596B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116050935A (en) * | 2023-02-24 | 2023-05-02 | 北京师范大学 | Method and device for determining information of biodiversity priority protection area |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104573405A (en) * | 2014-12-22 | 2015-04-29 | 中国科学院计算机网络信息中心 | Phylogenetic tree rebuilding method for building sub trees on basis of big trees |
| CN106844688A (en) * | 2017-01-23 | 2017-06-13 | 环境保护部南京环境科学研究所 | The plant habitat analyzed based on space overlapping and GAP protects red line demarcation method |
| WO2019087373A1 (en) * | 2017-11-02 | 2019-05-09 | 富士通株式会社 | Animal habitat value assessment system, animal habitat value assessment method, and animal habitat value assessment program |
| CN110852532A (en) * | 2019-11-26 | 2020-02-28 | 浙江大学 | Natural protected area system integration and boundary delineation method |
-
2020
- 2020-10-12 CN CN202011084995.6A patent/CN112418596B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104573405A (en) * | 2014-12-22 | 2015-04-29 | 中国科学院计算机网络信息中心 | Phylogenetic tree rebuilding method for building sub trees on basis of big trees |
| CN106844688A (en) * | 2017-01-23 | 2017-06-13 | 环境保护部南京环境科学研究所 | The plant habitat analyzed based on space overlapping and GAP protects red line demarcation method |
| WO2019087373A1 (en) * | 2017-11-02 | 2019-05-09 | 富士通株式会社 | Animal habitat value assessment system, animal habitat value assessment method, and animal habitat value assessment program |
| CN110852532A (en) * | 2019-11-26 | 2020-02-28 | 浙江大学 | Natural protected area system integration and boundary delineation method |
Non-Patent Citations (1)
| Title |
|---|
| 慈秀芹等: "系统发育多样性在植物区系研究与生物多样性保护中的应用", 《生物多样性》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116050935A (en) * | 2023-02-24 | 2023-05-02 | 北京师范大学 | Method and device for determining information of biodiversity priority protection area |
| CN116050935B (en) * | 2023-02-24 | 2024-03-15 | 北京师范大学 | Method and device for determining information of biodiversity priority protection area |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112418596B (en) | 2022-10-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Schaefer et al. | Evaluation of land cover change and agricultural protection sites: a GIS and remote sensing approach for Ho Chi Minh City, Vietnam | |
| Hu et al. | A bird's view of new conservation hotspots in China | |
| Ghobadi et al. | Land suitability assessment (LSA) for aquaculture site selection via an integrated GIS-DANP multi-criteria method; a case study of lorestan province, Iran | |
| Aertsen et al. | Comparison and ranking of different modelling techniques for prediction of site index in Mediterranean mountain forests | |
| Tang et al. | Spatio‐temporal urban landscape change analysis using the Markov chain model and a modified genetic algorithm | |
| Chandan et al. | Analysing spatial patterns and trend of future urban expansion using SLEUTH | |
| Faith et al. | Practical application of biodiversity surrogates and percentage targets for conservation in Papua New Guinea | |
| Behera et al. | Rapid assessment of biological richness in a part of Eastern Himalaya: an integrated three-tier approach | |
| Wali et al. | Development of a land suitability model for saffron (Crocus sativus L.) cultivation in Khost Province of Afghanistan using GIS and AHP techniques | |
| Li et al. | Assessing regional sustainability: the case of land use and land cover change in the middle Yiluo catchment of the Yellow River basin, China | |
| Johnston et al. | A unifying approach for evaluating the condition of wetland plant communities and identifying related stressors | |
| Heumann et al. | Land suitability modeling using a geographic socio-environmental niche-based approach: A case study from northeastern Thailand | |
| Gao et al. | Exploring detailed urban-rural development under intersecting population growth and food production scenarios: Trajectories for China’s most populous agricultural province to 2030 | |
| Lanz et al. | Sampling design and estimation procedures | |
| Yohannes et al. | Integration of remote sensing, GIS and MCDM for land capability classification in Andit Tid Watershed, Ethiopia | |
| Huang et al. | Lizard species richness patterns in China and its environmental associations | |
| Zampieri et al. | Simple statistical approach for computing land cover types and potential natural vegetation | |
| CN112418596B (en) | Biodiversity protection planning method | |
| Shapira et al. | Combining analytical hierarchy process and agglomerative hierarchical clustering in search of expert consensus in green corridors development management | |
| Liu et al. | Urban green infrastructure affects bird biodiversity in the coastal megalopolis region of Shenzhen city | |
| Kalichkin et al. | GIS-based multi-criteria analysis of the suitability of Western Siberian forest-steppe lands | |
| Bodesmo et al. | The role of environmental and socio-demographic indicators in the analysis of land use changes in a protected area of the Natura 2000 Network: the case study of Lake Trasimeno, Umbria, Central Italy | |
| Liu et al. | Security patterns and resistance surface model in urban development: Case study of Sanshui, China | |
| Flores-Tolentino et al. | Biogeographic regionalization by spatial and environmental components: Numerical proposal | |
| Leimgruber et al. | Rural areas between regional needs and global challenges |
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 |