A method for updating the urban entity area based on the ArcGIS platform

Abstract

The present invention provides a method for updating the urban entity area based on the ArcGIS platform. The land use patch layer is obtained according to the survey data, the initial urban area is extracted based on the urban-village attribute code, a buffer zone with a preset distance is set to obtain the expanded space, and the land object patches that meet the conditions in the intersecting area between the expanded space and the land use patch layer are screened to obtain the first patches to be included, which are added to the urban entity area of the previous iteration to obtain the current urban entity area. The buffer zone setting to patch addition steps are iteratively performed until the preset number of iterations is reached to obtain the initial urban entity area, which is compared with the urban entity area of the previous cycle to obtain the second patches to be included, and added to the initial urban entity area; the linear features are trimmed, it is judged whether the added patches meet the inclusion rules, and boundary verification is carried out to obtain the target urban entity area. The present invention can achieve the standardized and rapid update of the urban entity area, improving the update efficiency and accuracy.

Description

Technical Field

[0001] The present invention relates to the field of geographic information technology, and in particular to a method for updating a physical area of ​​an urban area based on an ArcGIS platform. Background Art

[0002] With the rapid development of my country's economy and society and the in-depth promotion of new urbanization, the urban spatial structure continues to evolve, and the dynamics of urban physical areas are becoming increasingly prominent. The precise definition and timely update of physical areas are not only crucial for the scientific formulation and effective implementation of national land space planning, but also directly affect the work in many fields such as land use management, urban construction, environmental protection, and public service resource allocation. However, the existing methods and technologies for updating urban physical areas face a series of severe challenges, which restrict their application effect and efficiency in practice.

[0003] The existing technical problems include: first, data lag and difficulty in updating. The update of urban physical areas usually requires the aggregation of data from multiple data sources, but the update cycles of multi-source data vary greatly, and the acquisition and update of some key data may be delayed. In addition, the imperfect data sharing mechanism leads to the lag of information update and the complexity of integration. The lag of data makes it difficult for planning managers to capture the actual changes of physical areas in a timely manner, which in turn affects the timeliness and accuracy of decision-making, while the difficulty of data integration increases the labor cost and reduces the efficiency of updating. Second, inconsistent implementation of standards and cumbersome operations. Due to individual differences in understanding and implementing standards, limitations of technical means and complexity of business processes in the actual implementation process, inconsistent implementation of standards is prone to occur. In addition, traditional urban physical area update methods are mostly based on a large number of manual operations, which are time-consuming and labor-intensive, and are easily affected by human factors, resulting in the difficulty in ensuring the accuracy and consistency of the update results. Third, the update speed conflicts with the needs of planning management. Modern planning management should have the ability to monitor the changes of urban physical areas in real time in order to adjust planning schemes, optimize resource allocation and respond to emergencies in a timely manner. However, current urban physical area renewal technologies generally have problems such as slow update speed and untimely response, which cannot meet the needs of refined and dynamic planning management. This lag makes planning decisions often based on outdated information, which not only affects the timeliness of decision-making, but also may lead to adverse consequences such as resource misallocation and policy failure.

[0004] Therefore, there is an urgent need for a method for updating urban physical areas that can achieve standardized and rapid updates of urban physical areas and improve the update speed and accuracy. Summary of the invention

[0005] Based on this, it is necessary to provide a method for updating urban entity regions based on the ArcGIS platform for the above-mentioned technical problems, so as to achieve the standardized and rapid update of urban entity regions, reduce manual operations, and improve work efficiency.

[0006] A method for updating urban entity regions based on the ArcGIS platform includes the following steps:

[0007] S1. Obtain a land use patch layer according to the survey data, and extract the patches in the land use patch layer according to the urban-village attribute code to obtain the initial urban range;

[0008] S2. Set a buffer zone with a preset distance based on the initial urban range to obtain an expanded space;

[0009] S3. Screen the land use patches in the intersection area of the land use patch layer and the expanded space according to the land use patch attributes to obtain the first patches to be included;

[0010] S4. Add the first patches to be included to the initial urban range or the urban entity region of the previous iteration to obtain the current urban entity region. Iteratively perform steps S2 to S4 until the preset number of iterations is reached to obtain the initial urban entity region;

[0011] S5. Obtain the urban entity region of the previous cycle, compare it with the initial urban entity region to obtain the second patches to be included, and add them to the initial urban entity region;

[0012] S6. Trim the linear features, determine whether the first patches to be included and the second patches to be included meet the inclusion rules, and conduct boundary verification to obtain the target urban entity region.

[0013] In one embodiment, step S1 includes: finding the survey data according to the national land use change survey database and extracting the land use patch layer in the survey data; screening the urban patches in the land use patch layer according to the urban-village attribute code to obtain the initial urban range; fusing all the screened urban patches to obtain the first layer.

