Plateau traffic corridor ecological zoning method, device, equipment and mapping system

Abstract

The invention discloses a plateau traffic corridor ecological zoning method, device and equipment and a mapping system. According to the scheme, meteorological station data, digital elevation model data and land utilization type data in a research area are obtained based on a vector range of the research area of the ecological division of the plateau traffic corridor. And obtaining the surface undulation degree and the ecosystem type of the research area based on the data, providing a meteorological station missing air temperature data interpolation method, and interpolating the missing air temperature. And determining and making three zoning indexes of a plateau climate zone type, a plateau landform type and a road area ecological type of the research area by combining a preset mapping relation based on the temperature data, the surface undulation degree and the ecological system type after interpolation, and performing spatial coupling analysis on the three indexes to obtain a coupling result. Performing multi-index-plaque clustering analysis based on the coupling result to obtain a clustering result; and performing ecological zoning on the research area based on a clustering result. And refined division of the plateau traffic corridor ecology is realized.

Description

Technical Field

[0001] The present application relates to the technical field of ecological zoning of transportation corridors, and in particular to a method, device, equipment and mapping system for ecological zoning of plateau transportation corridors. Background Art

[0002] Ecological zoning is one of the important contents of ecosystem services, ecological risks and ecological sensitivity assessment, and is also a reference for ecological restoration construction. At present, ecological zoning mainly serves urban space development, rural green construction, etc. Most of the division methods are based on land use types, combined with local topographic features and ecological development plans, to formulate reasonable ecological zoning.

[0003] In the high-altitude and cold regions of the plateau, natural and human factors are intertwined and complex, the ecological environment is harsh and changeable, and the construction of plateau road projects inevitably affects the surrounding ecological environment, causing certain difficulties in the ecological zoning of plateau transportation corridors. Therefore, the study of ecological zoning of transportation corridors in high-altitude and cold regions is difficult, and a clear systematic evaluation index system has not yet been formed; and the ecological differences of transportation corridors in high-altitude and cold regions of the plateau are large, and the traditional ecological zoning method can no longer meet the needs of refined transportation construction, operation, management and ecological protection. It is urgent to propose a new ecological zoning method for plateau transportation corridors. Summary of the invention

[0004] Based on the above problems, this application provides a method, device, equipment and mapping system for the ecological zoning of plateau transportation corridors, with the aim of proposing a systematic evaluation index system suitable for the ecological zoning of plateau transportation corridors, thereby realizing the refined zoning of the ecology of plateau transportation corridors.

[0005] The embodiments of the present application disclose the following technical solutions:

[0006] The first aspect of the present application provides a method for ecological zoning of plateau transportation corridors, the method comprising:

[0007] Acquire meteorological station data, digital elevation model data, and land use type data within the plateau transportation corridor ecological zoning study area based on the vector range of the plateau transportation corridor ecological zoning study area;

[0008] Interpolating the temperature data of the missing period in the meteorological station data based on the existing temperature data in the meteorological station data to obtain interpolated temperature data;

[0009] Based on the digital elevation model data, the surface undulation within the plateau transportation corridor ecological zoning study area is obtained;

[0010] According to the preset mapping relationship between temperature and plateau climate zone type, the plateau climate zone type of the plateau transportation corridor ecological zoning study area is determined based on the interpolated temperature data; according to the preset mapping relationship between surface undulation and plateau landform type, the plateau landform type of the plateau transportation corridor ecological zoning study area is determined based on the surface undulation in the plateau transportation corridor ecological zoning study area; according to the preset mapping relationship between land use type and road ecological type, the road ecological type of the plateau transportation corridor ecological zoning study area is determined based on the land use type data in the plateau transportation corridor ecological zoning study area;

[0011] The plateau climate zone type, plateau landform type and road ecological type of the plateau transportation corridor ecological zoning study area are used as three indicators to perform multi-indicator spatial coupling to obtain a coupling result; the coupling result includes multiple patches, each patch corresponds to an ecological zoning type, and the ecological zoning type is a type expressed by a combination of the three indicators;

[0012] Performing multi-index-plaque clustering based on the coupling result to obtain a clustering result;

[0013] Based on the clustering results, ecological zoning is carried out for the plateau transportation corridor ecological zoning study area.

[0014] In an optional implementation, the meteorological station data includes a station number, a station name, a station accuracy, a station latitude, and a temperature value observed at the station;

[0015] The interpolating the temperature data of the missing period in the meteorological station data based on the existing temperature data in the meteorological station data to obtain the interpolated temperature data includes:

[0016] The temperature data of the missing period in the meteorological station data is averaged and assigned by the existing temperature data of the first three and last three periods of the same length as the missing period, and the calculation formula is:

[0017]

[0018] Among them, TV i represents the average temperature value of the i-th period; i-3, i-2, i-1, i+1, i+2, i+3 represent the i-3rd period, i-2nd period, i-1st period, i+1st period, i+2nd period, i+3rd period respectively; TV i-3 , TV i-2 , TV i-1 , TV i+1 , TV i+2 , TV i+3 They represent the average temperature values ​​of the corresponding time periods.

[0019] In an optional implementation, the plateau climate zone type, plateau landform type and road ecological type of the plateau transportation corridor ecological zoning study area are used as three indicators to perform multi-indicator spatial coupling to obtain coupling results, including:

[0020] The plateau climate zone type is used as the first indicator, the plateau landform type is used as the second indicator, and the road ecological type is used as the third indicator;

[0021] Obtain the number of plaques a, b and c respectively included in the first index, the second index and the third index; a, b and c are all positive integers;

[0022] The first index, the second index and the third index are spatially superimposed to obtain a total index, wherein the total index includes a×b×c patches and is a coupling result.

[0023] In an optional implementation, the performing multi-index-plaque clustering based on the coupling result includes:

[0024] Based on the coupling result, the centroids of patches of N types of ecological zoning are selected as cluster centers respectively; N is an integer greater than or equal to 1 and less than the total number of ecological zoning types in the coupling result, and the ecological zoning type of the cluster center is the top N ecological zoning types in terms of the number of patches and the proportion of patch area in the coupling result;

[0025] The ecological zoning types of the remaining patches are determined based on the distances between the centroids of the remaining patches and the centers of each cluster.

