An adaptive scale ecological product value accounting system

Abstract

The application discloses a self-adaptive scale ecological product value accounting system, relates to the technical field of specific regional ecological product value accounting, and solves the problem of insufficient self-adaptive degree of accounting data source and accounting scale in the prior art.The application comprises a multi-source data base module, an ecological product type automatic discrimination module, an ecological product value integrated accounting module, a multi-scale self-adaptive aggregation and interaction module and a map visualization and special topic display module.The system takes a map spot as the minimum accounting unit, automatically discriminates the ecological product type according to the map spot land class code, and completes the index accounting including material supply, regulation service and cultural service in a unified framework;based on the spatial topological attribution relationship between the map spot and the management units at all levels, the accounting result of any specified spatial unit such as administrative village, township, county, single project and the like is automatically aggregated from bottom to top in stages and output on demand, and the accounting result is stored in a vector attribute form.

Description

Technical Field

[0001] This invention relates to the field of ecological product value accounting technology in specific regions, and particularly to an adaptive-scale ecological product value accounting system. Background Technology

[0002] Ecological products refer to the collective products and services provided by ecosystems for human production and life, mainly including three categories: material supply, regulatory services, and cultural services. Ecological product valuation (also known as Gross Ecosystem Product (GEP) valuation) aims to measure, in a monetized manner, the total economic value of the final products and services provided by an ecosystem in a specific region for human well-being within a certain period. It has significant theoretical and practical value for establishing and improving the mechanism for realizing the value of ecological products and promoting ecological civilization.

[0003] Currently, various technical methods have been developed both domestically and internationally in the field of ecological product value accounting, mainly including biophysical modeling, equivalent factoring, emergy analysis, and statistical reporting. Regarding tools and platforms, internationally, assessment tools such as InVEST, ARIES, SolVES, and SWAT are available; domestically, operational accounting systems such as the Shenzhen GEP automated accounting platform and the Lishui GEP accounting system have been established. In terms of standards and specifications, the United Nations has published the "Environmental Economic Accounting System – Ecosystem Accounting" (SEEA-EA, 2021), and domestically, normative documents such as the "Specifications for Gross Ecosystem Product Accounting (Trial)" (2022) and the "Technical Guidelines for Gross Ecosystem Product (GEP) Accounting" (2020) have been issued. Various provinces and cities have also published more than 60 local technical specifications.

[0004] However, existing technologies generally use fixed administrative regions or fixed-resolution raster pixels as the accounting unit. The processing flow for multi-source heterogeneous data is cumbersome and non-standardized, requiring different accounting personnel to repeatedly perform basic tasks such as data collection, format conversion, and spatial registration, resulting in a large amount of repetitive work. When management decisions require accounting results at different spatial scales, data must be reorganized, parameters configured, and the accounting process run again for each target scale, multiplying the workload. Furthermore, the smallest accounting unit in existing methods is mainly county-level and above administrative regions. While biophysical modeling can output rasterized results, there is a lack of correspondence between raster pixels and actual land parcels in land management, making it unsuitable for refined management based on land parcels, such as ecological compensation and carbon trading. Existing multi-level accounting systems conduct accounting independently at each level, potentially using different data sources or parameter sets, leading to a lack of strict mathematical logic between results at each level and difficulty in ensuring vertical consistency. Therefore, in terms of presenting accounting results, existing methods mostly present results in the form of statistical reports or regional aggregates, lacking a precise correspondence between accounting results and geographic space. Even if some methods can output spatialized raster layers, they are limited to the physical quantity distribution of a single service type and cannot overlay and display the monetization value distribution of multiple types of ecological products on the same spatial base map. Most existing accounting platforms are customized for specific regions, and the accounting methods are highly coupled with specific regions, resulting in a lack of horizontal comparability of accounting results between different regions.

[0005] Therefore, there is an urgent need for an ecological product value accounting system that can be compatible with multi-source heterogeneous data, cover major ecological products, and support micro-level map patch accounting and multi-scale adaptive aggregation. Summary of the Invention

[0006] In view of the above-mentioned shortcomings of the existing technology, the purpose of this invention is to provide an adaptive scale ecological product value accounting system to solve the problems of insufficient data source and insufficient degree of adaptability of accounting scale in the existing technology.

[0007] This invention provides an adaptive scale ecological product value accounting system, which enables scale aggregation from unified fine units to multi-level management units through a standardized multi-source data foundation and an integrated value accounting system that covers all indicators and can be executed on micro-spatial units, thereby achieving flexible response to different management scales and customized ranges.

[0008] To achieve the above objectives, the present invention adopts the following technical solution: An adaptive-scale ecological product value accounting system includes: The multi-source data base module serves as the system's data entry point, responsible for collecting, formatting, unifying coordinates, and classifying various heterogeneous multi-source data required for accounting, forming a standardized data foundation. This data includes at least land survey data, remote sensing imagery data, meteorological and hydrological monitoring data, socio-economic statistics, and points of interest. Simultaneously, this module receives user-imported or specified geographic boundary vectors of the study area and uses these boundaries as constraints to decompose the accounting area into the system's unified smallest computational unit, namely, a micro-spatial computational unit, which is a vector patch. The output of this module provides standardized data services for all subsequent modules. The multi-source data base module also includes a patch preprocessing unit, used for topological checking, gap repair, and spatial overlay analysis of the acquired patch vector data to ensure no overlap or gaps between patches, forming a working base map that completely covers the accounting area spatially.

[0009] An automated ecological product type identification module is included. This module stores the mapping relationship between land categories and ecological product types. The rule base associates land category name attributes with ecological product types, enabling different land category patches to automatically match different combinations of ecological product types, and then derives accounting indicators based on the ecological product types. The ecological product types cover three major categories: material supply, regulatory services, and cultural services. Among them, the material supply category includes agricultural products, forestry products, livestock products, aquatic products, and water resource supply; the regulatory services category includes water conservation, soil retention, flood control, air purification, water purification, carbon sequestration, oxygen release, local climate regulation, provision of negative oxygen ions, and species conservation; and the cultural services category includes tourism and health care, and leisure and recreation. In summary, the automated ecological product type identification module uses the aforementioned land survey data as its core input. Based on preset land category-ecological product type identification rules, it automatically identifies and converts the land category name attribute of each micro-spatial calculation unit into the corresponding ecosystem type (such as farmland ecosystem, forest ecosystem, wetland ecosystem, etc.), forming a set of calculation units carrying ecosystem type attributes, which serves as the spatial basis for subsequent value accounting. Furthermore, the automated ecological product type identification module adopts a modular indicator management architecture. The accounting items in the identification rule base are organized in a configurable manner, supporting the addition, modification, and deactivation of accounting items without altering the core system logic.

