CN114611866A

CN114611866A - Visualization presentation method and device for total production value of plot ecosystem

Info

Publication number: CN114611866A
Application number: CN202210076147.3A
Authority: CN
Inventors: 周祖煜; 陈煜人; 李兆伟; 张澎彬; 白博文; 林波; 莫志敏; 张�浩; 李天齐; 刘俊
Original assignee: Hangzhou Lingjian Digital Agricultural Technology Co ltd
Current assignee: Hangzhou Lingjian Digital Agricultural Technology Co ltd
Priority date: 2022-01-23
Filing date: 2022-01-23
Publication date: 2022-06-10

Abstract

The application provides a visual presentation method and device for a total production value of a plot ecosystem, which relate to the technical field of remote sensing and comprise the following steps: selecting a fixed land, making an ecological product list of the fixed land and acquiring related data to obtain first data; calculating the value quantity of the first data according to the ecological product list according to a visualization method to obtain second data; and according to the second data, carrying out superposition calculation on the value quantity of the material product, the value quantity of the adjusting service and the value quantity of the cultural service through a grid calculator to obtain a total value image of the ecological system production of the fixed plot. The easy acquirement nature of this scheme utilization remote sensing data, visual and minimum statistical range depend on the characteristics of image in space, error, the data threshold that reduces GEP account that bring when can avoid discerning ecosystem improve the real-time, and the formation is visual image, and the condition that can more audio-visual analysis GEP through this image carries out ecological civilization construction and guides.

Description

Visualization presentation method and device for total production value of plot ecosystem

Technical Field

The invention belongs to the technical field of remote sensing, and particularly relates to a visual presentation method and device for a total production value of a plot ecosystem.

Background

The natural ecosystem not only provides abundant ecological products and services for human beings, but also has huge ecological benefits, and simultaneously, the ecological value can be converted into economic benefits to benefit the human beings. Natural ecosystem such as forest, lake, meadow, marsh, river and ocean has not only provided food, drinking water, timber, medicine etc. more provides the necessary ecosystem service function of human existence and development, including adjust the weather, release oxygen, fixed carbon dioxide, keep in good health the water source, the regulation is flooded, prevent wind and fix sand etc. has contributed huge economic benefits to human welfare. The ecosystem production total value (GEP) means that in a certain area, the ecosystem provides the sum of final products and services for human beings and the economic value of the final products and services, and is the total currency value contributed by the ecosystem of the certain area to human welfare.

In the existing ecosystem production total value (GEP) calculation method, the following problems exist: (1) the type of the ecological system needs to be identified firstly, the identification precision can directly influence the final calculation result of the total production value of the ecological system, and interference factors influencing the precision of the calculation result exist; (2) data for performing accounting relates to a plurality of government mechanisms, including but not limited to departments such as water conservancy, natural resources, ecological environment, forestry and the like, and authorities such as homeland survey, forestry survey results and the like, and has the problem that the accounting data is difficult to obtain; (3) the form data is finally obtained by calculation according to the standard in the technical specification file of the accounting, but the problems of complex calculation mode, easy error and high updating difficulty exist; (4) the calculated total production value of the ecological system is updated with frequency generally taking one year as a time unit, so that certain time lag exists, and data have the problem of lag; (5) in addition, the prior art cannot realize the visualization of the total production value of the ecosystem, and cannot be used for analyzing the GEP condition in a certain area and the factors influencing the GEP change more finely.

Disclosure of Invention

The invention provides a visual presentation method and device for a total production value of a land ecosystem, and aims to solve the problems that interference factors influencing the accuracy of a calculation result exist, the calculation data are difficult to obtain, the calculation mode is complex, errors are easy to make, the updating difficulty is high, the calculation data have hysteresis, the visualization of the total production value of the ecosystem cannot be realized, and the calculation data cannot be used for analyzing the GEP condition in a certain area and influencing the GEP change factors more precisely.

In order to achieve the above object, the present application adopts the following technical solutions, including:

selecting a fixed land, formulating an ecological product list of the fixed land and acquiring related data to obtain first data, wherein the related data comprises remote sensing data, basic geographic data, literature data, statistical yearbook and meteorological data;

calculating the value quantity of the first data according to the ecological product list by a visualization method to obtain second data, wherein the ecological product list comprises material products, adjusting services and cultural services;

and according to the second data, carrying out superposition calculation on the value quantity of the material product, the value quantity of the adjusting service and the value quantity of the cultural service through a grid calculator to obtain a total value image of the ecological system production of the fixed plot.

Preferably, the calculating the value amount of the first data according to the ecological product list according to a visualization method to obtain second data includes:

calculating the value quantity of the material product and/or the cultural service according to a coefficient method and the first data to obtain a first value quantity, wherein the visualization method comprises the coefficient method, an image substitution method and a weighted distribution method;

and calculating the value quantity of the adjusting service according to a visualization method and the first data to obtain a second value quantity, wherein the adjusting service comprises water source conservation, soil conservation, carbon fixation, oxygen supply, atmosphere purification, water quality purification and climate adjustment.

Preferably, the calculating the value amount of the adjustment service according to the visualization method and the first data to obtain a second value amount includes:

replacing the data of water source conservation, soil conservation, carbon fixation, oxygen supply and climate regulation by using a remote sensing data product according to the image substitution method, and processing and calculating the image of the remote sensing data product by using a grid calculator to obtain a first value image;

carrying out numerical value distribution on the atmospheric purification data according to remote sensing image pixel points of the fixed plot according to a weighted distribution method, and processing and calculating the image subjected to numerical value distribution according to a grid calculator to obtain a second value quantity image;

and calculating the value quantity of water purification according to the first data by a coefficient method to obtain a water purification value quantity coefficient, assigning the water purification value quantity coefficient in a vector area of the fixed plot to obtain a first vector area, and performing vector-to-grid processing to obtain a third value quantity image.

Preferably, the obtaining of the total value image of the ecosystem production of the fixed plot by performing superposition calculation on the value quantity of the material product, the value quantity of the adjustment service and the value quantity of the culture service through a grid calculator according to the second data includes:

performing grid calculation on the first value quantity image, the second value quantity image and the third value quantity image, and summing to obtain a value quantity image of the adjusting service;

and superposing the value quantity image of the adjusting service and the first value quantity according to a grid calculator to obtain a total value image of the ecosystem production.

A visual presentation device of a total value of a land ecosystem production, comprising:

a parcel data acquisition module: the system comprises a data acquisition module, a data processing module and a data processing module, wherein the data acquisition module is used for selecting a fixed plot, formulating an ecological product list of the fixed plot and acquiring related data to obtain first data, and the related data comprises remote sensing data, basic geographic data, literature data, statistical yearbook and meteorological data;

ecological product value volume calculation module: the system comprises a data acquisition module, a data processing module and a data processing module, wherein the data acquisition module is used for acquiring first data and second data, and the first data is calculated according to a visual method and value quantity according to an ecological product list to obtain the second data, and the ecological product list comprises material products, adjusting services and cultural services;

the ecosystem production total value image generation module: and the ecological system production total value image of the fixed plot is obtained by performing superposition calculation on the value quantity of the material product, the value quantity of the adjusting service and the value quantity of the cultural service through a grid calculator according to second data.

Preferably, the ecological product value calculation module includes:

a first value amount calculation module: the system comprises a visualization method, a coefficient method, an image substitution method and a weighted distribution method, wherein the visualization method is used for calculating the value quantity of the material product and/or the cultural service according to the coefficient method and the first data to obtain a first value quantity;

a second value amount calculation module: and the system is used for calculating the value quantity of the adjusting service according to a visualization method and the first data to obtain a second value quantity, wherein the adjusting service comprises water source conservation, soil conservation, carbon fixation, oxygen supply, atmosphere purification, water quality purification and climate adjustment.

Preferably, the second value calculation module includes:

the first image calculation module: the system comprises a remote sensing data product, a grid calculator and a data processing and calculating unit, wherein the remote sensing data product is used for replacing data of water source conservation, soil conservation, carbon fixation, oxygen supply and climate regulation according to the image substitution method, and the image of the remote sensing data product is processed and calculated by the grid calculator to obtain a first value image;

the second image calculation module: the system is used for carrying out numerical value distribution on the atmospheric purification data according to remote sensing image pixel points of the fixed plot according to a weighted distribution method, and processing and calculating the image subjected to numerical value distribution according to a grid calculator to obtain a second value quantity image;

an image third calculation module: and the system is used for calculating the value quantity of water purification according to the first data by a coefficient method to obtain a water purification value quantity coefficient, assigning the water purification value quantity coefficient in a vector area of the fixed plot to obtain a first vector area, and performing vector-to-grid processing to obtain a third value quantity image.

Preferably, the ecosystem production total value image generation module includes:

a first value image generation module: the value quantity image used for summing the first value quantity image, the second value quantity image and the third value quantity image after grid calculation to obtain the value quantity image of the adjusting service;

a second value image generation module: and the grid calculator is used for superposing the value quantity image of the adjusting service and the first value quantity to obtain a total value image of the ecosystem production.

A visual presentation apparatus of a total value produced by a land ecosystem, comprising a memory for storing one or more computer instructions and a processor, wherein the one or more computer instructions are executed by the processor to implement a visual presentation method of a total value produced by a land ecosystem as claimed in any one of the above.

A computer-readable storage medium storing a computer program which, when executed by a computer, implements a method for visually presenting a total production value of a land ecosystem as claimed in any one of the above.

The invention has the following technical effects:

(1) before the scheme is formulated, the reason that different ecosystems need to be distinguished in the original standard is considered and analyzed, and the fact that parameters such as evaporation capacity and vegetation coverage factors of the different ecosystems are obviously different is found, and statistical calculation needs to be carried out according to the different ecosystems, so that in the scheme, remote sensing data inversion products are applied to replace the parameters, namely when data are visually processed, an image substitution method is adopted, each pixel point has different calculation results, the problem that the different types of ecosystems use different experience parameters can be solved, the preprocessing process of recognizing the ecosystems is avoided, factors interfering with final calculation results are reduced, and the precision of the total value calculation results produced by the ecosystems is improved;

(2) in the technical scheme, data sources comprise statistical data, water resource bulletin, sentinel No. 2 multispectral data, MODIS remote sensing inversion product data, regional temperature change data and the like, the data are government-published data, global shared data or easily-acquired weather forecast data, the data can be acquired more conveniently, and the situations that the calculation operation is complex due to difficulty in data acquisition or the calculation result is inaccurate due to the fact that some data cannot be acquired do not exist;

(3) the method utilizes the characteristic that the remote sensing image can be visualized in space, assigns the statistical values of different types of parameters on specific pixel points by a visualization method such as a coefficient method, an image substitution method, a weighted distribution method and the like, enables the statistical data which is originally only a single value to be implemented on each pixel point to obtain a final visualized total value image produced by the ecological system, presents the GEP calculation result of a certain area in an image form, is more visual than a form, and is easier to compare different land blocks, and can be used for analyzing specific reasons influencing the total value produced by the ecological system because the image is obtained by superposing numerical values corresponding to various indexes related to the ecological system, and can play a role in dealing with global climate change and realizing the guidance of ecological environment and sustainable development of human society;

(4) according to the technical scheme, after data of indexes related to the ecosystem are replaced by a remote sensing data source, some accounting data can be refined, the remote sensing data can generate monthly data, 8-day synthesis data, daily data and the like according to different satellites, time resolution is improved, time change intervals of the data are shortened, the effect of improving updating frequency is achieved, the updating frequency of GEP accounting is improved, the total value image produced by the ecosystem finally generated by the scheme has real-time performance, and value quantity calculation errors in the presented image cannot be caused by data lag;

(5) according to the technical scheme, the remote sensing image is introduced to calculate the GEP, the parameters can be determined to each pixel point, the calculation scale at the moment is closely related to the spatial resolution of the remote sensing product, the precision of the calculation scale is far higher than the statistical result, GEP accounting can be performed on a certain small range, factors influencing the GEP can be analyzed, the large-range GEP can be calculated only, the small-range GEP can be calculated, and the application range of GEP calculation is widened.

