A public service facility service area division method, device, equipment and medium

Abstract

The application discloses a public service facility service area division method, device, equipment and medium, which obtains address information of all public service facilities in a to-be-divided service area in a research range, and constructs a public service facility space point according to the address information of the public service facilities. Then, navigation schemes of all resident travel OD pairs are obtained respectively, and the navigation schemes of all resident travel OD pairs associated with the same public service facility space point are integrated to construct a travel cost grid corresponding to the public service facility space point. Finally, the travel cost grids of public service facility space points of the same level are fused to obtain the service area range of each public service facility level, so that the city public service facility service area can be scientifically and reasonably demarcated.

Description

Technical Field

[0001] This invention relates to the field of urban planning technology, and in particular to a method, apparatus, equipment and medium for dividing service areas of public service facilities. Background Technology

[0002] There is a considerable amount of research in China on delineating the service areas and service scope of urban public service facilities. Traditional methods include buffer zone delineation, Thiessen polygon delineation, network analysis and identification, and two-step moving search.

[0003] The buffer zone delineation method uses the public service facility as the center and its service radius to define a circular buffer zone, which serves as the service area of ​​the public service facility. Due to its simplicity and convenience, this method is commonly used for delineating urban public service facility service areas. However, it does not consider the accessibility of the road network or physical barriers, and the delineation results often do not match the actual travel distances of residents.

[0004] The Thiessen polygon delineation method uses the principle of minimizing the Euclidean distance to nearby public service facilities to define the service areas of these facilities. This method only spatially allocates public service facilities without considering whether the defined service areas meet the requirements of public service facility planning.

[0005] Network analysis and identification methods construct road network traffic models, consider the travel costs of residents using public service facilities according to road traffic theory, and define costs based on a certain threshold. However, this method requires significant time to construct road network traffic patterns and suffers from issues such as the model information not fully reflecting reality and timeliness.

[0006] The two-step mobile search method determines the supply-demand ratio by moving back and forth between public service facility supply points and resident demand points, and identifies the service area of ​​public service facilities by setting travel time thresholds. This method involves a huge computational burden when dealing with multiple public service facility supply points and multiple resident demand points using Cartesian products, making it unsuitable for large-scale public service facility planning.

[0007] The above methods all define the scope of public service areas based on theoretical spatial distance and coverage. However, they have shortcomings such as incomplete consideration of factors, long modeling time, and overly complex calculations. Furthermore, they fail to consider the actual distance perception of users of public service facilities due to actual road traffic conditions, which can easily lead to discrepancies between theoretical and actual distances. This results in situations where public service facilities are too large, making them inconvenient to use, or too small, leaving facilities idle, thus causing an unreasonable layout of urban public service facilities. Summary of the Invention

[0008] This invention provides a method, apparatus, equipment, and medium for delineating service areas of public service facilities, which can more conveniently, scientifically, and rationally delineate service areas of urban public service facilities.

[0009] The first aspect of this invention provides a method for dividing service areas of public service facilities, comprising:

[0010] Obtain the address information of all public service facilities within the research area to be divided into service areas, and construct the spatial points of public service facilities based on the address information of the public service facilities;

[0011] The navigation schemes for each resident travel OD pair are obtained separately, and the navigation schemes for all resident travel OD pairs associated with the same public service facility spatial point are integrated to construct a travel cost grid corresponding to the public service facility spatial point; wherein, the resident travel OD pair is formed with the public service facility spatial point as the starting point and the sampling points established within the research scope as the ending point, and there are multiple sampling points established within the research scope.

[0012] By integrating the travel cost grids of public service facility locations of the same level, the service area range of each public service facility level is obtained.

[0013] A second aspect of the present invention provides a service area delineation device for public service facilities, comprising:

[0014] The spatial point construction module is used to obtain the address information of all public service facilities in the service area to be divided within the research scope, and to construct the spatial points of public service facilities based on the address information of the public service facilities;

[0015] The spatial point grid construction module is used to obtain the navigation schemes of each resident travel OD pair separately, and integrate the navigation schemes of all resident travel OD pairs associated with the same public service facility spatial point to construct the travel cost grid corresponding to the public service facility spatial point; wherein, the resident travel OD pair is formed with the public service facility spatial point as the starting point and the sampling points established within the research scope as the ending point, and there are multiple sampling points established within the research scope;