[0014] In one embodiment, step S2 includes: obtaining the preset distance, generating a buffer zone in the first layer based on the preset distance to obtain a buffer zone layer; performing an intersection-inverse operation on the buffer zone layer and the first layer to obtain the expanded space.

[0015] In one embodiment, step S3 includes: based on the attributes of the land feature patches, screening out the patches to be confirmed in the land type patch layer that intersect with the expanded space; removing the patches that have been included in the initial urban range from the patches to be confirmed to obtain the first patches to be included, and the first patches to be included carry the marked value of the current iteration number.

[0016] In one embodiment, step S4 includes: at the initial iteration, using the append tool to add the first patches to be included to the initial urban range to obtain the urban entity area; at non-initial iterations, using the append tool to add the first patches to be included to the urban entity area obtained in the previous iteration to obtain the current urban entity area; repeating the operations of buffer setting, obtaining the first patches to be included, and adding patches until the preset iteration number is reached to obtain the initial urban entity area.

[0017] In one embodiment, step S5 includes: importing the urban entity area of the previous cycle, and screening out the identical patches that are exactly the same as the initial urban entity area; based on the identical patches, performing inverse selection in the urban entity area of the previous cycle to obtain the different patches that are different from the initial urban entity area, denoted as the third layer; setting the target layer as the land type patch layer and the source layer as the third layer, and using the spatial selection method including the source layer elements to screen out the patches in the land type patch layer that match the boundary of the third layer according to the land type patch attributes to obtain the fourth layer; performing field addition operations on the third layer and the fourth layer respectively, calculating the area and the longitude and latitude of the centroid of each patch, and adding a unique identifier to the fourth layer; calculating the overlapping part of the third layer and the fourth layer, and fusing to obtain the layer of interest, adding fields to the layer of interest for obtaining the longitude and latitude of the centroid of the patch and calculating the patch area; matching the layer of interest with the fourth layer to obtain the corresponding patch groups in the two layers; calculating the centroid distance between the two patches in the patch group, and the area proportion of the patch in the layer of interest in the corresponding patch in the fourth layer and storing them; screening out the patches in the fourth layer with the area proportion greater than the preset threshold, the centroid distance less than the preset distance threshold, and the land type name being the same as the corresponding patch in the third layer to obtain the second patches to be included, and adding them to the initial urban entity area.

[0018] In one embodiment, step S5 includes: importing the urban entity area of the previous cycle, and screening out the same patches that are exactly the same as the initial urban entity area; based on the same patches, performing an inverse selection in the urban entity area of the previous cycle to obtain different patches that are different from the initial urban entity area, denoted as the third layer; setting the target layer as the land type patch layer and the source layer as the third layer, and using the spatial selection method within the range of the source layer features to screen out the patches that match the boundary of the third layer in the land type patch layer according to the land type patch attributes to obtain the fifth layer; performing field addition operations on the third layer and the fifth layer respectively, calculating the area and the longitude and latitude of the centroid of each patch, and adding a unique identifier to the fifth layer; calculating the overlapping part of the third layer and the fifth layer, and fusing to obtain the layer of interest, adding fields to the layer of interest for obtaining the longitude and latitude of the centroid of the patch and calculating the area of the patch; matching the layer of interest with the fifth layer to obtain a corresponding patch group in the two layers; calculating the centroid distance between the two patches in the patch group, and the area ratio of the patch in the layer of interest in the corresponding patch in the fifth layer and storing it; screening out the patches in the third layer where the area ratio is greater than the preset threshold, the centroid distance is less than the preset distance threshold, and the land type name is the same as the corresponding patch in the fifth layer to obtain the second patches to be included, and adding them to the initial urban entity area.

[0019] In one embodiment, the calculation formula for the longitude and latitude of the centroid is:

[0020]

[0021]

[0022] In the formula, D represents the area where the patch is located, ∫∫Ddxdy and ∫∫Dy dx dy are the double integrals of calculating the x coordinate and y coordinate of the patch area D respectively, representing the first moment of the patch about the y-axis and x-axis; and the denominator ∫∫Ddxdy is the area of the patch.

[0023] In one embodiment, step S6 includes: using a clipping tool to trim the linear features extending outside the concentrated contiguous area; determining whether the first patches to be included and the second patches to be included meet the inclusion rules, discarding them if they do not meet the rules, and retaining them if they meet the rules; comparing the urban entity area with the control line, marking the patches intersecting the control line, and performing boundary verification, and discarding or retaining the marked patches according to the boundary verification results to output the target urban entity area.