[0010] Integrated Ecological Product Value Accounting Module: This module is the core computing engine of the system. Based on the ecosystem type and the type of ecological product to be accounted for in each micro-spatial computing unit, it calls the corresponding dedicated spatial accounting models in parallel or sequentially. Each spatial accounting model only calls the relevant parameters stored in the multi-source data base module, independently completing various value calculations. The calculation results are attached to the corresponding micro-spatial computing unit in units of map patches or pixels. This module covers 17 indicators in three major categories: material supply, regulatory services, and cultural services. Specifically, it includes specialized accounting sub-modules for agricultural products, forestry products, livestock products, fishery products, water resource supply, flood control, air purification, carbon sequestration, oxygen release, water purification, climate regulation, water conservation, soil retention, negative oxygen ions, species conservation, tourism and health care, and leisure and recreation.

[0011] The multi-scale adaptive aggregation and interaction module receives various value data attached to micro-spatial calculation units from the integrated ecological product value accounting module. This module supports spatial aggregation at different scales based on geographical or administrative boundaries, achieving "one-time fine-grained calculation, multi-scale output on demand." Specifically, based on any field specified in the user-defined spatial aggregation layer, it categorizes and summarizes various value quantities of micro-units into spatial units corresponding to those values, thus achieving adaptive aggregation output based on different administrative levels (such as administrative villages, townships, counties (districts, cities)), project scope, or demand dimensions. Simultaneously, this module provides statistical analysis functions in chart form, offering multi-dimensional visualization of the accounting results at each level and for each indicator.

[0012] Map Visualization and Thematic Display Module: As the system's visualization output, this module uses a geographic base map as a medium to spatially render and overlay the value of each micro-spatial calculation unit. It provides functions such as layer management, spatial query, thematic color configuration, and interactive element information viewing, supporting users in geospatial exploration and thematic display of the accounting results.

[0013] Preferably, the land survey data includes land category names and planting attributes.

[0014] Preferably, the integrated accounting module for ecological product value includes an agricultural product value accounting submodule. This submodule uses a set of micro-spatial calculation units determined to be farmland ecosystems as the basic accounting object. First, through the agricultural product crop type mapping unit, based on the land use name and planting attributes of the plots, it performs secondary subdivision according to predefined rules and maps them to specific agricultural product crop types. Subsequently, the plot-level value parameter quantification calculation unit obtains the theoretical yield per unit area Y, regional production correction coefficient C, and product market price P corresponding to each mapped plot for that crop type, and assigns an area adjustment coefficient A according to the crop type, and calculates the agricultural product value V_plot carried by the plot. The calculation formula is: V_plot=(S×A)×Y×C×P, where S is the area of ​​the plot.

[0015] Preferably, the integrated accounting module for ecological product value includes a forest product value accounting submodule. This submodule uses a set of micro-spatial calculation units identified as forest ecosystems as the basic accounting object. First, through the forest product type mapping unit, based on the land use name of the map patch, it performs secondary subdivision according to predefined rules, mapping it to specific forest product types, including at least timber and bamboo. Subsequently, the map patch-level value parameter quantification calculation unit obtains the key accounting parameters corresponding to each mapped map patch. For timber map patches, V_plot = S × Y_a × P_a is calculated, and for bamboo map patches, V_plot = S × Y_b × P_b is calculated, where Y_a is the timber harvesting volume per unit area, P_a is the timber unit price, Y_b is the bamboo harvesting volume per unit area, and P_b is the bamboo unit price.

[0016] Preferably, the integrated ecological product value accounting module includes a livestock product value accounting submodule. This submodule is based on a set of micro-spatial calculation units identified as farmland ecosystems, and only includes plots with land use designation as facility agricultural land within the scope of livestock product value accounting. For each selected plot, the unit area output Y_l and product unit price P_l are obtained, and the plot area S is adjusted using a preset area adjustment coefficient K, and V_plot=(S×K)×Y_l×P_l is calculated.

[0017] Preferably, the integrated ecological product value accounting module includes a fishery product value accounting submodule. This submodule uses a set of micro-spatial calculation units identified as wetland ecosystems as its foundation, and only includes plots with land use names of aquaculture ponds or adjustable aquaculture ponds within the scope of fishery product value accounting. For each selected plot, the output per unit area Y_f and the product unit price P_f are obtained, and V_plot = S × Y_f × P_f is calculated, where S is the area of ​​the plot.

[0018] Preferably, the integrated accounting module for ecological product value includes a tourism and health care value accounting submodule. This submodule includes: a spatial data association unit, used to spatially connect and overlay the land parcel vector layer with the tourism and health care point of interest (POI) database to form a parcel-POI association dataset; an POI value parameter matching unit, used to read a preset tourism and health care POI type-value parameter mapping table, which defines the unit value P_type corresponding to each type of POI; and a tourism and health care value calculation unit, used to, for each parcel containing POIs, count the number N_type of each type of POI within that parcel, calculate V_type = N_type × P_type, and sum the values ​​of all types to obtain the total tourism and health care value V_tourism of that parcel.

[0019] Preferably, the integrated ecological product value accounting module includes a recreational value accounting submodule. This submodule filters map patches based on a preset list of recreational ecosystem types, obtains the area S of each filtered map patch, and matches and obtains the unit area benchmark land price P_base of the spatial location of the map patch from the land benchmark price spatial database, and calculates V_recreation=S×P_base.

[0020] Preferably, the multi-scale adaptive aggregation and interaction module is configured to: classify and summarize the various values ​​of the micro-spatial calculation units calculated by the integrated accounting module for ecological product value into the spatial units corresponding to the values ​​of the field according to any field in the spatial aggregation layer specified by the user, thereby realizing adaptive aggregation output according to different administrative levels, project scopes or demand dimensions.