Drawings

FIG. 1 is a flowchart of a method for visually presenting a total production value of a plot ecosystem according to an embodiment of the present invention

FIG. 2 is a schematic diagram of a selected evaluation plot in an embodiment of the present invention

FIG. 3 is a schematic view of the discharge amount of pollutants from municipal sewage in the embodiment of the invention

FIG. 4 is a schematic diagram of the excretion coefficients of livestock and poultry recommended by Ministry of environmental protection in the embodiment of the present invention

FIG. 5 is a schematic diagram of a water quality purification value coefficient according to an embodiment of the present invention

FIG. 6 is a schematic diagram of an ecosystem production total value image in an embodiment of the present invention

FIG. 7 is a schematic structural diagram of a visual presentation device for realizing a total production value of a plot ecosystem according to an embodiment of the present invention

FIG. 8 is a schematic structural diagram of an ecological product value calculation module 200 in an apparatus for visually presenting a total production value of a plot ecosystem according to an embodiment of the present invention

FIG. 9 is a schematic structural diagram of a second value calculation module 220 in a device for visually presenting a total production value of a plot ecosystem according to an embodiment of the present invention

FIG. 10 is a schematic structural diagram of an ecosystem production total value image generation module 300 in a device for visually presenting a total value of a production of a plot ecosystem according to an embodiment of the present invention

FIG. 11 is a schematic view of an electronic device in a device for visually presenting a total production value of a plot ecosystem according to an embodiment of the present invention

Detailed Description

The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, rather than all embodiments, and all other embodiments obtained by those skilled in the art without any inventive work based on the embodiments of the present invention belong to the protection scope of the present invention.

The terms "first," "second," and the like in the claims and in the description of the application are used for distinguishing between similar elements and not necessarily for describing a particular sequential or chronological order, it being understood that the terms so used are interchangeable under appropriate circumstances and are merely used to describe a different manner of distinguishing between similar elements in the embodiments of the application and that the terms "comprising" and "having," and any variations thereof, are intended to cover a non-exclusive inclusion such that a process, method, system, article, or apparatus that comprises a list of elements is not necessarily limited to those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus.

Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs, and the terms used herein in the specification of the present application are for the purpose of describing particular embodiments only and are not intended to limit the present application.

Example 1

As shown in fig. 1, a method for visually presenting a total production value of a land ecosystem includes the following steps:

s11, selecting a fixed land, formulating an ecological product list of the fixed land and acquiring related data to obtain first data, wherein the related data comprises remote sensing data, basic geographic data, literature data, statistical yearbook and meteorological data;

s12, calculating the value quantity of the first data according to the ecological product list by a visualization method to obtain second data, wherein the ecological product list comprises material products, adjusting services and cultural services;

and S13, according to the second data, carrying out superposition calculation on the value quantity of the material product, the value quantity of the adjusting service and the value quantity of the culture service through a grid calculator to obtain an ecosystem production total value image of the fixed plot.

In this embodiment, a taizhou city-road-bridge area is taken as an example for explanation, and specifically includes the following contents:

1. determining an evaluation range area:

here, the checking work is performed by taking a taizhou city-road-bridge area as an example, and therefore, the range is determined as a road-bridge area, and the range is shown in fig. 2.

2. Compiling an ecological product list:

the method has the advantages that the difficulty degree of data acquisition and the minimum number of evaluation items in the national standard are considered, and the list is divided into three dimensions, namely material products, adjustment services and culture services; among the material products, agricultural products, forestry products, animal husbandry products, fishery products, fresh water resources; the regulation service comprises water source conservation, soil maintenance, carbon fixation, oxygen supply, atmosphere purification, water quality purification and climate regulation; in cultural services, including leisure travel.

3. Data collection and download, i.e. "get relevant data":

the collection of data comprises basic geographic data, water resource bulletin, literature data, meteorological data, statistical yearbook and remote sensing data, namely 'first data';

wherein the basic geographic data typically uses network public data; water resource communiques are commonly available from the water bureau; the literature data can be obtained from the learning network and the accounting standard; meteorological data is typically available from meteorological stations; the yearbook is usually obtained from the bureau of statistics; the remote sensing data has more sources, such as GEE platform, ESA official network, NASA official network and the like.

4. Introduction to visualization methods:

different methods are used for performing visual transformation on different accounting projects.

4.1 coefficient method (coefficient distribution method):

the coefficient method is used for carrying out accounting work aiming at four items of agriculture, forestry, animal husbandry and fishery products, fresh water resources, leisure travel and water quality purification. Data using the coefficient method is often closely related to human activities and cannot be obtained from remote sensing channels.

Firstly, obtaining various branch output values from a table for counting the total output values of agriculture, forestry, animal husbandry and fishery in each county and city area in the yearbook, obtaining the local fresh water resource usage amount from a water resource gazette, and obtaining the total tourism income of the area according to the international and domestic tourism conditions of each county and city area in the statistics yearbook; secondly, area statistics is carried out by the evaluation range vector of the research area; and finally, dividing the value quantity by the area to obtain the value quantity of a certain product in unit area.

The unit values of agriculture, forestry, pasture, fishing, fresh water and tourism are respectively calculated, namely the corresponding value quantity provided by each square meter of land is respectively as follows:

average value of agricultural product value: 3.4313 yuan/square meter

Average value of forestry product value: 0.0145 yuan/sq m

Average value of the animal husbandry product value: 0.0618 yuan/sq m

Average value of fishery product value: 4.0077 yuan/square meter

Mean value of fresh water resource value: 0.4971 yuan/square meter

Ecological tourist value mean value: 25.065 yuan/square meter

Different from other three items, the water quality purification item can be refined to a more accurate street level, and similarly, the unit area pollution amount of three pollution amount indexes, namely COD (chemical oxygen demand), NH3-N (ammonia nitrogen) and TP (total phosphorus), is calculated, so that the unit area water quality purification value mean value used by each street is calculated. The average value of the water purification value per unit area of each street is not listed here, and the average value of the water purification value in the road and bridge area is only listed as follows:

average value of water purification value: 0.2230 yuan per square meter.

4.2 image substitution:

the image substitution method refers to a step of substituting statistical data or site measurement data with a remote sensing image product. For example, the evapotranspiration amount is obtained by replacing MOD16A2.ET products in MODIS satellites; the vegetation coverage factor (factor C) was calculated NDVI using sentinel No. 2 multispectral data, inverted to the vegetation coverage factor by empirical formula.

The factors using the method in the technical scheme are as follows:

evapotranspiration amount: MOD16A2.ET

Vegetation coverage factor: sentinel-2

Net primary productivity: MOD17A2H. PsnNet

Potential annual evapotranspiration: MOD16A2.PET

Vegetation transpiration consumes energy: MOD16A2.LE

4.3 weighted distribution method:

the weighted distribution method is characterized in that remote sensing product data are used as weights, a statistical value in a large range is distributed to each pixel point, and GEP accounting in a smaller range is achieved. For example, the statistical yearbook only has the total SO2 emission of the Taizhou area, but the weighting distribution can be carried out according to the concentration distribution image in sentinel-5p, and the target area road bridge area is cut after the image is drawn.

The factors using the method in the technical scheme are as follows:

emission of SO 2: statistical yearbook, Sentinel-5PSO2OFFL

NOx emission amount: statistical yearbook, Sentinel-5PNO2OFFL

5. Determining a functional quantity and value quantity accounting method, namely 'calculating the value quantity of the first data according to the ecological product list according to a visualization method to obtain second data':

according to the acquirability of data, only the functional quantity calculation is carried out on the adjustment service, and the value quantity calculation is directly carried out on the material product and the culture service according to a coefficient method, namely farming, forestry, herding, fresh water resources and ecological tourism.

Average value of agricultural product value: 3.4313 yuan/square meter

Average value of forestry product value: 0.0145 yuan/sq m

Average value of the animal husbandry product value: 0.0618 yuan/sq m

Average value of fishery product value: 4.0077 yuan/square meter

Mean value of fresh water resource value: 0.4971 yuan/square meter

Ecological tourist value mean value: 25.065 yuan/square meter

The index value amounts are summarized as "first value amount".

5.1 conservation of Water resources

And calculating the water conservation functional quantity by using a water quantity balance equation. The water balance equation refers to that water keeps mass conservation in the ecological system in a certain space and time, namely the water conservation quantity of the ecological system is the difference value between rainfall input and storm runoff and the water consumption quantity of the ecological system.

Qwr＝(Pi-Ri-ETi)*10^-3

In the formula: qwr conservation capacity (m) of water source³Pi is runoff yield (mm/a), Ri is surface runoff (mm/a), and ETi is evapotranspiration (mm/a);

from Taizhou Water resources bulletin 2019The flow and rainfall capacity of the bridge area is 1988.2mm/a, the surface runoff capacity is 1294.412mm/a, the evapotranspiration capacity is MOD16A2.ET (image substitution method), namely, a water source conservation functional quantity image of the bridge area is obtained through calculation, after the functional quantity is obtained, the value quantity calculation needs to be carried out according to a value quantity calculation method, and the formula is as follows: v_wr＝Q_wr×C_weIn the formula: v_wrValue (Yuan/a) for conservation of Water resources, Q_wrFor the total water source conservation (m) in the accounting area³/a)，C_weEngineering cost and maintenance cost (yuan/m) of unit storage capacity of reservoir³)；

The long pond is located at the upstream of the Yongning river of the Zaojiang tributary, and the dam site of the long pond is 22 kilometers away from the yellow rock urban area between the long pond mountain and the Fuhu mountain in the yellow rock area. The total investment of reservoir engineering is 2837 ten thousand yuan, the construction is carried out in 10 months in 1958, and the construction is completed in 12 months in 1964; the total investment of reservoir danger eliminating and strengthening engineering is 17761.58 ten thousand yuan, which is started in 10 and 16 months in 2002 and completed in 10 and 15 months in 2004 (completion acceptance at 3 months end in 2005).

Total index of substrate prices is determined by 'Taizhou city statistical yearbook 2020(2019 data)' 11-2 main years, 542.1 in 2002, 830.9 in 2019 and 1.5327 in magnification. Since the total price index of the fixed matrix only dates back to 1978, the total production value of 2.74 million yuan in 1957, 5134.05 million yuan in 2019 and the multiplying power of 1873.7409 can be obtained by estimating the labor cost in 1958 according to the total production value of 1-9 major years. Therefore, the construction cost and the maintenance cost converted to 2019 are 534.3026 billion yuan, the construction cost and the maintenance cost per cubic meter are 72.99 yuan/cubic meter, and a water source conservation value image is obtained after the calculation is finished.

5.2 soil conservation

The soil preservation function is that the ecological system (such as forest, grassland, etc.) protects the soil, reduces rainfall erosion power, increases soil corrosion resistance, reduces soil loss and preserves the soil through all levels such as canopy layer, litter, root system, etc. Q_srR × K × L × S × (1-C), wherein: q_srFor the total soil conservation (t/a), R is rainfall erosive power factor (MJ.mm/(hm)²H.a)), K is a soil erodible factor (t.hm)²·h/(hm²mJ. mm)), L is the slope lengthThe factor S is a gradient factor, and C is a vegetation coverage factor;

the R factor is calculated from the rainfall of many years, the K factor is calculated from the slope length and the like, and the L, S factor is calculated from DEM elevation data. The four items are dynamic monitoring result data which can be directly obtained. The C factor uses the multispectral data (image substitution method) of sentinel No. 2, and refers to a formula given in a paper research for predicting soil erosion amount of small watershed by applying a USLE model and a geographic information system IDRISI:

C＝1,C＝0

C＝0.6508-0.3436lgc,0＜C＜78.3％

C＝0,C＞78.3％

wherein C is the vegetation coverage, and this value is calculated using the pixel dichotomy, namely:

VFC＝(NDVI-NDVIsoil)/(NDVIveg-NDVIsoil)

in the formula, VFC is vegetation coverage, NDVI is a normalized vegetation index value of a current pixel point, NDVIsoil is a normalized vegetation index value of bare soil in a current area, and NDVIvegg is a normalized vegetation index value of a position with dense vegetation coverage in the current area. The normalized vegetation index (NDVI) is calculated by the following formula:

NDVI＝(NIR-RED)/(NIR+RED)

in the formula, NIR is the reflectivity of a near infrared band, and RED is the reflectivity of a RED light band.

The formula in the paper needs to be modified, i.e. the threshold for C is changed to: 0< c < 0.0964%; 0.0964% < c < 78.3%; c > 78.3%. Note that when C is substituted into the formula, no percentile is taken. When the calculated C factor has a negative value, the negative value is reassigned to 0.

After the total soil retention amount is calculated, obtaining a grid image of the functional soil retention amount, and calculating the value amount according to the functional amount, wherein the formula is as follows: v_sr＝λ×(Q_sr/[ rho ]). times.C, wherein: v_srFor maintaining the total value (yuan/a) of the soil, lambda is the silt deposition coefficient, Q_srFor the total soil retention (t/a), ρ is the soil bulk weight (t/m)³) C is the cost (Yuan/m) of the reservoir unit capacity dredging project³)；

And (3) taking the sediment deposition coefficient as 1, taking the soil volume weight as 1.5, taking the desilting engineering cost of the unit reservoir capacity of the reservoir as 26.27, and calculating to obtain a soil retention value image.

5.3 carbon fixation

The carbon fixation function refers to a function that a natural ecosystem absorbs carbon dioxide (CO2) in the atmosphere to synthesize organic matters and fixes carbon in plants or soil. The function is beneficial to reducing the concentration of carbon dioxide in the atmosphere and slowing down the greenhouse effect. The carbon fixation function of the ecological system has important significance for reducing emission reduction pressure. Here, the net ecosystem productivity method is used for estimation. Net Ecosystem Productivity (NEP) is an important scientific index for quantitative analysis of ecosystem carbon source/sink, and the carbon fixation amount of ecosystem can be measured by NEP. NEP is widely used in carbon cycle studies, and can be obtained by subtracting iso-aerobic respiration from Net Primary Productivity (NPP) or by conversion according to the relative conversion coefficient of NPP and NEP, and then the quality of the carbon dioxide fixed by the land ecosystem is measured.

In the formula:

for total carbon dioxide fixation (t. CO) in land ecosystem₂/a)，

Is CO₂The ratio of molecular weight to C, i.e., 44/12, NEP is net ecosystem productivity (t.C/a); NEP ═ NPP-RS, where NPP is net primary productivity (t · C/a) and RS is carbon consumption by soil respiration (t · C/a);

the NPP data used the psnent layer product (image substitution) in MOD17A2H, which is the net primary productivity layer. And accumulating the original data all year round, multiplying the accumulated data by a proportionality coefficient of 0.0001, and outputting a net primary productivity map layer. The value is expressed in kg C/(m 2 a), which is the number of kilograms of carbon immobilized in one year per square meter. The soil foundation respiration is inverted, and the concrete formula is as follows:

namely, it is

In the formula, A_qFor respiration of soil foundation, PPT

R_H＝NPP

The annual precipitation is, PET is the annual potential evapotranspiration, b is a temperature sensitive constant factor, and the value is 0.0306;

wherein the NPP is a PsnNet layer in MOD17A2 product, PPT data is 1988.2mm, PET data is a PET (image substitution method) layer in MOD16A2 product, b value is 0.0306, and T is annual average temperature is 16.9 ℃. And calculating to obtain a soil respiration map layer with the unit of kg C/(m ^2 a).

And after the carbon fixed functional quantity map layer is calculated by integrating the data, the carbon fixed value quantity map layer is calculated according to the carbon trading price of 23.72 yuan/ton.

In the formula: v_cfFor the total carbon sequestration value (Yuan/a),

the total fixed carbon dioxide amount (t/a),

is the carbon trade price (dollar/t).

5.4 oxygen supply

The oxygen release function of the ecosystem refers to the function of the plants to release oxygen in the process of photosynthesis. The function has important significance for maintaining the stability of oxygen in the atmosphere and improving the living environment.

In the formula: q_OPReleasing oxygen (tO) for ecosystem₂/a)，

Namely CO₂Conversion to O₂The coefficient of (a) is calculated,

carbon sequestration for land ecosystem (tCO)₂/a)；

The oxygen supply functional quantity is calculated by a carbon fixed functional quantity map layer (image substitution method) in 5.3, and an oxygen supply value quantity image is calculated according to the medical oxygen generation price of 1200 yuan/ton after the calculation is finished.

V_OP＝Q_OP×C_oIn the formula: v_OPQ is the oxygen release value (yuan/a) of the ecosystem_OPFor the oxygen release (tO) of the ecosystem₂/a)，C_oOxygen prices (yuan/t) for medical treatment.

5.5 atmospheric purification (weighted distribution method)

The purification amount of the atmospheric pollutants is used as an accounting index, the indexes of pollutants such as sulfur dioxide, nitrogen oxide and the like are selected according to the regulation of environmental air quality control projects in GB3095-2012, and when the environmental air quality is equal to or superior to the national second-level, the purification amount of the atmospheric pollutants is the pollutant emission amount, namely, a pollutant mass balance model is used.

According to annual book data counted in Taizhou city, the emission of SO2 is 6034 tons, and the emission of NOx is 9567 tons. In order to visualize the statistical data, the visualization was performed by using the SO2 and NO2 concentration distribution images in Sentinel-5 pOFFL. And distributing the total emission by taking the pixel values of all positions in the Taizhou city as weights, cutting a road and bridge area, and outputting to obtain the SO2 functional quantity of the atmospheric purification and the NO2 functional quantity of the atmospheric purification, wherein the unit of the SO2 functional quantity and the NO2 functional quantity is kg/m ^ 2.

The formula for calculating the value of the atmosphere purification is as follows:

in the formula: v_apFor the total value (Yuan/a), Q of air purification_apiIs the purification amount (t/a) of the atmospheric pollutants in the ith, C_iThe treatment cost (Yuan/t) of the ith type of atmospheric pollutants, and n is the number of types of atmospheric pollutants in the accounting region.

