Arcpy-based land parcel boundary point identification and optimization sorting method

Abstract

The invention relates to an Arcpy-based land parcel boundary point identification and optimization sorting method, which comprises the following steps of: processing land parcel data by utilizing ArcGIS software to generate a land parcel boundary point data set, calculating a minimum enclosing rectangle of a land parcel and taking the minimum enclosing rectangle as a land parcel range, and determining a starting point and a sorting basis of land parcel boundary points according to the land parcel range. And traversing the land parcel boundary point data set to start from the determined land parcel boundary point starting point, sorting and numbering all boundary points of the land parcel to obtain a sorted and numbered boundary point sequence, and finally processing the sorted and numbered boundary point sequence to obtain complete boundary point result data of the land parcel. Therefore, the boundary point identification and sorting of the land parcel can be completed only by using one ArcGIS software, and the efficiency and accuracy of land parcel boundary point identification and optimized sorting are improved.

Description

Technical Field

[0001] The present invention relates to the field of surveying and mapping, and in particular to a method for identifying and optimizing land boundary points based on Arcpy. Background Art

[0002] Boundary points are turning points of parcels or property boundaries, that is, inflection points. The method of sorting boundary points is to take the parcel as the unit, starting from the northwest corner in a clockwise direction, starting with the first boundary point numbered J1, first number the outer boundary points of the parcel within the parcel, and then continue to number and sort the boundary points of the internal parcels according to the sequence number of the outer boundary points.

[0003] In the actual surveying and mapping process, manual identification and sorting methods and CASS software identification and sorting methods are usually used to identify the boundary points of the plot, and then these identified boundary points are numbered and sorted.

[0004] However, the existing methods for identifying and sorting parcel boundary points have shortcomings: although the manual identification and sorting method can achieve high-precision point identification and sorting based on intuitive judgment when facing regular-shaped parcels, its operation process is cumbersome and time-consuming, and it is difficult to adapt to the needs of some complex parcels and large-scale work tasks; although the CASS software identification and sorting method can automatically identify a large amount of parcel data, the accuracy of its automatic identification is often difficult to guarantee, and it often needs to rely on subsequent manual review and manual adjustment, which undoubtedly increases the workload and reduces the overall identification and sorting efficiency. Summary of the invention

[0005] The technical problem to be solved by the present invention is to provide an Arcpy-based method for identifying and optimizing the land boundary points, which can automatically, quickly and accurately identify and sort the land boundary points.

[0006] The technical solution adopted by the present invention to solve the above technical problems is: a method for identifying and optimizing the land parcel boundary points based on Arcpy, which is characterized by comprising the following steps:

[0007] Step 1, based on the ArcPy library of ArcGIS software, and using ArcGIS software to process the land parcel data, to generate a land parcel boundary point data set; wherein the land parcel boundary point data set is formed by processing all the boundary point data of the land parcel;

[0008] Step 2, calculate the minimum enclosing rectangle of the plot, and use the minimum enclosing rectangle as the plot range;

[0009] Step 3, according to the determined land parcel scope, determine the starting point and sorting basis of the land parcel boundary points;

[0010] Step 4, traversing the land parcel boundary point data set, starting from the determined land parcel boundary point starting point, sorting and numbering all the land parcel boundary points to obtain a sorted and numbered boundary point sequence;

[0011] Step 5, process the obtained sequence of boundary points after sorting and numbering to obtain the complete boundary point result data of the land parcel.

[0012] Improved, in the Arcpy-based method for identifying and optimizing the sorting of parcel boundary points, in step 3, the starting point of the minimum enclosing rectangle is first calculated using the azimuth function, and then the first boundary point in the parcel boundary point numbering sorting is calculated based on the parcel boundary point starting point function and the starting point of the minimum enclosing rectangle, and finally the boundary point sorting function is used to determine the positions and numbering order of all boundary points of the minimum enclosing rectangle.

[0013] Furthermore, in the Arcpy-based method for identifying and optimizing the parcel boundary points, the minimum enclosing rectangle function of the GIS software is used to create the minimum enclosing rectangle of the parcel, so as to determine the first boundary point of the parcel boundary point numbering sequence;