[0021] Compared with the prior art, the beneficial effects of the present invention are as follows: Beneficial effect one: It achieves scale-adaptive accounting capabilities at the patch level and multi-level adaptive aggregation, solving the problems of fixed accounting scale and coarse accuracy in existing technologies. Specifically: Existing ecological product value accounting technologies primarily use county-level and above administrative regions as the smallest accounting unit. While biophysical modeling methods (such as InVEST) can output rasterized results, raster pixels, as homogeneous regular grid units, lack a corresponding relationship with the map patches characterized by ownership, use, and land type in land management practice. This makes it difficult to support the refined management needs of ecological compensation, carbon trading, and other operations based on land parcels. Furthermore, once the accounting scale is determined, it cannot be flexibly adjusted, requiring repeated reconstruction of data and models for different scales. This invention uses land survey map patches as the smallest accounting unit, assigning each map patch clear spatial boundaries, area, land type, and ownership attributes, thus reducing the accounting accuracy to the land parcel level. This allows the accounting results to directly serve refined management scenarios such as land rights registration, ecological compensation fund allocation, and carbon trading. Based on this, the system establishes the spatial topological affiliation relationship between map features and administrative divisions at all levels. Only one map feature-level calculation is required to automatically aggregate from bottom to top along the path of "map feature → village → township → county" and output the calculation results of any level as needed. The system achieves adaptive switching of scale by performing one calculation and multiple-level aggregation, which significantly improves the accuracy of calculation while greatly reducing the workload of repeated calculations at multiple scales.

[0022] Benefit 2: It achieves a comprehensive spatial representation of the accounting results, making up for the shortcomings of existing technologies in that the accounting results lack spatial correlation and are difficult to support refined spatial decision-making, as detailed below: Existing accounting methods mostly present results in statistical reports or regional aggregate data, lacking a precise correspondence with geographic space. While specialized tools like InVEST can output spatial raster layers, they are limited to the physical quantity distribution of a single service type, unable to overlay and display the monetization value of multiple types of ecological products on the same spatial base map. Decision-makers struggle to intuitively identify the spatial differentiation patterns of value. This invention stores the value of various ecological products for each map patch in vector attribute form, forming a spatial database of "one record per map patch, one field per type." This supports overlay display and querying on GIS maps using thematic layers, categorized by product type, value level, administrative region, and other dimensions. This gives the accounting results spatial characteristics that allow for mapping, traceability, and comparison. Decision-makers can intuitively identify the spatial distribution patterns and high-value hotspots of ecological product value, providing precise geographic information support for spatial decisions such as identifying priority ecological protection areas, differentiating ecological compensation standards, and controlling land use. Meanwhile, since the results are generated based on the same set of spatial base maps and a unified accounting indicator system, the accounting results of different regions are naturally comparable, laying a data foundation for cross-regional comparative analysis and differentiated policy formulation. Attached Figure Description

[0023] To more clearly illustrate the specific embodiments of the present invention or the technical solutions in the prior art, the accompanying drawings used in the description of the specific embodiments or the prior art will be briefly introduced below. In all the drawings, similar elements or parts are generally identified by similar reference numerals. In the drawings, the elements or parts are not necessarily drawn to scale.

[0024] Figure 1 This is a schematic diagram of the overall architecture and data flow of the system of the present invention; Figure 2 Example diagram of indicator calculation method (taking material supply category - agricultural products as an example); Figure 3 This is a schematic diagram of the system modules of the present invention—an automated identification module for ecological product types; Figure 4 This is a schematic diagram of the module functions of the system of the present invention—the multi-scale adaptive aggregation and interaction module; Figure 5 This is a schematic diagram of the module functions of the system of the present invention—the map visualization and thematic display module.

[0025] Figure 6 This is a partial interface display of the system of the present invention. Figure 1 Multi-source data base; Figure 7 This is a partial interface display diagram of the system of the present invention—adaptive aggregation and interaction; Figure 8 This is a partial interface display of the system of the present invention—map visualization and thematic display. Detailed Implementation

[0026] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0027] It should be understood that, when used in this specification and the appended claims, the terms "comprising" and "including" indicate the presence of the described features, integrals, steps, operations, elements and / or components, but do not exclude the presence or addition of one or more other features, integrals, steps, operations, elements, components and / or collections thereof.

[0028] It should also be understood that the terminology used in this specification is for the purpose of describing particular embodiments only and is not intended to limit the invention. As used in this specification and the appended claims, the singular forms “a,” “an,” and “the” are intended to include the plural forms unless the context clearly indicates otherwise.

[0029] It should also be further understood that the term "and / or" as used in this specification and the appended claims refers to any combination of one or more of the associated listed items and all possible combinations, and includes such combinations.

[0030] I. Overall System Workflow In one embodiment of the present invention, an adaptive-scale ecological product value accounting system is provided, such as... Figure 1 As shown, the adaptive-scale ecological product value accounting system provided by this invention includes five functional modules, which are arranged in the following order according to data flow: multi-source data foundation module, automated ecological product type identification module, integrated ecological product value accounting module, multi-scale adaptive aggregation and interaction module, and map visualization and thematic display module. The overall workflow of the system is as follows: Step S1: The multi-source data base module acquires and integrates basic data. After determining the accounting area, this module systematically acquires various multi-source data required for accounting, including land survey data as the spatial carrier, and remote sensing image data, meteorological and hydrological monitoring data, socio-economic statistics, and point-of-interest data to support the calculation of various indicators. All data undergoes coordinate unification, format conversion, and classification for storage, forming a standardized multi-source data base. Simultaneously, this module uses the map features in the land survey data as basic spatial units to construct a map feature base map that fully covers the accounting area spatially, and establishes spatial topological affiliation relationships between each map feature and administrative divisions at the village, township, and county levels, forming a multi-level spatial affiliation index. For some adjustment service indicators that need to be calculated at the raster pixel scale, this module simultaneously prepares raster calculation units of corresponding resolutions.

[0031] Step S2: The automated ecological product type identification module identifies the calculable types. This module reads the land use code of each map patch and, based on a preset land use-ecological product type identification rule library, automatically identifies the calculable ecological product types and their combinations for each map patch, generating a map patch-accounting item mapping table. The identification rule library links the land use code system of land survey map patches with the accounting indicator system, enabling map patches of different land use types to automatically match different combinations of accounting indicators.