The treatment price of SO2 is 2000 yuan/ton, and the treatment price of NOx is 2518.25 yuan/ton. And calculating to obtain an atmospheric purification value image.

5.6 Water purification (coefficient method)

And when the pollutant emission concentration does not exceed the environmental function standard limit value of the surface water area, calculating the purification amount of various pollutants by the ecological system in the area through calculating the emission amount according to the pollutant mass balance model.

In the formula, Q_wpTotal amount of pollutants (Kg/a), Q_eiThe total amount (Kg/a) of the i-th pollutants, Q_aiThe total discharge amount (Kg/a) of the i-th pollutants in the area, Q_diIs the exit volume (Kg/a) of the i-th pollutants, Q_siAmount of i-type pollutant (Kg/a) treated by sewage treatment plant, i being pollutant class, i ═ 1, 2_aiIncluding rural life (W)_n) City life (W)_t) Agricultural non-point source pollution (W)_m) Cultivation pollution (W)_a)；

Wherein, the amount of the pollutants entering and exiting is not considered, the industrial pollution is not considered, and only three items of COD and NH3-N, TP are counted in the pollutant category. The non-point source pollution and the breeding pollution of rural life, urban life and agriculture are calculated.

Rural domestic sewage discharge is determined by referring to strong data of pollution source investigation sources determined by the national environmental protection department, and the average sewage discharge amount of people is 80 liters/person/day, the COD discharge amount is 16.4 grams/person/day, the ammonia nitrogen discharge amount is 4.0 grams/person/day, the total nitrogen discharge amount is 5.0 grams/person/day, and the total phosphorus discharge amount is 0.44 grams/person/day.

The pollutant discharge of urban domestic sewage is determined by the handbook of pollution discharge coefficient of urban domestic resource for the second national pollution source survey, as shown in fig. 3.

The number of breeding in the breeding pollution is obtained by a statistical yearbook, and is calculated by using a livestock and poultry excretion coefficient recommended by the ministry of environmental protection, and the pollution data obtained by calculation is distributed to all villages and towns according to the area, wherein the livestock and poultry excretion coefficient is shown in fig. 4.

The monitoring range of agricultural non-point source pollution is respectively calculated according to the land utilization type in the dynamic monitoring result, and the functional quantity is calculated according to a 'handbook for general survey of agricultural pollution source fertilizer loss coefficients for national pollution sources for the first time'. The simplified version calculation coefficient is obtained through weighted calculation, and the unit is kilogram/mu, as shown in table 1:

	ammonia nitrogen	Total phosphorus
			Paddy field	0.37	0.0415
Watering the ground	0.017	0.047
			Dry land	0.07	0.0455
(Orchard)	0.029	0.033
			Tea garden	0.029	0.033

TABLE 1

Summing agricultural non-point source pollution, rural living pollution, urban living pollution and cultivation pollution to obtain a water quality purification function quantity coefficient, and calculating the water quality purification value quantity coefficient according to the treatment cost of COD (chemical oxygen demand) 8000 yuan/ton, the treatment cost of NH 3-N9572.92 yuan/ton and the treatment cost of TP 10000 yuan/ton. The coefficients are shown in fig. 5. Because the street coefficients are different, in order to facilitate calculation, the water quality purification coefficient is assigned to a corresponding vector area in combination with an administrative division, the area is a 'first vector area', a tool 'vector is used for converting a vector into a grid', grid data, namely a water quality purification value image, is obtained, and is kept consistent with the value data form of 5.1-5.7, and calculation is facilitated.

5.7 climate control (image substitution method)

The ecological system climate regulation service refers to an ecological function that the ecological system absorbs solar energy through vegetation transpiration and water surface evaporation, reduces air temperature, increases air humidity and improves the comfort degree of human living environment. The energy consumed in the evaporation process of the ecological system is selected as the evaluation index of the climate regulation service of the ecological system. The method comprises the following steps: the energy consumed by vegetation transpiration and the energy consumed by water surface evaporation.

E_tt＝E_pt+E_weIn the formula: e_ttTotal energy (KWh/a), E consumed for transpiration of the ecosystem_ptEnergy consumed for transpiration of vegetation (KWh/a), E_weEnergy expended for water surface evaporation (KWh/a);

the energy consumed by vegetation transpiration is obtained from a remote sensing product image of MOD16A2.LE (image substitution method), and the energy consumed by water surface evaporation is calculated according to the following formula.

E_we＝E_w×q×ρ×10³/3600+E_wX y, wherein: e_wIs the evaporation capacity (m) of water surface³A) and rho is the density of water (1 g/cm)³) Q is the latent heat of vaporization, i.e., the amount of heat (J/g) required to evaporate 1 gram of water, and y is the power consumed by the humidifier to convert 1 cubic meter of water to steam (KWh/m)³)；

Wherein, the evaporation capacity of the water surface is calculated by MOD16A2.ET (image substitution method), the density of the water is 1 g/cubic centimeter, the latent heat of volatilization of the water at 25 ℃ is 44KJ/mol, and the power consumption of the humidifier for converting 1 cubic meter of water into steam is 900 kilowatt-hours.

And after the calculation is finished, obtaining a climate regulation function quantity map layer, and calculating a grid image of the climate regulation value quantity according to the residential electricity price of 0.54 yuan/kilowatt hour.

In summary, the five value quantity images of the water source conservation value quantity image, the soil conservation value quantity image, the carbon fixation value quantity image layer, the oxygen supply value quantity image and the grid image of the climate regulation value quantity, which are respectively obtained in the water source conservation, the soil conservation 5.2, the carbon fixation 5.4, the oxygen supply 5.4 and the climate regulation 5.7, are summarized, namely, the first value quantity image; 5.5 obtaining an atmospheric purification value image, namely a second value image in atmospheric purification; 5.6 water quality purification value quantity image in water quality purification, namely 'third value quantity image'; summarizing the first valence quantity image, the second valence quantity image and the third valence quantity image to obtain total valence quality, namely 'second valence quantity'; 5, summarizing the value mean values of agriculture, forestry, animal husbandry, fishery, fresh water resources and ecological tourism, namely a first value amount; the first and second value amounts are summed, i.e., "second data".