[0024] Compared with the prior art, the advantages and beneficial effects of the present invention are as follows: obtaining a land type patch layer based on survey data, extracting the patches therein based on the urban and rural village attribute codes to obtain the initial urban area; setting a buffer zone with a preset distance based on the initial urban area to obtain an expanded space, screening out the land type patches in the land type patch layer that intersect with the expanded space according to the land object patch attributes to obtain the first patches to be incorporated; adding the first patches to be incorporated to the initial urban area or the urban entity area of the previous iteration to obtain the current urban entity area, and iteratively performing the steps of buffer zone setting to patch addition until the preset number of iterations is reached to obtain the initial urban entity area; obtaining the urban entity area of the previous cycle, comparing it with the initial urban entity area to obtain the second patches to be incorporated, and adding them to the initial urban entity area; trimming linear features, determining whether the first patches to be incorporated and the second patches to be incorporated meet the incorporation rules, and performing boundary verification to obtain the target urban entity area. By integrating the powerful functions of the ArcGIS platform and the flexibility of Python programming, the rapid update of the standardization of the urban entity area is realized, reducing the manual patch selection operation, effectively solving the problems of data lag, inconsistent standard execution, cumbersome operation, and mismatch between the update speed and the planning management requirements, providing a real-time and accurate urban entity area update solution for the field of national land space planning, and strongly supporting the refined and dynamic planning management requirements under the background of new urbanization. BRIEF DESCRIPTION OF THE DRAWINGS

[0025] Figure 1 It is a schematic flowchart of a method for updating the urban entity area based on the ArcGIS platform in an embodiment;

[0026] Figure 2 It is a flowchart of a method for updating the urban entity area based on the ArcGIS platform in an embodiment;

[0027] Figure 3 It is a distribution map of the initial urban area of a certain city in an embodiment;

[0028] Figure 4 It is the patches participating in the iteration in the initial urban area of a certain city in an embodiment;

[0029] Figure 5 It is the fusion result of the patches participating in the iteration in the initial urban area of a certain city in an embodiment;

[0030] Figure 6 It is the buffer zone generated from the fusion result of the patches participating in the iteration in the initial urban area of a certain city in an embodiment;

[0031] Figure 7 It is the expanded space during the iteration process of determining the urban area of a certain city in an embodiment;

[0032] Figure 8 Schematic diagram of the first patch to be incorporated in an embodiment;

[0033] Figure 9 Result after iteration in an embodiment;

[0034] Figure 10 Second patch to be incorporated determined with reference to the data of the previous cycle in an embodiment;

[0035] Figure 11 Result of the target urban entity area after correction and verification in an embodiment. Detailed implementation manners

[0036] Before describing the detailed implementation manners of the present invention, the overall concept of the present invention is described as follows:

[0037] The present invention is mainly developed for the update process of urban entity areas. At present, there are problems such as data lag, difficult update, inconsistent standard execution, and slow update speed in the update of urban entity areas.

[0038] Therefore, the present invention proposes a method for updating urban entity areas based on the ArcGIS platform. The land use patch layer is obtained according to the survey data, and the patches therein are extracted based on the urban-village attribute code to obtain the initial urban range; a buffer zone with a preset distance is set based on the initial urban range to obtain an expanded space, and the land use patches in the intersecting area between the expanded space and the land use patch layer are screened according to the land use patch attributes to obtain the first patch to be incorporated; the first patch to be incorporated is added to the initial urban range or the urban entity area of the previous iteration to obtain the current urban entity area, and the steps of buffer zone setting to patch addition are iteratively performed until the preset number of iterations is reached to obtain the initial urban entity area; the urban entity area of the previous cycle is obtained, compared with the initial urban entity area to obtain the second patch to be incorporated, and added to the initial urban entity area; the linear features are trimmed, it is judged whether the first patch to be incorporated and the second patch to be incorporated meet the incorporation rules, and boundary verification is performed to obtain the target urban entity area. By integrating the ArcGIS platform and Python programming, the real-time automatic update of the urban entity area is realized, meeting the requirements of unified standards and fast update speed, greatly improving the update efficiency, so as to adapt to the needs of urban planning and management.

[0039] After introducing the overall concept of the present invention, in order to make the purpose, technical solution and advantages of the present invention clearer, the present invention will be further described in detail below through specific implementation manners in conjunction with the drawings. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention.

[0040] In one embodiment, as Figure 1 and Figure 2As shown in the figure, a method for updating the urban entity area based on the ArcGIS platform is provided, including the following steps:

[0041] Step S101: Obtain the land use patch layer according to the survey data, and extract the patches in the land use patch layer according to the urban-village attribute code to obtain the initial urban area.

[0042] Specifically, in the ArcGIS platform, through the survey data of the urban area, the land use patch layer is extracted, and the corresponding patches are extracted from the land use patch layer according to the urban-village attribute code. Based on the extracted patches, the initial urban area is obtained. The ArcGIS platform has a powerful data management function, which can achieve seamless docking with various geographical data sources and realize the unified storage, organization, query, and maintenance of data.

[0043] Furthermore, by combining the ArcGIS platform and Python programming, an automated data update script is developed to grab the latest data from each data source regularly or on demand. Through preset data cleaning and conversion rules, the consistency and integrity of the data format are ensured. Through integrated and automated data management, the speed and accuracy of data update are greatly improved, the manual processing links are reduced, and the timeliness of urban entity area information is ensured.