[0014] Use the database creation function of the GIS software to process the sequence of boundary points after sorting and numbering to obtain the complete boundary point result data of the land parcel.

[0015] Further improvement, in the Arcpy-based method for identifying and optimizing the parcel boundary points, the starting point of the minimum enclosing rectangle is calculated as follows:

[0016] Obtain the maximum value of the ordinate and the minimum value of the abscissa of the minimum circumscribed rectangle;

[0017] Determine the longitudinal position point corresponding to the maximum value of the longitudinal coordinate and the transverse position point corresponding to the minimum value of the transverse coordinate respectively;

[0018] Calculate the azimuth angle between the longitudinal position point and the transverse position point;

[0019] Make judgments based on the calculated azimuth:

[0020] When the azimuth angle is greater than 45°, the determined longitudinal position point is selected as the starting point of the minimum circumscribed rectangle; otherwise, the determined transverse position point is selected as the starting point of the minimum circumscribed rectangle.

[0021] Furthermore, in the Arcpy-based method for identifying and optimizing the sorting of land parcel boundary points, the first boundary point of the land parcel number sorting is calculated as follows:

[0022] Traverse the land parcel boundary point data set, and calculate the Euclidean distance value from each boundary point in the land parcel boundary point data set to the starting point of the minimum circumscribed rectangle corresponding to the land parcel;

[0023] The boundary point corresponding to the minimum Euclidean distance value among all the obtained Euclidean distance values ​​is taken as the first boundary point in the parcel number sorting.

[0024] Further improved, in this invention, the Arcpy-based method for identifying and optimizing the parcel boundary points also includes:

[0025] Step a1, sorting the boundary points clockwise according to the first boundary point sorted by the parcel number and the original sequence number of the first boundary point in the boundary points;

[0026] Step a2, determine whether there is a situation where the outer ring parcels are enclosed in the inner ring parcels:

[0027] If it exists, go to step a3; otherwise, sort and number the boundary points on the land in a clockwise order according to the order of the land numbering;

[0028] Step a3, calculate the Euclidean distance from the minimum circumscribed rectangle to each boundary point on the outer ring parcel and each boundary point on the inner ring parcel, and use the boundary point located on the outer ring parcel and corresponding to the minimum Euclidean distance value as the starting point for the parcel numbering, use the boundary point located on the inner ring parcel and corresponding to the minimum Euclidean distance value as the starting point for the inner ring parcel numbering, use the next number of the last boundary point on the outer ring parcel as the starting number for the inner ring parcel numbering, and then number each boundary point on the inner ring parcel in a counterclockwise order.

[0029] Furthermore, in the Arcpy-based method for identifying and optimizing the sorting of land parcel boundary points, in step 5, the processing of the sorted and numbered boundary point sequence includes performing a turning point conversion operation on the land parcel; wherein the turning point conversion operation includes the following steps:

[0030] Use ArcGIS software to convert the parcel vertices into a turning point feature class;

[0031] Add the boundary point result layer and the converted turning point feature class to the map document;

[0032] Connect the boundary point result feature class and the obtained turning point feature class according to the polygon_id field, and export the connected data as result data to obtain the boundary point result data, so that the boundary point result data inherits the attributes of the parcel data;

[0033] ArcGIS software was used to calculate the X and Y coordinates of the boundary point results data;

[0034] Add the boundary point results data to the map document.

[0035] Further improvement, in the Arcpy-based method for identifying and optimizing the parcel boundary points, the complete boundary point result data of the parcel includes the boundary point coordinates, the boundary point sequence number and the information of the parcel where the boundary point is located.

[0036] Compared with the prior art, the advantages of the present invention are: the Arcpy-based method for identifying and optimizing the sorting of land parcel boundary points of the present invention uses ArcGIS software to process land parcel data to generate a land parcel boundary point data set, calculates the minimum enclosing rectangle of the land parcel and uses the minimum enclosing rectangle as the land parcel range, and then determines the starting point and sorting basis of the land parcel boundary points according to the land parcel range, and traverses the land parcel boundary point data set starting from the determined land parcel boundary point starting point, sorts and numbers all the land parcel boundary points to obtain a sorted and numbered boundary point sequence, and finally processes the obtained sorted and numbered boundary point sequence to obtain the complete boundary point result data of the land parcel. In this way, only one ArcGIS software is needed to complete the identification and sorting of the land parcel boundary points, thereby improving the efficiency and accuracy of the identification and optimization sorting of the land parcel boundary points.

