Land resource planning method and system based on cadastral information

Abstract

This application provides a land resource planning method and system based on cadastral information, relating to the field of territorial spatial planning technology. This application obtains spatial patches of the target land parcel, extracts land category identifiers and area control indicators, and establishes a mapping relationship between the area control indicators and the spatial patches. It determines the corresponding land use assessment parameters based on the land category identifiers; responds to spatial replacement requests, extracts source land and target land, and generates a transfer path; along the transfer path, performs a matching test based on the mapping relationship and the replacement area, combined with the land use assessment parameters, and generates a land planning scheme based on the verification results and the vector direction of the transfer path. This enables path-based management of all elements and spatial graphics, solves the problem of unrecorded accounting for land use balance and land consolidation indicators, and improves the accuracy of territorial spatial use control and the scientific nature of planning schemes.

Description

Technical Field

[0001] This application relates to the field of land spatial planning technology, and in particular to a land resource planning method and system based on cadastral information. Background Technology

[0002] With the establishment and improvement of my country's territorial spatial planning system, territorial spatial use control has become an inevitable requirement for improving the spatial governance system and promoting the modernization of governance capabilities. Relevant Chinese documents have clearly defined detailed planning as the sole legal basis for implementing and managing territorial spatial use control. The scope of detailed planning has expanded from the original urban planning areas to the entire territorial space, and the controlled elements have extended from construction land to all elements. Against this backdrop, how to achieve precise control with full coverage and comprehensive coordination of all elements in detailed planning has become a key technical issue that urgently needs to be addressed in the field of territorial spatial planning.

[0003] In existing technologies, implementing land use control through detailed planning faces several prominent problems. For example, at the technical level, the lack of cadastral information exposes a weakness in the legal basis of planning outcomes. Planning and cadastral management are the responsibilities of different departments. In the past, planners were accustomed to preparing plans at the physical planning level, without sufficient consideration for the socio-economic attributes of land parcels. They often focused only on the future development rights of the land, ignoring the current state of the land, resulting in problems such as misallocation of development rights and frequent planning revisions. In terms of control measures, existing control plans lack indicator-based constraint mechanisms. Although there are overall scale constraints for construction land, cultivated land, and forest land, etc., these issues persist. However, the flexibility in converting between land uses is relatively large, and spatial replacement is often carried out through methods such as balancing land use with compensation and linking increases and decreases in urban and rural construction land. When optimizing spatial layout, the land use and indicators involved in balancing land use with compensation and linking increases and decreases cannot be clearly seen in the planning map, resulting in indicators being recorded in vain, making it difficult to effectively guide land transfer and weakening the control role of planning. In addition, the control objectives of the three types of land use—construction land, agricultural land, and ecological land—are different, but the existing control system lacks differentiated control measures, and the relationship and dynamic balance between the conversion between different land use types are difficult to be effectively reflected in a unified planning platform.

[0004] Therefore, there is an urgent need for a land resource planning method that can deeply integrate cadastral information into the planning process, realize the management of all elements and spatial graphic paths, and support differentiated land use control, so as to overcome the shortcomings of existing technologies such as the disconnect between indicators and space and the single control method. Summary of the Invention

[0005] The purpose of this application is to provide a land resource planning method and system based on cadastral information to solve the problems of poor effectiveness of land resource planning caused by the disconnect between cadastral information and planning, the lack of clear accounting for land occupation and compensation balance and land use linkage indicators, and insufficient accuracy of land use control in the existing technology.

[0006] To address the aforementioned technical problems, firstly, this application provides a land resource planning method based on cadastral information, applied to a spatial use control platform. The method includes:

[0007] Obtain spatial map features of the target land parcel, extract the land use identifiers and area control indicators of the spatial map features, and establish a mapping relationship between the area control indicators and the corresponding spatial map features;

[0008] Based on the land use classification, corresponding land use assessment parameters are determined. Specifically, when the land use classification is construction land, the land use assessment parameter is the utilization intensity parameter; when the land use classification is agricultural land, the land use assessment parameter is the land occupation and compensation balance parameter; and when the land use classification is ecological land, the land use assessment parameter is the red line intrusion parameter.

[0009] In response to a spatial replacement request for the target parcel, the source land and target land corresponding to the spatial replacement request are extracted, and a flow path from the source land to the target land is generated, wherein both the source land and the target land are located within the target parcel;

[0010] Obtain the replacement area of ​​the spatial replacement request, perform matching verification along the transfer path based on the mapping relationship and the replacement area, combined with the land use assessment parameters, and generate the corresponding land planning scheme based on the verification results and the vector direction of the transfer path.

[0011] Optionally, along the transfer path, a matching verification is performed based on the mapping relationship and the replacement area, combined with the land use assessment parameters, and a corresponding land planning scheme is generated based on the verification results and the vector direction of the transfer path, including:

[0012] If the land use assessment parameters meet the corresponding preset compliance conditions, and an incremental map patch with an area consistent with the replacement area is retrieved for the spatial map patch corresponding to the target land use, then the matching is determined to be successful;

[0013] The mapping relationship is updated based on the coordinate data of the incremental map features to generate a land planning scheme;

[0014] Alternatively, if the land use assessment parameters do not meet the preset compliance conditions, or if the incremental map patch is not found, the matching is deemed to have failed.

[0015] In the event of a failed match, a search guide is extracted based on the land use category identifier of the source land, and alternative land parcels that meet the preset compliance conditions are searched around the preset anchor point along the search guide.

[0016] The spare land parcels are compensated to the spatial parcels corresponding to the target land use, and the mapping relationship is updated to generate a land planning scheme.

[0017] Optionally, for the spatial patch corresponding to the target land use, an incremental patch with an area consistent with the replacement area is retrieved, including:

[0018] Based on the spatial permutation request, generate the target boundary coordinates of the spatial patch corresponding to the target land after spatial permutation;

[0019] By comparing the target boundary coordinates with the initial boundary coordinates of the spatial patch corresponding to the target land before spatial permutation, a set of candidate patches generated by the boundary expansion of the spatial patch is extracted, and the expansion area of ​​each candidate patch is calculated.

[0020] Calculate the first area difference between the expanded area and the replacement area;

[0021] When the first area difference is less than the preset area tolerance, the corresponding candidate patch is determined as an incremental patch.

[0022] Optionally, a search guide is extracted based on the land use category identifier of the source land, and alternative land parcels that meet the preset compliance conditions are searched around preset anchor points along the search guide, including:

[0023] The corresponding preset anchor point is determined based on the land use type identifier of the source land, and the vector direction of the flow path is generated based on the spatial control rules matched by the land use type identifier.

[0024] Using the centroid of the preset anchor point as the reference point, a fan-shaped retrieval area is constructed by expanding the preset angle range along the vector direction;

[0025] Extract the set of unoccupied free patches within the sector-shaped search area;

[0026] From the set of idle patches, select idle patches that meet the preset compliance conditions and whose area is not less than the second area difference, and use them as backup patches. The second area difference is the difference between the replacement area and the area of ​​the retrieved incremental patches.

[0027] Optionally, the spare land parcels are compensated to the spatial parcels corresponding to the target land use, and the mapping relationship is updated to generate a land planning scheme, including:

[0028] Obtain the backup boundary coordinates of the backup patch, and concatenate the backup boundary coordinates with the target boundary coordinates to obtain the target extended boundary;