[0044] Python programming can call the rich API interfaces provided by the ArcGIS platform to convert the rules and standards in the urban area determination procedure into executable algorithm codes. By writing Python scripts, a full-chain automated processing from data preparation, boundary recognition, feature extraction to result verification is realized, ensuring that each step strictly follows the industry standard of the "Urban Area Determination Procedure" TD / T 1064-2021. Reducing the deviation of standard execution caused by human factors, the automated update process can not only greatly improve work efficiency and reduce human errors, but also ensure the standardization and consistency of the update results.

[0045] Among them, step S101 includes: finding the survey data according to the national land change survey database and extracting the land use patch layer in the survey data; screening the urban patches in the land use patch layer according to the urban-village attribute code to obtain the initial urban area; fusing all the screened urban patches to obtain the first layer.

[0046] Specifically, in the ArcGIS platform, connect to the national land change survey database through the data source interface, obtain the survey data based on the national land change survey database, and extract the land use patch layer in the survey data through the ArcGIS platform. This layer contains at least fields such as land use name, land use code, patch area, and urban-village attribute code; screen the urban patches in the land use patch layer according to the urban-village attribute code (for example, 201, 201A) to obtain the urban patches and use them as the initial urban area.Figure 3 As shown, determine the urban patches participating in the iteration within the initial urban area range, such as Figure 4 As shown, fuse all the obtained urban patches participating in the iteration to obtain the first layer, such as Figure 5 As shown.

[0047] Step S102: Set a buffer zone with a preset distance based on the initial urban area range to obtain an expanded space.

[0048] Specifically, set the preset distance in the ArcGIS platform or set the default preset distance through Python programming. For example, set the preset distance to 100m according to the "Regulations for Determining Urban Area Range". Since the range of a town usually expands on the basis of the previous year, after obtaining the initial urban area range, obtain the preset distance, set a buffer zone with the preset distance outside the initial urban area range to obtain a fused buffer zone, and based on the initial urban area range and the buffer zone, obtain the expanded space of the initial urban area range, so that the obtained expanded space is more in line with the actual expansion situation of the urban area.

[0049] Among them, step S102 includes: obtaining the preset distance, generating a buffer zone in the first layer based on the preset distance to obtain a buffer zone layer; performing an intersection negation operation on the buffer zone layer and the first layer to obtain the expanded space.

[0050] Specifically, obtain the buffer zone preset distance in the system, generate a buffer zone layer after fusing the preset distance in the first layer according to the preset distance, such as Figure 6 As shown, the buffer distance with the preset distance set to 100m; use the intersection negation tool to obtain the non-intersecting part of the buffer zone layer and the first layer to obtain the expanded space; or use the erase tool to erase the buffer zone layer with the initial urban area range to obtain an expanded ring-shaped area, which is the expanded space, such as Figure 7 As shown.

[0051] Step S103: Screen the feature patches in the intersection area of the feature patch layer and the expanded space according to the feature patch attributes to obtain the first patches to be included.

[0052] Specifically, based on the ArcGIS platform, judge the feature patch attributes of the feature patches in the intersection area of the feature patch layer and the expanded space, screen out the feature patches that meet the requirements in the feature patch layer to obtain the first patches to be included, which are used to be added to the initial urban area range to improve the accuracy of the urban entity area.

[0053] Step S103 includes: screening out the patches to be confirmed in the intersection area with the expanded space in the feature patch layer based on the feature patch attributes; removing the patches that have been included in the initial urban area range from the patches to be confirmed to obtain the first patches to be included, and the first patches to be included carry the marked value of the current iteration number.

[0054] Specifically, based on the attributes of the feature patches, a spatial selection method is adopted to screen the patches to be confirmed in the land use patch layer that intersect with the expanded space, and remove the feature patches that have been included in the initial urban area in the previous steps to obtain the first patches to be included. As shown in Figure 8 the figure; in order to obtain more accurate first patches to be included, multiple iterative screenings can be performed. The iterative process follows industry standards, and the current iteration number is marked in the specified field. The iteration number is updated and recorded by carrying the marked value through the first patches to be included.

[0055] Step S104: Add the first patches to be included to the initial urban area or the urban entity area of the previous iteration to obtain the current urban entity area. Iteratively perform steps S102 to S104 until the preset iteration number is reached to obtain the initial urban entity area.

[0056] Specifically, in the ArcGIS platform, add the first patches to be included to the layer iterated in the previous step. If the marked value carried in the first patches to be included is detected as 1, indicating the first iteration, add the first patches to be included to the initial urban area to obtain the current urban entity area; according to the preset iteration number of the ArcGIS platform or Python programming, repeat the operations of steps S102 to S104 to obtain the initial urban entity area. As shown in Figure 9 the figure.