[0016] The service area delineation module is used to merge the travel cost grids of public service facility spatial points of the same level to obtain the service area of ​​each public service facility level.

[0017] A third aspect of the present invention provides a terminal device, the terminal device including a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, wherein the processor executes the computer program to implement the public service facility service area division method provided in the first aspect above.

[0018] A fourth aspect of the present invention provides a storage medium comprising a stored computer program, wherein, when the computer program is executed, it controls the device where the storage medium is located to perform the public service facility service area division method provided in the first aspect above.

[0019] Compared with the prior art, the public service facility service area division method, apparatus, equipment and storage medium provided by the present invention have the following beneficial effects:

[0020] This invention acquires the address information of all public service facilities within the research area to be delineated as service zones, and constructs spatial points for these facilities based on this address information. Then, it acquires navigation schemes for each resident's origin-destination (OD) pair, and integrates the navigation schemes for all resident OD pairs associated with the same public service facility spatial point to construct a travel cost grid corresponding to that public service facility spatial point. Finally, it merges the travel cost grids of public service facility spatial points of the same level to obtain the service area range for each public service facility level. This solution, by providing route navigation for resident travel plans related to public service facilities, can assess the characteristics of resident travel costs in real time based on the current road network structure and traffic volume, thereby accurately reflecting the road capacity and actual service coverage during the public service facility delineation period. This allows for the scientific and rational delineation of urban public service facility service areas. Attached Figure Description

[0021] Figure 1 This is a flowchart illustrating an embodiment of the public service facility service area division method provided by the present invention;

[0022] Figure 2 This is an example diagram of a spatial point dataset of medical public service facilities provided in an embodiment of the present invention;

[0023] Figure 3 This is an example diagram of a Python script for obtaining the spatial location of public service facilities provided in an embodiment of the present invention;

[0024] Figure 4 This is an example diagram of a Python script for converting the spatial location coordinates of public service facilities, provided in an embodiment of the present invention.

[0025] Figure 5 This is an example diagram of a resident travel navigation scheme for reaching Meizhou People's Hospital provided in an embodiment of the present invention;

[0026] Figure 6 This is an example image of the travel cost raster dataset of Meizhou People's Hospital provided in an embodiment of the present invention;

[0027] Figure 7This is an example image of the Meizhou City Central Urban Area Travel Cost Sampling Point Dataset provided in this embodiment of the invention;

[0028] Figure 8 This is an example diagram of a Python script for obtaining travel cost data by calling the map route planning API, provided in an embodiment of the present invention.

[0029] Figure 9 This is an example image of the raster dataset of travel costs for municipal-level medical public service facilities in Meizhou City, provided in an embodiment of the present invention.

[0030] Figure 10 This is an example diagram of a district-level medical public service facility service area in the central urban area of ​​Meizhou City, provided by an embodiment of the present invention. Detailed Implementation

[0031] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0032] See Figure 1 , Figure 1 This is a flowchart illustrating an embodiment of the public service facility service area division method provided by the present invention.

[0033] The method for dividing service areas of public service facilities provided in this embodiment of the invention includes steps S11 to S13:

[0034] Step S11: Obtain the address information of all public service facilities in the service area to be divided within the research scope, and construct the spatial points of public service facilities based on the address information of the public service facilities.

[0035] Step S12: Obtain the navigation schemes for each resident travel OD pair, and integrate the navigation schemes for all resident travel OD pairs associated with the same public service facility spatial point to construct a travel cost grid corresponding to the public service facility spatial point; wherein, the resident travel OD pair is formed with the public service facility spatial point as the starting point and the sampling points established within the research scope as the ending point, and there are multiple sampling points established within the research scope.

[0036] Step S13: Merge the travel cost grids of public service facility spatial points of the same level to obtain the service area range of each public service facility level.

[0037] Based on the technical solution provided by this invention, the invention obtains the address information of all public service facilities within the research area to be divided into service zones, and constructs spatial points of public service facilities according to the address information. Then, it obtains the navigation schemes for each resident travel origin-destination (OD) pair, and integrates the navigation schemes of all resident travel OD pairs associated with the same public service facility spatial point to construct a travel cost grid corresponding to the public service facility spatial point. Finally, it merges the travel cost grids of public service facility spatial points of the same level to obtain the service area range of each public service facility level. This invention, by providing route navigation for resident travel plans related to public service facilities, can assess the characteristics of resident travel costs in real time based on the current road network structure and traffic volume, thereby accurately reflecting the road capacity and actual service coverage range during the public service facility delineation period, and thus enabling the scientific and reasonable delineation of urban public service facility service areas.