[0029] Adjust the original area control indicators associated with the target land in the mapping relationship according to the second area difference to obtain the target area control indicators, and re-associate the target area control indicators with the corresponding spatial patches to update the mapping relationship;

[0030] Based on the updated mapping relationship and the target extended boundary, a planning scheme map containing the adjusted plot outline is output on the spatial use control platform as a land planning scheme.

[0031] Optionally, the area control indicators include the land occupation and compensation balance indicator, the increase and decrease linkage indicator, or the ecological red line indicator;

[0032] Establishing the mapping relationship between the area control indicators and the corresponding spatial patches includes:

[0033] When the land type is identified as agricultural land, agricultural land patches with agricultural land attributes and whose area matches the occupation-compensation balance index are selected from the spatial patches, and the occupation-compensation balance index is associated with the boundary node of the agricultural land patch to establish a first mapping association.

[0034] When the land type is identified as construction land, construction land parcels with construction land attributes and whose areas match the increase-decrease linkage index are selected from the spatial parcels. The increase-decrease linkage index is then associated with the boundary nodes of the construction land parcels to establish a second mapping association.

[0035] When the land use type is identified as ecological land, ecological land parcels with ecological land use attributes and whose area matches the ecological red line index are selected from the spatial parcels. The ecological red line index is then associated with the boundary nodes of the ecological land parcels to establish a third mapping association.

[0036] Merge the first mapping association, the second mapping association, and the third mapping association to obtain the mapping relationship.

[0037] Optionally, when the land type is identified as construction land, the utilization intensity parameters are determined, including:

[0038] Extract the three-dimensional vertex coordinates of the planned building model, and calculate the first volume of the planned building model based on the three-dimensional vertex coordinates;

[0039] Obtain preset three-dimensional control lines and calculate the second volume inside the three-dimensional constraint boundary formed by the three-dimensional control lines;

[0040] Calculate the volume difference between the second volume and the first volume, and determine the ratio of the volume difference to the second volume as the strength parameter.

[0041] Optionally, when the land type is identified as ecological land, the parameters for red line intrusion are determined, including:

[0042] Extract the first spatial region within the vector boundary of the target patch, and the second spatial region defined by the ecological protection red line;

[0043] Extract the overlapping area between the first spatial region and the second spatial region;

[0044] Calculate the projected area of ​​the overlapping region and divide the projected area by the projected area of ​​the first spatial region to obtain the overlap ratio;

[0045] The numerical combination of the projected area and the overlap ratio is determined as the red line intrusion parameter.

[0046] Optionally, the source land and target land corresponding to the spatial replacement request are extracted, and a transfer path from the source land to the target land is generated, including:

[0047] The spatial permutation request is parsed to extract the source spatial coordinates and source land class attributes of the source land, as well as the target spatial coordinates and target land class attributes of the target land.

[0048] The source space coordinates and the target space coordinates are used as path nodes, and connecting line segments with vector directions are constructed between the path nodes;

[0049] The attribute change identifier, which represents the transformation from the source land type attribute to the target land type attribute, is recorded on the connecting line segment to obtain the flow path.

[0050] Secondly, this application provides a land resource planning system based on cadastral information, applied to a spatial use control platform. The system includes:

[0051] A module is established to acquire spatial patches of the target land parcel, extract the land type identifier and area control index of the spatial patches, and establish a mapping relationship between the area control index and the corresponding spatial patches;

[0052] The determination module is used to determine the corresponding land use assessment parameters based on the land use identifier, wherein when the land use identifier is construction land, the land use assessment parameter is the utilization intensity parameter; when the land use identifier is agricultural land, the land use assessment parameter is the land occupation and compensation balance parameter; and when the land use identifier is ecological land, the land use assessment parameter is the red line intrusion parameter.

[0053] The extraction module is used to respond to a spatial replacement request for the target parcel, extract the source land and target land corresponding to the spatial replacement request, and generate a flow path from the source land to the target land, wherein both the source land and the target land are located within the target parcel;

[0054] The generation module is used to obtain the replacement area of ​​the spatial replacement request, perform matching verification along the transfer path according to the mapping relationship and the replacement area, combined with the land use assessment parameters, and generate a corresponding land planning scheme based on the verification result and the vector direction of the transfer path.

[0055] The land resource planning method based on cadastral information provided in this application has the following beneficial effects: By introducing cadastral information across the entire region at the initial stage of planning, a mapping relationship between area control indicators and spatial patches is established, enabling the land occupation and compensation balance indicators and land increase / decrease linkage indicators to be immediately implemented on the planning map and find their corresponding spaces. This achieves path-following and real-time observable management of all-element indicators and spatial graphics, avoiding the problem of indicators being suspended and creating a cumulative cyclical effect. Differentiated land use assessment parameters are designed for construction land and ecological land respectively. Construction land uses utilization intensity parameters to assess the degree of intensive and efficient utilization, agricultural land uses land occupation and compensation balance parameters to ensure the bottom line of arable land protection, and ecological land uses red line intrusion parameters to assess the strength of protection, realizing differentiated quantitative expression of control objectives for different land use types. Through the transfer path mechanism, a vectorized tracking system from source land to target land is established in the spatial replacement process. Combined with matching verification and land planning scheme generation, this ensures the accurate transmission of indicators and the rational allocation of space in the spatial replacement process, improving the scientificity and operability of land planning schemes. Attached Figure Description

[0056] To more clearly illustrate the technical solutions of the embodiments of this application or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0057] Figure 1 This is an overall flowchart of a land resource planning method based on cadastral information disclosed in this application;

[0058] Figure 2 This is a schematic diagram showing the location of the centroid coordinates of an agricultural patch disclosed in this application;

[0059] Figure 3 This is a schematic diagram of the structure of a land resource planning system based on cadastral information disclosed in this application. Detailed Implementation

[0060] To enable those skilled in the art to better understand the present application, the present application will be further described in detail below with reference to the accompanying drawings and specific embodiments. Obviously, the described embodiments are merely some embodiments of the present application, and not all embodiments. Based on the embodiments in this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0061] Under the integrated control of land use, detailed planning for implementing land use control requires both spatial layout and the transmission and spatial allocation of binding indicators. Due to varying control intensities, ecological protection red lines and permanent basic farmland indicators are relatively easy to control. However, while construction land, arable land, and forest land have overall scale constraints, the flexibility in converting between land uses is greater. Spatial substitution is often achieved through methods such as land occupation-compensation balance and linking increases and decreases in urban and rural construction land, increasing the difficulty of control. In existing technologies, these involved plots are scattered and fragmented. When optimizing spatial layout, the lack of comprehensive indicators and spatial path-based graphical management methods leads to indicators being recorded in vain, ultimately resulting in low effectiveness of land resource planning.

[0062] The core of this application is to provide a land resource planning method based on cadastral information, applied to a spatial use control platform. A flowchart of one specific implementation method is shown below. Figure 1 As shown, the method includes:

[0063] S101. Obtain the spatial map of the target land parcel, extract the land type identifier and area control index of the spatial map, and establish the mapping relationship between the area control index and the corresponding spatial map.

[0064] In this embodiment, the spatial use control platform refers to a planning management system that integrates data from the national land space basic information platform. This platform can connect to the real estate registration database, land survey database, and planning results database.