[0032] Step S3: The integrated accounting module for ecological product value performs integrated accounting on a per-plot basis. Based on the plot-account mapping table, this module uses plots as the smallest accounting unit, calling the corresponding accounting methods and parameter sets for each plot to complete the physical quantity and value accounting of various ecological products, generating detailed plot-level accounting results. For regulation service indicators that need to be calculated at the raster pixel scale, this module performs the accounting on the corresponding raster calculation unit and maps the results to plots. Each accounting model calls the relevant parameters and InVEST model stored in Step S1 to form various value calculation results.

[0033] Step S4: The multi-scale adaptive aggregation and interaction module aggregates multi-scale results on demand. Based on the spatial topological affiliation established in Step S1, this module automatically summarizes the map patch-level calculation results from bottom to top along the topological path of "map patch → village → township → county". The aggregation process satisfies the additivity condition, that is, the calculation result of any superior administrative division is equal to the sum of the calculation results of all its subordinate units. Users can independently select the target level of aggregation (such as administrative village, township, county (district, city)) and spatial scope (such as a single project, custom area), and obtain calculation results of any level and any area on demand, realizing the adaptive capability of "one-time fine calculation, multi-scale on-demand output".

[0034] Step S5: The map visualization and thematic display module performs spatial storage and visualization. This module writes the accounting results of various ecological products of each map patch into the map patch spatial database in the form of vector attributes (i.e., it links the accounting results of various ecological products to the map patch spatial database), forming a spatial accounting result database. On the GIS map, this module supports displaying the spatial distribution pattern of the value of various ecological products in thematic layers, providing functions such as layer management, spatial query, thematic color configuration, and interactive element information viewing, and supports generating statistical charts and exportable summary reports of accounting results at each level.

[0035] II. Specific Implementation of the Multi-Source Data Base Module The multi-source data foundation module serves as the system's data infrastructure, primarily managing the following five types of data: (1) Land survey data. Based on the results of the Third National Land Survey, vector data of map features in the accounting area were obtained. Each map feature is a vector surface feature with a closed boundary, and its attribute table records information such as map feature number, land use code, area, and ownership unit.

[0036] (2) Remote sensing data. This includes land cover raster data, vegetation types and coverage, etc., which are used for various physical quantity calculations.

[0037] (3) Meteorological and climate data. This includes meteorological data such as the multi-year average temperature, precipitation, and evapotranspiration of each meteorological station in the accounting area, which are used for the physical quantity accounting of regulation services such as water conservation and climate regulation.

[0038] (4) Socioeconomic data. This includes market price data for various agricultural, forestry, livestock and aquatic products, data on alternative costs such as water supply prices and carbon trading prices, and economic statistics required for cultural service accounting, such as benchmark land prices and total tourism revenue.

[0039] (5) Other data. This includes soil type and physicochemical properties, POI data, nature reserves, water resources bulletin data, water quality parameters, etc., used for physical quantity accounting of various regulation services.

[0040] This module establishes the spatial attribution relationship between each map patch and its corresponding administrative division at the village, township, county, and other levels. The specific methods are as follows: (a) Spatial overlay analysis is performed between the map patch vector layer and the administrative division vector layer; (b) For map patches that fall entirely within the scope of a certain administrative division, a one-to-one attribution relationship is directly established; (c) For map patches that cross the boundaries of administrative divisions, they are divided or allocated proportionally according to the area proportion of the map patch in different administrative divisions to ensure that the calculation result of each map patch can be uniquely or proportionally assigned to the corresponding administrative division; (d) A multi-level spatial attribution index table of "map patch number - village code - township code - county code" is formed and stored in the system database.

[0041] III. Specific Implementation of the Automated Identification Module for Ecological Product Types The core of this module is the construction of a rule base for distinguishing land use types and ecological product types. It links the land use type codes of map patches with the accounting indicator system of this system, and pre-determines the possible ecological product types and corresponding accounting methods for each land use type.

[0042] After the system starts the accounting task, the automated discrimination process of this module is as follows: (a) read the land use code of each map patch one by one; (b) query the list of calculable ecological product types corresponding to the land use code of the map patch in the discrimination rule base; (c) generate a list of accounting subjects for the map patch, including the name, category (material / regulation / culture) and value accounting method identifier of each accounting subject (some ecological products also include the physical quantity accounting method identifier); (d) summarize the discrimination results of all map patches to form a "map patch-accounting subject mapping table" for the accounting area, which serves as the input for the subsequent integrated accounting module. The entire discrimination process does not require manual intervention for each map patch, realizing automated batch discrimination across the entire area.

[0043] IV. Specific Implementation of the Integrated Accounting Module for Ecological Product Value The integrated accounting module obtains the types of ecological products that need to be accounted for for each land category based on the land category-ecological product type discrimination rule. Then, based on the land category attribute corresponding to the map patch, it performs the value accounting of all ecological products corresponding to the map patch one by one. Some products also include physical quantity accounting.

[0044] It should be known that an ecosystem can be identified by its land use name, and the type of ecological product can also be identified by a "one-to-many" or "many-to-one" relationship.

[0045] For example, the mapping rules are represented as follows: Table 1 Example of Mapping Rules It should be understood that the above is only an example excerpt to more clearly illustrate the mapping relationship of this solution, and is not intended to limit the technology of this invention. In actual applications, there are many more types of mapping relationships, which will not be elaborated here.

[0046] 4.1 Calculation of Material Supply Indicators The following uses agricultural products, forestry products, livestock products, and fishery products as examples to illustrate the method of calculating the value of material supply indicators at the map-level. Figure 2 Taking agricultural products as an example, this paper demonstrates an example of indicator calculation methods.

[0047] 1. Calculation of Material Supply Indicators This method targets the micro-scale value allocation of material supply indicators (including agricultural products, forestry products, livestock products, fishery products, and water resources supply). Its core idea is to conduct refined accounting based on the spatial differences in resource endowment and production capacity.

[0048] (1) Agricultural product value accounting The data is calculated using land parcels identified as farmland ecosystems. First, based on the land use type and planting attributes of the parcels, predefined rules are applied for secondary subdivision, mapping them to specific agricultural crop types. If the land type name is "paddy field", then the accounting category of the map patch is determined as "rice". If the land type is "dry land" and the planting attribute is "planting grain crops", then it is determined to be "corn". If the land type is "dry land" and the planting attribute is "planting non-grain crops", then it is determined to be "sugar crops"; If the land type is "irrigated land", "facility agricultural land", or "dry land" and the planting attribute is "grain and non-grain rotation", then it is included in the accounting scope of the "vegetable" category.