[0026] In an optional implementation, determining the ecological zoning type of the remaining patches according to the distances between the centroids of the remaining patches and the centers of each cluster includes:

[0027] If the centroid of a patch among the remaining patches has the shortest straight-line distance to a cluster center, the ecological zoning type of the patch is the ecological zoning type of the cluster center.

[0028] In an optional implementation, the mapping relationship between the preset temperature and the plateau climate zone type is to obtain five plateau climate zone types by dividing them into -4°C, -2°C, 4°C, and 6°C as critical values, and the five plateau climate zone types include high cold, cold, cold temperature, cool temperature, and warm;

[0029] The preset mapping relationship between the surface relief and the plateau landform type is that 20m, 30m, 200m, 500m, 1000m and 2500m are used as the critical values ​​to obtain six plateau landform types, and the six plateau landform types include plains, terraces, hills, medium mountains, high mountains and extremely high mountains;

[0030] The mapping relationship between the preset land use type and the road ecological type is to reclassify the land use type to obtain eight road ecological types, and the eight road ecological types include urban ecological type, cultivated ecological type, mixed agricultural and pastoral ecological type, mixed agricultural and forestry ecological type, grassland ecological type, woodland ecological type, water wetland ecological type and arid desert ecological type.

[0031] In an optional implementation, the reclassification of land use types includes:

[0032] Reclassify cities, towns, and industrial and mining areas into urban ecological types;

[0033] Reclassify paddy fields, dry fields, and oases into cultivated ecological types;

[0034] reclassification of agro-pastoral and pastoral farming into mixed agro-pastoral ecological types;

[0035] reclassification of agroforestry and forestry as mixed agroforestry ecological types;

[0036] Typical grass, patchy grass, and sparse grass were reclassified into pasture ecotypes;

[0037] reclassified dense forest, semi-dense forest, open forest and shrub forest into woodland ecotypes;

[0038] Reclassify water bodies, swamps, beaches, ice and snow, and frozen deserts into water wetland ecotypes;

[0039] Deserts, deserts and bare rock gobi were reclassified as arid desert ecological types.

[0040] The second aspect of the present application provides a plateau transportation corridor ecological zoning device, which includes:

[0041] A data acquisition module, used to acquire meteorological station data, digital elevation model data and land use type data within the plateau transportation corridor ecological zoning study area based on the vector range of the plateau transportation corridor ecological zoning study area;

[0042] An interpolation module, used to interpolate the temperature data of the missing period in the meteorological station data based on the existing temperature data in the meteorological station data to obtain the interpolated temperature data;

[0043] A surface relief determination module, used for obtaining the surface relief in the plateau transportation corridor ecological zoning study area based on the digital elevation model data;

[0044] A type determination module, for determining the plateau climate zone type of the plateau transportation corridor ecological zoning study area based on the interpolated temperature data according to a preset mapping relationship between temperature and plateau climate zone type; determining the plateau landform type of the plateau transportation corridor ecological zoning study area based on the surface undulation in the plateau transportation corridor ecological zoning study area according to a preset mapping relationship between surface undulation and plateau landform type; determining the road area ecological type of the plateau transportation corridor ecological zoning study area based on the land use type data in the plateau transportation corridor ecological zoning study area according to a preset mapping relationship between land use type and road area ecological type;

[0045] A coupling module is used to perform multi-indicator spatial coupling on the plateau climate zone type, plateau landform type and road ecological type of the plateau transportation corridor ecological zoning study area as three indicators to obtain a coupling result; the coupling result includes multiple patches, each patch corresponds to an ecological zoning type, and the ecological zoning type is a type expressed by a combination of the three indicators;

[0046] A clustering module, used for performing multi-index-plaque clustering based on the coupling result to obtain a clustering result;

[0047] The ecological zoning module is used to carry out ecological zoning for the plateau transportation corridor ecological zoning study area according to the clustering results.

[0048] The third aspect of the present application provides a plateau transportation corridor ecological zoning device, the device comprising: a processor and a memory:

[0049] The memory is used to store program code and transmit the program code to the processor;

[0050] The processor is used to execute the steps of the plateau transportation corridor ecological zoning method introduced in any implementation method of the first aspect according to the instructions in the program code.

[0051] The fourth aspect of the present application provides a plateau transportation corridor ecological zoning mapping system, the mapping system comprising: a plateau transportation corridor ecological zoning visualization mapping information library, a plateau transportation corridor ecological zoning thematic map mapping component, and the plateau transportation corridor ecological zoning device introduced in the third aspect;

[0052] The plateau transportation corridor ecological zoning visualization mapping information base and the plateau transportation corridor ecological zoning thematic map mapping component are respectively connected to the processor in the plateau transportation corridor ecological zoning device;

[0053] The plateau transportation corridor ecological zoning visualization mapping information database includes a symbol library, a color library and an annotation library;

[0054] The plateau transportation corridor ecological zoning thematic map mapping component is used to receive the plateau transportation corridor ecological zoning results of the processor; based on the plateau transportation corridor ecological zoning results, the plateau transportation corridor ecological zoning visualization mapping information library is used to produce a plateau transportation corridor ecological zoning thematic map; the plateau transportation corridor ecological zoning thematic map includes a scale, a compass and a self-made legend.

[0055] Compared with the prior art, this application has the following beneficial effects:

[0056] In the technical solution of the present application, firstly, the meteorological station data, digital elevation model data and land use type data within the plateau transportation corridor ecological zoning study area are obtained based on the vector range of the plateau transportation corridor ecological zoning study area.

[0057] Secondly, the temperature data of the missing period in the meteorological station data is interpolated based on the existing temperature data in the meteorological station data to obtain the interpolated temperature data. The surface undulation in the plateau transportation corridor ecological zoning study area is obtained based on the digital elevation model data.

[0058] Again, according to the preset mapping relationship, the plateau climate zone type, plateau landform type and road ecological type of the plateau transportation corridor ecological zoning study area are determined based on the interpolated temperature data, surface undulation and land use type data.

[0059] Then, the plateau climate zone type, plateau landform type and road ecological type in the plateau transportation corridor ecological zoning study area were used as three indicators for multi-indicator spatial coupling to obtain the coupling results; the coupling results include multiple patches, each patch corresponds to an ecological zoning type, and the ecological zoning type is the type expressed by the combination of three indicators.