6. Ecosystem production total accounting

The air purification value, the climate regulation value, the water source conservation value, the water quality purification value, the carbon fixation value, the soil retention value and the oxygen supply value are obtained from 5.1 to 5.7, and the total number of the 7 raster images is 7, and the unit is Yuan/sq meter.

Performing grid calculation and summation on the 7 images to obtain a grid map of the adjusting service value, namely a grid calculator executes the steps of 'atmospheric purification value + climate adjusting value + water source conservation value + water quality purification value + carbon fixed value + soil retention value + oxygen supply value', and then 'grid map of the adjusting service value', namely 'performing grid calculation on the first value image, the second value image and the third value image, and then summing to obtain a value image of the adjusting service'; and then, superposing the average value of the unit area value of agriculture, forestry, animal husbandry, ecological tourism and fresh water resources, namely, executing a grid calculator to obtain a grid image of the total value of the ecological system production, namely, superposing the value image of the adjustment service and the first value according to the grid calculator to obtain the grid image of the total value of the ecological system production, wherein the grid image of the adjustment service is obtained by executing a grid image of the adjustment service value +3.4313+0.0145+0.0618+4.0077+0.4971+ 25.065'. The process of the above 6 is summarized, that is, the value amount of the material product, the value amount of the adjustment service and the value amount of the culture service are superposed and calculated by a grid calculator according to the second data to obtain the total value image of the ecosystem production of the fixed land parcel. In the formula, the numerical values have the following meanings:

average value of agricultural product value: 3.4313 yuan/square meter

Mean value of forestry product values: 0.0145 yuan/square meter

Average value of the animal husbandry product value: 0.0618 yuan/sq m

Average value of fishery product value: 4.0077 yuan/square meter

Mean value of fresh water resource value: 0.4971 yuan/square meter

Ecological tourist value mean value: 25.065 yuan/square meter

As shown in fig. 6, for the total value result chart of ecosystem production, the GEP height per unit area is displayed in a visual form, and the spatial resolution is 10 meters.

The embodiment has the following beneficial effects:

Example 2

As shown in fig. 7, a visual presentation device for total production value of a land ecosystem includes:

the parcel data acquisition module 100: the system comprises a data acquisition module, a data processing module and a data processing module, wherein the data acquisition module is used for selecting a fixed plot, formulating an ecological product list of the fixed plot and acquiring related data to obtain first data, and the related data comprises remote sensing data, basic geographic data, literature data, statistical yearbook and meteorological data;

ecological product value amount calculation module 200: the system comprises a data acquisition module, a data processing module and a data processing module, wherein the data acquisition module is used for acquiring first data and second data, and the first data is calculated according to a visual method and value quantity according to an ecological product list to obtain the second data, and the ecological product list comprises material products, adjusting services and cultural services;

ecosystem production total value image generation module 300: and the ecological system production total value image of the fixed plot is obtained by performing superposition calculation on the value quantity of the material product, the value quantity of the adjusting service and the value quantity of the cultural service through a grid calculator according to second data.

In one embodiment of the above apparatus, in the plot data acquisition module 100, a fixed plot is selected, an ecological product list of the fixed plot is prepared and related data is acquired, so as to obtain first data, the related data comprises remote sensing data, basic geographic data, literature data, statistical yearbook and meteorological data, in an ecological product value calculation module 200, calculating the value of the first data according to the ecological product list by a visualization method to obtain second data, wherein the ecological product list comprises material products, adjustment services and cultural services, in the ecosystem production total value image generating module 300, the value amount of the material product, the value amount of the adjustment service and the value amount of the culture service are superposed and calculated through a grid calculator according to the second data, and the ecosystem production total value image of the fixed land parcel is obtained.

Example 3

As shown in fig. 8, an ecological product value calculating module 200 in a visualization presentation apparatus for total production value of a land ecosystem includes:

the first value amount calculation module 210: the system comprises a visualization method, a coefficient method, an image substitution method and a weighted distribution method, wherein the visualization method is used for calculating the value quantity of the material product and/or the cultural service according to the coefficient method and the first data to obtain a first value quantity;

the second value amount calculation module 220: and the system is used for calculating the value quantity of the adjusting service according to a visualization method and the first data to obtain a second value quantity, wherein the adjusting service comprises water source conservation, soil conservation, carbon fixation, oxygen supply, atmosphere purification, water quality purification and climate adjustment.

In one embodiment of the above apparatus, the first value calculation module 210 calculates the value of the material product and/or the cultural service according to a coefficient method and the first data to obtain a first value, the visualization method includes the coefficient method, an image substitution method and a weighted distribution method, and the second value calculation module 220 calculates the value of the adjustment service according to the visualization method and the first data to obtain a second value, and the adjustment service includes water conservation, soil conservation, carbon fixation, oxygen supply, atmospheric purification, water purification and climate adjustment.

Example 4

As shown in fig. 9, a second value amount calculating module 220 in a device for visually presenting a total production value of a land ecosystem includes:

the image first calculation module 221: the system comprises a remote sensing data product, a grid calculator and a data processing and calculating unit, wherein the remote sensing data product is used for replacing data of water source conservation, soil conservation, carbon fixation, oxygen supply and climate regulation according to the image substitution method, and the image of the remote sensing data product is processed and calculated by the grid calculator to obtain a first value image;

the image second calculation module 222: the system is used for carrying out numerical value distribution on the atmospheric purification data according to remote sensing image pixel points of the fixed plot according to a weighted distribution method, and processing and calculating the image subjected to numerical value distribution according to a grid calculator to obtain a second value quantity image;

the image third calculation module 223: and the system is used for calculating the value quantity of water purification according to the first data by a coefficient method to obtain a water purification value quantity coefficient, assigning the water purification value quantity coefficient in a vector area of the fixed plot to obtain a first vector area, and performing vector-to-grid processing to obtain a third value quantity image.