[0065] The target parcel refers to a specific plot of land that needs to undergo land use planning or spatial replacement in this planning process; the spatial plot refers to all registered or surveyed land plots within the target parcel, and each spatial plot contains attributes such as vector boundary coordinates, land category code, ownership information, and area data.

[0066] Land use type identification refers to the land use type code assigned to the target land parcel according to the national land and sea use classification standard, such as construction land, agricultural land, ecological land, etc.; area control indicators refer to the binding area control values ​​transmitted from the higher-level plan to the local plan. These indicators include the land occupation and compensation balance indicator, the increase and decrease linkage indicator, and the ecological red line indicator.

[0067] It should be understood that the above mapping relationship can reflect cadastral information.

[0068] It should also be understood that the area control indicators include land occupation and compensation balance indicators, increase and decrease linkage indicators, or ecological red line indicators; different land categories correspond to different area control indicators.

[0069] Based on this, the mapping relationship establishment process includes three branch paths. With the support of the cadastral data of the whole region, the land use identifier of the spatial map pieces is first identified, and then a differentiated mapping strategy is adopted according to different land use types. Specifically, the specific process of establishing the mapping relationship between the area control index and the spatial map pieces is as follows: When the land use identifier is agricultural land, agricultural map pieces with agricultural land attributes and whose area matches the occupation and compensation balance index are selected from the spatial map pieces. The occupation and compensation balance index is associated with the boundary node of the agricultural map piece to establish the first mapping relationship.

[0070] Among them, the "balance between cultivated land occupation and compensation" index refers to the binding value that requires the replenishment of an equal quantity and quality of cultivated land after construction occupies cultivated land; the boundary node refers to the coordinate point on the vector boundary of the spatial patch. By associating the index with specific boundary nodes, the index is precisely anchored in space. The specific mapping and association process is as follows:

[0071] In the spatial map, all map patches with the land use code as the filtering condition are retrieved to form a candidate set of agricultural land patches, where each map patch contains its vector boundary coordinate set. and area value The numerical value of the land occupation balance index Izb is compared with the area value of each agricultural land patch in the candidate set. A comparison can be performed. Specifically, an area matching algorithm can be used, for example: calculating the absolute value of the difference between the value of the occupation-compensation balance index and the area of ​​each map patch. When the difference is less than the preset area tolerance, the area of ​​the patch is determined to match the occupancy balance index.

[0072] It should be noted that the area tolerance can be set according to the actual surveying accuracy, for example, it can be set to 1% of the value of the occupation-compensation balance index. If the area of ​​a single patch cannot meet the matching requirements, a combination matching strategy is adopted, that is, several adjacent patches are selected from the candidate set, their areas are calculated, and it is determined whether the combined area matches the value of the occupation-compensation balance index.

[0073] For successfully matched agricultural land parcels, the occupation-compensation balance index is associated with the boundary nodes of those parcels. Boundary nodes are coordinate points on the vector boundary of the spatial land parcel. The specific association method is as follows: in the vector boundary coordinate set of the matched land parcel... In this process, the centroid coordinates of the patch are selected as the primary anchor node, and the coordinates of the vertex farthest from the centroid on the vector boundary are selected as the secondary anchor node. The value of the occupancy balance index, the index type code, and the associated timestamp are written as attribute information into the primary and secondary anchor nodes to establish the first mapping association. Through this dual-node anchoring method, even if the patch undergoes minor boundary adjustments, the validity of the mapping relationship can still be verified through the secondary anchor node. First Mapping Association The specific formula can be formula (1):

[0074] ={ , Izb, (xa, ya), (xf, yf), type='ZB', timestamp} (1)

[0075] in, Let be each map patch in the candidate set of agricultural map patches; Izb is the value of the occupation-compensation balance index; (xa, ya) is the agricultural map patch. The centroid coordinates are the geometric center positions of all vertices of the vector boundary of the patch; (xf, yf) represents the agricultural patch. The coordinates of the vertex furthest from the centroid on the vector boundary are the coordinates of the auxiliary anchor node; type='ZB' is the index type code, where ZB is the abbreviation of the first letter of the pinyin for "occupation and supplementation", indicating that the mapping association belongs to the occupation and supplementation balance type; timestamp is the timestamp for establishing the mapping association, used to record the specific time when the index and spatial patch association were created.

[0076] Suppose that an agricultural plot exists within a certain coding unit. ,like Figure 2 As shown, its vector boundary consists of 8 coordinate nodes, and the coordinates of each node are as follows: (500176.0, 3045694.5) (500156.9, 3045712.6) (500134.0, 3045739.3) (500104.8, 3045719.1) (500068.9, 3045697.6) (500074.8, 3045660.0) (500099.5, 3045623.2) (500151.2, 3045636.6), with centroid coordinates of (500123.5, 3045678.2), the figure shows an agricultural patch within a certain tubular unit. The spatial relationship between the vector boundary, coordinate nodes, and centroid coordinates of the agricultural patch. The vector boundary consists of 8 coordinate nodes. to Connecting the nodes end to end sequentially forms a closed, irregular polygon. The positions of each coordinate node are marked with solid dots. The position marked with a pentagram is the centroid coordinate (500123.5, 3045678.2), calculated from the above eight vertices using the polygon centroid formula, and serves as the primary anchor node in the first mapping association. The position marked with a triangle is the vertex on the vector boundary farthest from the centroid. The vertex coordinates serve as the auxiliary anchor node in the first mapping association; the dashed line represents the distance from the centroid coordinates to the farthest vertex. The connecting lines are used to visually demonstrate the spatial distance relationship between the auxiliary anchor node and the main anchor node. The value of the land occupation balance index transmitted from the higher-level planning is 4980 square meters, and the area tolerance is set to 50 square meters. Therefore, |5000-4980|=20<50, indicating a successful match. The land occupation balance index value of 4980 and the index type code 'ZB' are associated with the centroid coordinates (500123.5, 3045678.2) of the map patch and the coordinates of the vertex farthest from the centroid coordinates on the vector boundary, completing the establishment of the first mapping association.

[0077] When the local land use category is identified as construction land, construction land parcels with construction land attributes and whose area matches the increase / decrease linkage index are selected from the spatial parcels. The increase / decrease linkage index is then associated with the boundary node of the construction land parcel to establish a second mapping association.

[0078] It should be understood that the process of establishing the mapping association is similar to that for agricultural land, but the filtering criteria are changed to construction land attributes. Second mapping association. The specific formula can be shown in formula (2):

[0079] ={ , Izj, (xb, yb), (xe, ye), type='ZJ', timestamp} (2)

[0080] in, Let I be each map patch in the candidate set for construction; let Izj be the value of the increase / decrease linkage index; and let (xb, yb) be the map patch to be constructed. The centroid coordinates; (xe, ye) represent the construction patch. The coordinates of the vertex furthest from the centroid on the vector boundary, i.e., the coordinates of the auxiliary anchor node; type='ZJ' is the index type code, where ZJ is the abbreviation of the first letter of the pinyin for increase and decrease, indicating that the mapping relationship belongs to the increase and decrease hook type.

[0081] When the land use category is identified as ecological land, the ecological red line indicator is extracted from the area control indicators. Ecological patches with ecological land use attributes and areas matching the ecological red line indicator are selected from the spatial patches. The ecological red line indicator is then associated with the boundary node of the ecological patch to establish a third mapping association. The ecological red line indicator refers to the area constraint value that needs strict protection within the ecological protection red line. The specific formula for the third mapping association M3 can be shown in formula (3):