[0060] Then, multi-index-plaque clustering is performed based on the coupling results to obtain the clustering results.

[0061] Finally, based on the clustering results, ecological zoning was carried out for the plateau transportation corridor ecological zoning study area.

[0062] It can be seen that in the technical solution of this application, by comprehensively considering multi-dimensional indicator data such as meteorology, topography and land use type, a refined division of the ecological environment of the region is achieved. Based on the pre-set clear preset relationship as a systematic evaluation indicator system, and by interpolating the actual temperature data, it is possible to more accurately judge the plateau climate zone type, landform type and road ecological type of the study area. Through multi-indicator spatial coupling and cluster analysis, a scientifically based ecological zoning result is obtained to meet the needs of refined transportation construction, operation, management and ecological protection. BRIEF DESCRIPTION OF THE DRAWINGS

[0063] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative labor.

[0064] Figure 1 A flow chart of a plateau transportation corridor ecological zoning method provided in an embodiment of the present application;

[0065] Figure 2 A schematic diagram of multi-index spatial coupling provided in an embodiment of the present application;

[0066] Figure 3 A schematic diagram of multi-indicator-plaque clustering provided in an embodiment of the present application;

[0067] Figure 4 A schematic diagram of the structure of a plateau transportation corridor ecological zoning device provided in an embodiment of the present application;

[0068] Figure 5 A schematic diagram of a thematic map provided in an embodiment of the present application;

[0069] Figure 6 A plateau transportation corridor ecological zoning thematic map provided in the embodiment of the present application;

[0070] Figure 7 Schematic diagram of some ecological zoning types provided for the embodiments of this application. DETAILED DESCRIPTION

[0071] As described above, due to the great ecological differences in transportation corridors in plateau and cold regions, the current traditional ecological zoning methods can no longer meet the needs of refined transportation construction, operation, management and ecological protection.

[0072] After research, the inventor proposed a method, device, equipment and mapping system for ecological zoning of plateau transportation corridors.

[0073] Firstly, based on the vector range of the plateau transportation corridor ecological zoning study area, the meteorological station data, digital elevation model data and land use type data in the plateau transportation corridor ecological zoning study area are obtained. Secondly, based on the existing temperature data in the meteorological station data, the temperature data of the missing period in the meteorological station data are interpolated to obtain the interpolated temperature data. The surface undulation in the plateau transportation corridor ecological zoning study area is obtained based on the digital elevation model data. Again, according to the preset mapping relationship, the plateau climate zone type, plateau landform type and road area ecological type of the plateau transportation corridor ecological zoning study area are determined based on the interpolated temperature data, surface undulation and land use type data. Then, the plateau climate zone type, plateau landform type and road area ecological type of the plateau transportation corridor ecological zoning study area are used as three indicators for multi-indicator spatial coupling to obtain the coupling result; the coupling result includes multiple patches, each patch corresponds to an ecological zoning type, and the ecological zoning type is the type expressed by the combination of three indicators. Then, multi-indicator-patch clustering is performed based on the coupling result to obtain the clustering result. Finally, based on the clustering result, ecological zoning is performed for the plateau transportation corridor ecological zoning study area.

[0074] In order to enable those skilled in the art to better understand the solution of the present application, the technical solution in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0075] See also Figure 1 , which is a flow chart of a plateau transportation corridor ecological zoning method provided in an embodiment of the present application. Figure 1 As shown, the method comprises the following steps:

[0076] S101. Obtain meteorological station data, digital elevation model data, and land use type data within the plateau transportation corridor ecological zoning study area based on the vector range of the plateau transportation corridor ecological zoning study area.

[0077] In the embodiment of the present application, the vector range of the plateau transportation corridor ecological zoning study area refers to the specific geographical boundaries delineated in the plateau area for the purpose of studying the impact of transportation corridors (such as railways, highways, etc.) on the ecological environment and conducting ecological zoning, and these boundaries are expressed in the form of vector data. In an example implementation, the vector range of the plateau transportation corridor ecological zoning study area is in shapefile format.

[0078] In the embodiment of the present application, a weather station is an observation point set on the surface of the earth or in the air to collect weather data. Weather station data refers to relevant information from a weather observation station. In the embodiment of the present application, the weather station data is in Excel table format.

[0079] Digital Elevation Model (DEM) data refers to a digital model that describes the elevation information of the earth's surface. It usually exists in the form of raster data or vector data and is used for terrain analysis and geomorphological research.

[0080] Land use type data refers to data that describes different land use types on the earth's surface (such as towns, townships, agriculture and forestry, etc.), which are used for land resource management and urban planning.

[0081] In an example implementation, according to the vector range of the study area, meteorological station data can be obtained from institutions such as the National Meteorological Science Data Center or the National Meteorological Information Center. DEM data can be obtained from institutions such as the Geospatial Data Cloud and the Resource and Environmental Science Data Center of the Chinese Academy of Sciences. Land use type data can be obtained from institutions such as the National Geographic Information Resource Directory Service System and the Geospatial Data Cloud.

[0082] S102. Interpolating the temperature data of the missing time period in the meteorological station data based on the existing temperature data in the meteorological station data to obtain interpolated temperature data.

[0083] In the embodiment of the present application, the reason why there is missing temperature data in the weather station data may be equipment failure, data transmission problem, observation condition limitation, human factors and other reasons. In order to make the temperature data more accurate, the temperature data of the missing period in the weather station data is interpolated based on the existing temperature data in the weather station data to obtain the interpolated temperature data.

[0084] S103. Obtain the surface undulation within the plateau transportation corridor ecological zoning study area based on digital elevation model data.

[0085] In the implementation of this application, terrain analysis is performed based on digital elevation model data to obtain the relief degree of land surface (RDLS). The relief degree of land surface is the difference between the maximum and minimum elevation values ​​of all pixels in the neighborhood. The calculation formula is as follows:

[0086] RDLS=MaxDEM-MinDEM

[0087] in, RDLS is the surface relief, MaxDEM , MinDEMThey are the maximum and minimum elevation values ​​of all pixels in the neighborhood, respectively.