[0057] Among them, step S104 includes: in the initial iteration, use the append tool to add the first patches to be included to the initial urban area to obtain the urban entity area; in non-initial iterations, use the append tool to add the first patches to be included to the urban entity area obtained in the previous iteration to obtain the current urban entity area; repeat the operations of buffer setting, obtaining the first patches to be included, and adding patches until the preset iteration number is reached to obtain the initial urban entity area.

[0058] Specifically, detect the current iteration number according to the marked value carried by the first patches to be included. If the marked value is 1, it indicates the initial iteration. Use the append tool in the ArcGIS platform to add the first patches to be included to the initial urban area to obtain the corresponding urban entity area; if the marked value is any value greater than 1, it indicates a non-initial iteration. Use the append tool in the ArcGIS platform to add the first patches to be included to the urban entity area obtained in the previous iteration to obtain the current urban entity area. Repeat the operations of iterative steps S102 to S104 until the preset iteration number is reached to obtain the initial urban entity area, thereby realizing the accurate acquisition of the patches of the urban entity area.

[0059] Step S105: Obtain the urban entity area of the previous cycle, compare it with the initial urban entity area to obtain the second to-be-incorporated patches, and add them to the initial urban entity area.

[0060] Specifically, in order to avoid patch omission in the above steps, resulting in inaccurate urban entity areas obtained finally, the data interface in the ArcGIS platform can be used to obtain the urban entity area of the previous cycle, such as the urban entity area results of a certain city in the previous year, and compare it with the obtained initial urban entity area to screen out the second to-be-incorporated patches, as Figure 10 shown, and add the second to-be-incorporated patches to the initial urban entity area, so as to obtain a more accurate urban entity area. In addition, directly using the urban entity area of the previous cycle can simplify the process of obtaining the urban entity area, improving work efficiency while ensuring the accuracy of the obtained area.

[0061] Among them, step S105 includes: Import the urban entity area of the previous cycle, and screen out the identical patches that are exactly the same as the initial urban entity area; Based on the identical patches, perform an inverse selection in the urban entity area of the previous cycle to obtain the different patches that are different from the initial urban entity area, denoted as the third layer; Set the target layer as the land use patch layer, and the source layer as the third layer. Use the spatial selection method that includes the elements of the source layer to screen out the patches that match the boundary of the third layer in the land use patch layer according to the land use patch attributes to obtain the fourth layer; Perform field addition operations on the third layer and the fourth layer respectively, calculate the area and centroid longitude and latitude of each patch, and add a unique identifier to the fourth layer; Calculate the overlapping part of the third layer and the fourth layer, and fuse to obtain the layer of interest. Add fields to the layer of interest to obtain the centroid longitude and latitude of the patches and calculate the patch area; Match the layer of interest with the fourth layer to obtain the corresponding patch groups in the two layers; Calculate the centroid distance between the two patches in the patch group, as well as the area proportion of the patch in the layer of interest in the corresponding patch in the fourth layer and store it; Screen out the patches in the fourth layer whose area proportion is greater than the preset threshold, the centroid distance is less than the preset distance threshold, and the land use name is the same as the corresponding patch in the third layer to obtain the second to-be-incorporated patches, and add them to the initial urban entity area.

[0062] Specifically, as Figure 2As shown in the figure, when comparing the urban entity area of the previous cycle with the initial urban entity area, based on the ArcGIS platform, import the results of the urban entity area of the previous year, denoted as stdy_old, and import the urban entity area completed in this year's iteration, denoted as stdy_new; use the Select by Location tool to screen out the patches that are exactly the same in stdy_old and stdy_new; deselect in stdy_old, that is, select the patches that are inconsistent with stdy_new, denoted as the layer patches, which is the third layer; use the Select by Location tool, set the target layer to the DLTB layer, the source layer to the patches layer, and the spatial selection method of the target layer features to "Contain source layer features". At this time, it means that the area of the patches in the land use type patches layer is greater than or equal to the area of the corresponding patches in the urban entity area of the previous cycle. The selected patches are named DLTB_include to obtain the fourth layer.

[0063] Add fields to the patches layer and the DLTB_include layer respectively to calculate the area of each feature (the area is denoted as the area_include field) and the longitude and latitude of the centroid; add the uni_code_inlude field to the DLTB_include layer and assign a unique non-repeating marker value; use the Intersect tool to calculate the overlapping part of the patches layer and the DLTB_include layer; use the Dissolve tool, set the dissolve field to the uni_code_inlude field, and the resulting layer is denoted as the intersect layer. Add fields to this layer to obtain the longitude and latitude of the centroid of the patches and add the area_intersect field to calculate the feature area.

