Urban edge development measurement method and system based on continuous slice projection technology

Abstract

The application provides a city edge development measurement method and system based on a continuous slice projection technology, the method comprises the following steps: acquiring land use data in a set space range; processing the land use data by using a continuous slice technology, expanding outward to obtain a buffer zone, and determining edge land use data; preprocessing the edge land use data; performing reverse buffer zone projection to obtain edge development length; calculating total single-edge edge development length and double-edge edge length; calculating edge rates in the set space range; and comparing edge development in different space ranges according to the edge rates in the different space ranges. The application can accurately point to a city edge development research target, accurately measure edge degree by projection, intersection area determination and reduction of repeated calculation of overlapping areas in edge line length calculation, and accurately compare city edge degrees of different areas after unified calculation of edge rates.

Description

Technical Field

[0001] This invention belongs to the field of land and space planning data analysis, and specifically relates to a method and system for measuring urban development along the edge based on continuous slice projection technology. Background Technology

[0002] Urban development along city boundaries refers to the land use characteristics resulting from the development of adjacent urban administrative districts within their border areas. Due to differences in resource endowments, comparative advantages, development goals, and management measures, different cities exhibit diverse land use characteristics on both sides of the boundary line. Particularly in developed cities, this development along boundaries can lead to continuous stretches of construction land, causing urban problems such as ecological damage and conflicts arising from NIMBY (Not In My Backyard) infrastructure development. Therefore, the phenomenon of urban development along city boundaries is increasingly attracting the attention of local governments and is an important aspect of current regional planning and regional collaboration.

[0003] Currently, there is a lack of effective technical methods for measuring and monitoring development along city boundaries. The accuracy and universality of the analysis are insufficient, making it difficult to conduct scientific comparisons between different regions to ensure the effective implementation of regional coordination in territorial spatial planning. Existing methods mainly rely on intuitive visual representations such as satellite remote sensing images and land use distribution maps, or use statistics on land use scale within the administrative boundaries adjacent to city limits to roughly analyze land use along city boundaries. These traditional methods suffer from problems such as imprecise target targeting and imprecise measurement of the degree of boundary overlap, making it difficult to accurately measure the degree of boundary overlap and conduct scientific comparisons between different regions. This hinders the scientific formulation of refined planning policies and the effective implementation of spatial control measures. Summary of the Invention

[0004] One of the objectives of this invention is to provide a method for measuring urban edge development based on continuous slice projection technology. This invention can accurately target the research objectives of urban edge development, and through projection, determination of intersecting areas, and reduction of redundant calculations in edge line length calculations for overlapping areas, it can accurately measure the degree of edge development. After uniformly calculating the edge development rate, it can accurately compare the degree of urban edge development in different areas, which is conducive to the scientific formulation of refined planning policies and the effective implementation of spatial control measures.

[0005] This invention provides the following technical solution: a method for measuring urban edge development based on continuous slice projection technology, wherein the method processes land use data based on continuous slice projection technology, and the method includes the following steps: S1: Obtain land use data within a defined spatial range; S2: The land use data within the set spatial range obtained in step S1 is processed using continuous slicing technology, and a buffer zone is obtained by expanding outward. The edge development technology is then used to determine the edge land use data. S3: Preprocess the edge land use data determined in step S2, and project the buffer established in step S2 onto the preprocessed edge land use data to obtain the edge line of the land parcel. S4: Perform reverse buffer projection to obtain the edge development length; S5: Based on the edge development length obtained in step S4, calculate the total edge development length SRL on one side and the edge development length DRL on both sides; S6: Calculate the edge-fitting rate LR within the set spatial range, and compare the edge-fitting development in different spatial ranges based on the edge-fitting rate LR within different spatial ranges.

[0006] Furthermore, step S2 includes the following steps: S2.1: Using the urban boundary line of the study area where the land use data obtained in step S1 is located as the benchmark, and using b meters as the slicing step unit, gradually expand outward to establish a buffer zone, and establish a buffer zone of a meters on both sides respectively; where 0 < b ≤ 10, 500 < a ≤ 1000; S2.2: Obtain land use data within a meter range on both sides of the urban boundary line of the study area, and define all land use data within this range as edge-bound land use data.