[0088] S104. According to the preset mapping relationship between temperature and plateau climate zone type, the plateau climate zone type of the plateau transportation corridor ecological zoning study area is determined based on the interpolated temperature data; according to the preset mapping relationship between surface undulation and plateau landform type, the plateau landform type of the plateau transportation corridor ecological zoning study area is determined based on the surface undulation in the plateau transportation corridor ecological zoning study area; according to the preset mapping relationship between land use type and road ecological type, the road ecological type of the plateau transportation corridor ecological zoning study area is determined based on the land use type data in the plateau transportation corridor ecological zoning study area.

[0089] In the embodiment of the present application, the plateau climate zone type refers to a unique climate zone formed in the plateau area due to factors such as altitude, geographical location, and atmospheric circulation.

[0090] Plateau landform type refers to the unique landform formed in the plateau area due to factors such as geological structure, climate change and erosion.

[0091] Road ecological type refers to the ecosystem type and ecological environment characteristics on both sides of the highway and within a certain range along the road.

[0092] In the embodiment of the present application, determining the plateau climate zone type, plateau landform type and road ecological type of the plateau transportation corridor ecological zoning study area can provide an important reference for subsequent ecological zoning.

[0093] S105. The plateau climate zone type, plateau landform type and road ecological type in the plateau transportation corridor ecological zoning study area are used as three indicators for multi-indicator spatial coupling to obtain the coupling results.

[0094] In an embodiment of the present application, the coupling result includes multiple patches, each patch corresponds to an ecological zoning type, and the ecological zoning type is a type expressed by a combination of three indicators: plateau climate zone type, plateau landform type and road ecological type of the plateau transportation corridor ecological zoning study area.

[0095] S106. Perform multi-index-plaque clustering based on the coupling result to obtain a clustering result.

[0096] Since the coupling results contain a large number of patches, each patch corresponds to an ecological zoning type. These patches may be dispersed in space, making the data very complicated. Therefore, through multi-indicator-patch clustering, these patches are classified according to certain rules, thereby reducing the complexity of the data and facilitating more reasonable ecological zoning in the future.

[0097] S107. Based on the clustering results, ecological zoning is carried out for the plateau transportation corridor ecological zoning study area.

[0098] In an embodiment of the present application, ecological zoning is performed on the plateau transportation corridor ecological zoning study area according to the ecological zoning types corresponding to the patches in the clustering results after clustering.

[0099] In the embodiment of the present application, the temperature data of the missing period in the meteorological station data is interpolated based on the existing temperature data in the meteorological station data, which ensures the quality of the temperature data and makes the subsequent determination of the plateau climate zone type more accurate. In view of the specificity of the ecological characteristics of the plateau transportation corridor, the multi-index spatial coupling analysis and multi-index-patch clustering analysis methods are integrated to make the ecological zoning of the plateau transportation corridor more efficient, while overcoming the problems of the traditional zoning method being too general and insufficient in adaptability.

[0100] Optionally, the meteorological station data includes station number, station name, station accuracy, station latitude, and temperature value observed at the station;

[0101] Optionally, the temperature data of the missing period in the meteorological station data in the above step S102 is averaged and assigned by the existing temperature data of the first three and last three periods of the same length as the missing period, and the calculation formula is:

[0102]

[0103] Among them, TV i represents the average temperature value of the i-th period; i-3, i-2, i-1, i+1, i+2, i+3 represent the i-3rd period, i-2nd period, i-1st period, i+1st period, i+2nd period, i+3rd period respectively; TV i-3 , TV i-2 , TV i-1 , TV i+1 , TV i+2 , TV i+3 They represent the average temperature values ​​of the corresponding time periods.

[0104] In an embodiment of the present application, when the temperature data of a certain day is missing, the temperature data of the missing day is averaged and assigned by the temperature data of the three days before and after; when the temperature data of a certain week is missing, the temperature data of the missing week is averaged and assigned by the temperature data of the three weeks before and after; for the temperature data of other missing periods of different lengths, the interpolation method is similar. For example, when the temperature data of the i-th day is missing, the average temperature value data of the i-3rd day, the i-2nd day, the i-1st day, the i+1st day, the i+2nd day and the i+3th day are already in the meteorological station data, which are 23℃, 22℃, 24℃, 23℃, 25℃ and 24℃ respectively. According to the calculation formula, the average temperature value of the i-th day is (23+22+24+23+25+24) / 6=23.5℃, that is, the missing temperature data of the i-th day is interpolated to 23.5℃.

[0105] Optionally, the above step S105 includes:

[0106] The plateau climate zone type is taken as the first indicator, the plateau landform type is taken as the second indicator, and the road area ecological type is taken as the third indicator;

[0107] Obtain the number of plaques a, b and c respectively included in the first index, the second index and the third index; a, b and c are all positive integers;

[0108] The first, second and third indicators are spatially superimposed to obtain the total indicator, which contains a×b×c patches and is the coupling result.

[0109] In the embodiment of the present application, multi-indicator spatial coupling is performed based on the aforementioned determined plateau climate zone type, plateau landform type and road ecological type as the first indicator, the second indicator and the third indicator respectively.

[0110] See also Figure 2 , Figure 2 A schematic diagram of multi-index spatial coupling provided in an embodiment of the present application. Figure 2 For example, the first index (8 patches), the second index (4 patches) and the third index (6 patches) of the same range are spatially coupled to obtain the total index, which has a total of 8×4×6=192 patches. Among them, 1-1 to 1-8 represent the 8 patches in the first index, 2-1 to 2-4 represent the 4 patches in the second index, and 3-1 to 3-6 represent the 6 patches in the third index. For the total index, Figure 2 Only 1 to 4 are marked as examples of patches in the total index. In practice, Figure 2 The total index shown has a total of 192 patches.

[0111] In an example implementation, the steps of coupling multiple index spaces are as follows:

[0112] Step 1: Define the first indicator (layer 1), the second indicator (layer 2), the third indicator (layer 3) and the total indicator (new layer), which respectively indicate the plateau climate zone layer, the plateau landform type layer, the road ecological type layer and the preliminary layer of plateau transportation corridor zoning.