[0038] In an alternative implementation, prior to step S11, "obtaining the address information of all public service facilities within the research area to be divided into service zones," the method further includes:

[0039] Obtain descriptive text data of public service facilities within the research area to be divided into service areas, and extract the address information of the public service facilities from the descriptive text data to obtain the address information of the public service facilities.

[0040] For example, descriptive text data for public service facilities includes the name, address, and level information of the public service facility. The address should include complete information such as city (prefecture), district (county), town (street), and house number. The level information should include information that can be distinguished into four levels according to the "Urban Public Service Facilities Planning Standard GB50442 (Revised) (Draft for Comments)": city level, district (county) level, street (town) level, and community (village) level. For example, a city-level public service facility is represented as P. city ={p c1 p c2 p c3 ... p cm District-level public service facilities are represented by P. county ={p co1 p co2 p co3 ... p con Township-level public service facilities are represented by P. town ={p t1 p t2 p t3 ... p ti Community-level public service facilities are represented by P. village ={p v1 p v2 pv3 ... p vj Each type of public service facility includes m, n, i, and j public service facility points p.

[0041] In one optional implementation, step S11, "constructing spatial points of public service facilities based on the address information of the public service facilities," specifically includes:

[0042] Based on the address information of the public service facilities, spatial location data of each public service facility is obtained through a map API interface;

[0043] The latitude and longitude in the spatial location data of the public service facilities are converted into geographic coordinate system data that conforms to the WGS84 standard to obtain the converted spatial location of each public service facility.

[0044] Using the XY conversion function of ArcGIS software, spatial points of public service facilities are constructed based on the converted spatial locations.

[0045] For example, a Python script can be written to convert the latitude and longitude of the spatial location data of the above public service facilities into geographic coordinate system data that conforms to the WGS84 standard, so as to obtain the converted spatial location and ensure that there is no address offset between the spatial points of each public service facility.

[0046] In an optional implementation, the "navigation scheme for each resident's travel OD pair" in S12 is obtained in the following manner:

[0047] Obtain the research area shapefile and use ArcGIS software's fishing net creation tool to create several sampling points within the research area;

[0048] Starting from the public service facility locations and ending at the sampling points, OD pairs for resident travel are formed.

[0049] The map route planning API is invoked, and navigation is performed on each of the aforementioned resident travel OD pairs according to the navigation transportation mode corresponding to the preset level of public service facility spatial points, so as to obtain the navigation scheme for each of the aforementioned resident travel OD pairs.

[0050] In this embodiment of the invention, the navigation and transportation mode corresponding to the level of the public service facility spatial point can be set as: municipal-level, district (county)-level public service facility P city P county Using driving navigation mode, public service facilities at the street (town) level and community (village) level P town P villageUsing a pedestrian navigation mode, the route planning functions of Gaode Maps, Baidu Maps, or Tencent Maps are invoked during implementation. The aforementioned resident travel origin-destination (OD) pairs are sequentially input, and the corresponding navigation mode is selected to obtain navigation schemes for each resident travel OD pair. The navigation scheme for each resident travel OD pair includes the longitude (lon), latitude (lat), and OD pair travel time duration information of the sampling point, such as the location of municipal public service facilities (p). c1 Resident travel navigation information collection r c1 The following table shows (where s1 is the sampling point):

[0051] lon lat duration <![CDATA[s1 longitude]]> <![CDATA[s1 latitude]]> <![CDATA[from s1 to p c1 travel time]]> <![CDATA[s2 longitude]]> <![CDATA[s2 latitude]]> <![CDATA[From s2 to p c1 travel time]]> …… …… …… <![CDATA[s n Longitude <![CDATA[s n Latitude <![CDATA[From s n to p c1 travel time]]>

[0052] It is worth noting that, unlike the traditional method of delineating service areas for public service facilities, this embodiment of the invention uses network data transmission technology to obtain route planning and travel scheme data provided by map APIs such as Gaode Maps. This allows for real-time assessment of residents' travel costs based on the current road network structure and traffic volume, accurately reflecting the road capacity and actual service coverage of public service facilities during the designated work period, thus exhibiting strong timeliness.

[0053] Furthermore, traditional methods estimate residents' travel costs by constructing road network models. The road network data used to construct these models is typically obtained through manual vectorization or internet sharing, resulting in insufficient data availability and reliability. Moreover, this indirect estimation of travel costs does not reflect actual travel time, leading to poor data timeliness. In contrast, this invention utilizes map navigation data. Various internet map service providers offer free API interfaces, allowing access to current travel cost data for any region in China, ensuring good data availability and timeliness.

[0054] Furthermore, traditional road network segmentation methods suffer from data quality variations due to differences in regional economic levels and surveying technology, leading to tedious and time-consuming quality verification processes such as road network data topology checks. This invention utilizes travel cost data from map navigation data, providing unbiased coverage across China, saving time on data verification and road network modeling, and offering ease of operation and good portability.

[0055] In an optional implementation, S12, "integrating navigation schemes for all resident travel OD pairs associated with the same public service facility spatial point and constructing a travel cost grid corresponding to the public service facility spatial point," specifically includes:

[0056] By integrating the navigation schemes of all resident travel origin-destination (OD) pairs associated with the same public service facility spatial point, a navigation scheme for all resident travel corresponding to each public service facility spatial point is obtained. The resident travel navigation scheme includes OD pair travel time information.

[0057] Using the XY conversion function of ArcGIS software, all resident travel navigation schemes corresponding to the spatial points of the public service facilities are converted into travel cost points, and the travel cost points of the spatial points of the public service facilities are converted into individual travel cost grids using Kriging interpolation; wherein, each element in the travel cost points represents the OD pair travel time corresponding to each resident travel navigation scheme.

[0058] For example, if sampling point s1 reaches the municipal public service facility p c1 The resident travel navigation plan includes information such as travel mode, travel route, travel cost, and travel time. Travel time can be used as a travel cost indicator, and the travel time can be assigned to the duration field attribute value of sampling point s1. Then, the municipal public service facility p... c1 The travel cost lattice is [d1, d2, ..., d n ], d n Represented as sampling point s n Reaching municipal public service facilities p c1 The travel time is calculated; finally, Kriging interpolation is performed based on the duration field attribute values ​​of all sampling points to obtain the city-level public service facility p. c1 The travel cost grid is for municipal public service facilities g c1 .

[0059] In one optional implementation, S13, "integrating the travel cost grid of public service facility spatial points of the same level to obtain the service area range of each public service facility level," specifically includes:

[0060] According to the level of public service facilities, multiple public service facility spatial point travel cost grids of the same level are embedded to obtain the travel cost grids corresponding to each level of public service facilities.

[0061] Import the travel cost raster corresponding to each level of public service facility into ArcGIS software, and use ArcGIS software’s contour tool to generate contour lines with preset contour time intervals.

[0062] Based on the preset time interval contour lines corresponding to each level of public service facilities, select the contour lines corresponding to each level of public service facilities from all the generated contour lines as the service area of ​​each level of public service facilities.

[0063] For example, municipal public service facilities p c1 The travel cost grid is g c1 Each level of public service facility P forms a travel cost raster dataset G, such as a cost raster dataset G consisting of m municipal public service facilities.city ={g c1 g c2 g c3 ... g cm}, will g c1 g c2 g c3 ... g cm By tessellation, a grid of travel costs from any point in the space within the study area to the nearest municipal public service facility is obtained.

[0064] In this embodiment of the invention, the time interval contour lines set for each level of public service facilities can be N times the set contour time interval. For example, if the contour time interval is set to 5 minutes, that is, contour lines are generated every 5 minutes, then for municipal public service facilities, 30-minute contour lines are set; for district-level public service facilities, 20-minute contour lines are set; for street-level public service facilities, 15-minute contour lines are set; and for community-level public services, 10-minute contour lines are set.

[0065] In one implementation, the tiling operation mode used when tiling multiple public service facility spatial point travel cost grids of the same level is MINIMUM.

[0066] Accordingly, embodiments of the present invention also provide a public service facility service area division device, comprising:

[0067] The spatial point construction module is used to obtain the address information of all public service facilities in the service area to be divided within the research scope, and to construct the spatial points of public service facilities based on the address information of the public service facilities;

[0068] The spatial point grid construction module is used to obtain the navigation schemes of each resident travel OD pair separately, and integrate the navigation schemes of all resident travel OD pairs associated with the same public service facility spatial point to construct the travel cost grid corresponding to the public service facility spatial point; wherein, the resident travel OD pair is formed with the public service facility spatial point as the starting point and the sampling points established within the research scope as the ending point, and there are multiple sampling points established within the research scope;

[0069] The service area delineation module is used to merge the travel cost grids of public service facility spatial points of the same level to obtain the service area of ​​each public service facility level.

[0070] In one alternative implementation, the navigation scheme for each resident's travel origin-destination (OD) pair is obtained as follows:

[0071] Obtain the research area shapefile and use ArcGIS software's fishing net creation tool to create several sampling points within the research area;

[0072] Starting from the public service facility locations and ending at the sampling points, OD pairs for resident travel are formed.

[0073] The map route planning API is invoked, and navigation is performed on each of the aforementioned resident travel OD pairs according to the navigation transportation mode corresponding to the preset level of public service facility spatial points, so as to obtain the navigation scheme for each of the aforementioned resident travel OD pairs.

[0074] In one optional implementation, the navigation scheme that integrates all resident travel origin-destination pairs associated with the same public service facility spatial point, and constructs a travel cost grid corresponding to the public service facility spatial point, specifically includes:

[0075] By integrating the navigation schemes of all resident travel origin-destination (OD) pairs associated with the same public service facility spatial point, a navigation scheme for all resident travel corresponding to each public service facility spatial point is obtained. The resident travel navigation scheme includes OD pair travel time information.

[0076] Using the XY conversion function of ArcGIS software, all resident travel navigation schemes corresponding to the spatial points of the public service facilities are converted into a travel cost matrix. Then, Kriging interpolation is used to convert the travel cost matrix of the spatial points of the public service facilities into a single travel cost grid. Each element in the travel cost matrix represents the OD pair travel time corresponding to each resident travel navigation scheme.

[0077] In one optional implementation, the travel cost grid of merging public service facility spatial points of the same level to obtain the service area range of each public service facility level specifically includes:

[0078] According to the level of public service facilities, multiple public service facility spatial point travel cost grids of the same level are embedded to obtain the travel cost grids corresponding to each level of public service facilities.

[0079] Import the travel cost raster corresponding to each level of public service facility into ArcGIS software, and use ArcGIS software’s contour tool to generate contour lines with preset contour time intervals.

[0080] Based on the preset time interval contour lines corresponding to each level of public service facilities, select the contour lines corresponding to each level of public service facilities from all the generated contour lines as the service area of ​​each level of public service facilities.

[0081] In one optional implementation, the tiling operation mode used for tiling the multiple public service facility spatial point travel cost grids of the same level is MINIMUM.

[0082] In one optional implementation, before obtaining the address information of all public service facilities within the research area to be divided into service zones, the method further includes:

[0083] Obtain descriptive text data of public service facilities within the research area to be divided into service areas, and extract the address information of the public service facilities from the descriptive text data to obtain the address information of the public service facilities.

[0084] In one optional implementation, constructing the spatial points of the public service facilities based on their address information specifically includes:

[0085] Based on the address information of the public service facilities, spatial location data of each public service facility is obtained through a map API interface;

[0086] The latitude and longitude in the spatial location data of the public service facilities are converted into geographic coordinate system data that conforms to the WGS84 standard to obtain the converted spatial location of each public service facility.

[0087] Using the XY conversion function of ArcGIS software, spatial points of public service facilities are constructed based on the converted spatial locations.

[0088] It should be noted that the public service facility service area division device provided in this embodiment of the invention is used to execute all the processes and steps of the public service facility service area division method in the above embodiment. The working principles and effects of the two are one-to-one, and will not be elaborated further here.

[0089] It should be noted that the device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs.

[0090] A third aspect of the present invention provides a terminal device, including a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor. When the processor executes the computer program, it implements the public service facility service area division method provided in the above embodiments, such as S11 to S13.

[0091] The terminal device may be a desktop computer, laptop, handheld computer, or cloud server, etc. The terminal device may include, but is not limited to, a processor and a memory.

[0092] The processor can be a Central Processing Unit (CPU), or other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor can be a microprocessor or any conventional processor. The processor is the control center of the terminal device, connecting all parts of the terminal device via various interfaces and lines.

[0093] The memory can be used to store the computer programs and / or modules. The processor, by running or executing the computer programs and / or modules stored in the memory and calling the data stored in the memory, realizes various functions of the public service facility service area division device / terminal equipment. The memory may mainly include a program storage area and a data storage area. The program storage area may store the operating system, at least one application program required for a function (such as sound playback function, image playback function, etc.), etc.; the data storage area may store data created based on the use of the mobile phone (such as audio data, phonebook, etc.). In addition, the memory may include high-speed random access memory, and may also include non-volatile memory, such as hard disk, memory, plug-in hard disk, smart media card (SMC), secure digital (SD) card, flash card, at least one disk storage device, flash memory device, or other volatile solid-state storage device.

[0094] Wherein, if the modules / units integrated by the public service facility service area division device / terminal equipment are implemented as software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, all or part of the processes in the methods of the above embodiments of the present invention can also be implemented by a computer program instructing related hardware. The computer program can be stored in a computer-readable storage medium, and when executed by a processor, it can implement the steps of the various method embodiments described above. The computer program includes computer program code, which can be in the form of source code, object code, executable files, or certain intermediate forms, etc.

[0095] A fourth aspect of this invention provides a storage medium comprising a stored computer program, wherein, when the computer program is executed, it controls the device containing the storage medium to perform a public service facility service area division method as provided in the above embodiments, for example... Figure 1 S11 to S13.

[0096] The storage medium may include: any entity or device capable of carrying the computer program code, recording media, USB flash drive, portable hard drive, magnetic disk, optical disk, computer memory, read-only memory (ROM), random access memory (RAM), electrical carrier signal, telecommunication signal, and software distribution medium, etc.

[0097] The following description uses the classification of medical public service facilities within Meizhou City as an example to illustrate this application.

[0098] Step S1: Collect descriptive text data of medical public service facilities within Meizhou City. The descriptive text data of medical service facilities includes the name, address, and level information of the public service facilities. The address includes complete information such as city (prefecture), district (county), town (street), and house number. The level information includes information that can be distinguished into four levels in accordance with the "Urban Public Service Facilities Planning Standard GB50442 (Revised) (Draft for Comments)"—city level, district (county) level, street (town) level, and community (village) level, as shown in Table 1.

[0099] Table 1

[0100] name address level Meizhou People's Hospital No. 63 Huangtang Road, Meijiang District, Meizhou City City level Meizhou City Meixian District Traditional Chinese Medicine Hospital No. 52, Xianzi Avenue, Meixian District, Meizhou City District level Chengbei Town Health Center No. 155 Huanshi North Road, Meijiang District, Meizhou City Street level Wuliting Health Station No. 54, Wuliting Road, Meijiang District, Meizhou City Community level …… …… ……

[0101] Step S2, construct a spatial point dataset of public service facilities, such as Figure 2 The process of constructing a spatial point dataset for public service facilities includes the following steps:

[0102] Step 201: Using the address information in the descriptive text data of public service facilities from step S1, write a Python script (e.g., ...). Figure 3 Request the map API interface (provided by Gaode Map, Baidu Map and Tencent Map) to obtain the spatial location data of various medical public service facilities in Meizhou City. The spatial location data of each public service facility includes latitude and longitude, name and level information, as shown in Table 2.

[0103] Table 2

[0104] longitude latitude name level 116.103612 24.311194 Meizhou People's Hospital City level 116.097703 24.279291 Meizhou City Meixian District Traditional Chinese Medicine Hospital District level 116.101157 24.327378 Chengbei Town Health Center Street level 116.112166 24.317673 Wuliting Health Station Community level …… …… …… ……

[0105] Step 202: By writing a Python script (such as...) Figure 4 The latitude and longitude coordinates of the spatial location data of the above public service facilities are converted into geographic coordinate system data that conforms to the WGS84 standard to form spatial location data without geographic offset.

[0106] Step 203: Using ArcGIS software's XY conversion function, construct a public service facility spatial point dataset from the above spatial location data without geographical offset, such as... Figure 2 ;

[0107] Step S3: Obtain resident travel navigation information for each medical public service facility location, such as... Figure 5 Construct a raster dataset of travel costs, such as Figure 6 The process of constructing a raster dataset of travel costs specifically includes the following steps:

[0108] Step 301: Obtain the shapefile file for the central urban area of ​​Meizhou City. Using the fishing net creation tool in ArcGIS software, create a 500*500m sampling point dataset within the study area, such as... Figure 7 ;

[0109] Step 302: Starting from the sampling points and ending with the spatial point data of public service facilities in step S2, form OD pairs for resident travel, as shown in Table 3:

[0110] Table 3

[0111] Longitude of point O Latitude of point O Longitude of point D Latitude of point D 116.113612 24.311194 116.103612 24.311194 116.107703 24.279291 116.103612 24.311194 116.091157 24.327378 116.103612 24.311194 116.082166 24.317673 116.103612 24.311194 …… …… …… ……

[0112] Step 303: By writing a Python script (such as...) Figure 8Call the map route planning API (provided by Gaode Map, Baidu Map, and Tencent Map) and, during off-peak traffic hours, sequentially input the OD pairs of residents' travel from the previous step. Among them, the city-level and district-level public service facilities adopt driving navigation mode, and the street-level and community-level public service facilities adopt walking navigation mode. Collect residents' travel navigation plan information for each public service facility point. The residents' travel navigation plan information should include the latitude and longitude of the sampling point and travel time information, as shown in Table 4.

[0113] Longitude of point O Latitude of point O duration 116.113612 24.311194 14377 116.107703 24.279291 7593 116.091157 24.327378 9491 116.082166 24.317673 11290 …… …… ……

[0114] Step 304: Using the XY conversion function of ArcGIS software, the resident travel navigation information of the above-mentioned public service facility points is converted into a travel cost matrix. The values ​​of the travel cost matrix are travel times. Then, using Kriging interpolation, the travel cost matrix is ​​converted into a travel cost raster for individual public service facility spatial points, such as... Figure 6 ;

[0115] Step S4: Integrate the raster data of travel costs for spatial points of public service facilities of the same level, such as... Figure 9 Generate the service area of ​​public service facilities, such as Figure 10 The process of generating the service area boundaries of public service facilities specifically includes the following steps:

[0116] Step 401: Using the single public service facility travel cost grid in step S3, multiple public service facility travel cost grids of the same level are mosaicked according to the public service facility level. Overlapping grids adopt the minimum value of each layer of grids (i.e. MINIMUM mosaicking operation mode) to obtain the travel cost grids of the nearest municipal, district, street and community-level medical public service facilities in the central urban area of ​​Meizhou City.

[0117] Step 402: Using the contour tool in ArcGIS software, import the travel cost raster of medical public service facilities at all levels and generate contour lines at 5-minute intervals;

[0118] Step 403: Select the 30-minute isoline of municipal public service facilities, the 20-minute isoline of district-level public service facilities, the 15-minute isoline of street-level public service facilities, and the 10-minute isoline of community-level public service facilities as the service areas of municipal, district, street, and community-level public service facilities, respectively. Figure 10 .

[0119] The above description represents the preferred embodiments of the present invention. It should be noted that those skilled in the art can make various improvements and modifications without departing from the principles of the present invention, and these improvements and modifications are also considered to be within the scope of protection of the present invention.

Claims

1. A public service facility service area division method characterized by, include: Obtain the address information of all public service facilities within the research area to be divided into service areas, and construct the spatial points of public service facilities based on the address information of the public service facilities; The navigation schemes for each resident travel OD pair are obtained separately, and the navigation schemes for all resident travel OD pairs associated with the same public service facility spatial point are integrated to construct a travel cost grid corresponding to the public service facility spatial point; wherein, the resident travel OD pair is formed with the public service facility spatial point as the starting point and the sampling points established within the research scope as the ending point, and there are multiple sampling points established within the research scope. By integrating the travel cost grid of public service facility spatial points of the same level, the service area of ​​each public service facility level is obtained, specifically including: According to the level of public service facilities, multiple public service facility spatial point travel cost grids of the same level are mosaicked to obtain the travel cost grid corresponding to each public service facility level; the mosaicking operation mode used for multiple public service facility spatial point travel cost grids of the same level is MINIMUM. Import the travel cost raster corresponding to each level of public service facility into ArcGIS software, and use ArcGIS software’s contour tool to generate contour lines with preset contour time intervals. Based on the preset time interval contour lines corresponding to each level of public service facilities, select the contour lines corresponding to each level of public service facilities from all the generated contour lines as the service area of ​​each level of public service facilities.

2. The method for dividing service areas of public service facilities as described in claim 1, characterized in that, The navigation schemes for each resident's travel origin-destination (OD) pair are obtained in the following way: Obtain the research area shapefile and use ArcGIS software's fishing net creation tool to create several sampling points within the research area; Starting from the public service facility locations and ending at the sampling points, OD pairs for resident travel are formed. The map route planning API is invoked, and navigation is performed on each of the aforementioned resident travel OD pairs according to the navigation transportation mode corresponding to the preset level of public service facility spatial points, so as to obtain the navigation scheme for each of the aforementioned resident travel OD pairs.

3. The method for dividing service areas of public service facilities as described in claim 1, characterized in that, The navigation scheme that integrates all resident travel origin-destination pairs associated with the same public service facility spatial point, and constructs a travel cost grid corresponding to the public service facility spatial point, specifically includes: By integrating the navigation schemes of all resident travel origin-destination (OD) pairs associated with the same public service facility spatial point, a navigation scheme for all resident travel corresponding to each public service facility spatial point is obtained. The resident travel navigation scheme includes OD pair travel time information. Using the XY conversion function of ArcGIS software, all resident travel navigation schemes corresponding to the spatial points of the public service facilities are converted into a travel cost matrix. Then, Kriging interpolation is used to convert the travel cost matrix of the spatial points of the public service facilities into a single travel cost grid. Each element in the travel cost matrix represents the OD pair travel time corresponding to each resident travel navigation scheme.

4. The method for dividing service areas of public service facilities as described in claim 1, characterized in that, Before obtaining the address information of all public service facilities in the service area to be divided within the research scope, the method further includes: Obtain descriptive text data of public service facilities within the research area to be divided into service areas, and extract the address information of the public service facilities from the descriptive text data to obtain the address information of the public service facilities.

5. The method for dividing service areas of public service facilities as described in claim 1, characterized in that, The step of constructing spatial points for public service facilities based on their address information specifically includes: Based on the address information of the public service facilities, spatial location data of each public service facility is obtained through a map API interface; The latitude and longitude in the spatial location data of the public service facilities are converted into geographic coordinate system data that conforms to the WGS84 standard to obtain the converted spatial location of each public service facility. Using the XY conversion function of ArcGIS software, spatial points of public service facilities are constructed based on the converted spatial locations.

6. A service area delineation device for public service facilities, characterized in that, include: The spatial point construction module is used to obtain the address information of all public service facilities in the service area to be divided within the research scope, and to construct the spatial points of public service facilities based on the address information of the public service facilities; The spatial point grid construction module is used to obtain the navigation schemes of each resident travel OD pair separately, and integrate the navigation schemes of all resident travel OD pairs associated with the same public service facility spatial point to construct the travel cost grid corresponding to the public service facility spatial point; wherein, the resident travel OD pair is formed with the public service facility spatial point as the starting point and the sampling points established within the research scope as the ending point, and there are multiple sampling points established within the research scope; The service area delineation module is used to merge the travel cost grids of public service facility spatial points of the same level to obtain the service area range for each public service facility level, specifically including: According to the level of public service facilities, multiple public service facility spatial point travel cost grids of the same level are mosaicked to obtain the travel cost grid corresponding to each public service facility level; the mosaicking operation mode used for multiple public service facility spatial point travel cost grids of the same level is MINIMUM. Import the travel cost raster corresponding to each level of public service facility into ArcGIS software, and use ArcGIS software’s contour tool to generate contour lines with preset contour time intervals. Based on the preset time interval contour lines corresponding to each level of public service facilities, select the contour lines corresponding to each level of public service facilities from all the generated contour lines as the service area of ​​each level of public service facilities.

7. A terminal device, characterized in that, It includes a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, wherein the processor executes the computer program to implement the public service facility service area division method as described in any one of claims 1 to 5.

8. A storage medium, characterized in that, The storage medium includes a stored computer program, wherein, when the computer program is executed, it controls the device where the storage medium is located to perform the public service facility service area division method as described in any one of claims 1 to 5.

Images