[0007] Furthermore, step S3 includes the following steps: S3.1: Traverse the map within the a-meter buffer using the data of the plot images, determine the shape of the plots based on the plot images, find plots that differ greatly from the city boundary line, or have irregular shapes and are projected onto the development length of the boundary line, and statistically record such plots. S3.2: For the land use data of the plots statistically recorded in step S3.1, the surface elements are modified according to the irregularity type of the plots based on the merging and surface segmentation technology in GIS; S3.3: Extract the parcel boundary lines of the parcel image after it has been modified by the polygon features within the buffer.

[0008] Furthermore, in step S3.2, based on merging and surface segmentation techniques in GIS, the surface features are modified according to the irregularity type of the land parcel, including the following steps: S3.21: For a set of seed points within the buffer zone of a meters, form a seed point set G, and define the line connecting any two points p and q in the seed point set G as the equidistant line for region division. : ;in, , Let x and y be the x and y coordinates of point p in the GIS system, respectively. , These are the x and y coordinates of point q in the GIS system, respectively. S3.22: Following the method in step S3.21, traverse any two points in the seed point set G within the buffer zone of a meters to form multiple equidistant lines, dividing the land space within the buffer zone of a meters into multiple polygonal regions; S3.23: For any two polygonal regions A and B among multiple polygonal regions, calculate the intersecting and overlapping areas between them. and the merged area between the two. : ; ; in, Modify the topological area of ​​the face features in polygon region A; Modify the topological area of ​​the face features in polygon region B; S3.24: Using the method described in step S3.23, traverse the multiple polygonal regions obtained in step S3.22, and calculate the buffer S where the boundaries of the different length buffers intersect: ; The number of combinations calculated in step S3.23 is the number of combinations of any two polygon regions from the plurality of polygon regions formed in step S3.22. M represents the total number of polygonal regions formed in step S3.22. This is to calculate the number of combinations of selecting any two polygons from M polygon regions and calculating their areas.

[0009] Furthermore, in step S3.23, the topological area of ​​the polygonal region A is modified by calculating the face features. Modify the topological area of ​​the polygon region B's face features. The calculation formulas are as follows: ; ; in, This represents the total number of seed points within polygonal region A; , Let x and y be the x and y coordinates of the i-th seed point within polygonal region A, respectively. , Let x and y be the x and y coordinates of the (i+1)th seed point within polygonal region A, respectively. , Let x and y be the x and y coordinates of the first seed point within polygonal region A, respectively. , These are the first two terms within polygonal region A. The x and y coordinates of each seed point; This represents the total number of seed points within the polygonal region B. , Let x and y be the x and y coordinates of the j-th seed point within polygonal region B, respectively. , Let x and y be the x and y coordinates of the (j+1)th seed point within polygon region B, respectively. and are the x and y coordinates of the first seed point within polygon region B, respectively. These are the first two terms within polygonal region B. The x and y coordinates of each seed point.

[0010] Furthermore, the calculation method for the plot boundary line in step S3.3 is to overlay and intersect the buffer zone of a meters with the preprocessed land use data to obtain the tangent of the buffer zone S that intersects the boundaries of buffer zones of different lengths, and extract the length of the longest buffer zone tangent, which is recorded as the plot boundary line.

[0011] Furthermore, the method for performing reverse buffer projection in step S4 is as follows: using the edge line of each plot, a reverse buffer projection is performed on the city boundary line; the buffer distance of the reverse buffer is set to a+200 meters, and the length of the intersection between the reverse buffer and the city boundary line is the corresponding edge development length.

[0012] Furthermore, the formulas for calculating the total development length SRL of a single-sided edge and the length DRL of both sides in step S5 are as follows: ; ; in, It is the edge development length obtained by the nth reverse buffer, K is the total number of reverse buffers; n=1,2,…,K; G is the total number of overlapping regions in the reverse buffers, G<K, Let G be the edge development length of the m-th overlapping region in the reverse buffer; m = 1, 2, ..., G.

[0013] Furthermore, the formula for calculating the edge-fitting ratio LR in step S6 is as follows: .