[0113] Step 2: Assume that the first indicator has 8 patches, the second indicator has 4 patches, and the third indicator has 6 patches. Perform spatial superposition of the indicators in pairs (such as spatial superposition of the first indicator and the second indicator), that is, multiply them in pairs (such as 8×4).

[0114] Step 3: Repeat step 2 in sequence to finally obtain the superposition results of multiple indicators, which is the coupling result.

[0115] In the embodiment of the present application, due to the presence of fragmented patches in the above-mentioned multi-index spatial coupling results, based on the ecological and environmental characteristics of the given plateau transportation corridor ecological zoning study area, combined with the road ecological type, plateau climate zone type and plateau landform type, ecological zoning is performed according to the zoning principle of "consistency within the region, difference between intervals, taking the primary and discarding the secondary, and merging similar ones". This zoning principle is specifically reflected through the subsequent multi-index-patch clustering.

[0116] Optionally, the above step S106 includes:

[0117] Based on the coupling results, the centroids of patches of N ecological zoning types are selected as cluster centers; N is an integer greater than or equal to 1 and less than the total number of ecological zoning types in the coupling results. The ecological zoning type of the cluster center is the ecological zoning type with the top N patches in terms of the number of patches and the proportion of patch area in the coupling results.

[0118] The ecological zoning types of the remaining patches are determined based on the distances between the centroids of the remaining patches and the centers of each cluster.

[0119] In the embodiment of the present application, multi-index-plaque clustering is performed on the coupling results, wherein the value of N is not specifically limited. Figure 3 A schematic diagram of multi-index-plaque clustering provided in an embodiment of the present application, where N is 4, Figure 3 The multi-index-patch clustering shown selects patches of ecological zoning types A, B, C, and D as cluster centers, and aggregates different patches into one cluster center due to their similarities.

[0120] In an exemplary implementation, when N is 9, the multi-index-plaque clustering steps are as follows:

[0121] Step 1: Based on the superposition results of multiple indicators, select several cluster centers (centroids), for example, select the patches of ecological zoning types A, B, C, D, E, F, G, H and I as cluster centers, that is, ecological zoning types A, B, C, D, E, F, G, H and I are the top 9 ecological zoning types in terms of patch number and patch area ratio in the above coupling results;

[0122] Step 2: Perform initial classification based on the distance between the remaining broken patch centers (centroids) and each cluster center (centroid), that is, the straight-line distance between two points. The calculation formula is as follows:

[0123]

[0124] In the formula, the absolute value of O1O2 is the distance between the center (center of mass) of a certain patch among the remaining patches to be clustered and the center (center of mass) of a certain cluster. O1 represents a certain cluster center (center of mass), and O2 represents a certain patch center (center of mass) among the remaining patches to be clustered. x1 , O y1 ) is the longitude and latitude coordinates of the cluster center (centroid), (O x2 , O y2 ) are the latitude and longitude coordinates of the center (centroid) of the patch;

[0125] Step 3: Obtain the distance between a patch and each cluster center in turn. If the patch is closest to a cluster center, the ecological zoning type of the patch is the ecological zoning type to which the cluster center belongs. The calculation formula is as follows:

[0126] D min ={d A , d E ,d C ,d D , ...}

[0127] Where D min is the minimum distance from the center (center of mass) of the patch to the center (center of mass) of each cluster, d A d B d C d D ,…represent the straight-line distance from the patch to the cluster centers (centroids) of A, B, C, D,… respectively.

[0128] Optionally, the mapping relationship between the preset temperature and the plateau climate zone type in the above step S104 is to obtain five plateau climate zone types with -4°C, -2°C, 4°C and 6°C as critical values. The five plateau climate zone types include high cold, cold, cold and warm, cool and warm. For specific mapping relationships, see Table 1.

[0129] Table 1 Classification of plateau climate zones

[0130] Serial number Climate zone scope 1 Alpine TV≤-4℃ 2 cold -4℃<TV≤-2℃ 3 cold temperature -2℃<TV≤4℃ 4 Cool 4℃<TV≤6℃ 5 warmth TV>6℃

[0131] The preset mapping relationship between surface relief and plateau landform type is 20m, 30m, 200m, 500m, 1000m and 2500m as the critical values ​​to obtain six plateau landform types. The six plateau landform types include plains, terraces, hills, medium mountains, high mountains and extremely high mountains. The specific mapping relationship is shown in Table 2.

[0132] Table 2 Classification of plateau landform types

[0133]

[0134]

[0135] The preset mapping relationship between land use types and road ecological types is to reclassify land use types to obtain eight road ecological types, including urban ecological type, cultivated ecological type, mixed agricultural and pastoral ecological type, mixed agricultural and forestry ecological type, grassland ecological type, woodland ecological type, water wetland ecological type and arid desert ecological type.

[0136] Optionally, reclassify land use types including:

[0137] Reclassify cities, towns, and industrial and mining areas into urban ecological types;

[0138] Reclassify paddy fields, dry fields, and oases into cultivated ecological types;

[0139] reclassification of agro-pastoral and pastoral farming into mixed agro-pastoral ecological types;

[0140] reclassification of agroforestry and forestry as mixed agroforestry ecological types;

[0141] Typical grass, patchy grass, and sparse grass were reclassified into pasture ecotypes;

[0142] reclassified dense forest, semi-dense forest, open forest and shrub forest into woodland ecotypes;

[0143] Reclassify water bodies, swamps, beaches, ice and snow, and frozen deserts into water wetland ecotypes;

[0144] Deserts, deserts and bare rock gobi were reclassified as arid desert ecological types.

[0145] The specific mapping relationship between land use type and road ecological type is shown in Table 3.

[0146] Table 3 Classification of road area ecological types

[0147]

[0148]

[0149] In an embodiment of the present application, plateau climate zone type, plateau landform type and road ecological type are used as three indicators to perform multi-indicator spatial coupling to obtain a coupling result. The ecological zoning type corresponding to the patch is a combination of indicator categories of plateau climate zone type, plateau landform type and road ecological type. A total of 240 ecological zoning types are divided, such as the high-altitude and medium-mountain urban ecological zone, the cold-temperate plain farming ecological zone, the warm extremely high mountain arid desert ecological zone, etc. Table 4 shows the zoning schemes for some ecological zoning types.

