CN113469512A - Natural protected area determining method and device, electronic device and storage medium - Google Patents

Natural protected area determining method and device, electronic device and storage medium Download PDF

Info

Publication number
CN113469512A
CN113469512A CN202110693013.1A CN202110693013A CN113469512A CN 113469512 A CN113469512 A CN 113469512A CN 202110693013 A CN202110693013 A CN 202110693013A CN 113469512 A CN113469512 A CN 113469512A
Authority
CN
China
Prior art keywords
data
land
natural
habitat
analyzed
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
Application number
CN202110693013.1A
Other languages
Chinese (zh)
Other versions
CN113469512B (en
Inventor
单樑
唐文魁
岳隽
陈君丽
袁富林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Urban Planning And Design Institute Co ltd
Original Assignee
URBAN PLANNING & DESIGN INSTITUTE OF SHENZHEN (UPDIS)
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by URBAN PLANNING & DESIGN INSTITUTE OF SHENZHEN (UPDIS) filed Critical URBAN PLANNING & DESIGN INSTITUTE OF SHENZHEN (UPDIS)
Priority to CN202110693013.1A priority Critical patent/CN113469512B/en
Publication of CN113469512A publication Critical patent/CN113469512A/en
Application granted granted Critical
Publication of CN113469512B publication Critical patent/CN113469512B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06395—Quality analysis or management
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • 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—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/10—Services
    • G06Q50/26—Government or public services
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05—Geographic models
    • 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
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Human Resources & Organizations (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Software Systems (AREA)
  • Strategic Management (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Tourism & Hospitality (AREA)
  • Economics (AREA)
  • Geometry (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Business, Economics & Management (AREA)
  • Evolutionary Computation (AREA)
  • Marketing (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Graphics (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Remote Sensing (AREA)
  • Medical Informatics (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a natural protected area defining method, a natural protected area defining device, electronic equipment and a storage medium, and relates to the technical field of ecological planning, wherein the natural protected area defining method comprises the following steps: acquiring land data to be analyzed; extracting land species distribution data from land data to be analyzed, and extracting threat source data from the land data to be analyzed; calculating habitat quality evaluation data by using a preset habitat quality model, soil species distribution data and threat source data; calculating natural protection area data by using a preset intelligent agent evolution model and the habitat quality evaluation data; and generating the natural protection area according to the natural protection area data and the habitat quality evaluation data. The natural protection area dividing method can simulate, predict, optimize and display the ecological suitability evaluation result to indicate the species to divide the natural protection area in the space gathering process, and makes up for the functional deficiency of the conventional cellular automata modeling in the simulation and optimization of the complex geographic space-time process.

Description

Natural protected area determining method and device, electronic device and storage medium
Technical Field
The present invention relates to the field of ecological planning technologies, and in particular, to a method and an apparatus for defining a natural protected area, an electronic device, and a storage medium.
Background
Cellular automata is a common modeling mode for spatial pattern evolution, but the description of the driving force of the geographic environment is single and fixed, the cellular automata is easily influenced by historical rule inertia driving, and the subjectivity of describing spatial evolution participation main bodies is lacked.
Disclosure of Invention
The present invention is directed to solving at least one of the problems of the prior art. Therefore, the embodiment of the invention provides a natural protection area dividing method which can divide the natural protection area by simulating, predicting, optimizing and displaying the ecological suitability evaluation result in the space gathering process of indication species, and makes up for the functional deficiency of the conventional cellular automata modeling in the simulation and optimization of the complex geographical time-space process.
The natural reserve determination method according to the first aspect of the present invention includes:
acquiring land data to be analyzed;
extracting land species distribution data from the land data to be analyzed, and extracting threat source data from the land data to be analyzed;
calculating birth environment quality evaluation data by using a preset habitat quality model, the soil species distribution data and the threat source data;
calculating natural protection area data by using a preset intelligent agent evolution model and the habitat quality evaluation data;
and generating a natural protection area according to the natural protection area data and the habitat quality evaluation data.
According to the natural reserve determination method of the first aspect of the present invention, at least the following advantageous effects are obtained: firstly, land data to be analyzed is obtained, land species distribution data are extracted from the land data to be analyzed, threat source data are extracted from the land data to be analyzed, a preset habitat quality model is used for calculating habitat quality evaluation data for the soil species distribution data and the threat source data, a preset intelligent agent evolution model and the habitat quality evaluation data are used for calculating natural protection area data, a natural protection area is generated according to the natural protection area data and the habitat quality evaluation data, and the space gathering process of the indicated species can be simulated, predicted, optimized and displayed based on the ecological suitability evaluation result, the gathering area of the indicated species is rapidly displayed in the map, natural protection area division is realized, and the functional deficiency of the conventional cellular automata modeling in the simulation and optimization of the complex geographic time-space process is overcome.
According to some embodiments of the invention, the obtaining land data to be analyzed comprises: acquiring land utilization data and DEM elevation data; and performing correlation preprocessing on the land utilization data and the DEM elevation data to obtain the land data to be analyzed.
According to some embodiments of the invention, the extracting of the land feature distribution data from the land data to be analyzed comprises: and carrying out species distribution extraction on the land data to be analyzed to obtain the land species distribution data.
According to some embodiments of the invention, the extracting threat source data from the land data to be analyzed comprises: drawing up threat factor data according to the land data to be analyzed; acquiring a threat factor type corresponding to the threat factor data; and carrying out assignment processing on the threat factor data according to the threat factor type to obtain the threat source data.
According to some embodiments of the invention, the calculating of the habitat quality evaluation data using the preset habitat quality model, the soil species distribution data and the threat source data comprises: substituting the soil species distribution data and the threat source data into the preset habitat quality model to calculate a habitat degradation index and a habitat quality evaluation index; and obtaining the habitat quality evaluation data according to the habitat degradation index and the habitat quality evaluation index.
According to some embodiments of the present invention, the calculating using the preset agent evolution model and the habitat quality evaluation data to obtain natural reserve data includes: acquiring ecological protection target parameters and utility functions; and performing regional optimization calculation according to the preset intelligent agent evolution model, the ecological protection target parameters and the utility function to obtain the natural protection area data.
According to some embodiments of the invention, the generating a natural reserve based on the natural reserve data and the habitat quality evaluation data comprises: generating natural protection area defining data according to the natural protection area data and the habitat quality evaluation data; and loading data according to the natural protection area to obtain the natural protection area.
According to a second aspect of the present invention, a natural preservation area delineation apparatus comprises:
the acquisition module is used for acquiring land data to be analyzed;
the extraction module is used for extracting land species distribution data from the land data to be analyzed and extracting threat source data from the land data to be analyzed;
the habitat quality evaluation module is used for calculating habitat quality evaluation data by utilizing a preset habitat quality model, the soil species distribution data and the threat source data;
the intelligent agent evolution module is used for calculating natural protection area data by utilizing a preset intelligent agent evolution model and the habitat quality evaluation data;
and the demarcation module is used for generating a natural protection area according to the natural protection area data and the habitat quality evaluation data.
The natural protection area defining device according to the embodiment of the second aspect of the present invention has at least the following beneficial effects: by executing the natural protected area determining method provided by the embodiment of the first aspect of the invention, the aggregation process of the indication species in the space can be simulated, predicted, optimized and displayed based on the ecological suitability evaluation result, the aggregation area of the indication species can be rapidly displayed in a map, the natural protected area can be determined, and the functional deficiency of the conventional cellular automaton modeling in the process of simulating and optimizing the complex geographical space-time process can be overcome.
An electronic device according to an embodiment of the third aspect of the invention includes: at least one processor, and a memory communicatively coupled to the at least one processor; wherein the memory stores instructions for execution by the at least one processor to cause the at least one processor, when executing the instructions, to implement the natural protected area determination method of the first aspect.
According to the electronic device of the embodiment of the third aspect of the invention, at least the following beneficial effects are achieved: by executing the natural protected area determining method provided by the embodiment of the first aspect of the invention, the aggregation process of the indicated species in the space can be simulated, predicted, optimized and displayed based on the ecological suitability evaluation result, the aggregation area of the indicated species can be rapidly displayed in a map, the natural protected area determination is realized, and the functional deficiency of the conventional cellular automata modeling in the simulation and optimization of the complex geographic space-time process is made up.
A computer-readable storage medium according to an embodiment of a fourth aspect of the present invention stores computer-executable instructions for causing a computer to perform the natural protected area determination method of the first aspect.
The computer-readable storage medium according to the fourth aspect of the present invention has at least the following advantageous effects: by executing the natural protected area determining method provided by the embodiment of the first aspect of the invention, the aggregation process of the indication species in the space can be simulated, predicted, optimized and displayed based on the ecological suitability evaluation result, the aggregation area of the indication species can be rapidly displayed in a map, the natural protected area can be determined, and the functional deficiency of the conventional cellular automaton modeling in the process of simulating and optimizing the complex geographical space-time process can be overcome.
Additional aspects and advantages of the invention will be set forth in part in the description which follows and, in part, will be obvious from the description, or may be learned by practice of the invention.
Drawings
The above and/or additional aspects and advantages of the present invention will become apparent and readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings of which:
fig. 1 is a flowchart illustrating a natural protected area determining method according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of algorithm optimization of a pre-defined agent evolution model according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of a spatial data set according to an embodiment of the present invention;
FIG. 4 is a schematic structural diagram of a natural reserve delimiting apparatus according to an embodiment of the present invention;
fig. 5 is a functional block diagram of an electronic device according to an embodiment of the invention.
Reference numerals:
the system comprises an acquisition module 400, an extraction module 410, a habitat quality evaluation module 420, an agent evolution module 430, a demarcation module 440, a processor 500, a memory 510, a data transmission module 520, a camera 530 and a display screen 540.
Detailed Description
Reference will now be made in detail to embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein like reference numerals refer to the same or similar elements or elements having the same or similar function throughout. The embodiments described below with reference to the accompanying drawings are illustrative only for the purpose of explaining the present invention, and are not to be construed as limiting the present invention.
In the description of the present invention, unless otherwise explicitly defined, terms such as arrangement, installation, connection and the like should be broadly construed, and those skilled in the art can reasonably determine the specific meanings of the above terms in the present invention in combination with the detailed contents of the technical solutions.
First, several terms referred to in the present application are resolved:
1. DEM elevation data: DEM elevation data includes two components: ASTER GDEM30 meter resolution elevation data and SRTM90 meter resolution elevation data. ASTER GDEM data is from NASA, and the data coverage is all land areas between 83 degrees north latitude and 83 degrees south latitude, and the time range is around 2000 years; SRTM data is derived from CIAT, covers all land areas between 60 degrees north latitude and 60 degrees south latitude, and is in the time range of about 2000 years.
2. Habitat Quality (InVEST) Habitat Quality model: the Habitat Quality model (Habitut Quality) is one of the models of the INVEST, and the INVEST model is a comprehensive evaluation model for performing ecosystem service and balance. The evaluation of the Habitat Quality is mainly to use a Habitat Quality model (Habitut Quality) in an INVEST model, wherein the model is mainly to determine a series of data such as threat factors and sensitivity factors and then to process the data in a GIS.
3. And GIS: the Geographic Information System or Geo-Information System, GIS, Geographic Information System, sometimes also referred to as a "Geographic Information System". It is a specific and very important spatial information system. The system is a technical system for collecting, storing, managing, operating, analyzing, displaying and describing relevant geographic distribution data in the whole or partial space of the earth surface layer (including the atmospheric layer) under the support of a computer hardware and software system.
4. CSV: Comma-Separated Values, CSV, Comma Separated value file format. Sometimes also referred to as character separation values, because the separation character may not be a comma. The file stores table data (numbers and text) in plain text form
Referring to fig. 1, a natural reserve determination method according to an embodiment of a first aspect of the present invention includes:
and S100, acquiring land data to be analyzed.
The land data to be analyzed can comprise land utilization data and DEM elevation data of a research area, and the data type is raster data. The land use data may include data such as a type of land cover or a type of land use for the area of study. Optionally, the land utilization data and the DEM elevation data may be acquired, and data association and preprocessing are performed on the two data as analysis basic data, so that land data to be analyzed is obtained.
And S110, extracting land species distribution data from the land data to be analyzed, and extracting threat source data from the land data to be analyzed.
The land species distribution data may be, among other things, specific information about the distribution of organisms (including animals and plants) on the land over the area of study. Species distribution is an important component of biodiversity data, has important significance on biodiversity situation research and biophysical research, is also important data for discovering and evaluating biodiversity hot spot regions, and plays a decision support role in hot spot region protection, natural protection region planning and the like, so that species distribution extraction needs to be carried out on land data to be analyzed to obtain land species distribution data which is used as analysis basic data. The threat source data can be determined according to the drawn threat factor type, for example, the threat factor data is drawn according to the road network, population, GDP and other data and the land data to be analyzed, the drawn threat factor type is obtained, and the type data with the threat is processed to obtain the threat source data.
And step S120, calculating the birth environment quality evaluation data by using a preset habitat quality model, soil species distribution data and threat source data.
The preset habitat quality model can be a preset habitat quality evaluation mathematical model; the habitat quality evaluation data may be data of the quality of the biological habitat in the evaluation area, including a habitat degradation index, a habitat quality evaluation index, and the like. Optionally, the preset Habitat Quality model is set as a Habitut Quality (InVEST) Habitat Quality model, and the Habitat Quality of the area can be quantitatively evaluated according to the land coverage/utilization type of the research area. The Habitat degradation index and the Habitat Quality evaluation index can be calculated and obtained based on a Habitut Quality (InVEST) Habitat Quality model by using the land species distribution data and the threat source data extracted in the step S110, and the data represent the Habitat Quality change trend of land use change, so that the Habitat Quality evaluation data are obtained.
And step S130, calculating natural preservation area data by using a preset intelligent agent evolution model and the habitat quality evaluation data.
The preset intelligent agent evolution model can be a pre-compiled improved group intelligent agent evolution mathematical model, and the basic idea of the group intelligent agent is to perform optimization solution by simulating the predation and migration process of the biological group, namely searching for the shortest path for acquiring food and migration; the natural reserve data may include: total area of the protected area, habitat quality, compactness index of space, ecological suitability and the like. Optionally, in the preset agent evolution model provided in the embodiment, the biological habit of the regional species may be simulated by constructing the swarm agents, and the model is used for finding an optimized path in the grid image of the research area, where the key for constructing the swarm agent model is to establish an optimized pheromone mechanism for spatial layout optimization. Through the compiled preset intelligent evolution model, the habitat quality evaluation data in the step S120 are used, factors such as ecological suitability and compactness are considered, ecological protection target parameters are determined, a utility function is introduced, regional optimization calculation is completed, and final natural protection area data are obtained.
And step S140, generating a natural preservation area according to the natural preservation area data and the habitat quality evaluation data.
Optionally, natural protection area delimiting data may be generated according to the natural protection area data and the habitat quality evaluation data to obtain a final analysis result, and the final analysis result is loaded into a map for browsing to form a final natural protection area.
The natural protection area dividing method comprises the steps of firstly, obtaining land data to be analyzed, then extracting land species distribution data from the land data to be analyzed, extracting threat source data from the land data to be analyzed, then utilizing a preset habitat quality model to calculate habitat quality evaluation data for the soil species distribution data and the threat source data, utilizing a preset intelligent body evolution model and the habitat quality evaluation data to calculate natural protection area data, finally generating a natural protection area according to the natural protection area data and the habitat quality evaluation data, simulating, predicting, optimizing and displaying the indication species in a space gathering process based on an ecological suitability evaluation result, the gathering area of the indicated species is rapidly displayed in the map, natural protection area division is realized, and the functional deficiency of the conventional cellular automata modeling in the simulation and optimization of the complex geographic space-time process is overcome.
In some embodiments of the invention, obtaining land data to be analyzed comprises:
and acquiring land utilization data and DEM elevation data.
And performing correlation preprocessing on the land utilization data and the DEM elevation data to obtain land data to be analyzed.
Optionally, the land cover/utilization type of the research area may be obtained to obtain land utilization data; the DEM elevation data is used as common data in the GIS, the resolution ratio and the source of the data are various, and the DEM elevation data can be acquired from various platforms at home and abroad, for example, the DEM elevation data can be acquired from platforms such as a national earth system scientific data sharing platform and a geographic space data cloud. The data type of the land use data and the DEM elevation data is raster data.
Optionally, the land utilization data and the DEM elevation data may be subjected to correlation preprocessing, for example, GIS analysis including buffer area analysis, intersection analysis, terrain analysis, distance analysis, geostatistical analysis, etc., on the land utilization data and the DEM elevation data, so as to obtain land data to be analyzed. By acquiring the land utilization data and the DEM elevation data and performing correlation preprocessing on the land utilization data and the DEM elevation data, the land data to be analyzed serving as an analysis basis can be obtained, and the dividing precision of the natural protected area is improved.
In some embodiments of the invention, extracting land species distribution data from land data to be analyzed comprises:
and carrying out species distribution extraction on the land data to be analyzed to obtain land species distribution data.
Optionally, the land data to be analyzed may be analyzed, for example: land species distribution data can be quickly extracted according to GPS points of species distribution, or the land species distribution data can be extracted by analyzing data to be analyzed according to an algorithm model of species distribution, so that the land species distribution data can be obtained and used as analysis basic data.
In some embodiments of the invention, extracting threat source data from the analyzed land data comprises:
drawing up threat factor data according to the land data to be analyzed;
acquiring a threat factor type corresponding to the threat factor data;
and assigning the threat factor data according to the type of the threat factor to obtain threat source data.
Wherein the threat factor data may include: regional human activity data (including population distribution, GDP industry distribution and the like), traffic network vector data (including railways, expressways and urban traffic roads) and the like. Optionally, the threat factor data may be set according to requirements, for example: and (3) drawing up threat factor data in a CSV (China general analysis) format by combining the data of the land to be analyzed, the road network, the population, the GDP (graphics data processing) and the like, sorting each item of data through literature or scoring by an obtained expert, and determining parameters such as the influence distance, the weight, the attenuation linearity and the like of the threat factor.
Wherein the threat factor type may be a category for threat factor data. Optionally, the threat factor type may be drawn up according to the threat factor data, for example, the population distribution data is drawn up as regional human activity data type, the railway traffic road data is drawn up as traffic network vector data type, and the like.
Optionally, data of each threat source may be acquired according to the proposed threat factor type, specifically: and assigning the type data with the threat to be 1, and assigning other types of data to be 0, so that the threat source data can be processed. Threat source data are obtained through threat factor data and threat factor type processing, authenticity and effectiveness of the statistic threat source data can be improved, and scientificity of natural protection area division is enhanced.
In some embodiments of the present invention, the calculating of the habitat quality evaluation data using the preset habitat quality model, the soil species distribution data, and the threat source data includes:
and substituting the soil species distribution data and the threat source data into a preset habitat quality model, and calculating a habitat degradation index and a habitat quality evaluation index.
And obtaining the habitat quality evaluation data according to the habitat degradation index and the habitat quality evaluation index.
Optionally, the preset Habitat Quality model is set as a Habitut Quality (InVEST) Habitat Quality model, and the Habitut Quality (InVEST) Habitat Quality model can quantitatively evaluate the Habitat Quality of the area according to the land covering/utilization type of the research area. Setting land species distribution data to comprise land covering/using habitat types, threat source data to comprise parameters such as ecological stress factors, sensitivity of the habitat types to the stress factors, maximum influence distance of the stress factors and the like, and inputting parameters such as land covering/using habitat types, sensitivity of the ecological stress factors, the sensitivity of the habitat types to the stress factors, the maximum influence distance of the stress factors and the like of a research area into a preset habitat quality model to obtain normalized (with the value range of 0-1) habitat degradation degrees and habitat quality raster images, namely habitat degradation indexes and habitat quality evaluation indexes. Wherein, the habitat degradation degree is the interference intensity of the stress factor to the habitat, the habitat quality is the suitability degree of the habitat, and the calculation formula is shown as the following formula:
Figure BDA0003126898490000091
wherein S isxjThe habitat quality score of the grid unit x in the land covering/utilizing type j is represented by a value range of [0, 1%]The habitat quality score is 1 when the habitat suitability is the highest, and the score is 0 when the habitat suitability is the lowest; hjAssigning an initial habitat suitability value when the land covering/utilizing type j is not disturbed by the stress factor, representing the supporting capability of the land covering type as the habitat, and taking the value of [0, 1]；DxjThe habitat degradation degree of the grid cells x in the land covering/utilizing type j after being disturbed by the stress factors is obtained; k is a half-saturation constant, typically set to 1/2 which is the maximum value of the habitat degradation degree; z is a normalized model default parameter with a constant of 2.5. In the above formula, the calculation formula of the habitat degradation degree is shown as the following formula (ii):
Figure BDA0003126898490000092
in the formula: r is a stress factor; r is the number of the stress factors; y is a single grid unit in the grid layer of the threat factor r; y isrThe total number of grids of the threat factor r image layer; w is arRepresenting the weight of the stress factor r; r isyThe grid unit y is used for judging whether the grid unit y is the source of the stress factor r; i.e. irxyA distance function between the habitat and the threat factor; beta is axIs the level of accessibility of the coercion factor to grid cell x; sjrThe relative sensitivity of type j to the stress factor r for land cover/use. r isy、irxy、βx、SjrAll values of (1) are [0, 1 ]]. Further, irxyThe distance function reflects the variation of the stress factor with distance when the stress factor disturbs the habitat grid cell x, and generally, the farther the habitat grid cell x is away from the stress factor, the less the habitat grid cell x is affected by the disturbance. Aiming at different stress factors, the following two distance attenuation formulas (c) and (d) can be respectively used for calculation according to actual conditions:
Figure BDA0003126898490000101
Figure BDA0003126898490000102
wherein, formula (c) is linear attenuation, formula (iv) is exponential attenuation, and in the formula: dxyIs the linear distance, d, between the environmental grid cell x and the grid cell yrmaxThe maximum stress distance of the stress factor r. And (4) calculating the habitat degradation index and the habitat quality evaluation index according to the formulas I to IV so as to obtain habitat quality evaluation data, wherein the data represent the habitat quality change trend of land use change. Substituting data such as land utilization data, threat factors and threat source data into a preset habitat quality model for grid calculation, and calculating the quality index of each pixel of the land utilization grid data to obtainHabitat quality evaluation data, which can be evaluated based on the habitat suitability of the regional indicator species.
In some embodiments of the present invention, the calculating using the preset agent evolution model and the habitat quality evaluation data to obtain the natural preservation area data includes:
acquiring ecological protection target parameters and a utility function;
and performing regional optimization calculation according to a preset intelligent agent evolution model, ecological protection target parameters and a utility function to obtain natural protection area data.
Optionally, as shown in fig. 2, fig. 2 is a schematic diagram of algorithm optimization of a preset agent evolution model according to an embodiment of the present invention. In the swarm agent algorithm of the embodiment, the swarm agents are constructed to simulate the biological habits of the regional species for finding the optimal path in the grid image of the research area. The key for constructing the swarm intelligence model is to establish an optimized pheromone mechanism for space layout optimization. The detailed steps are as follows:
first, the probability that the kth agent selects a grid cell to form a natural-protected area at time t may be defined as follows:
Figure BDA0003126898490000111
in the formula:
Figure BDA0003126898490000115
representing the selection probability of the agent k to the grid cell i at the time t; allowedkRepresenting the set of grid cells that agent k can visit next tour; alpha is pheromone enlightenment factor and reflects the concentration tau of the pheromone of the intelligent body in the process of movingi(t) the effect of; beta is an expected information heuristic factor used for representing the heuristic function eta of the intelligent body in the process of movementi(t) the function of the compound (a). Heuristic function etai(t) directing agent movement by including habitat quality on grid cell i, agent tendency to select movement clusters to grid cells with high habitat quality scores, this characteristic substantially corresponding toBiological habits of general species, thereby forming a relatively reasonable biological protection pattern. Heuristic function etai(t) may be represented by the following formula (i):
Figure BDA0003126898490000112
the formula is as follows: seiThe habitat quality score is the grid unit i;
Figure BDA0003126898490000113
which is the sum of the habitat quality scores for all grid cells in the study area, is calculated in step 120.
Secondly, the compactness index is introduced, so that the layout of the protection area is kept reasonable in form, management and control are facilitated, and the fragmentation of the land utilization pattern is avoided. And integrating the habitat quality and compactness target through the weight, constructing a benefit target function, and determining the benefit target function as a total benefit value of the spatial layout scheme. The compactness index is calculated from the total area of the protected scene composed of all selected elements and its perimeter. The indicator is defined according to the following ratio function (c), (b):
Figure BDA0003126898490000114
Uutility＝weS+wcP ⑧
in the formulas (III), (IV), (V), (IV), (V), (IV), (V), (IV), (V; u shapeutilityExpressing the benefit objective function, S and P being compactness index values of the environmental quality and space, weAnd wcRespectively, corresponding weights. Predicting an optimal protection mode UutilityThe highest values should be generated at the average overall ecological suitability and compactness index.
Finally, a benefit objective function is introduced into the pheromone concentration tauiIn the calculation of (t), the formula is as follows:
Figure BDA0003126898490000121
wherein d (x) is the distance from the protective unit x to the central cell i. The distance variable d (x) reflects the neighborhood impact of the grid cells in the protection area, and if the neighborhood cells of a grid cell are already listed in the protection area, then the grid cell is more likely to be included in the protection area. The size of the neighborhood window is typically chosen to be 3 x 3, 5 x 5 or 7 x 7 window depending on the application. By means of a preset intelligent agent evolution model shown in fig. 2, ecological suitability, compactness and other factors are considered, ecological protection target parameters are established, a utility function is introduced, regional optimization calculation is completed, and final natural protection area data are obtained, interaction between a main body and a regional ecological environment can be effectively simulated, interpretation modes of spatial simulation on geographic process evolution are further enriched, and the process of indicating species in spatial aggregation is simulated, predicted, optimized and displayed based on ecological suitability evaluation results, so that data are automatically subjected to iterative aggregation.
In some embodiments of the invention, generating the natural protected area from the natural protected area data and the habitat quality evaluation data comprises:
generating natural protection area delimiting data according to the natural protection area data and the habitat quality evaluation data;
and dividing data loading according to the natural protection area to obtain the natural protection area.
Optionally, the natural preservation area defining data includes habitat quality evaluation data (including a habitat degradation index and a habitat quality evaluation index) and final natural preservation area data. The grid images of the habitat quality evaluation and the habitat degradation degree evaluation of the research area are intermediate generation data, and the natural reserve area is determined as a final analysis result. As shown in fig. 3, the natural reserve area data may include, but is not limited to, the spatial data set shown in fig. 3, the data of land cover/utilization of the research area as shown in fig. 3, regional human activity data (including population distribution, GDP industry distribution), DEM digital elevation data, traffic network vector data (including railway, expressway and urban traffic road), and the like. After the nature protected area demarcating data is generated, it can be loaded on a map for browsing. The method is characterized in that the space gathering process of the indicating species is simulated, predicted, optimized and displayed based on the ecological suitability evaluation result, the gathering area of the indicating species is rapidly displayed in a map, the species protection area is defined, and the method can be used as a land space simulation and optimization platform to make up for the functional deficiency of the conventional cellular automata modeling in the simulation and optimization of the complex geographic space-time process.
Referring to fig. 4, a natural preservation area delineation apparatus according to an embodiment of the second aspect of the present invention includes:
an acquisition module 400, configured to acquire land data to be analyzed;
the extraction module 410 is used for extracting land species distribution data from land data to be analyzed and extracting threat source data from the land data to be analyzed;
the habitat quality evaluation module 420 is configured to calculate habitat quality evaluation data by using a preset habitat quality model, soil species distribution data and threat source data;
the agent evolution module 430 is used for calculating natural protection area data by using a preset agent evolution model and the habitat quality evaluation data;
and the delimiting module 440 is used for generating the natural protection area according to the natural protection area data and the habitat quality evaluation data.
By executing the natural protected area defining method of the embodiment of the first aspect of the invention, the natural protected area defining device can simulate, predict, optimize and display the aggregation process of the indication species in the space based on the ecological suitability evaluation result, quickly display the aggregation area of the indication species in the map, realize the definition of the natural protected area and make up for the functional deficiency of the conventional cellular automata modeling in the simulation and optimization of the complex geographic space-time process.
Referring to fig. 5, an embodiment of the third aspect of the present invention further provides a functional module diagram of an electronic device, including: at least one processor 500, and a memory 510 communicatively coupled to the at least one processor 500; and the system also comprises a data transmission module 520, a camera 530 and a display screen 540.
Wherein the processor 500 is adapted to perform the natural protected area determination method in the first aspect embodiment by calling a computer program stored in the memory 510.
The data transmission module 520 is connected to the processor 500, and is configured to implement data interaction between the data transmission module 520 and the processor 500.
The cameras 530 may include front cameras and rear cameras. Generally, a front camera is disposed at a front panel of a terminal, and a rear camera is disposed at a rear surface of the terminal. In some embodiments, the number of the rear cameras is at least two, and each of the rear cameras is any one of a main camera, a depth-of-field camera, a wide-angle camera and a telephoto camera, so that the main camera and the depth-of-field camera are fused to realize a background blurring function, the main camera and the wide-angle camera are fused to realize panoramic shooting and a VR (Virtual Reality) shooting function or other fusion shooting functions. In some embodiments, camera 530 may also include a flash. The flash lamp can be a single-color temperature flash lamp or a double-color temperature flash lamp. The double-color-temperature flash lamp is a combination of a warm-light flash lamp and a cold-light flash lamp, and can be used for light compensation at different color temperatures.
The display 540 may be used to display information entered by the user or provided to the user. The Display screen 540 may include a Display panel, and optionally, the Display panel may be configured in the form of a Liquid Crystal Display (LCD), an Organic Light-Emitting Diode (OLED), or the like. Further, the touch panel may cover the display panel, and when the touch panel detects a touch operation thereon or nearby, the touch panel transmits the touch operation to the processor 500 to determine the type of the touch event, and then the processor 500 provides a corresponding visual output on the display panel according to the type of the touch event. In some embodiments, the touch panel may be integrated with the display panel to implement input and output functions.
The memory, as a non-transitory storage medium, may be used to store a non-transitory software program and a non-transitory computer-executable program, such as the natural protected area determination method in the embodiment of the first aspect of the present invention. The processor implements the natural protected area determination method in the first embodiment by executing the non-transitory software program and the instructions stored in the memory.
The memory may include a program storage area and a data storage area, wherein the program storage area may store an operating system, an application program required for at least one function; the storage data area may store a natural protection area determination method executed in the embodiment of the first aspect. Further, the memory may include high speed random access memory, and may also include non-transitory memory, such as at least one magnetic disk storage device, flash memory device, or other non-transitory solid state storage device. In some embodiments, the memory optionally includes a memory located remotely from the processor, and these remote memories may be connected to the terminal over a network. Examples of such networks include, but are not limited to, the internet, intranets, local area networks, mobile communication networks, and combinations thereof.
Non-transitory software programs and instructions required to implement the natural protection zone determination method in the first aspect embodiment described above are stored in a memory and, when executed by one or more processors, perform the natural protection zone determination method in the first aspect embodiment described above.
Embodiments of the fourth aspect of the present invention also provide a computer-readable storage medium storing computer-executable instructions for: the natural guard area determination method in the first aspect embodiment is performed.
In some embodiments, the storage medium stores computer-executable instructions, which are executed by one or more control processors, for example, by one of the processors in the electronic device of the third aspect, and may cause the one or more processors to perform the natural protection zoning method of the first aspect.
The embodiments of the present invention have been described in detail with reference to the accompanying drawings, but the present invention is not limited to the above embodiments, and various changes can be made within the knowledge of those skilled in the art without departing from the gist of the present invention.
The above described embodiments of the apparatus are merely illustrative, wherein the units illustrated as separate components may or may not be physically separate, may be located in one place, or may be distributed over a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the embodiment.
It will be understood by those of ordinary skill in the art that all or some of the steps, systems, and methods disclosed above may be implemented as software, firmware, hardware, and suitable combinations thereof. Some or all of the physical components may be implemented as software executed by a processor, such as a central processing unit, digital signal processor, or microprocessor, or as hardware, or as an integrated circuit, such as an application specific integrated circuit. Such software may be distributed on computer readable media, which may include computer storage media (or non-transitory media) and communication media (or transitory media). As is well known to those of ordinary skill in the art, the term computer storage media includes volatile and nonvolatile, removable and non-removable media implemented in any method or technology for storage of information such as computer readable instructions, data structures, program modules or other data. Computer storage media includes, but is not limited to, RAM, ROM, EEPROM, flash memory or other memory technology, CD-ROM, Digital Versatile Disks (DVD) or other optical disk storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to store the desired information and which can accessed by a computer. In addition, communication media typically embodies computer readable instructions, data structures, program modules or other data in a modulated data signal such as a carrier wave or other transport mechanism and includes any information delivery media as known to those skilled in the art.
In the description herein, references to the description of the terms "one embodiment," "some embodiments," "an illustrative embodiment," "an example," "a specific example," or "some examples" or the like mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, the schematic representations of the terms used above do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
While embodiments of the invention have been shown and described, it will be understood by those of ordinary skill in the art that: various changes, modifications, substitutions and alterations can be made to the embodiments without departing from the principles and spirit of the invention, the scope of which is defined by the claims and their equivalents.

Claims (10)

1. A natural reserve determination method is characterized by comprising:
acquiring land data to be analyzed;
extracting land species distribution data from the land data to be analyzed, and extracting threat source data from the land data to be analyzed;
calculating habitat quality evaluation data by using a preset habitat quality model, the soil species distribution data and the threat source data;
calculating natural protection area data by using a preset intelligent agent evolution model and the habitat quality evaluation data;
and generating a natural protection area according to the natural protection area data and the habitat quality evaluation data.
2. The method of claim 1, wherein the obtaining land data to be analyzed comprises:
acquiring land utilization data and DEM elevation data;
and performing correlation preprocessing on the land utilization data and the DEM elevation data to obtain the land data to be analyzed.
3. The method of claim 2, wherein said extracting land species distribution data from said land data to be analyzed comprises:
and carrying out species distribution extraction on the land data to be analyzed to obtain the land species distribution data.
4. The method of claim 3, wherein said extracting threat source data from said land data to be analyzed comprises:
drawing up threat factor data according to the land data to be analyzed;
acquiring a threat factor type corresponding to the threat factor data;
and assigning the threat factor data according to the threat factor type to obtain the threat source data.
5. The method of claim 4, wherein the computing of habitat quality evaluation data using a preset habitat quality model, the soil species distribution data and the threat source data comprises:
substituting the soil species distribution data and the threat source data into the preset habitat quality model to calculate a habitat degradation index and a habitat quality evaluation index;
and obtaining the habitat quality evaluation data according to the habitat degradation index and the habitat quality evaluation index.
6. The method of claim 1, wherein the calculating using the pre-set agent evolution model and the habitat quality evaluation data to obtain natural reserve data comprises:
acquiring ecological protection target parameters and a utility function;
and performing regional optimization calculation according to the preset intelligent agent evolution model, the ecological protection target parameters and the utility function to obtain the natural protection area data.
7. The method of claim 1, wherein generating a physical protective area from the physical protective area data and the habitat quality evaluation data comprises:
generating natural protection area delimiting data according to the natural protection area data and the habitat quality evaluation data;
and loading data according to the natural protection area to obtain the natural protection area.
8. Nature protects district demarcation device, its characterized in that includes:
the acquisition module is used for acquiring land data to be analyzed;
the extraction module is used for extracting land species distribution data from the land data to be analyzed and extracting threat source data from the land data to be analyzed;
the habitat quality evaluation module is used for calculating habitat quality evaluation data by utilizing a preset habitat quality model, the soil species distribution data and the threat source data;
the intelligent agent evolution module is used for calculating natural protection area data by utilizing a preset intelligent agent evolution model and the habitat quality evaluation data;
and the demarcation module is used for generating a natural protection area according to the natural protection area data and the habitat quality evaluation data.
9. An electronic device, comprising:
at least one processor, and,
a memory communicatively coupled to the at least one processor; wherein the content of the first and second substances,
the memory stores instructions for execution by the at least one processor to cause the at least one processor, when executing the instructions, to implement the natural protected area determination method of any of claims 1 to 7.
10. A computer-readable storage medium having stored thereon computer-executable instructions for causing a computer to perform the natural protected area determination method as defined in any one of claims 1 to 7.
CN202110693013.1A 2021-06-22 2021-06-22 Natural protected area determining method and device, electronic equipment and storage medium Active CN113469512B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110693013.1A CN113469512B (en) 2021-06-22 2021-06-22 Natural protected area determining method and device, electronic equipment and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110693013.1A CN113469512B (en) 2021-06-22 2021-06-22 Natural protected area determining method and device, electronic equipment and storage medium

Publications (2)

Publication Number Publication Date
CN113469512A true CN113469512A (en) 2021-10-01
CN113469512B CN113469512B (en) 2022-08-30

Family

ID=77869145

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110693013.1A Active CN113469512B (en) 2021-06-22 2021-06-22 Natural protected area determining method and device, electronic equipment and storage medium

Country Status (1)

Country Link
CN (1) CN113469512B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116049342A (en) * 2023-03-24 2023-05-02 清华大学 Habitat quality monitoring method and device

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106294289A (en) * 2016-08-25 2017-01-04 环境保护部南京环境科学研究所 A kind of ecological red line demarcation method protecting animal
CN107977711A (en) * 2017-12-21 2018-05-01 武汉大学 A kind of multi-Agent Genetic Algorithm towards " three lines " collaboration optimization
CN110334864A (en) * 2019-07-03 2019-10-15 洛阳市规划建筑设计研究院有限公司 A method of more rule unification Urban-rural space subregion based on GIS delimited
CN110852532A (en) * 2019-11-26 2020-02-28 浙江大学 Natural protected area system integration and boundary delineation method
CN112528359A (en) * 2020-08-06 2021-03-19 湖北民族大学 City expansion simulation method based on bargaining model and ant colony optimization algorithm

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106294289A (en) * 2016-08-25 2017-01-04 环境保护部南京环境科学研究所 A kind of ecological red line demarcation method protecting animal
CN107977711A (en) * 2017-12-21 2018-05-01 武汉大学 A kind of multi-Agent Genetic Algorithm towards " three lines " collaboration optimization
CN110334864A (en) * 2019-07-03 2019-10-15 洛阳市规划建筑设计研究院有限公司 A method of more rule unification Urban-rural space subregion based on GIS delimited
CN110852532A (en) * 2019-11-26 2020-02-28 浙江大学 Natural protected area system integration and boundary delineation method
CN112528359A (en) * 2020-08-06 2021-03-19 湖北民族大学 City expansion simulation method based on bargaining model and ant colony optimization algorithm

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
张学儒 等: ""基于土地利用格局重建的区域生境质量时空变化分析"", 《地理学报》 *
王 尧: ""基于景观格局和蚁群算法的横山区农村居民点布局优化研究"", 《中国优秀博硕士学位论文全文数据库（硕士） 经济与管理科学辑》 *
黄木易 等: ""基于InVEST 模型的皖西大别山区生境质量时空演化及景观格局分析"", 《生态学报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116049342A (en) * 2023-03-24 2023-05-02 清华大学 Habitat quality monitoring method and device

Also Published As

Publication number Publication date
CN113469512B (en) 2022-08-30

Similar Documents

Publication Publication Date Title
Mustafa et al. Coupling agent-based, cellular automata and logistic regression into a hybrid urban expansion model (HUEM)
Terribile et al. A Web-based spatial decision supporting system for land management and soil conservation
CN109001091B (en) Satellite remote-sensing monitoring method, device and the computer-readable medium of atmosphere pollution
Macdonald et al. Modelling space use and dispersal of mammals in real landscapes: a tool for conservation
CN108351959A (en) Image Acquisition task system based on variation detection
JP2019533180A (en) Change detection based map interface update system
Thies et al. Projecting land-use and land-cover changes in a tropical mountain forest of Southern Ecuador
Rhodes et al. The relative value of field survey and remote sensing for biodiversity assessment
Adams et al. Distribution, demography and dispersal model of spatial spread of invasive plant populations with limited data
Guo et al. Integrating airborne lidar and satellite imagery to model habitat connectivity dynamics for spatial conservation prioritization
Gómez et al. Modelling desert locust presences using 32-year soil moisture data on a large-scale
Hand et al. Assessing multi-taxa sensitivity to the human footprint, habitat fragmentation and loss by exploring alternative scenarios of dispersal ability and population size: a simulation approach
Schulz et al. A spatial application of a marten habitat model
WO2022009739A1 (en) Vegetation management system and vegetation management method
Monsimet et al. Explicit integration of dispersal-related metrics improves predictions of SDM in predatory arthropods
CN113469512B (en) Natural protected area determining method and device, electronic equipment and storage medium
Borana et al. Prediction of land cover changes of Jodhpur city using cellular automata Markov modelling techniques
Bourdouxhe et al. Comparison of methods to model species habitat networks for decision-making in nature conservation: The case of the wildcat in southern Belgium
Lamichhane et al. Land use land cover (LULC) change projection in Kathmandu valley using the clue-s model
Ma et al. Mapping potential wilderness in China with location-based services data
CN104318319A (en) Multi-agent urban expansion simulation method based on land acquisition process competition result
Suárez-Seoane et al. Habitat-partitioning improves regional distribution models in multi-habitat species: A case study with the European bilberry
Fogarty et al. Serious game modeling of caribou behavior across Lake Huron using cultural algorithms and influence maps
Greenspan et al. Large felid habitat connectivity in the transboundary Dawna-Tanintharyi landscape of Myanmar and Thailand
Sullivan et al. Changes in habitat associations during range expansion: disentangling the effects of climate and residence time

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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Dan Liang

Inventor after: Li Qijun

Inventor after: Yue Juan

Inventor after: Tang Wenkui

Inventor after: Chen Junli

Inventor after: Yuan Fulin

Inventor before: Dan Liang

Inventor before: Tang Wenkui

Inventor before: Yue Juan

Inventor before: Chen Junli

Inventor before: Yuan Fulin

GR01 Patent grant
GR01 Patent grant
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 10 / F, Jianyi building, 3 Zhenxing Road, Futian District, Shenzhen, Guangdong 518028

Patentee after: Shenzhen Urban Planning and Design Institute Co.,Ltd.

Address before: 10 / F, Jianyi building, 3 Zhenxing Road, Futian District, Shenzhen, Guangdong 518028

Patentee before: URBAN PLANNING & DESIGN INSTITUTE OF SHENZHEN (UPDIS)