[0014] The present invention also provides an urban edge development measurement system based on continuous slice projection technology using the method described above, characterized in that the system includes a land use data acquisition module, an edge land use data calculation module, a plot edge line calculation module, a reverse buffer projection module, an edge development length calculation module, and an edge development comparison module. The land use data acquisition module is used to acquire land use data within a set spatial range; The edge-adjacent land use data calculation module is used to process the land use data within a set spatial range obtained by the land use data acquisition module using continuous slicing technology, expand outward to obtain a buffer zone, and further use edge-adjacent development technology to determine the edge-adjacent land use data. The plot boundary line calculation module is used to preprocess the boundary land use data determined by the boundary land use data calculation module, and project the buffer established by the boundary land use data calculation module onto the preprocessed boundary land use data to obtain the plot boundary line. The reverse buffer projection module is used to perform reverse buffer projection to obtain the edge development length; The edge development length calculation module is used to calculate the total edge development length SRL on one side and the edge development length DRL on both sides based on the edge development length obtained by the reverse buffer projection module. The edge-fitting development comparison module is used to calculate the edge-fitting rate LR within the set spatial range, and compare the edge-fitting development in different spatial ranges based on the edge-fitting rate LR in different spatial ranges.

[0015] The present invention has the following beneficial effects: 1. This invention can accurately target the research goals of urban edge development, and accurately measure the degree of edge integration by projection, determining the intersection area, and reducing the repeated calculation of overlapping areas in the edge line length calculation. After uniformly calculating the edge integration rate, it can accurately compare the degree of urban edge integration in different areas.

[0016] 2. Based on the effective acquisition of continuous slice projection technology, a method for measuring the degree of urban edge development using continuous slice projection is designed. For single-sided edge development, continuous slice and linear projection methods are used to accurately calculate the linear length of single-sided edge development. For double-sided edge development, the calculation results of single-sided edge development are used, and the linear length of double-sided edge development is accurately calculated by using the reverse projection method on both sides.

[0017] 3. The method of the present invention has a stable and reliable data acquisition source that can be dynamically updated. The overall algorithm has high efficiency and accuracy in calculation, and can achieve consistency in results.

[0018] Other features and advantages of the invention will be set forth in the description which follows, and will be apparent in part from the description, or may be learned by practicing the invention. The objects and other advantages of the invention may be realized and obtained by means of the structures particularly pointed out in the written description, claims, and drawings.

[0019] The technical solution of the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. Attached Figure Description

[0020] The accompanying drawings are provided to further illustrate the invention and form part of the specification. They are used in conjunction with embodiments of the invention to explain the invention and do not constitute a limitation thereof. In the drawings: Figure 1 This is a schematic diagram of a method for measuring urban development along the edge based on continuous slice projection technology in an embodiment of the present invention; Figure 2 This is an example of the changes in irregular plots before and after preprocessing in an embodiment of the present invention; Figure 3 This is a comparison image before and after the extraction of the land parcel boundary line in an embodiment of the present invention; Figure 4 This is a schematic diagram of the reverse buffer zone projection on one side of the land parcel in an embodiment of the present invention; Figure 5 This is a schematic diagram of the double-sided reverse buffer zone projection of the land parcel in an embodiment of the present invention; Figure 6 This is a schematic diagram of an urban edge-fitting development measurement system based on continuous slice projection technology in an embodiment of the present invention; Figure 7 This is a schematic diagram of the development boundary of towns bordering Shanghai, measured using the method of this invention. Detailed Implementation

[0021] The preferred embodiments of the present invention will be described below with reference to the accompanying drawings. It should be understood that the preferred embodiments described herein are for illustration and explanation only and are not intended to limit the present invention.

[0022] This invention provides a method for measuring urban development along city edges based on continuous slice projection technology, such as... Figure 1 As shown, the method provided by the present invention includes the following steps: This invention provides a method for measuring urban development along city edges based on continuous slicing projection technology. The method processes land use data using continuous slicing projection technology, such as... Figure 1 As shown, the method includes the following steps: S1: Obtain land use data within a defined spatial range; Generally, land use data in a land information system is obtained at a specific point in time within a specific spatial area to be studied, based on research needs, through methods such as current status surveys and satellite remote sensing interpretation. Land use data includes aerial photographs of the land obtained by land remote sensing, land elevation, land area, and latitude and longitude. S2: The land use data within the set spatial range obtained in step S1 is processed using continuous slicing technology, and a buffer zone is obtained by expanding outward. The edge development technology is then used to determine the edge land use data. S3: Preprocess the edge land use data determined in step S2, and project the buffer established in step S2 onto the preprocessed edge land use data to obtain the edge line of the plot. S4: Perform reverse buffer projection to obtain the edge development length; S5: Based on the edge development length obtained in step S4, calculate the total edge development length SRL for one side and the edge development length DRL for both sides; S6: Calculate the edge-fitting rate LR within the set spatial range, and compare the edge-fitting development in different spatial ranges based on the edge-fitting rate LR in different spatial ranges.