[0082] ={ , Ist, (xc, yc), (xd, yd), type='ST', timestamp} (3)

[0083] in, For each ecological patch candidate set, Ist is the value of the ecological red line index; (xc, yc) represents the ecological patch. The centroid coordinates; (xd, yd) represent the ecological patch. The coordinates of the vertex furthest from the centroid on the vector boundary are the coordinates of the auxiliary anchor node; type='ST' is the index type code, where ST is the abbreviation of the first letter of the pinyin for "ecology", indicating that the mapping relationship belongs to the ecological red line type.

[0084] Finally, merge the first mapping association. Second mapping association and third mapping association To obtain the complete mapping relationship M, M= ⋃ ⋃ Furthermore, the merged mapping relationships are stored in a key-value pair data structure, where the key is the unique code of the map patch, and the value is the associated indicator information and anchor node information. Simultaneously, on the planning base map of the spatial use control platform, each mapping relationship is visually displayed using different colored markers. Specifically, the land occupation balance indicator is marked with a green circle at the centroid of the corresponding map patch, the increase / decrease linkage indicator is marked with an orange triangle at the centroid of the corresponding map patch, and the ecological red line indicator is marked with a red diamond at the centroid of the corresponding map patch, facilitating a clear view of the spatial distribution of indicators for planning personnel. The rules defining the mapping relationships are shown in Table 1.

[0085] Table 1 Mapping Rules for Area Control Indicators

[0086] Therefore, all land types have been managed using graphic data. The land occupation and compensation balance indicators, the increase and decrease linkage indicators, and the ecological red line indicators can be implemented in real time on the planning map and find their corresponding spaces, thus achieving the management goal of traceable paths and real-time observability.

[0087] S102. Determine the corresponding land use assessment parameters based on the land use type identifier, wherein when the land use type identifier is construction land, the land use assessment parameter is the utilization intensity parameter; when the land use type identifier is agricultural land, the land use assessment parameter is the land occupation and compensation balance parameter; and when the land use type identifier is ecological land, the land use assessment parameter is the red line intrusion parameter.

[0088] In this embodiment, land use assessment parameters are numerical indicators used to quantitatively assess the development and utilization status or protection status of the target land parcel. Different assessment methods are used for different types of land use. This differentiated assessment mechanism stems from the differences in the control objectives for different land use types in the national land use control. Specifically: construction land emphasizes intensive and efficient use, requiring an assessment of its spatial utilization degree; agricultural land emphasizes farmland protection, requiring an assessment of the implementation of the balance between land occupation and compensation; and ecological land emphasizes protection and restoration, requiring an assessment of the degree to which it is covered by the protection red line.

[0089] As a preferred embodiment, when the land type is identified as construction land, the utilization intensity parameter is determined by a three-dimensional spatial analysis method. The utilization intensity parameter represents the development intensity of the three-dimensional space of the target land parcel, and the sum of the utilization intensity parameter and the utilization rate of the three-dimensional space of the target land parcel is 1. Specifically, firstly, the three-dimensional vertex coordinates of the planning building model are extracted. The planning building model refers to the three-dimensional digital representation of the proposed building on the target land parcel in urban design or detailed planning. The model contains information such as the building's plan outline, building height, and floor plan. The three-dimensional vertex coordinates refer to the coordinate values ​​of each vertex of the planning building model in the three-dimensional coordinate system. Based on the three-dimensional vertex coordinates, the first volume of the planning building model is calculated using the three-dimensional convex hull decomposition method. The specific process is as follows: Firstly, the planning building model is decomposed into several horizontal cross-sectional layers according to the floors. The height of each cross-sectional layer is equal to the floor height h. Then, the area of ​​each cross-sectional layer is calculated using the shoelace formula. Finally, the product of the area of ​​each layer and the floor height is accumulated to obtain the first volume. The first volume is the physical space occupied by the planning building model. Its calculation formula is shown in formula (4).

[0090] =sum( ), i=1,2,...,N (4)

[0091] in, The first volume of the planned building model is N, where N is the total number of floors in the building. Let be the planar area of ​​the i-th layer cross section. The height of the i-th layer. The area of ​​each cross-sectional layer. The shoelace formula is used for calculation. This formula is a commonly used algorithm for calculating the area of ​​simple polygons, especially for polygons with n vertices {( , ), ( , ), ...,( , The area of ​​a simple polygon formed by )} is calculated using the formula shown in formula (5):

[0092] , j=1,2,...,n (5)

[0093] Where, when j=n, ( , )=( , ), that is, the two ends are connected to form a closed polygon.

[0094] Then, a pre-set three-dimensional control line is obtained. This three-dimensional control line refers to the three-dimensional boundary used in urban planning and management to control the spatial form of buildings. This control line can form a comprehensive envelope surface in three-dimensional space, which is formed by constraints such as building height limits, building setbacks, and sunlight spacing. After offsetting the setback distance d inward from the land boundary coordinates of the target plot, the buildable area polygon is obtained. Then, the building height limit is used as the basis for the calculation. As a height constraint, the second volume inside the three-dimensional constraint boundary formed by the three-dimensional control lines is calculated. The second volume is the overall spatial capacity allowed by the three-dimensional control lines, and its calculation formula is shown in formula (6):

[0095] = (6)

[0096] in, For the second volume, The area of ​​the polygon within the buildable area is calculated from the polygon obtained by offsetting the coordinates of the target land parcel's boundary line inward by a setback distance d. This refers to the building height limit, which is the upper limit constraint value of the three-dimensional control line on the building height. When the target land parcel has more complex constraints such as sunlight spacing, the calculation of the second volume needs to be modified according to the solar analysis results for the building height limit by region. In this case, the calculation of the second volume is carried out by accumulating the data from each region.

[0097] Calculate the volume difference between the second volume and the first volume. The specific formula is shown in formula (7):

[0098] = (7)

[0099] Volume difference This reflects the unused three-dimensional space capacity of the target land parcel under the current planning scheme. The ratio of the volume difference to the second volume is then calculated, and this ratio is determined as the utilization intensity parameter. The specific formula is as follows: = = The utilization intensity parameter ranges from 0 to 1. A larger value indicates more remaining usable space and a lower current utilization intensity; a smaller value indicates more efficient space utilization and a higher degree of intensification.

[0100] Assuming a land parcel has a boundary line enclosing an area of ​​8,000 square meters and a setback distance of 5 meters, what is the area of ​​the buildable area after the setback? =6800 square meters, building height limit =60 meters. Therefore, the second volume... =6800 × 60 = 408000 cubic meters. The land parcel is planned to house a 20-story building, with each floor 3 meters high and each floor having a cross-sectional area of ​​2000 square meters. Therefore, the first volume... =2000 × 3 × 20 = 120000 cubic meters. (Volume difference) =408000 - 120000 = 288000 cubic meters. Using strength parameters. = ≈0.706. This result indicates that the utilization rate of the three-dimensional space of the target site is approximately 29.4%. If the plan is adjusted to a 40-story building, then... =2000×3×40=240000 cubic meters = ≈0.412, therefore, the utilization rate increases to approximately 58.8%.