[0064] Match the uni_code_inlude field in the intersect layer with the uni_code_inlude field in the DLTB_include layer to obtain the corresponding patches in the original intersect layer and the DLTB_include layer one by one; calculate area_intersect / area_include to obtain the area ratio, and save the calculation result to the specified field in the DLTB_include layer; calculate the centroid distance of the corresponding patches in the intersect layer and the DLTB_include layer, and save the calculation result to the specified field in the DLTB_include layer; according to the relationship between the preset threshold and the area ratio, the relationship between the centroid distance and the preset distance threshold, and the land use type name, screen out the patches in the DLTB_include layer where the area ratio is greater than the preset threshold, the centroid distance is less than the preset distance threshold, and the land use type name is the same as the corresponding patches in the third layer to obtain the second to-be-incorporated patches a that meet the requirements and incorporate them into the initial urban entity area.

[0065] Among them, step S5 includes: importing the urban entity area of the previous cycle, and screening out the same patches that are exactly the same as the initial urban entity area; performing an inverse selection in the urban entity area of the previous cycle based on the same patches to obtain different patches that are different from the initial urban entity area, denoted as the third layer; setting the target layer as the land use type patch layer, setting the source layer as the third layer, and using the spatial selection method within the range of the source layer features to screen out the patches that match the boundary of the third layer in the land use type patch layer according to the land use type patch attributes to obtain the fifth layer; performing field addition operations on the third layer and the fifth layer respectively, calculating the area and the longitude and latitude of the centroid of each patch, and adding a unique identifier to the fifth layer; calculating the overlapping part of the third layer and the fifth layer, and fusing to obtain the layer of interest, adding fields to the layer of interest to obtain the longitude and latitude of the centroid of the patch and calculate the area of the patch; matching the layer of interest with the fifth layer to obtain a corresponding patch group in the two layers; calculating the centroid distance between the two patches in the patch group, and the area ratio of the patch in the layer of interest in the corresponding patch in the fifth layer and storing it; screening out the patches in the third layer with an area ratio greater than the preset threshold, a centroid distance less than the preset distance threshold, and the land use type name being the same as the corresponding patch in the fifth layer to obtain the second patches to be included, and adding them to the initial urban entity area.

[0066] Specifically, using the Select by Location tool, setting the target layer as the DLTB layer, the source layer as the patches layer, and setting the spatial selection method of the target layer features as "within the range of the source layer features". At this time, it means that the area of the patch representing the land use type patch layer is smaller than the area of the corresponding patch in the urban entity area of the previous cycle, and the screened patches are named DLTB_included to obtain the fifth layer.

[0067] Since the area relationship between the urban entity areas of this cycle and the previous cycle has changed, it is necessary to replace the DLTB_include layer in the foregoing steps with the patches layer, replace the patches layer with the DLTB_included layer, repeat the foregoing steps of field addition, obtaining the layer of interest, patch group matching, area and centroid distance calculation, and patch screening operations, and record the screened result as the second patches to be included b, and include them in the initial urban entity area to achieve accurate acquisition of the patches of the urban entity area.

[0068] Among them, the calculation formula for the longitude and latitude of the centroid is:

[0069]

[0070]

[0071] In the formula, D represents the area where the patch is located. ∫∫Ddxdy and ∫∫Dy dx dy are the double integrals of the x - coordinate and y - coordinate calculated for the patch area D, representing the first moment of the patch about the y - axis and x - axis respectively; while the denominator ∫∫Ddxdy is the area of the patch.

[0072] Specifically, the longitude and latitude of the centroid of the patch are calculated according to the above formula, so that the centroid distance between two patches in the patch group can be calculated, facilitating patch screening based on the centroid distance.

[0073] Step S106: Trim the linear feature, determine whether the first patch to be included and the second patch to be included meet the inclusion rules, and conduct a boundary verification to obtain the target urban entity area.

[0074] Specifically, after adding the second patch to be included, trim the linear feature of the resulting layer, determine whether the added first patch to be included and the second patch to be included meet the inclusion rules, and conduct a boundary verification of the urban area to obtain the target urban entity area. As Figure 11 shown, the above steps are implemented using the ArcGIS platform and Python programming technology, reducing manual operations, improving work efficiency, and enabling rapid update of the urban entity area.

[0075] Among them, step S106 includes: using a clipping tool to trim the linear feature; determining whether the first patch to be included and the second patch to be included meet the inclusion rules. If they do not meet the rules, discard them; if they meet the rules, retain them; compare the urban entity area with the control line, mark the patches that intersect the control line, and conduct a boundary verification. Discard or retain the marked patches according to the boundary verification results, and output to obtain the target urban entity area.

[0076] Specifically, in the ArcGIS platform, trim the linear feature of the urban entity area where the second patch to be included is added, such as whether there are connection conditions on both sides of rivers and railways; the intersection situation can be checked and recorded by defining a function to determine whether the first patch to be included and the second patch to be included meet the inclusion rules, such as whether the patch position conforms to the actual situation and whether the patch appears repeatedly, etc. If they do not meet the rules, discard them; if they meet the rules, retain them; compare the urban entity area with the control line, mark the patches that intersect the control line, and conduct a boundary verification. Discard or retain the marked patches according to the boundary verification results, such as retaining the marked patches within the boundary and discarding those outside the boundary, etc., output to obtain the target urban entity area, and save it to the set path.