[0037] Secondly, the Arcpy-based method for identifying and optimizing the sorting of land parcel boundary points of the invention obtains the maximum longitudinal coordinate and the minimum transverse coordinate of the minimum circumscribed rectangle, determines the longitudinal position point corresponding to the maximum longitudinal coordinate and the transverse position point corresponding to the minimum transverse coordinate, and then calculates the azimuth between the longitudinal position point and the transverse position point, so as to accurately determine the starting point of the land parcel; moreover, the method for optimizing the sorting of land parcel boundary points accurately sorts and rearranges the boundary points of the land parcel according to the starting point, Euclidean distance and other conditions, especially when there is an outer ring land parcel nested in an inner ring land parcel, the boundary point lists of the outer ring land parcel and the inner ring land parcel can be processed more reasonably. BRIEF DESCRIPTION OF THE DRAWINGS

[0038] Figure 1 It is a flow chart of a method for identifying and optimizing parcel boundary points based on Arcpy in an embodiment of the present invention;

[0039] Figure 2 It is a schematic diagram of the entire land parcel containing regular parcels and irregular parcels;

[0040] Figure 3 for Figure 2 A schematic diagram of generating the minimum enclosing rectangle for the entire plot;

[0041] Figure 4 for Figure 2 The schematic diagram shows the generation of the first boundary point (J1) corresponding to each plot in the overall plot. DETAILED DESCRIPTION

[0042] The present invention is further described in detail below with reference to the accompanying drawings.

[0043] This embodiment provides a method for identifying and optimizing the parcel boundary points based on Arcpy. Figure 1 As shown, the Arcpy-based method for identifying and optimizing the land parcel boundary points of this embodiment includes the following steps:

[0044] Step 1, based on the ArcPy library of ArcGIS software, and using ArcGIS software to process the land parcel data, to generate a land parcel boundary point data set; wherein the land parcel boundary point data set is formed by processing all the boundary point data of the land parcel;

[0045] Step 2, calculate the minimum enclosing rectangle of the plot, and use the minimum enclosing rectangle as the plot range;

[0046] Step 3, according to the determined land parcel scope, determine the starting point and sorting basis of the land parcel boundary points;

[0047] In this embodiment, the starting point of the minimum enclosing rectangle is first calculated by the azimuth function, and then the first boundary point (J1) of the parcel boundary point numbering sequence is calculated according to the parcel boundary point starting point function and the starting point of the minimum enclosing rectangle, and finally the position and numbering sequence of all boundary points of the minimum enclosing rectangle are determined by the boundary point sorting function;

[0048] Step 4, traversing the land parcel boundary point data set, starting from the determined land parcel boundary point starting point, sorting and numbering all the land parcel boundary points to obtain a sorted and numbered boundary point sequence;

[0049] Step 5, processing the obtained sorted and numbered boundary point sequence to obtain complete boundary point result data of the land parcel; wherein the complete boundary point result data of the land parcel here includes the boundary point coordinates, boundary point sequence number and information of the land parcel where the boundary point is located.

[0050] Specifically in the above step 2 of the embodiment, the minimum enclosing rectangle function of the GIS software is used to create the minimum enclosing rectangle of the plot, so as to determine the first boundary point in the sorted order of the plot boundary points; then the database creation function of the GIS software is used to process the sequence of boundary points after sorting and numbering to obtain the complete boundary point result data of the plot.

[0051] It should be noted that, for the land parcel boundary point identification and optimization sorting method of this embodiment, in the above step 2, the starting point of the minimum circumscribed rectangle is calculated as follows:

[0052] Get the maximum value of the vertical coordinate and the minimum value of the horizontal coordinate of the minimum bounding rectangle;