[0150] Table 4 Ecological zoning scheme for transportation corridors

[0151]

[0152]

[0153]

[0154] In the embodiments of the present application, the mapping relationship between temperature and plateau climate zone type, the mapping relationship between surface undulation and plateau landform type, and the mapping relationship between land use type and road ecological type are preset, and the plateau transportation corridor ecological zoning thematic database is improved to provide a refined data basis for the plateau transportation corridor ecological zoning. Subsequently, based on the meteorological station data, digital elevation model data and land use type data combined with the aforementioned data basis, multi-indicator spatial coupling and multi-indicator-patch clustering are integrated to realize the ecological zoning of the plateau transportation corridor ecological zoning study area. It overcomes the problems of traditional reinforcement methods being too general and lacking in adaptability.

[0155] Based on the plateau transportation corridor ecological zoning method provided in the aforementioned embodiment, the present application also provides a plateau transportation corridor ecological zoning device accordingly. Figure 4 This is a schematic diagram of the structure of a plateau transportation corridor ecological zoning device provided in an embodiment of the present application. Figure 4 As shown, the ecological zoning device of the plateau transportation corridor includes:

[0156] The data acquisition module 401 is used to acquire meteorological station data, digital elevation model data and land use type data in the plateau transportation corridor ecological zoning study area based on the vector range of the plateau transportation corridor ecological zoning study area;

[0157] An interpolation module 402 is used to interpolate the temperature data of the missing period in the meteorological station data based on the existing temperature data in the meteorological station data to obtain the interpolated temperature data;

[0158] A surface relief determination module 403 is used to obtain the surface relief in the plateau transportation corridor ecological zoning study area based on digital elevation model data;

[0159] The type determination module 404 is used to determine the plateau climate zone type of the plateau transportation corridor ecological zoning study area based on the interpolated temperature data according to the preset mapping relationship between the temperature and the plateau climate zone type; determine the plateau landform type of the plateau transportation corridor ecological zoning study area based on the surface undulation in the plateau transportation corridor ecological zoning study area according to the preset mapping relationship between the surface undulation and the plateau landform type; determine the road area ecological type of the plateau transportation corridor ecological zoning study area based on the land use type data in the plateau transportation corridor ecological zoning study area according to the preset mapping relationship between the land use type and the road area ecological type;

[0160] The coupling module 405 is used to perform multi-indicator spatial coupling on the plateau climate zone type, plateau landform type and road ecological type of the plateau transportation corridor ecological zoning study area as three indicators to obtain a coupling result; the coupling result includes multiple patches, each patch corresponds to an ecological zoning type, and the ecological zoning type is a type expressed by a combination of the three indicators;

[0161] A clustering module 406 is used to perform multi-index-plaque clustering based on the coupling result to obtain a clustering result;

[0162] The ecological zoning module 407 is used to perform ecological zoning for the plateau transportation corridor ecological zoning study area according to the clustering results.

[0163] Optionally, the meteorological station data includes station number, station name, station accuracy, station latitude, and temperature value observed at the station;

[0164] Optionally, the interpolation module 402 is specifically used to: average and assign values ​​to the temperature data of the missing period in the meteorological station data by the existing temperature data of the first three and last three periods of the same length as the missing period, and the calculation formula is:

[0165]

[0166] Among them, TV i represents the average temperature value of the i-th period; i-3, i-2, i-1, i+1, i+2, i+3 represent the i-3rd period, i-2nd period, i-1st period, i+1st period, i+2nd period, i+3rd period respectively; TV i-3 , TV i-2 , TV i-1 , TV i+1 , TV i+2 , TV i+3 They represent the average temperature values ​​of the corresponding time periods.

[0167] Optionally, the coupling module 405 is specifically used to: use the plateau climate zone type as the first indicator, the plateau landform type as the second indicator, and the road ecological type as the third indicator;

[0168] Obtain the number of plaques a, b and c respectively included in the first index, the second index and the third index; a, b and c are all positive integers;

[0169] The first, second and third indicators are spatially superimposed to obtain the total indicator, which contains a×b×c patches and is the coupling result.

[0170] Optionally, the clustering module 406 is specifically used to: select the centroids of patches of N ecological zoning types as cluster centers based on the coupling results; N is an integer greater than or equal to 1 and less than the total number of ecological zoning types in the coupling results, and the ecological zoning type of the cluster center is the ecological zoning type with the top N plaque numbers and plaque area proportions in the coupling results;

[0171] The ecological zoning types of the remaining patches are determined based on the distances between the centroids of the remaining patches and the centers of each cluster.

[0172] Optionally, the clustering module 406 determines the ecological zoning types of the remaining patches according to the distances between the centroids of the remaining patches and the centers of each cluster, including:

[0173] If the centroid of a patch among the remaining patches has the shortest straight-line distance to a cluster center, the ecological zoning type of the patch is the ecological zoning type of the cluster center.

[0174] Optionally, the device also includes a mapping relationship storage module for storing the mapping relationship between preset temperature and plateau climate zone type, the mapping relationship between surface undulation and plateau landform type, and the mapping relationship between land use type and road ecological type required by the type determination module 404.

[0175] Optionally, the mapping relationship between the preset temperature and the plateau climate zone type stored in the mapping relationship storage module is divided into five plateau climate zone types with -4°C, -2°C, 4°C and 6°C as critical values, and the five plateau climate zone types include high cold, cold, cold temperature, cool temperature and warm;

[0176] The preset mapping relationship between the surface relief and the plateau landform type stored in the mapping relationship storage module is that 20m, 30m, 200m, 500m, 1000m and 2500m are used as the critical values ​​to obtain six plateau landform types, and the six plateau landform types include plains, terraces, hills, medium mountains, high mountains and extremely high mountains;

[0177] The mapping relationship between the preset land use types and road ecological types stored in the mapping relationship storage module is to reclassify the land use types to obtain eight road ecological types, including urban ecological type, cultivated ecological type, mixed agricultural and pastoral ecological type, mixed agricultural and forestry ecological type, grassland ecological type, woodland ecological type, water wetland ecological type and arid desert ecological type.

[0178] Optionally, the mapping relationship storage module includes a land use type reclassification module, which is used to:

[0179] Reclassify cities, towns, and industrial and mining areas into urban ecological types;