[0077] In this embodiment, a land use polygon layer is obtained based on the survey data, and polygons are extracted from it based on the urban-village attribute code to obtain the initial urban area; a buffer zone with a preset distance is set based on the initial urban area to obtain the expanded space, and land use polygons in the intersection area of the land use polygon layer and the expanded space are screened according to the land use polygon attributes to obtain the first polygons to be included; the first polygons to be included are added to the initial urban area or the urban entity area of the previous iteration to obtain the current urban entity area, and the buffer zone setting to polygon addition steps are iterated until the preset number of iterations is reached to obtain the initial urban entity area; the urban entity area of the previous cycle is obtained, compared with the initial urban entity area to obtain the second polygons to be included, and added to the initial urban entity area; linear features are trimmed, and it is determined whether the first polygons to be included and the second polygons to be included meet the inclusion rules, and boundary verification is performed to obtain the target urban entity area. By integrating the powerful functions of the ArcGIS platform and the flexibility of Python programming, the rapid update of the standardization of the urban entity area is realized, reducing the manual polygon selection operation, effectively solving problems such as data lag, inconsistent standard implementation, cumbersome operation, and mismatch between the update speed and the planning management requirements, providing a real-time and accurate urban entity area update solution for the field of territorial space planning, and strongly supporting the refined and dynamic planning management requirements under the background of new urbanization.

[0078] Those of ordinary skill in the art can understand that all or part of the processes of the methods in the above embodiments can be completed by instructing relevant hardware through a computer program. The program can be stored in a computer-readable storage medium. When the program is executed, it can include the processes of the embodiments of the above methods. Among them, the storage medium can be a magnetic disk, an optical disc, a read-only memory (ROM), or a random access memory (RAM), etc.

[0079] Obviously, those skilled in the art should understand that the above modules or steps of the present invention can be implemented by a general-purpose computing device. They can be concentrated on a single computing device or distributed on a network composed of multiple computing devices. Optionally, they can be implemented by program codes executable by the computing device, so that they can be stored in a computer storage medium (ROM / RAM, magnetic disk, optical disc) and executed by the computing device. And in some cases, the steps shown or described can be executed in a different order than here, or they can be separately made into individual integrated circuit modules, or multiple modules or steps among them can be made into a single integrated circuit module to implement. Therefore, the present invention is not limited to any specific combination of hardware and software.

[0080] The above content is a further detailed description of the present invention in combination with specific embodiments, and it cannot be determined that the specific implementation of the present invention is only limited to these descriptions. For those of ordinary skill in the technical field to which the present invention pertains, without departing from the concept of the present invention, several simple deductions or substitutions can still be made, and all should be regarded as belonging to the protection scope of the present invention.

Claims

1. A method for updating urban physical areas based on ArcGIS platform, characterized in that: The following steps are involved: S1. Obtain a land type map layer according to the survey data, and extract the spots in the land type map layer according to the town and village attribute codes to obtain the initial scope of the urban area; S2. Setting a buffer zone of a preset distance based on the initial range of the urban area to obtain an external expansion space; S3, screening the ground feature spots in the intersection area between the ground feature spot layer and the outer expansion space according to the ground feature spot attributes, and obtaining the first spots to be included; S4, adding the first image patch to be included to the initial range of the urban area or the physical area of ​​the urban area of ​​the last iteration to obtain the current physical area of ​​the urban area, iterating steps S2 to S4 until a preset number of iterations is reached to obtain the initial physical area of ​​the urban area; S5, obtaining the physical area of ​​the urban area in the previous cycle, comparing it with the initial physical area of ​​the urban area, obtaining the second map spot to be included, and adding it to the initial physical area of ​​the urban area; S6. Trim the linear features, determine whether the first to-be-included map spot and the second to-be-included map spot meet the inclusion rules, and perform boundary verification to obtain the physical area of ​​the target urban area.

2. The method for updating the physical area of ​​a city based on the ArcGIS platform according to claim 1, characterized in that: The step S1 comprises: Searching for survey data according to a land change survey database and extracting a land type map layer from the survey data; According to the town and village attribute codes, the urban area spots are screened in the land type map spot layer to obtain the initial scope of the urban area; All the screened urban areas are merged to obtain the first layer.

3. The method for updating the physical area of ​​a city based on the ArcGIS platform according to claim 2, characterized in that: The step S2 comprises: Acquire a preset distance, and generate a buffer zone in the first layer based on the preset distance to obtain a buffer zone layer; An intersection inversion operation is performed on the buffer layer and the first layer to obtain an external expansion space.