[0023] Step S2 includes the following steps: S2.1: Using the urban boundary line of the study area where the land use data obtained in step S1 is located as the benchmark, and using b meters as the slicing step unit, gradually expand outward to establish a buffer zone, and establish a buffer zone of a meters on both sides respectively; where 0 < b ≤ 10, 500 < a ≤ 1000; S2.2: Obtain land use data within a meter radius on both sides of the urban boundary line of the study area, and define all land use data within this range as edge-bound land use data.

[0024] In practical applications, land use data often contains irregularly shaped plots, exhibiting issues such as deviations from urban boundaries or gaps within the plots. These issues significantly impact the accuracy of edge-following development analysis based on land use data. Therefore, preprocessing and modification of land use data are necessary to improve the accuracy of edge-following development analysis. Thus, as another preferred embodiment of the present invention, step S3 includes the following steps: S3.1: Traverse the map within the a-meter buffer using the data of the plot images, determine the shape of the plots based on the plot images, find plots that differ greatly from the city boundary line, or have irregular shapes and are projected onto the development length of the boundary line, and statistically record such plots. S3.2: For the land use data of the plots statistically recorded in step S3.1, based on the merging and surface segmentation technology in GIS (Geographic Information System), the surface elements are modified according to the irregular type of the plots; S3.3: Extract the parcel boundary lines of the parcel image after it has been modified by the polygon features within the buffer.

[0025] like Figure 2 (a) is a schematic diagram of the irregular plots before preprocessing. Figure 2(b) is a schematic diagram of the land parcel after preprocessing using steps S3.1-S3.2 of the method of the present invention. The present invention preprocesses irregular land parcels without affecting the actual land parcel edge development length, and ensures that the intersection of the land parcel and the buffer zone boundary line can yield the land parcel edge development length value that is closest to the actual value.

[0026] In GIS, merging polygon features refers to combining multiple polygons (faces) into a new polygon. This process involves geometric and topological operations and generates new polygon features. Merging polygon features can be achieved through geometric set operations (such as union operations). Therefore, as another preferred embodiment of the present invention, step S3.2, based on merging and polygon segmentation techniques in GIS, modifies polygon features according to the irregularity type of the land parcel, including the following steps: S3.21: For a set of seed points within a buffer zone of length a meters, form a seed point set G. Define the line connecting any two points p and q in the seed point set G as the equidistant line for region division. : ;in, , Let x and y be the x and y coordinates of point p in the GIS system, respectively. , These are the x and y coordinates of point q in the GIS system, respectively. S3.22: Following the method in step S3.21, traverse any two points in the seed point set G within the buffer zone of a meters to form multiple equidistant lines, dividing the land space within the buffer zone of a meters into multiple polygonal regions; S3.23: For any two polygonal regions A and B among multiple polygonal regions, calculate the intersecting and overlapping areas between them. and the merged area between the two. : ; ; in, Modify the topological area of ​​the face features of polygon region A after steps S3.21-S3.22; Modify the topological area of ​​the polygon region B after steps S3.21-S3.22; S3.24: Using the method in step S3.23, traverse the multiple polygonal regions obtained in step S3.22, and calculate the buffer S where the boundaries of the different length buffers intersect. ; Let S be the number of combinations calculated in step S3.23 from any two polygon regions formed in step S3.22. M represents the total number of polygonal regions formed in step S3.22. This is to calculate the number of combinations of selecting any two polygons from M polygon regions and calculating their areas.

[0027] Steps S3.21-S3.22 completed the polygon segmentation of land use data within a 'a' meter buffer in the GIS system. Step S3.23 completed the polygon merging (overlapping merging) of multiple regions obtained from the polygon segmentation. This process merges multiple adjacent or intersecting polygon regions obtained from the polygon segmentation into a new polygon region, often used to simplify complex geographic data and reduce data redundancy. The primary goal of polygon segmentation and merging is to modify polygonal regions in land images within a Geographic Information System (GIS) for further analysis and processing.