[0180] Reclassify paddy fields, dry fields, and oases into cultivated ecological types;

[0181] reclassification of agro-pastoral and pastoral farming into mixed agro-pastoral ecological types;

[0182] reclassification of agroforestry and forestry as mixed agroforestry ecological types;

[0183] Typical grass, patchy grass, and sparse grass were reclassified into pasture ecotypes;

[0184] reclassified dense forest, semi-dense forest, open forest and shrub forest into woodland ecotypes;

[0185] Reclassify water bodies, swamps, beaches, ice and snow, and frozen deserts into water wetland ecotypes;

[0186] Deserts, deserts and bare rock gobi were reclassified as arid desert ecological types.

[0187] In addition, an embodiment of the present application also provides a plateau transportation corridor ecological zoning device, which includes a processor and a memory.

[0188] The memory is used to store computer program code and transmit the program code to the processor;

[0189] The processor is used to execute the steps of the plateau transportation corridor ecological zoning method introduced in any of the above method embodiments according to the instructions in the program code.

[0190] In addition, the embodiment of the present application also provides a plateau transportation corridor ecological zoning mapping system, the mapping system comprising: a plateau transportation corridor ecological zoning visualization mapping information library, a plateau transportation corridor ecological zoning thematic map mapping component, and a plateau transportation corridor ecological zoning device of the above device embodiment;

[0191] The plateau transportation corridor ecological zoning visualization mapping information base and the plateau transportation corridor ecological zoning thematic map mapping component are respectively connected to the processor in the plateau transportation corridor ecological zoning device;

[0192] The information database for visual mapping of ecological zoning of plateau transportation corridors includes symbol library, color library and annotation library;

[0193] The plateau transportation corridor ecological zoning thematic map mapping component is used to receive the plateau transportation corridor ecological zoning results from the processor; based on the plateau transportation corridor ecological zoning results, the plateau transportation corridor ecological zoning visualization mapping information database is used to produce the plateau transportation corridor ecological zoning thematic map; the plateau transportation corridor ecological zoning thematic map includes a scale, a compass and a self-made legend.

[0194] In an example implementation, a plateau transportation corridor ecological zoning thematic map is produced based on the plateau transportation corridor ecological zoning results and using the plateau transportation corridor ecological zoning visualization mapping information database. Figure 5 , Figure 6 ,in Figure 5 A schematic diagram of a thematic map provided in an embodiment of the present application, Figure 6 A thematic map of ecological zoning of plateau transportation corridors provided in an embodiment of the present application.

[0195] Scale bars, compasses and self-made legends were inserted into the thematic maps. Table 5 shows the self-made legends corresponding to some ecological zoning types.

[0196] Table 5

[0197]

[0198]

[0199] For details about the legend in Table 5, see Figure 7 , Figure 7 Schematic diagram of some ecological zoning types provided for the embodiments of this application.

[0200] In an embodiment of the present application, a plateau transportation corridor ecological zoning mapping system is provided, which realizes the collaborative processing of multi-source mapping data, on-demand customization of layouts, multi-dimensional visual expression, and standardized mapping output, thereby improving the level of automated and specialized mapping in the transportation and ecological and environmental protection industries, and can be widely used in transportation construction, operation, management, and road biodiversity protection.

[0201] It should be noted that each embodiment in this specification is described in a progressive manner, and the same and similar parts between the embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments. In particular, for the device, equipment and system embodiments, since they are basically similar to the method embodiments, the description is relatively simple, and the relevant parts can be referred to the partial description of the method embodiment. The device, equipment and system embodiments described above are merely schematic, wherein the units described as separate components may or may not be physically separated, and the components indicated as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. Ordinary technicians in this field can understand and implement it without paying creative labor.

[0202] The above is only a specific implementation of the present application, but the protection scope of the present application is not limited thereto. Any changes or substitutions that can be easily thought of by a person skilled in the art within the technical scope disclosed in the present application should be included in the protection scope of the present application. Therefore, the protection scope of the present application should be based on the protection scope of the claims.

Claims

1. A plateau transportation corridor ecological zoning method, characterized in that: The method comprises: Acquire meteorological station data, digital elevation model data, and land use type data within the plateau transportation corridor ecological zoning study area based on the vector range of the plateau transportation corridor ecological zoning study area; Interpolating the temperature data of the missing period in the meteorological station data based on the existing temperature data in the meteorological station data to obtain interpolated temperature data; Based on the digital elevation model data, the surface undulation within the plateau transportation corridor ecological zoning study area is obtained; According to the preset mapping relationship between temperature and plateau climate zone type, the plateau climate zone type of the plateau transportation corridor ecological zoning study area is determined based on the interpolated temperature data; according to the preset mapping relationship between surface undulation and plateau landform type, the plateau landform type of the plateau transportation corridor ecological zoning study area is determined based on the surface undulation in the plateau transportation corridor ecological zoning study area; according to the preset mapping relationship between land use type and road ecological type, the road ecological type of the plateau transportation corridor ecological zoning study area is determined based on the land use type data in the plateau transportation corridor ecological zoning study area; The plateau climate zone type, plateau landform type and road ecological type of the plateau transportation corridor ecological zoning study area are used as three indicators to perform multi-indicator spatial coupling to obtain a coupling result; the coupling result includes multiple patches, each patch corresponds to an ecological zoning type, and the ecological zoning type is a type expressed by a combination of the three indicators; Performing multi-index-plaque clustering based on the coupling result to obtain a clustering result; Based on the clustering results, ecological zoning is carried out for the plateau transportation corridor ecological zoning study area.

2. The method according to claim 1, characterized in that The meteorological station data includes station number, station name, station accuracy, station latitude, and temperature value observed at the station; The interpolating the temperature data of the missing period in the meteorological station data based on the existing temperature data in the meteorological station data to obtain the interpolated temperature data includes: The temperature data of the missing period in the meteorological station data is averaged and assigned by the existing temperature data of the first three and last three periods of the same length as the missing period, and the calculation formula is: Among them, TV i represents the average temperature value of the i-th period; i-3, i-2, i-1, i+1, i+2, i+3 represent the i-3rd period, i-2nd period, i-1st period, i+1st period, i+2nd period, i+3rd period respectively; TV i-3 , TV i-2 , TV i-1 , TV i+1 , TV i+2 , TV i+3 They represent the average temperature values ​​of the corresponding time periods.