[0053] Determine the longitudinal position point corresponding to the maximum value of the longitudinal coordinate and the transverse position point corresponding to the minimum value of the transverse coordinate respectively;

[0054] Calculate the azimuth angle between the longitudinal position point and the transverse position point;

[0055] Make judgments based on the calculated azimuth:

[0056] When the azimuth angle is greater than 45°, the determined longitudinal position point is selected as the starting point of the minimum circumscribed rectangle; otherwise, the determined transverse position point is selected as the starting point of the minimum circumscribed rectangle.

[0057] In addition, for the parcel number sorting, the calculation method of the first boundary point of the above parcel number sorting is as follows: traverse the parcel boundary point data set, and calculate the Euclidean distance value from each boundary point in the parcel boundary point data set to the starting point of the minimum circumscribed rectangle corresponding to the parcel; take the boundary point corresponding to the minimum Euclidean distance value among all the obtained Euclidean distance values ​​as the first boundary point of the parcel number sorting.

[0058] In addition, considering that the outer ring parcel may be nested within the inner ring parcel during the surveying and mapping process, the parcel boundary point identification and optimization sorting method of this embodiment has also made the following improvements, that is, it includes the following steps:

[0059] Step a1, sorting the boundary points clockwise according to the first boundary point sorted by the parcel number and the original sequence number of the first boundary point in the boundary points;

[0060] Step a2, determine whether there is a situation where the outer ring parcels are enclosed in the inner ring parcels:

[0061] If it exists, go to step a3; otherwise, sort and number the boundary points on the land in a clockwise order according to the order of the land numbering;

[0062] Step a3, calculate the Euclidean distance from the minimum circumscribed rectangle to each boundary point on the outer ring parcel and each boundary point on the inner ring parcel, and use the boundary point located on the outer ring parcel and corresponding to the minimum Euclidean distance value as the starting point for the parcel numbering, use the boundary point located on the inner ring parcel and corresponding to the minimum Euclidean distance value as the starting point for the inner ring parcel numbering, use the next number of the last boundary point on the outer ring parcel as the starting number for the inner ring parcel numbering, and then number each boundary point on the inner ring parcel in a counterclockwise order.

[0063] More specifically, with respect to the above step 5, the processing of the sequence of boundary points after sorting and numbering includes performing a turning point conversion operation on the parcel; wherein the turning point conversion operation includes the following steps:

[0064] Step b1, using ArcGIS software to convert the parcel vertex into a turning point feature class;

[0065] Step b2, adding the boundary point result layer and the turning point feature class obtained after conversion to the map document;

[0066] Step b3, connecting the boundary point result feature class and the obtained turning point feature class according to the polygon_id field, and exporting the connected data as result data to obtain the boundary point result data, so that the boundary point result data inherits the attributes of the parcel data;

[0067] Step b4, using ArcGIS software to calculate the X coordinate and Y coordinate of the boundary point result data respectively;

[0068] Step b5, adding the boundary point result data to the map document.

[0069] Although the preferred embodiments of the present invention are described in detail above, it should be clearly understood that various modifications and variations are possible for those skilled in the art. Any modification, equivalent substitution, improvement, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims

1. A method for identifying and optimizing parcel boundary points based on Arcpy, characterized in that: The following steps are involved: Step 1, based on the ArcPy library of ArcGIS software, and using ArcGIS software to process the land parcel data, to generate a land parcel boundary point data set; wherein the land parcel boundary point data set is formed by processing all the boundary point data of the land parcel; Step 2, calculate the minimum enclosing rectangle of the plot, and use the minimum enclosing rectangle as the plot range; Step 3, according to the determined land parcel scope, determine the starting point and sorting basis of the land parcel boundary points; Step 4, traversing the land parcel boundary point data set, starting from the determined land parcel boundary point starting point, sorting and numbering all the land parcel boundary points to obtain a sorted and numbered boundary point sequence; Step 5, process the obtained sequence of boundary points after sorting and numbering to obtain the complete boundary point result data of the land parcel.