[0028] Furthermore, in step S3.23, the modified topological area of ​​the polygon region A after steps S3.21-S3.22 is calculated. Modify the topological area of ​​the polygon region B's surface features after steps S3.21-S3.22. The calculation formulas are as follows: ; ; in, This represents the total number of seed points within polygonal region A; , Let x and y be the x and y coordinates of the i-th seed point within polygonal region A, respectively. , Let x and y be the x and y coordinates of the (i+1)th seed point within polygonal region A, respectively. , Let x and y be the x and y coordinates of the first seed point within polygonal region A, respectively. , These are the first two terms within polygonal region A. The x and y coordinates of each seed point; This represents the total number of seed points within the polygonal region B. , Let x and y be the x and y coordinates of the j-th seed point within polygonal region B, respectively. , Let x and y be the x and y coordinates of the (j+1)th seed point within polygon region B, respectively. , Here, x and y are the x and y coordinates of the first seed point within polygon region B, respectively. , These are the first two terms within polygonal region B. The x and y coordinates of each seed point.

[0029] Preferably, the calculation method for the plot boundary line in step S3.3 is to overlay and intersect the 'a' meter buffer zone with the preprocessed boundary land use data to obtain the tangent lines of the buffer zone S where different length buffer zones intersect, and extract the length of the longest buffer zone tangent line, which is recorded as the plot boundary line. Figure 3 As shown, this is for Figure 3 (a) After processing multiple polygonal regions, the curved boundary of the green plot is taken to obtain... Figure 3 (b) The longest tangent length of the yellow continuous buffer zone is taken as the edge line of the plot.

[0030] A reverse buffer is an operation that shrinks a feature inward, creating a smaller area than the original geometry, thus providing unique functionality and applications in GIS analysis. Reverse buffers can be used to define prohibited or excluded zones within a given area. For spatial planning and management, especially in environmental protection and urban planning, defining areas where construction is prohibited or development is restricted is crucial. Reverse buffers can help reduce the impact range of certain areas and clearly indicate areas where certain activities are not permitted. For example, for certain industrial facilities or pollution sources, it is necessary to establish buffer zones around them. Reverse buffers can help us determine areas sufficiently close to the pollution source to exclude contaminated areas.

[0031] For example, in urban planning and design, ecological protection zones are established around lakes, and reverse buffer zones are used to exclude land development activities within a certain distance, which can effectively preserve the natural environment around the water body. Therefore, in the process of calculating urban development along the edges, the use of reverse buffer zones can effectively ensure that the natural environment around the water body is not incorporated into the urban land, thereby preventing the natural environment from being polluted by urban activities.

[0032] Therefore, as another preferred embodiment of the present invention, the method for performing reverse buffer projection in step S4 is as follows: using the edge line of each plot, a reverse buffer projection is performed onto the city boundary line; to ensure that the buffer can fully cover the city boundary line, the buffer distance of the reverse buffer is set to a+200 meters, and the length of the intersection of the reverse buffer and the city boundary line is the corresponding edge development length. To ensure that the length of the intersection of the reverse buffer and the city development boundary line is not counted repeatedly, the corresponding reverse buffers need to be merged to ensure that the intersecting reverse buffers are merged into one polygon data. This process uses the fusion tool (dissolve) in GIS software to merge multiple polygons into one polygon, retaining only one attribute. The length of all reverse buffers on the city development boundary after fusion is the corresponding edge development length.

[0033] The formulas for calculating the total development length SRL for a single edge and the length DRL for both edges in step S5 are as follows: ; ; in, It is the edge development length obtained by the nth reverse buffer, K is the total number of reverse buffers; n=1,2,…,K; G is the total number of overlapping regions in the reverse buffers, G<K, Let G be the edge development length of the m-th overlapping region in the reverse buffer; m = 1, 2, ..., G.

[0034] like Figure 4 The diagram shows the calculation of the total length of the single-sided edge development after the reverse buffer projection obtained by projecting one side of the plot is merged and the overlapping area is removed. Figure 5 This diagram illustrates the calculation of the total length of single-sided edge development after projecting the reverse buffer projection obtained from the double-sided plots, merging them, removing overlapping areas, and further calculating the total length of single-sided edge development.