[0101] When the land use category is identified as agricultural land, the land occupation and replenishment balance parameters are determined. These parameters are quantitative representations of the occupation and replenishment of arable land during the spatial replacement process. Specifically, these parameters may include the quality grade and area of ​​the occupied arable land, as well as the quality grade and area of ​​the replenished arable land. For example, the compliance determination of the land occupation and replenishment balance parameters simultaneously satisfies two conditions: quantity balance requires that the area of ​​replenished arable land is not less than the area of ​​occupied arable land; quality balance requires that the quality grade of the replenished arable land is not lower than the quality grade of the occupied arable land.

[0102] When land is designated as ecological land, the parameters for ecological boundary intrusion are determined. Specifically, firstly, the first spatial region of the vector boundary envelope is extracted, which is the closed area of ​​the target land parcel on the plane projection; simultaneously, the second spatial region defined by the ecological protection red line is extracted, which is the space within which development and construction are prohibited by the ecological protection red line. The ecological protection red line refers to an area within the ecological space that has particularly important ecological functions and is subject to mandatory and strict protection. The area formed by this red line includes areas with important ecological functions such as water conservation, biodiversity maintenance, soil and water conservation, and windbreak and sand fixation.

[0103] Next, the overlapping region between the first and second spatial regions is extracted; that is, the spatial intersection of the two regions. For example, the Weiler-Atherton polygon clipping algorithm is used to calculate the geometric intersection of the two regions. The Weiler-Atherton algorithm is a classic algorithm in computational geometry used to solve for the intersection of two arbitrary polygons. Its basic principle is: traverse along the boundaries of the two polygons, switching paths at the intersection point based on the direction of entry or exit, ultimately tracing the complete boundary of the intersecting polygons. This algorithm can handle concave polygons and multiple disconnected intersection regions. The calculated overlapping region may contain one or more disconnected sub-polygons.

[0104] Then, the shoelace formula is used to calculate the projected area of ​​the overlapping region. And divide the projected area by the projected area of ​​the first spatial region. The overlap ratio rho is obtained as shown in formula (8):

[0105] rho= (8)