2. The Arcpy-based method for identifying and optimizing land parcel boundary points according to claim 1, characterized in that: In step 3, the starting point of the minimum enclosing rectangle is first calculated using the azimuth function, and then the first boundary point in the parcel boundary point numbering sequence is calculated based on the parcel boundary point starting point function and the starting point of the minimum enclosing rectangle. Finally, the boundary point sorting function is used to determine the position and numbering order of all boundary points of the minimum enclosing rectangle.

3. The Arcpy-based method for identifying and optimizing land parcel boundary points according to claim 2, characterized in that: Using the minimum enclosing rectangle function of the GIS software to create the minimum enclosing rectangle of the land parcel, so as to determine the first boundary point in the order of the boundary point numbers of the land parcel; Use the database creation function of the GIS software to process the sequence of boundary points after sorting and numbering to obtain the complete boundary point result data of the land parcel.

4. The Arcpy-based method for identifying and optimizing land parcel boundary points according to claim 2, characterized in that: The starting point of the minimum enclosing rectangle is calculated as follows: Obtain the maximum value of the ordinate and the minimum value of the abscissa of the minimum circumscribed rectangle; Determine the longitudinal position point corresponding to the maximum value of the longitudinal coordinate and the transverse position point corresponding to the minimum value of the transverse coordinate respectively; Calculate the azimuth angle between the longitudinal position point and the transverse position point; Make judgment based on the calculated azimuth: When the azimuth angle is greater than 45°, the determined longitudinal position point is selected as the starting point of the minimum circumscribed rectangle; otherwise, the determined transverse position point is selected as the starting point of the minimum circumscribed rectangle.

5. The Arcpy-based method for identifying and optimizing parcel boundary points according to claim 4, characterized in that: The first boundary point of the parcel numbering sequence is calculated as follows: Traverse the land parcel boundary point data set, and calculate the Euclidean distance value from each boundary point in the land parcel boundary point data set to the starting point of the minimum circumscribed rectangle corresponding to the land parcel; The boundary point corresponding to the minimum Euclidean distance value among all the obtained Euclidean distance values ​​is taken as the first boundary point in the parcel number sorting.

6. The Arcpy-based method for identifying and optimizing land parcel boundary points according to claim 5, characterized in that: Also includes: Step a1, sorting the boundary points clockwise according to the first boundary point sorted by the parcel number and the original sequence number of the first boundary point in the boundary points; Step a2, determine whether there is a situation where the outer ring parcels are enclosed in the inner ring parcels: If it exists, go to step a3; otherwise, sort and number the boundary points on the land in a clockwise order according to the order of the land numbering; Step a3, calculate the Euclidean distance from the minimum circumscribed rectangle to each boundary point on the outer ring parcel and each boundary point on the inner ring parcel, and use the boundary point located on the outer ring parcel and corresponding to the minimum Euclidean distance value as the starting point for the parcel numbering, use the boundary point located on the inner ring parcel and corresponding to the minimum Euclidean distance value as the starting point for the inner ring parcel numbering, use the next number of the last boundary point on the outer ring parcel as the starting number for the inner ring parcel numbering, and then number each boundary point on the inner ring parcel in a counterclockwise order.

7. The Arcpy-based method for identifying and optimizing land parcel boundary points according to claim 6, characterized in that: In step 5, the processing of the sequence of boundary points after sorting and numbering includes performing a turning point conversion operation on the parcel; wherein the turning point conversion operation includes the following steps: Use ArcGIS software to convert the parcel vertices into a turning point feature class; Add the boundary point result layer and the converted turning point feature class to the map document; Connect the boundary point result feature class and the obtained turning point feature class according to the polygon_id field, and export the connected data as result data to obtain the boundary point result data, so that the boundary point result data inherits the attributes of the parcel data; ArcGIS software was used to calculate the X and Y coordinates of the boundary point results data; Add the boundary point results data to the map document.

8. The Arcpy-based method for identifying and optimizing land parcel boundary points according to any one of claims 1 to 7, characterized in that: The complete boundary point result data of the plot includes the boundary point coordinates, the boundary point serial number and the information of the plot where the boundary point is located.

Images