[0035] The formula for calculating the edge contact ratio LR in step S6 is as follows: .

[0036] like Figure 6 As shown, the urban edge development measurement system based on continuous slice projection technology, which adopts the above method, includes a land use data acquisition module, an edge land use data calculation module, a plot edge line calculation module, a reverse buffer projection module, an edge development length calculation module, and an edge development comparison module. The land use data acquisition module is used to acquire land use data within a defined spatial range; The edge land use data calculation module is used to process the land use data within a set spatial range obtained by the land use data acquisition module using continuous slicing technology, expand outward to obtain a buffer zone, and further use edge development technology to determine the edge land use data; The plot boundary line calculation module is used to preprocess the boundary land use data determined by the boundary land use data calculation module, and project the buffer established by the boundary land use data calculation module and the preprocessed boundary land use data to obtain the plot boundary line. The reverse buffer projection module is used to perform reverse buffer projection to obtain the edge development length; The edge development length calculation module is used to calculate the total edge development length SRL on one side and the edge development length DRL on both sides based on the edge development length obtained by the reverse buffer projection module. The edge-fitting development comparison module is used to calculate the edge-fitting rate LR within a set spatial range, and compare the edge-fitting development in different spatial ranges based on the edge-fitting rate LR. like Figure 7 The diagram shown is a schematic representation of the development boundary of towns bordering Shanghai, measured using the method of this invention. It is evident that the method provided by this invention can accurately and uniformly measure the development boundary of cities bordering Shanghai.

[0037] It should be noted that the sequence numbers of the above embodiments of the present invention are merely for descriptive purposes and do not represent the superiority or inferiority of the embodiments. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, apparatus, article, or method that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, apparatus, article, or method. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, apparatus, article, or method that includes that element.

[0038] Through the above description of the embodiments, those skilled in the art can clearly understand that the methods of the above embodiments can be implemented by means of software plus necessary general-purpose hardware platforms. Of course, they can also be implemented by hardware, but in many cases the former is a better implementation method. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product is stored in a storage medium (such as ROM / RAM, magnetic disk, optical disk) as described above, and includes several instructions to cause a terminal device (which may be a mobile phone, computer, server, or network device, etc.) to execute the methods described in the various embodiments of the present invention.

[0039] Specifically, according to embodiments provided by the present invention, the facial recognition smart lock security management method provided by the present invention can be implemented as a computer software program. For example, a computer program product provided by the present invention includes a computer program carried on a computer-readable medium, the computer program containing program code for performing the methods shown in the flowchart. In such embodiments, the computer program can be downloaded and installed from a network via a communication device, or installed from a storage device, or installed from a ROM. When the computer program is executed by a processing device, it performs the functions defined in the methods of the embodiments of the present disclosure.

[0040] It should be noted that the computer-readable medium provided by this invention can be a computer-readable signal medium or a computer-readable storage medium, or any combination thereof. A computer-readable storage medium can be, for example,—but not limited to—an electrical, magnetic, optical, electromagnetic, infrared, semiconductor system, apparatus, or device, or any combination thereof. More specific examples of a computer-readable storage medium may include, but are not limited to: an electrical connection having one or more wires, a portable computer disk, a hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fiber, portable compact disk read-only memory (CD-ROM), optical storage device, magnetic storage device, or any suitable combination thereof. In this disclosure, a computer-readable storage medium can be any tangible medium containing or storing a program that can be used by or in conjunction with an instruction execution system, apparatus, or device. In this disclosure, a computer-readable signal medium can include a data signal propagated in baseband or as part of a carrier wave, carrying computer-readable program code. Such propagated data signals can take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination thereof. A computer-readable signal medium can be any computer-readable medium other than a computer-readable storage medium, which can send, propagate, or transmit a program for use by or in connection with an instruction execution system, apparatus, or device. The program code contained on the computer-readable medium can be transmitted using any suitable medium, including but not limited to: wires, optical fibers, RF (radio frequency), etc., or any suitable combination thereof.

[0041] The aforementioned computer-readable medium may be included in the aforementioned electronic device; or it may exist independently and not assembled into the electronic device.