[0106] Finally, the numerical combination of the projected area and the overlap ratio was determined as the red line intrusion parameter. =( (rho). The redline intrusion parameters contain information in two dimensions: projected area. It reflects the actual scale of the red line intrusion, and the overlap ratio rho reflects the relative impact of the red line intrusion on the land parcel.

[0107] Assuming a target land parcel is designated for ecological use, its vector boundary encloses an area of... =15000 square meters. The prohibited construction space designated by the ecological protection red line within this area partially overlaps with the target land parcel. Using the Weiler-Atherton algorithm, two disconnected overlapping sub-polygons were calculated, with areas of 1800 square meters and 1200 square meters respectively. Therefore, the projected area is... =1800 + 1200 = 3000 square meters. Overlap ratio rho = =0.2. Therefore, the red line intrusion parameter... =(3000,0.2) indicates that 20% of the land area is covered by the ecological protection red line, with an absolute intrusion area of ​​3000 square meters.

[0108] S103. In response to the spatial replacement request for the target parcel, extract the source land and target land corresponding to the spatial replacement request, and generate a flow path from the source land to the target land, wherein the source land and the target land are both located within the target parcel.

[0109] In this embodiment, a spatial replacement request refers to a land use type conversion application initiated by planning personnel or the spatial use management platform based on planning adjustment needs, such as converting a plot of land from agricultural land to construction land, or reclaiming construction land into agricultural land; the source land refers to the land that provides indicators in the spatial replacement, and the target land refers to the land that receives indicators in the spatial replacement. Both the source land and the target land are located within the target parcel, which ensures that the spatial replacement is completed within the same parcel and avoids the management complexity caused by cross-parcel replacement.

[0110] Parse the spatial permutation request, extracting the source spatial coordinates and source land class attributes of the source land, and the target spatial coordinates and target land class attributes of the target land. The source and target spatial coordinates are the centroid coordinates of the vector boundaries of the source and target land, respectively. The source and target land class attributes are the usage classification codes of the two land parcels before the permutation. The centroid coordinates are calculated using the polygon centroid formula, for a polygon with n vertices {( , ), ( , ), ...,( , The centroid coordinates of a simple polygon formed by )} are () , The calculation formulas are shown in formulas (9) and (10):

[0111] (9)

[0112] (10)

[0113] in, The x-coordinate of the polygon's centroid; y is the ordinate of the centroid of the polygon; A is the area of ​​the polygon; , Let x and y be the x and y coordinates of the i-th vertex of the polygon; , Let be the coordinates of the (i+1)th vertex of the polygon; n is the total number of vertices on the polygon's boundary, and the summation range is i=1,2,...,n. When i=n, ​​take ( , )=( , ).

[0114] Source space coordinates With target space coordinates As path nodes, connecting line segments with vector directions are constructed between path nodes. The mathematical expression for a connecting line segment is a directed line segment L = → The vector direction refers to the direction from the source land to the target land. The formula for calculating the vector direction D is shown in formula (11):

[0115] D=( - , - (11)

[0116] in,( , ) Using the topographic center coordinates of the target The x and y components; , (Source: Topographic center coordinates) The x-axis and y-axis components.

[0117] Normalize the direction vector D to obtain the unit direction vector, as shown in formula (12):

[0118] (12)

[0119] in, The Euclidean distance between the two centroid coordinates, in unit direction vector. It represents the spatial direction of the indicator flow from the source land to the target land, and this directional information is of guiding significance for subsequent retrieval of alternative land parcels around the target land.

[0120] The characterization will be derived from the source region class attributes. Convert to target land class attribute The attribute change identifier is recorded on the connecting line segment L, thus obtaining the transfer path. The attribute change identifier includes information such as the land category code before and after the conversion, the types of indicators involved, and the indicator values.

[0121] Assuming the source land is a piece of agricultural land, its centroid coordinates are... =(500100,3045500), Land Category Attribute =“Farmland”; The target land is a piece of construction land, and its centroid coordinates are... =(500800,3045900), Land Category Attribute =“Urban construction land”. Then the vector direction is vector D = (700, 400). = = ≈806.23, unit direction vector =(0.868,0.496). The indicator type is land occupation and compensation balance (due to the conversion of agricultural land to construction land), and the indicator value is 5000 square meters. The above information is encapsulated as a transfer path R, and visualized on the planning base map as a directed line segment with an arrow pointing from the source topographic center to the target topographic center. The direction of the arrow is the direction of indicator transfer.

[0122] S104. Obtain the replacement area of ​​the spatial replacement request, perform matching verification along the transfer path according to the mapping relationship and the replacement area, combined with the land use assessment parameters, and generate the corresponding land planning scheme based on the verification result and the vector direction of the transfer path.

[0123] In this embodiment, the replacement area refers to the land area explicitly required to undergo a land use conversion in the spatial replacement request. Determining the replacement area as the replacement area means that the area control index of the source land will be reduced by this value, while the target land needs to obtain corresponding spatial plots to accommodate this increase in the index.

[0124] The specific process of matching and verification is as follows: First, it is determined whether the land use assessment parameters meet the preset compliance conditions. The preset compliance conditions refer to the compliance judgment standards set according to different land use types. For construction land, the pre-set compliance condition can be set as follows: the utilization intensity parameter should not exceed a pre-set intensity threshold. For example, a utilization intensity parameter of no more than 0.3 means that the utilized space is no less than 70% of the three-dimensional space of the target land parcel, indicating that the land parcel has reached a high level of intensive utilization and meets the prerequisite for space replacement. For agricultural land, the pre-set compliance condition can be set as follows: the occupation-compensation balance parameter should meet the dual requirements of quantity and quality balance. For example, if the area of ​​cultivated land to be occupied is 5,000 square meters and the quality grade is grade three, then the area of ​​supplementary cultivated land should not be less than 5,000 square meters and the average quality grade of the supplementary cultivated land should not be lower than grade three. That is, it simultaneously meets the two conditions that the supplementary area is greater than or equal to the occupied area and the supplementary quality grade is not lower than the occupied quality grade. For ecological land, the pre-set compliance condition can be set as follows: the overlap ratio in the red line intrusion parameter does not exceed a pre-set ratio threshold. For example, if the pre-set ratio threshold is 0.05, and the area of ​​an ecological land parcel is 10,000 square meters, and its overlap area with the ecological protection red line is 400 square meters, then the overlap ratio is... =0.04. Since 0.04 is less than the threshold of 0.05, it is determined that the preset compliance conditions are met. If the overlapping area is 800 square meters, the overlap ratio is 0.08, which exceeds the threshold of 0.05, and it is determined that the preset compliance conditions are not met. The above thresholds are only illustrative examples, and specific values ​​can be set according to local management needs.

[0125] If the land use assessment parameters meet the preset compliance conditions, and an incremental map patch with an area consistent with the replacement area is found for the spatial map patch corresponding to the target land use, then the match is considered successful.

[0126] The incremental patch retrieval process is as follows: Based on the spatial permutation request, the target boundary coordinates of the spatial patch corresponding to the target land after spatial permutation are generated; the target boundary coordinates are compared with the initial boundary coordinates of the spatial patch corresponding to the target land before spatial permutation, and a set of candidate patches generated by the boundary expansion of the spatial patch is extracted, and the expansion area of ​​each candidate patch is calculated; the first area difference between the expansion area and the replacement area is calculated; when the first area difference is less than the preset area tolerance, the corresponding candidate patch is determined as an incremental patch. The preset area tolerance is an allowable deviation range set to adapt to the actual measurement accuracy and irregularity of the land parcel; for example, the preset area tolerance can be set to 2% of the replacement area.

[0127] Assume the replacement area is =5000 square meters, area tolerance ratio coefficient =0.02, then according to the area tolerance formula: = × =100 square meters. A candidate patch was extracted using the difference operation. The outer expansion area is calculated using the shoelace formula. =4950 square meters, then according to the formula for calculating the first area difference: , =|4950-5000|=50<100, then determine The incremental patch has been successfully matched.

[0128] If an incremental land parcel is successfully retrieved, the match is considered successful. The mapping relationship is then updated based on the coordinate data of the incremental land parcel, and a land planning scheme is generated. Updating the mapping relationship includes incorporating the incremental land parcel into the mapping relationship, adjusting the relevant area control indicator values, and recording the transfer path.

[0129] Conversely, if the land use assessment parameters do not meet the preset compliance conditions, or if no incremental land parcel with a matching area is found, the matching is deemed to have failed. In this case, the corresponding land planning scheme generation process is initiated: first, based on the replacement area... Calculate the second area difference ,Right now = - ,in This represents the area of ​​the incremental land parcels that have been retrieved. Then, the search direction is extracted based on the land use category identifier of the source land. The search direction refers to the spatial search strategy direction determined based on the land use category identifier of the source land. Different land use category identifiers correspond to different spatial control rules, thus generating different search priority directions.

[0130] Specifically, a pre-defined anchor point is determined based on the land use category identifier of the source land. This pre-defined anchor point is a spatial reference point pre-set in the spatial layout, and its location is related to the spatial control rules corresponding to the land use category identifier. For example, when the source land is agricultural land, according to the spatial control rules for agricultural land—prioritizing the search for supplementary arable land resources around permanent basic farmland protection areas—the pre-defined anchor point is set as the centroid of the permanent basic farmland protection area within the designated management unit. This ensures that the supplementary arable land is spatially close to the existing farmland protection area, maintaining the concentrated and contiguous spatial pattern of farmland. When the source land is construction land, according to the spatial control rules for construction land—prioritizing the search for replaceable land around concentrated urban construction areas—the pre-defined anchor point is set as the centroid of the concentrated urban construction area within the designated management unit. This ensures the compactness of the construction land layout and the intensive use of infrastructure.

[0131] The retrieval strategy is determined based on spatial control rules for land use identifier matching. Specifically, differentiated retrieval radii and preset half-angle parameters are set according to different land use identifiers. For example, the spatial control rules for agricultural land require the replenishment of arable land in concentrated and contiguous areas, so a smaller retrieval radius and a narrower preset half-angle are used to ensure that the spare plots are spatially adjacent to existing farmland. Conversely, the spatial control rules for construction land allow for greater spatial flexibility, so a larger retrieval radius and a wider preset half-angle can be used. Based on this, the corresponding vector direction is extracted from the flow path. =( , Using the centroid of the preset anchor point as the reference point, along the vector direction Expand the preset angle range, assuming the preset angle range is a preset half angle. Construct a sector-shaped search region. The mathematical definition of a sector-shaped search region is: with the centroid coordinates 0 of a preset anchor point = ( , (with the center as the center) radius, vector direction The angle bisector extends to both sides by a pre-defined half-angle. Construct a sector-shaped region. Among them, The retrieval radius can be set according to the spatial range of the braiding unit, for example, set to half the diagonal length of the braiding unit. Preset half-angle. It can be set according to spatial proximity requirements, for example, to 30 degrees.

[0132] A point P within the sector region = ( , The criteria for determining whether something is within the search range include a first criterion corresponding to distance and a second criterion corresponding to angle, where the first criterion is shown in formula (13) and the second criterion is shown in formula (14):

[0133] (13)

[0134] (14)

[0135] in, For vectors and The angle between the two chords is calculated using the dot product formula (15):

[0136] (15)

[0137] The design of the sector-shaped retrieval area takes into account the principle of spatial proximity, that is, prioritizing the retrieval of alternative map patches near the direction of indicator flow in order to reduce spatial fragmentation of the planning scheme.

[0138] Within the sector-shaped search area, a set of unoccupied vacant land parcels is extracted. From this set, vacant land parcels that meet preset compliance conditions and whose area is not less than the second area difference are selected as reserve land parcels. Then, the reserve boundary coordinates of these reserve land parcels are obtained. A polygon merging algorithm is used to concatenate the reserve boundary coordinates with the target boundary coordinates of the spatial parcels corresponding to the target land use, resulting in the target extended boundary. Simultaneously, the original area control indicators associated with the target land use in the mapping relationship are adjusted according to the second area difference to obtain the target area control indicators. These indicators are then re-associated with the corresponding spatial parcels to form an updated mapping relationship. Finally, based on the updated mapping relationship and the target extended boundary, a planning scheme map containing the adjusted land parcel outline is output on the spatial use control platform as the land planning scheme. Furthermore, the land planning scheme map uses different colors to distinguish the original and extended areas and marks the compensation source and transfer path of the second area difference.

[0139] Assuming replacement area =5000 square meters, after incremental patch retrieval, only an incremental patch with an area of ​​3200 square meters was found, then the second area difference value =5000-3200=1800 square meters. The target is represented by the topographic center O=(500800,3045900), and the direction vector is... =(0.868,0.496), =2000 meters, =30 degrees. Within the sector-shaped search area, a free patch with an area of ​​2100 square meters was found, 850 meters from the target topographic center, satisfying the condition 2100... The area requirement is 1800. The map patch is then connected to the target land use boundary, the mapping relationship is updated, and a land planning scheme is generated.

[0140] As a preferred embodiment, this application employs differentiated control strategies for different land use types when performing matching verification. For example, for construction land, a combination of land parcel control and planning permits is used. The matching verification focuses on whether the utilization intensity parameters meet the requirements for intensive and efficient use. Verification content includes whether indicators such as legal use, plot ratio, building height, building density, and supporting facilities are within compliance ranges. For agricultural land, a combination of land category control and indicator constraints is used. The matching verification focuses on whether the land occupation and compensation balance indicators meet the bottom-line requirements for arable land protection. Verification content includes boundary control of permanent basic farmland, the equivalence of area and quality in land occupation and compensation balance, and the balance of indicators in the land increase / decrease linkage mechanism. For ecological land, a combination of zoning control and access conditions is used. The matching verification focuses on whether the red line intrusion parameters meet the rigid constraints of ecological protection. Verification content includes the prohibition of development and construction within the ecological protection red line, the restricted construction requirements in the ecological control zone, and the access conditions of the restricted and prohibited catalogs.

[0141] As another preferred embodiment, before establishing the mapping relationship between area control indicators and spatial plots, the system can solidify and verify the cadastral information of the entire region. Under a market-oriented context, the cadastral management of construction land is relatively complete, while agricultural land and ecological land, due to the large number of land categories involved and historically inconsistent registration departments, currently lack a complete rights structure system. Therefore, the system needs to connect to the real estate registration database to obtain complete cadastral information for each plot, including ownership type, right of use type, acquisition method, land use classification, and ownership unit. For urban construction land, acquisition methods include contractual ownership, administrative allocation, and negotiated transfer, and the content of rights under different acquisition methods directly determines the ease or difficulty of changing the use. For village construction land, the transfer methods for homestead use rights and collective construction land use rights differ; collective commercial construction land can directly enter the market, while homesteads require prior requisition and then supply. For agricultural land, the intensity of control varies depending on the type, with permanent basic farmland being the most strictly controlled. For ecological land, it is necessary to define land categories, plot boundaries, and corresponding ownership entities in conjunction with ecological protection red lines. The comprehensive acquisition of the aforementioned cadastral information laid a data foundation for establishing accurate mapping relationships and scientific land use assessments.

[0142] As another preferred embodiment, this application adopts the management unit as the basic unit of management in the space use control platform. The management unit refers to a space management unit used for both organizing planning and implementing land use control. Sub-units can be established within the management unit according to different planning zones. Sub-units can be divided into four categories: strictly controlled, restricted development, suitable development, and conditional development. Elements with similar or comparable control objectives are grouped together and incorporated into the same space, constructing a "sub-unit—block—plot" control model. Within each sub-unit, population and building scale forecasts, public service facilities pre-control, and other demand forecasts and supply are conducted using the block as the calculation unit, ultimately refining the control indicators for land use development or protection for each plot. This hierarchical management approach achieves effective transmission from macro-control objectives to micro-plot control. The classification and control focus of each sub-unit are shown in Table 2.

[0143] Table 2 Sub-unit Classification and Key Control Points

[0144] Figure 3 This is a schematic diagram illustrating a specific implementation of a land resource planning system based on cadastral information, as provided in this application embodiment. (Refer to...) Figure 3 The system may include:

[0145] Module 31 is established to obtain spatial patches of the target land parcel, extract the land type identifier and area control index of the spatial patches, and establish a mapping relationship between the area control index and the corresponding spatial patches.

[0146] The determining module 32 is used to determine the corresponding land use assessment parameters based on the land type identifier, wherein when the land type identifier is construction land, the land use assessment parameter is the utilization intensity parameter; when the land type identifier is agricultural land, the land use assessment parameter is the land occupation and compensation balance parameter; and when the land type identifier is ecological land, the land use assessment parameter is the red line intrusion parameter.

[0147] Extraction module 33 is used to respond to a spatial replacement request for the target parcel, extract the source land and target land corresponding to the spatial replacement request, and generate a flow path from the source land to the target land, wherein the source land and the target land are both located within the target parcel.

[0148] The generation module 34 is used to obtain the replacement area of ​​the spatial replacement request, perform matching verification along the transfer path according to the mapping relationship and the replacement area, combined with the land use assessment parameters, and generate a corresponding land planning scheme based on the verification result and the vector direction of the transfer path.

[0149] The land resource planning system based on cadastral information in this application is used to implement the aforementioned land resource planning method based on cadastral information. Therefore, the specific implementation of the land resource planning system based on cadastral information can be found in the embodiment section of the land resource planning method based on cadastral information above. The specific implementation can be referred to the description of the corresponding embodiments, which will not be repeated here.

[0150] This application also provides an electronic device, comprising: a memory for storing a computer program; and a processor for executing the computer program to implement the steps of any of the above-described land resource planning methods based on cadastral information.

[0151] This application also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the steps of any of the above-described land resource planning methods based on cadastral information.

[0152] In one exemplary embodiment, the aforementioned computer-readable storage medium may include, but is not limited to, various media capable of storing computer programs, such as USB flash drives, read-only memory, random access memory, portable hard drives, magnetic disks, or optical disks.

[0153] The embodiments of this application also provide a computer program product, which includes a computer program that, when executed by a processor, implements the steps in any of the above embodiments of the land resource planning method based on cadastral information.

[0154] Those skilled in the art will further recognize that the units and algorithm steps of the various examples described in conjunction 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 and steps 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 implementation should not be considered beyond the scope of this application.

[0155] The foregoing has provided a detailed description of a land resource planning method and system based on cadastral information provided in this application. Specific examples have been used to illustrate the principles and implementation methods of this application. The descriptions of the embodiments above are merely for the purpose of helping to understand the method and its core ideas. It should be noted that those skilled in the art can make various improvements and modifications to this application without departing from its principles, and these improvements and modifications also fall within the protection scope of this application.

Claims

1. A land resource planning method based on cadastral information, applied to a spatial use control platform, characterized in that, include: Obtain spatial map features of the target land parcel, extract the land use identifiers and area control indicators of the spatial map features, and establish a mapping relationship between the area control indicators and the corresponding spatial map features; Based on the land use classification, corresponding land use assessment parameters are determined. Specifically, when the land use classification is construction land, the land use assessment parameter is a utilization intensity parameter, which represents the developable intensity of the target land parcel in three-dimensional space, and the sum of the utilization intensity parameter and the utilization rate of the target land parcel in three-dimensional space is 1. When the land use classification is agricultural land, the land use assessment parameter is an occupation-compensation balance parameter, which includes the area and quality grade of the occupied cultivated land during the spatial replacement process, as well as the area and quality grade of the supplemented cultivated land. When the land use classification is ecological land, the land use assessment parameter is a red line intrusion parameter, which includes the projected area of ​​the first spatial region within the vector boundary of the target patch, and the overlap ratio corresponding to the overlapping area between the second spatial region defined by the ecological protection red line and the first spatial region. In response to a spatial replacement request for the target parcel, the source land and target land corresponding to the spatial replacement request are extracted, and a flow path from the source land to the target land is generated, wherein both the source land and the target land are located within the target parcel; Obtain the replacement area of ​​the spatial replacement request, perform matching verification along the transfer path based on the mapping relationship and the replacement area, combined with the land use assessment parameters, and generate the corresponding land planning scheme based on the verification results and the vector direction of the transfer path.

2. The method according to claim 1, characterized in that, Along the transfer path, based on the mapping relationship and the replacement area, and in conjunction with the land use assessment parameters, a matching verification is performed, and a corresponding land planning scheme is generated based on the verification results and the vector direction of the transfer path, including: If the land use assessment parameters meet the corresponding preset compliance conditions, and an incremental map patch with an area consistent with the replacement area is retrieved for the spatial map patch corresponding to the target land use, then the matching is determined to be successful; The mapping relationship is updated based on the coordinate data of the incremental map features to generate a land planning scheme; Alternatively, if the land use assessment parameters do not meet the preset compliance conditions, or if the incremental map patch is not found, the matching is deemed to have failed. In the event of a failed match, a search guide is extracted based on the land use category identifier of the source land, and alternative land parcels that meet the preset compliance conditions are searched around the preset anchor point along the search guide. The spare land parcels are compensated to the spatial parcels corresponding to the target land use, and the mapping relationship is updated to generate a land planning scheme.

3. The method according to claim 2, characterized in that, For the spatial patch corresponding to the target land use, incremental patches with areas consistent with the replacement area are retrieved, including: Based on the spatial permutation request, generate the target boundary coordinates of the spatial patch corresponding to the target land after spatial permutation; By comparing the target boundary coordinates with the initial boundary coordinates of the spatial patch corresponding to the target land before spatial permutation, a set of candidate patches generated by the boundary expansion of the spatial patch is extracted, and the expansion area of ​​each candidate patch is calculated. Calculate the first area difference between the expanded area and the replacement area; When the first area difference is less than the preset area tolerance, the corresponding candidate patch is determined as an incremental patch.

4. The method according to claim 2, characterized in that, Based on the land use category identifier of the source land, a search guide is extracted, and along the search guide, alternative land parcels that meet the preset compliance conditions are searched around preset anchor points, including: The corresponding preset anchor point is determined based on the land use type identifier of the source land, and the vector direction of the flow path is generated based on the spatial control rules matched by the land use type identifier. Using the centroid of the preset anchor point as the reference point, a fan-shaped retrieval area is constructed by expanding the preset angle range along the vector direction; Extract the set of unoccupied free patches within the sector-shaped search area; From the set of idle patches, select idle patches that meet the preset compliance conditions and whose area is not less than the second area difference, and use them as backup patches. The second area difference is the difference between the replacement area and the area of ​​the retrieved incremental patches.

5. The method according to claim 2, characterized in that, The process involves compensating the spare land parcels to the spatial parcels corresponding to the target land use, updating the mapping relationship, and generating a land planning scheme, including: Obtain the backup boundary coordinates of the backup patch, and concatenate the backup boundary coordinates with the target boundary coordinates to obtain the target extended boundary; Adjust the original area control indicators associated with the target land in the mapping relationship according to the second area difference to obtain the target area control indicators, and re-associate the target area control indicators with the corresponding spatial patches to update the mapping relationship; Based on the updated mapping relationship and the target extended boundary, a planning scheme map containing the adjusted plot outline is output on the spatial use control platform as a land planning scheme.

6. The method according to claim 1, characterized in that, The area control indicators include the land occupation and compensation balance indicator, the increase and decrease linkage indicator, or the ecological red line indicator; Establishing the mapping relationship between the area control indicators and the corresponding spatial patches includes: When the land type is identified as agricultural land, agricultural land patches with agricultural land attributes and whose area matches the occupation-compensation balance index are selected from the spatial patches, and the occupation-compensation balance index is associated with the boundary node of the agricultural land patch to establish a first mapping association. When the land type is identified as construction land, construction land parcels with construction land attributes and whose areas match the increase-decrease linkage index are selected from the spatial parcels. The increase-decrease linkage index is then associated with the boundary nodes of the construction land parcels to establish a second mapping association. When the land use type is identified as ecological land, ecological land parcels with ecological land use attributes and whose area matches the ecological red line index are selected from the spatial parcels. The ecological red line index is then associated with the boundary nodes of the ecological land parcels to establish a third mapping association. Merge the first mapping association, the second mapping association, and the third mapping association to obtain the mapping relationship.

7. The method according to claim 1, characterized in that, When the land type is identified as construction land, the utilization intensity parameters are determined, including: Extract the three-dimensional vertex coordinates of the planned building model, and calculate the first volume of the planned building model based on the three-dimensional vertex coordinates; Obtain preset three-dimensional control lines and calculate the second volume inside the three-dimensional constraint boundary formed by the three-dimensional control lines; Calculate the volume difference between the second volume and the first volume, and determine the ratio of the volume difference to the second volume as the strength parameter.

8. The method according to claim 1, characterized in that, When the land type is identified as ecological land, the parameters for determining the red line intrusion are as follows: Extract the first spatial region within the vector boundary of the target patch, and the second spatial region defined by the ecological protection red line; Extract the overlapping area between the first spatial region and the second spatial region; Calculate the projected area of ​​the overlapping region and divide the projected area by the projected area of ​​the first spatial region to obtain the overlap ratio; The numerical combination of the projected area and the overlap ratio is determined as the red line intrusion parameter.

9. The method according to claim 1, characterized in that, Extract the source land and target land corresponding to the spatial replacement request, and generate a transfer path from the source land to the target land, including: The spatial permutation request is parsed to extract the source spatial coordinates and source land class attributes of the source land, as well as the target spatial coordinates and target land class attributes of the target land. The source space coordinates and the target space coordinates are used as path nodes, and connecting line segments with vector directions are constructed between the path nodes; The attribute change identifier, which represents the transformation from the source land type attribute to the target land type attribute, is recorded on the connecting line segment to obtain the flow path.

10. A land resource planning system based on cadastral information, applied to a spatial use control platform, for executing a land resource planning method based on cadastral information as described in any one of claims 1 to 9, characterized in that, include: A module is established to acquire spatial patches of the target land parcel, extract the land type identifier and area control index of the spatial patches, and establish a mapping relationship between the area control index and the corresponding spatial patches; The determination module is used to determine the corresponding land use assessment parameters based on the land use identifier, wherein when the land use identifier is construction land, the land use assessment parameter is the utilization intensity parameter; when the land use identifier is agricultural land, the land use assessment parameter is the land occupation and compensation balance parameter; and when the land use identifier is ecological land, the land use assessment parameter is the red line intrusion parameter. The extraction module is used to respond to a spatial replacement request for the target parcel, extract the source land and target land corresponding to the spatial replacement request, and generate a flow path from the source land to the target land, wherein both the source land and the target land are located within the target parcel; The generation module is used to obtain the replacement area of ​​the spatial replacement request, perform matching verification along the transfer path according to the mapping relationship and the replacement area, combined with the land use assessment parameters, and generate a corresponding land planning scheme based on the verification result and the vector direction of the transfer path.

Images