[0049] For each map patch with completed category mapping, obtain the key accounting parameters corresponding to that category, including: theoretical yield per unit area Y, regional production correction coefficient C, and product market price P. Simultaneously, assign an area adjustment coefficient A based on the specific land use attribute of the map patch (e.g., for a map patch mapped to "vegetables" and with the land use category "facility agricultural land," the coefficient A is set according to the proportion of the actual planting area within the facility; for general arable land, A is typically set to 1). The calculation formula is: V_plot=(S×A)×Y×C×P Where S is the area of ​​the patch.

[0050] (2) Valuation of forest products Map patches identified as forest ecosystems are used as the accounting objects. Sub-categories are then made based on the land use classification of the map patches: If the land use category is arbor forest, other forest, shrub forest, etc., the accounting type of the map patch shall be determined as "timber"; If the land parcel name belongs to bamboo forest land, the accounting type of the parcel shall be determined as "moso bamboo".

[0051] For timber patches, obtain the timber harvesting quantity per unit area Y_a and the timber price per unit area P_a; for bamboo patches, obtain the bamboo harvesting quantity per unit area Y_b and the bamboo price per unit area P_b. The calculation formulas are as follows: Timber plot: V_plot=S×Y_a×P_a Bamboo plot: V_plot=S×Y_b×P_b (3) Accounting for the value of livestock products The calculation focuses on land parcels identified as farmland ecosystems and designated as "facility agricultural land." The unit area output Y_l and product unit price P_l are obtained. A preset area adjustment coefficient K is used to recalculate the area of ​​the land parcels, reflecting the proportion of effective area actually used for livestock production within the facility agricultural land. The calculation formula is as follows: V_plot=(S×K)×Y_l×P_l (4) Value of fishery products: accounting for the value of fishery products The calculation object is a map patch identified as a wetland ecosystem and named "aquaculture pond" or "adjustable aquaculture pond". The output per unit area Y_f and the product price P_f are obtained. The calculation formula is: V_plot=S×Y_f×P_f.

[0052] (5) Water resource supply value accounting The calculation is performed using all map features as the accounting object. In practice, the AnnualWaterYield module of the InVEST model can be used to obtain the spatial distribution data of annual water yield. Then, the water yield is allocated according to the area proportion of domestic, agricultural, and industrial land within the map feature, multiplied by the corresponding water price per unit area, and summed to obtain the water resource supply value of that map feature. W_total=S×(WY / 1000) V_water=Σ[(W_total×A_i / A_total)×P_i] Where S is the area of ​​the map patch (m²), A_i is the area of ​​the i-th type of water use land, A_total is the total area of ​​the map patch, P_i is the corresponding water unit price (yuan / m³), and V_water is the value of water resource supply.