[0042] The aforementioned computer-readable medium carries one or more programs that, when executed by the electronic device, cause the electronic device to: acquire at least two Internet Protocol (IP) addresses; send a node evaluation request to a node evaluation device, including at least two IP addresses, wherein the node evaluation device selects an IP address from the at least two IP addresses and returns it; and receive the IP address returned by the node evaluation device; wherein the acquired IP address indicates an edge node in the content delivery network.

[0043] Alternatively, the aforementioned computer-readable medium carries one or more programs that, when executed by the electronic device, cause the electronic device to: receive a node evaluation request including at least two Internet Protocol (IP) addresses; select an IP address from the at least two IP addresses; and return the selected IP address; wherein the received IP address indicates an edge node in the content delivery network.

[0044] Computer program code for performing the operations of this disclosure can be written in one or more programming languages ​​or a combination thereof, including object-oriented programming languages ​​such as Java, Smalltalk, and C++, and conventional procedural programming languages ​​such as the "C" language or similar programming languages. The program code can be executed entirely on the user's computer, partially on the user's computer, as a standalone software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In cases involving remote computers, the remote computer can be connected to the user's computer via any type of network—including a local area network (LAN) or a wide area network (WAN)—or can be connected to an external computer (e.g., via the Internet using an Internet service provider).

[0045] The above description is merely a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A method for measuring urban development along the edge based on continuous slicing projection technology, wherein the method processes land use data based on continuous slicing projection technology, characterized in that... The method includes the following steps: S1: Obtain land use data within a defined spatial range; S2: The land use data within the set spatial range obtained in step S1 is processed using continuous slicing technology, and a buffer zone is obtained by expanding outward. The edge development technology is then used to determine the edge land use data. S3: Preprocess the edge land use data determined in step S2, and project the buffer established in step S2 onto the preprocessed edge land use data to obtain the edge line of the land parcel. S4: Perform reverse buffer projection to obtain the edge development length; S5: Based on the edge development length obtained in step S4, calculate the total edge development length SRL on one side and the edge development length DRL on both sides; S6: Calculate the edge-fitting rate LR within the set spatial range, and compare the edge-fitting development in different spatial ranges based on the edge-fitting rate LR within different spatial ranges.

2. The urban edge development measurement method based on continuous slice projection technology according to claim 1, characterized in that, Step S2 includes the following steps: S2.1: Using the urban boundary line of the study area where the land use data obtained in step S1 is located as the benchmark, and using b meters as the slicing step unit, gradually expand outward to establish a buffer zone, and establish a buffer zone of a meters on both sides respectively; where 0 < b ≤ 10, 500 < a ≤ 1000; S2.2: Obtain land use data within a meter range on both sides of the urban boundary line of the study area, and define all land use data within this range as edge-bound land use data.

3. The urban edge development measurement method based on continuous slice projection technology according to claim 2, characterized in that, Step S3 includes the following steps: S3.1: Traverse the map within the a-meter buffer using the data of the plot images, determine the shape of the plots based on the plot images, find plots that differ greatly from the city boundary line, or have irregular shapes and are projected onto the development length of the boundary line, and statistically record such plots. S3.2: For the land use data of the plots statistically recorded in step S3.1, the surface elements are modified according to the irregularity type of the plots based on the merging and surface segmentation technology in GIS; S3.3: Extract the parcel boundary lines of the parcel image after it has been modified by the polygon features within the buffer.

4. The urban edge development measurement method based on continuous slice projection technology according to claim 3, characterized in that, Step S3.2, based on merging and polygon segmentation techniques in GIS, modifies polygon features according to the irregularity type of the land parcel, including the following steps: S3.21: For a set of seed points within the buffer zone of a meters, form a seed point set G, and define the line connecting any two points p and q in the seed point set G as the equidistant line for region division. : ;in, , Let x and y be the x and y coordinates of point p in the GIS system, respectively. , These are the x and y coordinates of point q in the GIS system, respectively. S3.22: Following the method in step S3.21, traverse any two points in the seed point set G within the buffer zone of a meters to form multiple equidistant lines, dividing the land space within the buffer zone of a meters into multiple polygonal regions; S3.23: For any two polygonal regions A and B among multiple polygonal regions, calculate the intersecting and overlapping areas between them. and the merged area between the two. : ； ； in, Modify the topological area of ​​the face features in polygon region A; Modify the topological area of ​​the face features in polygon region B; S3.24: Using the method described in step S3.23, traverse the multiple polygonal regions obtained in step S3.22, and calculate the buffer S where the boundaries of the different length buffers intersect: ; The number of combinations calculated in step S3.23 is the number of combinations of any two polygon regions from the plurality of polygon regions formed in step S3.