One embodiment of the above apparatus is that, in the image first calculation module 221, the remote sensing data product is used to replace the data of the water source conservation, the soil conservation, the carbon fixation, the oxygen supply and the climate adjustment according to the image substitution method, and the image of the remote sensing data product is processed and calculated by the grid calculator to obtain a first value quantity image, in the image second calculation module 222, the data of the atmosphere purification is subjected to numerical distribution according to the remote sensing image pixel points of the fixed land blocks according to the weighted distribution method, and the image after the numerical distribution is processed and calculated according to the grid calculator to obtain a second value quantity image, in the image third calculation module 223, the value quantity of the water quality purification is calculated according to the first data and by the coefficient method to obtain a water quality purification value quantity coefficient, and assigning the water quality purification value coefficient in a vector area of the fixed plot to obtain a first vector area, and performing vector-to-grid processing to obtain a third value image.

Example 5

As shown in fig. 10, an ecosystem production total value image generating module 300 in a device for visually presenting a land ecosystem production total value includes:

the first value image generation module 310: the value quantity image used for summing the first value quantity image, the second value quantity image and the third value quantity image after grid calculation to obtain the value quantity image of the adjusting service;

the second value image generation module 320: and the grid calculator is used for superposing the value quantity image of the adjusting service and the first value quantity to obtain a total value image of the ecological system production.

One embodiment of the foregoing apparatus is that, in the first value image generating module 310, the first value image, the second value image and the third value image are subjected to grid calculation and then summed to obtain the value image of the adjustment service, and in the second value image generating module 320, the value image of the adjustment service and the first value are superimposed according to a grid calculator to obtain the total value image of the ecosystem production.

Example 6

As shown in fig. 11, an electronic device comprises a memory 601 and a processor 602, wherein the memory 601 is used for storing one or more computer instructions, and the one or more computer instructions are executed by the processor 602 to implement any one of the methods described above.

It can be clearly understood by those skilled in the art that, for convenience and brevity of description, the specific working process of the electronic device described above may refer to the corresponding process in the foregoing method embodiment, and is not described herein again.

A computer readable storage medium storing a computer program which, when executed, causes a computer to implement any of the methods as described above.

Illustratively, a computer program may be divided into one or more modules/units, one or more modules/units are stored in the memory 601 and executed by the processor 602, and the input interface 605 and the output interface 606 perform I/O interface transmission of data to complete the present invention, and one or more modules/units may be a series of computer program instruction segments for describing the execution of the computer program in a computer device, which can perform specific functions.

The computer device may be a desktop computer, a notebook, a palm computer, a cloud server, or other computing devices. The computer device may include, but is not limited to, the memory 601 and the processor 602, and those skilled in the art will appreciate that the present embodiment is only an example of the computer device, and does not constitute a limitation of the computer device, and may include more or less components, or combine some components, or different components, for example, the computer device may further include the input device 607, the network access device, the bus, and the like.

The processor 602 may be a Central Processing Unit (CPU), other general-purpose processor 602, a digital signal processor 602 (DSP), an Application Specific Integrated Circuit (ASIC), an off-the-shelf programmable gate array (FPGA) or other programmable logic device, discrete gate or transistor logic device, discrete hardware component, etc. The general purpose processor 602 may be a microprocessor 602 or the processor 602 may be any conventional processor 602 or the like.

The storage 601 may be an internal storage unit of the computer device, such as a hard disk or a memory of the computer device. The memory 601 may also be an external storage device of the computer device, such as a plug-in hard disk, a Smart Memory Card (SMC), a Secure Digital (SD) card, a flash memory card (FlashCard) and the like provided on the computer device, further, the memory 601 may also include both an internal storage unit and an external storage device of the computer device, the memory 601 is used for storing computer programs and other programs and data required by the computer device, the memory 601 may also be used for temporarily storing in the output device 608, and the aforementioned storage media include various media capable of storing program codes, such as a usb disk, a removable hard disk, a read only memory ROM603, a random access memory RAM604, a disk and an optical disk.

The above description is only an embodiment of the present invention, but the technical features of the present invention are not limited thereto, and any changes or modifications within the technical field of the present invention by those skilled in the art are covered by the claims of the present invention.

Claims

1. A visual presentation method for total production value of a land ecosystem is characterized by comprising the following steps:

2. The method for visually presenting total production value of a land ecosystem according to claim 1, wherein the calculating the value amount of the first data according to the ecological product list according to the visualization method to obtain the second data comprises:

3. The method as claimed in claim 2, wherein the step of calculating the value amount of the adjustment service according to the visualization method and the first data to obtain a second value amount comprises:

4. The method as claimed in claim 3, wherein the step of calculating the total value of the ecosystem production of the fixed plot by superimposing the value amount of the material product, the value amount of the adjustment service and the value amount of the cultural service by the grid calculator according to the second data comprises:

5. A visual presentation device for total production value of a land ecosystem, which is used for realizing the visual presentation method for total production value of a land ecosystem according to claim 1, and is characterized by comprising the following steps:

6. The apparatus for visually presenting total production value of a land ecosystem according to claim 5, wherein the ecological product value calculation module comprises:

7. The apparatus for visually presenting total production value of a land ecosystem according to claim 6, wherein the second value amount calculation module comprises:

the first image calculation module: the system is used for replacing the data of water source conservation, soil conservation, carbon fixation, oxygen supply and climate regulation by using a remote sensing data product according to the image substitution method, and processing and calculating the image of the remote sensing data product through a grid calculator to obtain a first value image;

the second image calculation module: the remote sensing image processing system is used for carrying out numerical value distribution on the data of the atmospheric purification according to remote sensing image pixel points of the fixed land parcel according to a weighted distribution method, and processing and calculating the image after the numerical value distribution according to a grid calculator to obtain a second value quantity image;

8. The apparatus for visually presenting total value of ecosystem production in a land according to claim 7, wherein the ecosystem total value of production image generation module comprises:

9. A visual presentation apparatus of a total value produced by a land ecosystem, comprising a memory and a processor, wherein the memory is used for storing one or more computer instructions, and the one or more computer instructions are executed by the processor to realize the visual presentation method of the total value produced by the land ecosystem according to any one of claims 1 to 4.

10. A computer-readable storage medium, in which a computer program is stored, which, when being executed by a computer, carries out a method for the visual presentation of a total value of a production of a land ecosystem according to one of claims 1 to 4.