[0053] 2. Calculation of service-related indicators It should be noted that the specific calculation models for regulating service indicators mentioned below are commonly used and mature models in the field. Provided that their physical mechanisms and accounting objectives remain consistent, these models can be replaced by other validated models of equivalent effectiveness. Their innovation lies in the systematic application and solidification of these models into the value spatialization allocation process at the micro-level. The following lists the physical quantity accounting methods for major regulating service accounting items: (1) Water source conservation The water balance equation is used to calculate the water conservation capacity based on the precipitation, evapotranspiration, and runoff coefficient of the area where the map patch is located, combined with the conservation coefficient corresponding to the patch area and vegetation type. The calculation formula is: Q_conservation capacity = A × (P - ET - R), where A is the patch area, P is the precipitation, ET is the evapotranspiration, and R is the surface runoff. The value is calculated using the shadow project method (based on the reservoir construction unit's storage capacity cost).

[0054] (2) Soil conservation The modified general soil loss equation (RUSLE) is used to calculate potential soil erosion and actual soil erosion; the difference between the two is the soil retention. In practice, the Sediment Delivery Ratio module of the InVEST model is used for calculation. This module calculates the potential soil erosion (RKLS) and actual soil erosion (USLE) for each spatial unit based on rainfall erosivity factor (R factor), soil erodibility factor (K factor), slope length and gradient factor (LS factor), cover management factor (C factor), and soil and water conservation measures factor (P factor). The difference between the two is the soil retention for that unit. The value is calculated using the substitution cost method.

[0055] (3) Flood storage The appropriate water storage model is matched based on the ecosystem type of the land parcel: for forests and grasslands, the water storage capacity is calculated based on the depth of storm runoff; for paddy fields, the water storage capacity is calculated based on the height of the paddy field ridges and the water depth; for reservoirs, the storage capacity coefficient method is used; for lakes, the area power law formula is used; and for swamps, the soil water storage capacity and surface water retention capacity are calculated separately. The value is calculated using the shadow project method (based on the investment per unit storage capacity of the flood control project).

[0056] (4) Air purification Based on the preset "list of ecosystem types with air purification functions", the map patches are selected, and the purification capacity parameters per unit area of ​​sulfur dioxide, nitrogen oxides and dust for each ecosystem type are matched from the built-in parameter library. The annual purification volume of each pollutant is calculated, and then the economic value is calculated based on the treatment cost per unit price of each pollutant.

[0057] (5) Water purification Based on the preset "list of water and wetland ecosystem types with water purification functions", map patches are selected, and the purification capacity parameters per unit area of ​​chemical oxygen demand, total nitrogen, and total phosphorus for each ecosystem type are matched from the built-in parameter library. The annual purification amount of each pollutant is calculated, and then the economic value is calculated based on the treatment cost per unit price of each pollutant. (6) Carbon fixation Based on the preset "list of ecosystem types participating in carbon sequestration accounting", map patches are selected. The aboveground carbon density of vegetation, the underground carbon density of vegetation, the soil carbon density, and the carbon density of dead organic matter are matched from the built-in ecosystem carbon density parameter table to calculate the total carbon storage of the map patches. Then, the carbon dioxide equivalent is calculated according to the molecular weight conversion relationship between carbon and carbon dioxide (12g of carbon is equivalent to 44g of carbon dioxide). Finally, the carbon tax method (based on the carbon trading market price) is used to calculate the economic value of carbon sequestration.

[0058] (7) Oxygen release Based on the total carbon reserves of the map patches calculated in the carbon sequestration accounting, and according to the fixed stoichiometric ratio of carbon dioxide fixation and oxygen release in photosynthesis (32g of oxygen is released for every 12g of carbon fixed), the mass of oxygen released is calculated, and then the market value method (based on the cost of industrial oxygen production) is used to calculate the economic value of oxygen release. (8) Local climate regulation Map patches are categorized into vegetation cover and water area types for separate calculations. For vegetation cover patches, the equivalent air conditioning energy consumption is calculated based on vegetation transpiration; for water area patches, the equivalent cooling and humidification energy consumption is calculated based on water surface evaporation. The climate regulation value is obtained by multiplying the equivalent energy consumption by the electricity price.

[0059] (9) Provides negative oxygen ions Based on the preset "list of ecosystem types that can generate negative oxygen ions", map patches are selected, and the average concentration and characteristic height of negative oxygen ions of the ecosystem type are matched from the built-in parameter library. The annual production of negative oxygen ions within the spatial volume of the map patch is calculated, and then the economic value is calculated using the preset negative oxygen ion unit price.

[0060] (10) Species conservation Spatial overlay analysis was performed on the map features and the vector layer of the nature reserve boundary to screen out the map features located within the nature reserve. The conservation value of the species was calculated by multiplying the area of ​​the map feature by the annual conservation cost per unit area.

[0061] 3. Calculation of Cultural Service Indicators (1) Value assessment of tourism and health care Using all land parcels as the calculation object, spatial connection and overlay analysis are performed between the land parcel vector layer and the tourism and wellness point of interest database. For each land parcel, all tourism and wellness points of interest and their types within its spatial range are recorded. A pre-defined tourism and wellness point of interest type-value parameter mapping table is read, which defines the unit value P_type corresponding to each type of point of interest. For each land parcel containing points of interest, the number N_type of each type of point of interest within that land parcel is counted, and V_type = N_type × P_type is calculated. The total tourism and wellness value V_tourism of the land parcel is obtained by summing the values ​​of all types. Land parcels without any points of interest are recorded as having a value of 0.

[0062] (2) Value assessment of leisure and recreation Map patches are selected based on a pre-defined list of ecosystem types with significant recreational functions. In this embodiment, the list includes at least urban green space ecosystems (such as parks and plazas) and coastal wetland ecosystems (such as mangroves and coastal mudflats). The area S of each selected map patch is obtained, and the benchmark land price per unit area P_base for the spatial location of the map patch is obtained from a pre-defined land benchmark price spatial database. The calculation formula is: V_recreation = S × P_base. Map patches not included in the list have a recreational value of 0.

[0063] 4. A unified framework for value accounting The general formula for value accounting is: V_i = Q_i × P_i × γ_i, where Q_i is the physical quantity of the i-th type of ecological product, P_i is the corresponding unit price, and γ_i is the natural contribution coefficient. For regulating services, γ_i = 1, meaning the value is entirely attributable to the natural ecosystem; for material products, γ_i < 1, meaning the contribution from human inputs needs to be deducted.

[0064] External tool results integration. For certain accounting items requiring specialized model calculations (such as the sediment transport ratio module and annual water production module of the InVEST model), the integrated accounting module provides standardized data access interfaces. After external professional tools complete their calculations using their own methods, they output spatialized calculation results in raster or vector format. The system reads the external results through the data access interface, spatially overlays the raster data with the vector base map of the map features, and uses area-weighted statistics to summarize the raster values ​​to each map feature, converting them into map feature-level data. The converted external results are written into the system's accounting results database, managed uniformly with the system's built-in accounting results, and participate in subsequent aggregation and visualization processes.

[0065] Data Structure of Accounting Results. After the accounting is completed, the system generates detailed accounting results data at the map patch level. Each record includes: map patch number, map patch area, map patch land type code, administrative division code at all levels to which it belongs, physical quantity and value fields of various ecological products, and the total GEP value of the map patch.

[0066] V. Specific Implementation of Multi-Scale Adaptive Aggregation and Interaction Module This module is based on the spatial attribution index table established by the multi-source data foundation module, and summarizes it level by level according to the following path: The first level of aggregation (map patch → village) sums the accounting values ​​of all map patches belonging to the same administrative village according to each accounting item, generating summary data of the physical quantity and value of various ecological products of the village, as well as the total GEP value of the village.

[0067] The second level of aggregation (village → township) involves summing the total accounting values ​​of all villages under the same township by category to generate township-level accounting summary data and township-level GEP total value.

[0068] The third level of aggregation (township → county) involves summing the total accounting values ​​of all townships under the same county to generate the county-level total accounting data and the county-level GEP total. This process can continue upwards to the city and provincial levels.

[0069] The above aggregation process strictly satisfies the additivity condition: the total GEP of a county is equal to the sum of the GEPs of all its subordinate townships, the sum of the GEPs of all its villages, and the sum of the GEPs of all its map features. The accounting results at any level can be traced back to the accounting details of its constituent map features.

[0070] This module supports the following interactive functions: (1) Target level selection. Users can choose to output accounting results at the village, township, county, or other levels. The system will automatically perform aggregate calculations at the corresponding depth based on the level selected by the user.

[0071] (2) Roll up and drill down. When users view the accounting results at a certain level, they can perform roll up operations (such as rolling up from the village level to the township level to view more macro-level summary data) and drill down operations (such as drilling down from the township level to the village level, and then from the village level to the map patch level to view the accounting details), so as to realize the browsing of multi-level data.

[0072] (3) Itemized statistics. The aggregated results support itemized statistics and comparisons by ecological product type (primary and secondary categories) and by land type (arable land / forest land / grassland / water area, etc.).

[0073] VI. Specific Implementation of the Map Visualization and Thematic Display Module Spatial storage. This module allows the calculation results to be linked to the spatial database of map features. In the attribute table of the map feature vector layer, an attribute field is set for the physical quantity or value of each type of ecological product, forming a data structure of "one map feature, one record; one type, one field", so that the spatial location of each map feature is precisely associated with its calculation results.

[0074] Thematic layer display. This module supports generating various types of thematic layers on GIS maps: (1) Thematic map by ecological product type. When a user selects a certain type of ecological product (such as water conservation), the system uses the value field of the product to classify and color-render the map patches, and generates a thematic map of the spatial distribution of the product's value, which intuitively displays the spatial distribution of high-value areas and low-value areas.

[0075] (2) Thematic maps by value level. The system classifies the GEP total value of each patch / village / township / county (district, city) into levels (such as high value area, relatively high value area, medium value area, relatively low value area, low value area), and renders them with different colors to form a GEP value level distribution map.

[0076] (3) Thematic map by administrative region. When a user selects a specific administrative region (such as a township), the system only displays the calculation results of the map features within that region, achieving a focused display of spatial scope.

[0077] (4) Multi-layer overlay. Users can overlay multiple thematic layers simultaneously (such as overlaying the water conservation value layer and the administrative division layer) to conduct comprehensive spatial comparative analysis.

[0078] Statistical Display. This module also supports generating statistical charts and summary reports of the accounting results of ecological products at various levels and of various types, including: bar charts ranking the total GEP value of each administrative division, pie charts showing the proportion of value composition of various ecological products, comparison charts of GEP contribution of different types of products, and exportable summary reports containing detailed data on physical quantity and value of each level and item.

[0079] In one embodiment of the present invention, a specific implementation method for engineering the above-described system is provided, including: First, using the Python language and its geospatial processing ecosystem (including the GDAL / OGR spatial data engine, InVEST model, etc.), an independent spatial accounting module is built for each ecosystem product indicator. Each module encapsulates the complete data reading, spatial clipping, parameter matching, value calculation and result output.

[0080] Subsequently, a backend service layer was built based on the FastAPI asynchronous web framework, which uniformly registered and encapsulated the aforementioned series of accounting modules and their scheduling relationships. The accounting parameters were centrally managed and dynamically obtained through the indicator registry and configuration repository. At the same time, an independent executable accounting tool with a graphical user interface was built using the Tauri cross-platform desktop application framework and the Vue frontend framework.

[0081] In summary, this tool guides users step-by-step through configuring data paths and setting parameters, and automatically invokes the background accounting and scheduling engine to complete all calculations, ultimately outputting standardized spatial results and statistical reports. This implementation ensures one-click processing, standardization, and repeatability of the results. Key code snippets of the implementation are provided below: 1. Study Area Import and Geometric Check: This module is responsible for converting uploaded spatial files into system-recognizable study area objects. The first part ensures the integrity of the spatial file components, while the second part extracts the study area extent and outputs the geometric type, bounding box, and coordinate check results, ensuring that subsequent modules all work based on a unified study area input. Code examples are as follows: 2. Spatial Clipping and Computability Determination: This module first converts the study area into a unified geometric object, then places the study area and the target dataset in the same coordinate system to perform spatial filtering and intersection clipping. Subsequently, it extracts indicator parameters from the intersecting patches and determines which samples are truly computable according to indicator rules. Its output is not a simple "whether data exists" result, but rather a computable coverage result that is closer to business computation. Specifically: 3. Asynchronous Calculation Tasks and Indicator Scheduling: This module is responsible for converting a calculation request into a traceable background task and continuously writing back the task progress during execution. It also handles indicator scheduling, routing different indicators from the same request to either the map patch calculation or the region model calculation path, and then uniformly summarizing the output results. Details are as follows: 4. InVEST Model Integration (Current Instance: Soil Conservation): This module demonstrates general InVEST integration capabilities: it constructs a working directory and unified analysis grid around the study area, performs clipping, reprojection, and alignment on multiple input rasters, then assembles `invest_args` to call the InVEST model interface, and uniformly reads back the results after model execution. The current instance is soil conservation; this integration method can be reused for other InVEST indicators. Details are as follows: Those skilled in the art will recognize that the units of the various examples described in connection with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of both. To clearly illustrate the interchangeability of hardware and software, the components of the various examples have been generally described in terms of functionality in the foregoing description. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementations should not be considered beyond the scope of the invention.

[0082] In the embodiments provided by the present invention, it should be understood that the division of units is only a logical functional division. In actual implementation, there may be other division methods, such as multiple units can be combined into one unit, one unit can be split into multiple units, or some features can be ignored.

[0083] Furthermore, the functional units in the various embodiments of the present invention can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. The integrated unit can be implemented in hardware or as a software functional unit.

[0084] If the integrated unit is implemented as a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or all or part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, read-only memory (ROM), random access memory (RAM), portable hard drives, magnetic disks, or optical disks.

[0085] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some or all of the technical features therein. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present invention, and they should all be covered within the scope of the claims and specification of the present invention.

Claims

1. An adaptive-scale ecological product value accounting system, characterized in that, include: Multi-source data foundation module: For multi-source heterogeneous input, it collects, converts, unifies coordinates, and classifies various types of data required for accounting, forming a standardized multi-source data base that can be directly called by subsequent modules; the multi-source foundation data includes basic information data of map features with map features as the basic spatial unit and map feature accounting-related data; the multi-source data foundation module is also used to establish the spatial topological affiliation relationship between each map feature and administrative divisions at all levels, forming a multi-level spatial affiliation index, and the basic information data of the map features includes at least the land category name attribute; The automated ecological product type identification module is connected to the multi-source data base module. It obtains the land type name attribute of the map patch and, according to the preset land type-ecological product type identification rule, determines the ecosystem type corresponding to the map patch, forming a set of map patches carrying ecosystem type attributes. The land type-ecological product type identification rule includes at least the ecosystem type corresponding to the land type name attribute. Ecological product value integrated accounting module: connected to the ecological product type automatic identification module, used to perform value accounting of various ecological products on a map patch as the smallest accounting unit, and generate map patch-level accounting result data; Multi-scale adaptive aggregation and interaction module: connected to the integrated accounting module for ecological product value and the multi-source data base module respectively, used to obtain the spatial topology attribution relationship and the map-level accounting result data, and classify and summarize the map-level accounting result data from bottom to top along the preset topological path or for a custom range according to different needs, and generate the final accounting result; Map visualization and thematic display module: connected to the multi-scale adaptive aggregation and interaction module and the integrated accounting module for ecological product value, used to spatially store the accounting results of each map patch in the form of vector attributes, visualize them on the map in the form of thematic layers, and visualize the final accounting results in the form of charts in multiple dimensions.

2. The adaptive-scale ecological product value accounting system according to claim 1, characterized in that, The basic information data of the map patch also includes spatial boundaries, area, land use code, and ownership information. The spatial boundary is the boundary of the corresponding range of the map patch in the geographical space. The area is the area within the spatial boundary. The land use code is the code corresponding to the land use name attribute. The ownership information is the ownership relationship between the range defined by the spatial boundary of the map patch and the administrative regions at all levels.

3. The system according to claim 1, characterized in that, The automated identification module for ecological product types adopts a modular indicator management architecture. The accounting items in the identification rule base are organized in a configurable manner, supporting the addition, modification, and deactivation of accounting items without modifying the core logic of the system.

4. The adaptive-scale ecological product value accounting system according to claim 1, characterized in that, The integrated accounting module for the value of ecological products includes a sub-module for agricultural product value accounting. This sub-module uses a set of map patches identified as farmland ecosystems as the basic accounting object, and calculates the agricultural product value V_plot for each map patch, including: Agricultural product crop type mapping unit: Based on the land use name and planting attribute of the map patch, it performs secondary subdivision by executing predefined rules and maps it to specific agricultural product crop types; The value parameter quantification calculation unit at the map patch level: For each map patch that has completed crop type mapping, obtain the key accounting parameters corresponding to the crop type, including at least: theoretical yield per unit area Y, regional production correction coefficient C, and product market price P. Assign an area adjustment coefficient A according to the specific agricultural crop type of the map patch, and calculate the agricultural product value V_plot carried by the map patch. The calculation formula is: V_plot=(S×A)×Y×C×P, where S is the area of ​​the map patch.

5. The adaptive-scale ecological product value accounting system according to claim 1, characterized in that, The integrated ecological product value accounting module includes a forest product value accounting submodule. This submodule uses a set of map patches identified as forest ecosystems as the basic accounting object, and calculates the forest product value V_plot for each map patch, including: Forest product type mapping unit: Based on the land use name of the map patch, perform secondary subdivision by executing predefined rules and map it to a specific forest product type, wherein the forest product type includes at least timber and bamboo; Map-level value parameter quantification calculation unit: For each map patch with completed type mapping, obtain the key accounting parameters corresponding to that type. For timber maps, the parameters include: timber harvesting volume per unit area Y_a and timber unit price P_a; for bamboo maps, the parameters include: bamboo harvesting volume per unit area Y_b and bamboo unit price P_b. The calculation formula is as follows: For wood plots: V_plot=S×Y_a×P_a For the bamboo patch: V_plot=S×Y_b×P_b Where S is the area of ​​the patch.

6. The adaptive-scale ecological product value accounting system according to claim 1, characterized in that, The integrated accounting module for ecological product value includes a sub-module for livestock product value accounting. This sub-module uses a set of plots identified as farmland ecosystems as its basis, and only includes plots with the land use designation "facility agricultural land" in the scope of livestock product value accounting. For each selected plot, it obtains the output per unit area Y_l and the unit price P_l of the product. It uses a preset area adjustment coefficient K to convert the area of ​​the plot to reflect the effective area ratio of facility agricultural land actually used for livestock production. The calculation formula is: V_plot=(S×K)×Y_l×P_l, where S is the area of ​​the plot.

7. The adaptive-scale ecological product value accounting system according to claim 1, characterized in that, The integrated accounting module for ecological product value includes a sub-module for fishery product value accounting. This sub-module uses a set of plots identified as wetland ecosystems as its basis, and only includes plots with land use names of aquaculture ponds or adjustable aquaculture ponds in the scope of fishery product value accounting. For each selected plot, the output per unit area Y_f and the product unit price P_f are obtained. The calculation formula is: V_plot=S×Y_f×P_f, where S is the area of ​​the plot.

8. The adaptive-scale ecological product value accounting system according to claim 1, characterized in that, The integrated accounting module for the value of ecological products includes a sub-module for tourism and health care value accounting. This sub-module calculates the value of tourism and health care using all micro-spatial calculation units as the accounting objects, including: Spatial data association unit: Spatial connection and overlay analysis is performed between the land parcel vector layer and the tourism and health care point of interest database to generate an attribute for each parcel to record all tourism and health care points of interest contained within its spatial range and their types, forming a parcel-point of interest association dataset; Interest Point Value Parameter Matching Unit: Reads a preset tourism and health care interest point type-value parameter mapping table, which defines the unit value P_type corresponding to each type of tourism and health care interest point; Tourism and Wellness Value Calculation Unit: For each map patch containing tourism and wellness points of interest, based on the map patch-point of interest association data, count the number N_type of each type of tourism and wellness point of interest within the map patch, and then calculate the value contributed by each type of tourism and wellness point of interest within the map patch according to the matched unit value. V_type = N_type × P_type Finally, the values ​​of all tourism and wellness interest types within the map patch are summed to obtain the total tourism and wellness value V_tourism of the map patch. For map patches that do not contain any tourism and wellness interest types, their value is recorded as 0.

9. The adaptive-scale ecological product value accounting system according to claim 1, characterized in that, The integrated accounting module for the value of ecological products includes a sub-module for calculating the value of leisure and recreation. This sub-module is used to calculate the value of leisure and recreation, and the calculation process is as follows: Based on the preset list of recreational ecosystem types, map patches are selected, and the area S of each selected map patch is obtained. At the same time, the benchmark land price per unit area P_base of the spatial location of the map patch is obtained from the land benchmark price spatial database. The recreational value V_recreation it carries is calculated using the following formula: V_recreation=(S×P_base). Map features not included in the aforementioned list have a recreational value of 0.

10. The system according to claim 1, characterized in that, The multi-scale adaptive aggregation and interaction module sums up step by step along the topological path of "map patch → village → township → county". The aggregation process satisfies the additivity condition: the calculation result of any superior administrative division is equal to the sum of the calculation results of all its subordinate units. The accounting results at each level can be traced back to the detailed accounting data of each map patch that makes up that level; moreover, the multi-scale adaptive aggregation and interaction module allows users to independently select the target level and spatial range for aggregation, and obtain accounting results of any level and any region as needed; and supports the statistical and comparative analysis of the aggregation results by ecological product type or land category.

Images