22. M represents the total number of polygonal regions formed in step S3.

22. This is to calculate the number of combinations of selecting any two polygons from M polygon regions and calculating their areas.

5. The urban edge development measurement method based on continuous slice projection technology according to claim 4, characterized in that, In step S3.23, the topological area of ​​the polygon region A is modified by calculating the face features. Modify the topological area of ​​the polygon region B's face features. The calculation formulas are as follows: ； ； in, This represents the total number of seed points within polygonal region A; , Let x and y be the x and y coordinates of the i-th seed point within polygonal region A, respectively. , Let x and y be the x and y coordinates of the (i+1)th seed point within polygonal region A, respectively. , Let x and y be the x and y coordinates of the first seed point within polygonal region A, respectively. , These are the first two terms within polygonal region A. The x and y coordinates of each seed point; This represents the total number of seed points within the polygonal region B. , Let x and y be the x and y coordinates of the j-th seed point within polygonal region B, respectively. , Let x and y be the x and y coordinates of the (j+1)th seed point within polygon region B, respectively. , Here, x and y are the x and y coordinates of the first seed point within polygon region B, respectively. , These are the first two terms within polygonal region B. The x and y coordinates of each seed point.

6. The urban edge development measurement method based on continuous slice projection technology according to claim 5, characterized in that, The calculation method for the edge line of the plot in step S3.3 is to overlay and intersect the buffer zone of a meter with the preprocessed edge land use data to obtain the tangent of the buffer zone S that intersects the boundaries of buffer zones of different lengths, and extract the length of the longest buffer zone tangent, which is recorded as the edge line of the plot.

7. The method for measuring urban development along land use boundaries based on land use zones according to claim 3, characterized in that, The method for performing reverse buffer projection in step S4 is as follows: using the edge line of each plot, a reverse buffer projection is made onto the city boundary line; the buffer distance of the reverse buffer is set to a+200 meters, and the length of the intersection between the reverse buffer and the city boundary line is the corresponding edge development length.

8. The urban edge development measurement method based on continuous slice projection technology according to claim 3, characterized in that, The formulas for calculating the total development length SRL of a single-sided edge and the length DRL of both sides in step S5 are as follows: ； ； in, It is the edge development length obtained by the nth reverse buffer, K is the total number of reverse buffers; n=1,2,…,K; G is the total number of overlapping regions in the reverse buffers, G<K, Let G be the edge development length of the m-th overlapping region in the reverse buffer; m = 1, 2, ..., G.

9. The urban edge development measurement method based on continuous slice projection technology according to claim 8, characterized in that, The formula for calculating the edge contact ratio LR in step S6 is as follows: 。 10. A city edge-fitting development measurement system based on continuous slice projection technology using any of the methods described in claims 1-9, characterized in that, The system includes a land use data acquisition module, a land use data calculation module, a land parcel edge line calculation module, a reverse buffer projection module, an edge development length calculation module, and an edge development comparison module. The land use data acquisition module is used to acquire land use data within a set spatial range; The edge-adjacent land use data calculation module is used to process the land use data within a set spatial range obtained by the land use data acquisition module using continuous slicing technology, expand outward to obtain a buffer zone, and further use edge-adjacent development technology to determine the edge-adjacent land use data. The plot boundary line calculation module is used to preprocess the boundary land use data determined by the boundary land use data calculation module, and project the buffer established by the boundary land use data calculation module onto the preprocessed boundary land use data to obtain the plot boundary line. The reverse buffer projection module is used to perform reverse buffer projection to obtain the edge development length; The edge development length calculation module is used to calculate the total edge development length SRL on one side and the edge development length DRL on both sides based on the edge development length obtained by the reverse buffer projection module. The edge-fitting development comparison module is used to calculate the edge-fitting rate LR within the set spatial range, and compare the edge-fitting development in different spatial ranges based on the edge-fitting rate LR in different spatial ranges.

Images