4. The method for updating the physical area of ​​a city based on the ArcGIS platform according to claim 1, characterized in that: The step S3 comprises: Based on the attributes of the ground feature spots, the spots to be confirmed that intersect the area of ​​the outer expansion space are screened out in the ground feature spot layer; The spots that have been included in the initial range of the urban area are removed from the spots to be confirmed to obtain a first spot to be included, wherein the first spot to be included carries a mark value of the current iteration number.

5. The method for updating the physical area of ​​a city based on the ArcGIS platform according to claim 1, characterized in that: The step S4 comprises: In the initial iteration, the first to-be-included map spot is added to the initial range of the urban area by using an appending tool to obtain a physical area of ​​the urban area; In non-initial iteration, the first to-be-included map patch is added to the urban physical area obtained in the previous iteration by using an appending tool to obtain the current urban physical area; Repeat the buffer zone setting, first map spot acquisition and map spot addition operations until the preset number of iterations is reached to obtain the initial urban physical area.

6. The method for updating the physical area of ​​a city based on the ArcGIS platform according to claim 1, characterized in that: The step S5 comprises: Importing the physical area of ​​the urban area in the previous cycle, and screening out the same spots that are completely the same as the physical area of ​​the initial urban area; Based on the same image spots, reverse selection is performed in the urban physical area of ​​the previous cycle to obtain different image spots different from the initial urban physical area, which are recorded as the third layer; The target layer is set as the land class spot layer, the source layer is set as the third layer, and a spatial selection method including source layer elements is used to select spots matching the boundaries of the third layer in the land class spot layer according to the land class spot attributes, so as to obtain a fourth layer; Performing field adding operations on the third layer and the fourth layer respectively, calculating the area and centroid longitude and latitude of each map spot, and adding a unique identifier to the fourth layer; Calculate the overlapped part of the third layer and the fourth layer, fuse them to obtain the layer of interest, and add a field to the layer of interest to obtain the centroid longitude and latitude of the patch and calculate the patch area; Matching the layer of interest with the fourth layer to obtain a one-to-one corresponding pattern group in the two layers; Calculate and store the centroid distance between two spots in the spot group and the area ratio of the spots in the layer of interest in the spots of the corresponding fourth layer; In the fourth layer, the spots whose area proportion is greater than a preset threshold, whose centroid distance is less than a preset distance threshold, and whose land type names are the same as those of the corresponding spots in the third layer are screened out to obtain the second spots to be included and added to the initial urban physical area.

7. The method for updating the physical area of ​​a city based on the ArcGIS platform according to claim 1, characterized in that: The step S5 comprises: Importing the physical area of ​​the urban area in the previous cycle, and screening out the same spots that are completely the same as the physical area of ​​the initial urban area; Based on the same image spots, reverse selection is performed in the urban physical area of ​​the previous cycle to obtain different image spots different from the initial urban physical area, which are recorded as the third layer; The target layer is set as the land class spot layer, the source layer is set as the third layer, and a spatial selection method within the source layer element range is used to select spots matching the boundaries of the third layer in the land class spot layer according to the land class spot attributes to obtain a fifth layer; Performing field adding operations on the third layer and the fifth layer respectively, calculating the area and centroid longitude and latitude of each map spot, and adding a unique identifier to the fifth layer; Calculate the overlapped part of the third layer and the fifth layer, fuse them to obtain the layer of interest, and add a field to the layer of interest to obtain the centroid longitude and latitude of the patch and calculate the patch area; Matching the layer of interest with the fifth layer to obtain a one-to-one corresponding pattern group in the two layers; Calculate and store the centroid distance between two spots in the spot group and the area ratio of the spots in the layer of interest in the spots in the corresponding fifth layer; In the third layer, the spots whose area proportion is greater than a preset threshold, whose centroid distance is less than a preset distance threshold, and whose land type names are the same as those of the corresponding spots in the fifth layer are screened out to obtain the second spots to be included and added to the initial urban physical area.

8. The method for updating the physical area of ​​an urban area based on the ArcGIS platform according to any one of claims 6 or 7, characterized in that: The calculation formula of the centroid longitude and latitude is: Where D represents the area where the patch is located, ∫∫ D xdxdy and ∫∫ D xdxdy is the quadratic integral of the x-coordinate and y-coordinate of the image area D, representing the first moment of the image about the y-axis and the x-axis; and the denominator ∫∫ D dxdy is the area of ​​the patch.

9. The method for updating the physical area of ​​a city based on the ArcGIS platform according to claim 1, characterized in that: The step S6 comprises: Use the cropping tool to trim linear features that extend outside the concentrated contiguous area; Determine whether the first to-be-included image spot and the second to-be-included image spot meet the inclusion rules, discard them if they do not meet the rules, and retain them if they meet the rules; The physical area of ​​the urban area is compared with the control line, the spots intersecting with the control line are marked, and a boundary check is performed. The marked spots are discarded or retained according to the boundary check results, and the physical area of ​​the target urban area is output.

Images