3. The method according to claim 1, characterized in that The plateau climate zone type, plateau landform type and road ecological type of the plateau transportation corridor ecological zoning study area are used as three indicators to perform multi-indicator spatial coupling to obtain coupling results, including: The plateau climate zone type is used as the first indicator, the plateau landform type is used as the second indicator, and the road ecological type is used as the third indicator; Obtain the number of plaques a, b and c respectively included in the first index, the second index and the third index; a, b and c are all positive integers; The first index, the second index and the third index are spatially superimposed to obtain a total index, wherein the total index includes a×b×c patches and is a coupling result.

4. The method according to claim 1, characterized in that: The performing multi-index-plaque clustering based on the coupling result includes: Based on the coupling result, the centroids of patches of N types of ecological zoning are selected as cluster centers respectively; N is an integer greater than or equal to 1 and less than the total number of ecological zoning types in the coupling result, and the ecological zoning type of the cluster center is the top N ecological zoning types in terms of the number of patches and the proportion of patch area in the coupling result; The ecological zoning types of the remaining patches are determined based on the distances between the centroids of the remaining patches and the centers of each cluster.

5. The method according to claim 4, characterized in that The ecological zoning types of the remaining patches are determined according to the distances between the centroids of the remaining patches and the centers of each cluster, including: If the centroid of a patch among the remaining patches has the shortest straight-line distance to a cluster center, the ecological zoning type of the patch is the ecological zoning type of the cluster center.

6. The method according to any one of claims 1 to 5, characterized in that: The mapping relationship between the preset temperature and the plateau climate zone type is to divide the five plateau climate zone types into -4°C, -2°C, 4°C, and 6°C as the critical values, and the five plateau climate zone types include high cold, cold, cold temperature, cool temperature, and warm; The preset mapping relationship between the surface relief and the plateau landform type is that 20m, 30m, 200m, 500m, 1000m and 2500m are used as the critical values ​​to obtain six plateau landform types, and the six plateau landform types include plains, terraces, hills, medium mountains, high mountains and extremely high mountains; The mapping relationship between the preset land use type and the road ecological type is to reclassify the land use type to obtain eight road ecological types, and the eight road ecological types include urban ecological type, cultivated ecological type, mixed agricultural and pastoral ecological type, mixed agricultural and forestry ecological type, grassland ecological type, woodland ecological type, water wetland ecological type and arid desert ecological type.

7. The method according to claim 6, characterized in that The reclassification of land use types includes: Reclassify cities, towns, and industrial and mining areas into urban ecological types; Reclassify paddy fields, dry fields, and oases into cultivated ecological types; reclassification of agro-pastoral and pastoral farming into mixed agro-pastoral ecological types; reclassification of agroforestry and forestry as mixed agroforestry ecological types; Typical grass, patchy grass, and sparse grass were reclassified into pasture ecotypes; reclassified dense forest, semi-dense forest, open forest and shrub forest into woodland ecotypes; Reclassify water bodies, swamps, beaches, ice and snow, and frozen deserts into water wetland ecotypes; Deserts, deserts and bare rock gobi were reclassified as arid desert ecological types.

8. A plateau transportation corridor ecological zoning device, characterized in that: The device comprises: A data acquisition module, used to acquire meteorological station data, digital elevation model data and land use type data within the plateau transportation corridor ecological zoning study area based on the vector range of the plateau transportation corridor ecological zoning study area; An interpolation module, used to interpolate the temperature data of the missing period in the meteorological station data based on the existing temperature data in the meteorological station data to obtain the interpolated temperature data; A surface relief determination module, used for obtaining the surface relief in the plateau transportation corridor ecological zoning study area based on the digital elevation model data; A type determination module, for determining the plateau climate zone type of the plateau transportation corridor ecological zoning study area based on the interpolated temperature data according to a preset mapping relationship between temperature and plateau climate zone type; determining the plateau landform type of the plateau transportation corridor ecological zoning study area based on the surface undulation in the plateau transportation corridor ecological zoning study area according to a preset mapping relationship between surface undulation and plateau landform type; determining the road area ecological type of the plateau transportation corridor ecological zoning study area based on the land use type data in the plateau transportation corridor ecological zoning study area according to a preset mapping relationship between land use type and road area ecological type; A coupling module is used to perform multi-indicator spatial coupling on the plateau climate zone type, plateau landform type and road ecological type of the plateau transportation corridor ecological zoning study area as three indicators to obtain a coupling result; the coupling result includes multiple patches, each patch corresponds to an ecological zoning type, and the ecological zoning type is a type expressed by a combination of the three indicators; A clustering module, used for performing multi-index-plaque clustering based on the coupling result to obtain a clustering result; The ecological zoning module is used to carry out ecological zoning for the plateau transportation corridor ecological zoning study area according to the clustering results.

9. A plateau transportation corridor ecological zoning device, characterized in that: include: Processor and memory: The memory is used to store program code and transmit the program code to the processor; The processor is used to execute the steps of the plateau transportation corridor ecological zoning method described in any one of claims 1 to 7 according to the instructions in the program code.

10. A plateau transportation corridor ecological zoning mapping system, characterized in that: The mapping system comprises: a plateau transportation corridor ecological zoning visualization mapping information base, a plateau transportation corridor ecological zoning thematic map mapping component, and the plateau transportation corridor ecological zoning device according to claim 9; The plateau transportation corridor ecological zoning visualization mapping information base and the plateau transportation corridor ecological zoning thematic map mapping component are respectively connected to the processor in the plateau transportation corridor ecological zoning device; The plateau transportation corridor ecological zoning visualization mapping information database includes a symbol library, a color library and an annotation library; The plateau transportation corridor ecological zoning thematic map mapping component is used to receive the plateau transportation corridor ecological zoning results of the processor; based on the plateau transportation corridor ecological zoning results, the plateau transportation corridor ecological zoning visualization mapping information library is used to produce a plateau transportation corridor ecological zoning thematic map; the plateau transportation corridor ecological zoning thematic map includes a scale, a compass and a self-made legend.

Images