An online interactive public service facility accessibility analysis method

Abstract

The application discloses an online interactive public service facility accessibility analysis method, relates to the technical fields of city planning, tourism planning, natural language processing and city traffic, and first encapsulates the improved grid cost accessibility calculation as an independent ArcToolbox in ArcGIS software.Secondly, the response function design of the server-side accessibility calculation is carried out.Thirdly, the function design of the client-side accessibility calculation parameter sending request is carried out.Meanwhile, the function design of the interactive source adjustment of the client-side is also realized.Finally, after the client-side obtains the result returned by the server-side, the result is parsed, drawn and counted, and the calculation conversion function of different travel modes is realized.The application can dynamically adjust the source of the grid accessibility analysis, batch-give the independent accessibility result of each source in the grid accessibility calculation and the accessibility result of different travel modes in the client-side, and directly converts the accessibility result of different travel modes in the client-side, and the like.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of urban planning, GIS, and particularly to an online interactive public service facility accessibility analysis method. BACKGROUND

[0002] Accessibility calculation is an important technology in the research of GIS algorithms, and is widely used in application fields such as urban planning, land management, and tourism development. Among them, grid accessibility calculation has the characteristics of rapidness, convenience, and the calculation results can reveal the overall characteristics of the study area, and has always been the focus of accessibility application. However, the current grid accessibility calculation lacks online calculation methods, and generally users operate on GIS software. This brings an inconvenient problem in real application. For example, when a research group conducts grid accessibility calculation on the same study area, if a good collaborative work mode cannot be formed, different users may need to operate the basic data or need to constantly adjust the spatial data according to the application requirements, and then perform calculation after data exchange between users.

[0003] Further, in many application researches, the source or cost of grid accessibility calculation needs to be dynamically adjusted, so the software system needs strong interactive operation, and can quickly feedback the calculation results according to the user-adjusted spatial data results. However, the existing software often ignores the dynamic interactive operation function design of the grid accessibility calculation tool. At the same time, in the application of urban planning, the source of grid accessibility calculation is often for multiple spatial facility points, such as calculating the overall 15-minute walk accessibility of multiple parks in the city. However, in actual planning analysis, after obtaining the overall accessibility calculation result, the accessibility of each park needs to be further calculated and counted to deeply analyze the coverage range of the park accessibility. SUMMARY

[0004] The technical problem to be solved by the present application is to overcome the shortcomings of the prior art and provide an online interactive public service facility accessibility analysis method. The present application can realize online interactive public service facility accessibility calculation without installing GIS software on the client side, and can dynamically adjust the source of grid accessibility calculation, batch give the independent accessibility results of each source in grid accessibility calculation, and directly convert the accessibility results of different travel modes on the client side. It has important technical method innovation significance for current public service facility accessibility analysis.

[0005] The present application adopts the following technical solutions to solve the above technical problems:

[0006] According to the online interactive public service facility accessibility analysis method provided by the present application, the following steps are included:

[0007] Step 1, improve the grid cost accessibility calculation, and encapsulate the improved grid cost accessibility calculation as an independent ArcToolbox tool in ArcGIS software; the specific steps are as follows:

[0008] Step 1.1, dynamically generating road cost grid data based on online map navigation data, and taking the road speed provided by the online map navigation as the attribute value of the road cost grid data;

[0009] The closed road network and the priority road through rule are added to the grid cost accessibility calculation;

[0010] Step 1.2, based on step 1.1, encapsulating the ArcToolbox tool of the grid cost accessibility calculation by using the API of the secondary development of ArcGIS;

[0011] Step 2, deploying the ArcToolbox tool in step 1 to the server side, and designing the accessibility calculation response function of the server side;

[0012] Step 3, designing the function of sending request of the client accessibility calculation parameters;

[0013] Step 4, designing the function of realizing interactive adjustment of the source of the client, and re-performing new grid cost accessibility calculation;

[0014] Step 5, after the client obtains the new grid cost accessibility calculation result of step 4 returned by the server side, the result is parsed, drawn and counted.

[0015] As a further optimization scheme of the online interactive public service facility accessibility analysis method, step 2 is specifically as follows:

[0016] Step 2.1, deploying the ArcToolbox tool in step 1 to the server side, and publishing it as a geographic processing service, and setting the geographic processing service as an asynchronous processing mode;

[0017] Step 2.2, the server side exposes the parameters of the geographic processing service in step 2.1 to the client.

[0018] As a further optimization scheme of the online interactive public service facility accessibility analysis method, step 3 is specifically as follows:

[0019] Step 3.1, encapsulating the parameters provided in step 2.2 into a whole request task package;

[0020] Step 3.2, submitting the encapsulated result to the server for calculation task.

[0021] As a further optimization of the online interactive public service facility accessibility analysis method described in the application, step 4 is specifically as follows:

[0022] Step 4.1, convert the source data of the pre-stored grid cost accessibility calculation on the server side into GeoJson format data;

[0023] Step 4.2, the client uses the provided function of editing GeoJson format data to edit the source data of the grid cost accessibility calculation, and saves the result as GeoJson format data;

[0024] Step 4.3, the client saves the edited GeoJson format data as new source data for grid cost accessibility calculation;

[0025] Step 4.4, the client re-sends the new source data for grid cost accessibility calculation obtained in step 4.3, and the server side also performs new grid cost accessibility calculation.

[0026] As a further optimization of the online interactive public service facility accessibility analysis method described in the application, step 5 is specifically as follows:

[0027] Step 5.1, the client parses, draws and counts the new grid cost accessibility calculation result returned by the server in step 4.4;

[0028] Step 5.2, at the same time, the new grid cost accessibility calculation result returned in step 4.4 is spatially overlaid and analyzed with residential land to judge the coverage of the accessibility service of residential land.

[0029] As a further optimization of the online interactive public service facility accessibility analysis method described in the application, the online interactive public service facility accessibility analysis method further includes: implementing a calculation conversion function for different travel modes: according to the conversion formula of the time required for walking and cycling, the accessibility coverage range under different traffic travel modes is counted.

[0030] As a further optimization of the online interactive public service facility accessibility analysis method described in the application, step 1 is based on the secondary development of ArcGIS 10.1~10.8 software version to package a custom grid cost accessibility calculation ArcToolbox tool.

[0031] As a further optimization of the online interactive public service facility accessibility analysis method described in the application, the method for dynamically generating road cost grid data based on online map navigation data in step 1.1 includes the following steps:

[0032] Step 1.1.1: Set up a batch list of OD point pairs for online map navigation queries. Each OD point pair includes the geospatial coordinates of the origin and destination.

[0033] Step 1.1.2: Use the travel API in the online map navigation interface to obtain the queried road segments and the speed of each OD point pair; at the same time, assign the obtained road segment speeds to the corresponding road segments.

[0034] Step 1.1.3: Convert the road segments obtained through online map navigation into road cost raster data. The value of the road cost raster is formed according to the speed assigned to each road segment.

[0035] As a further optimization of the online interactive public service facility accessibility analysis method described in this invention, when the server exposes the parameters of a geoprocessing service to the client in step 2.2, the parameters include: different types of road networks, water systems, study area scope, and source data for raster cost accessibility calculation.

[0036] As a further optimization of the online interactive public service facility accessibility analysis method described in this invention, the conversion formula for the time required for walking and cycling is that the speed of a bicycle is 5 times the walking speed.

[0037] Compared with the prior art, the present invention, employing the above technical solution, has the following technical effects:

[0038] (1) This invention provides an online interactive public service facility accessibility analysis method, which can eliminate the need to install GIS software, dynamically adjust the sources of raster accessibility analysis, batch provide independent accessibility results for each source in raster accessibility calculation, and directly convert accessibility results of different travel modes on the client side.

[0039] (2) Based on the principle of GIS raster accessibility spatial analysis, this invention proposes an online interactive public service facility accessibility analysis method. Attached Figure Description

[0040] Figure 1 This is a schematic diagram of the overall process of the present invention.

[0041] Figure 2 This is a schematic diagram of the overall architecture for implementing the present invention.

[0042] Figure 3 This is a schematic diagram of the cost grid mosaic sequence.

[0043] Figure 4 This is a schematic diagram of the accessibility analysis function panel.

[0044] Figure 5 is a client result drawing schematic diagram.

[0045] Figure 6 is a statistical analysis chart display schematic diagram.

[0046] Figure 7 is a calculation method result comparison schematic diagram.

[0047] Figure 8 is a reachability analysis result schematic diagram.

[0048] Figure 9 is a reachability analysis result schematic diagram after scheme adjustment. DETAILED DESCRIPTION

[0049] The technical solutions of the present application will be further described in detail below in combination with the drawings:

[0050] The present research provides an online interactive grid reachability calculation method, and the method has the following seven characteristics: (1) the client uses a WEB browser for operation, and no GIS software needs to be installed to realize the calculation of grid reachability. (2) The user does not need to calculate the cost layer of grid reachability separately, and the user only needs to input the data information for calculation, which greatly reduces the workload of the user in use. (3) The user can dynamically adjust the source of grid reachability calculation, and the new grid reachability calculation result can be directly fed back after adjustment. (4) Batch grid reachability results of each source in the calculation can be given. (5) The online map navigation road information is integrated into the road cost network to enhance the scientificity of the analysis result. (6) The time equivalent conversion calculation of different travel modes is carried out to provide a quick reference basis for different travel mode reachability analysis. (7) The scenario setting of closed road travel and priority road passing is considered in the grid cost calculation. The grid reachability in the present application refers to the grid cost reachability.

[0051] The method is divided into server side and client side two parts in the overall idea. The server side is responsible for the actual calculation task of grid reachability, and the result is fed back to the client. The client writes analysis parameters, sends a calculation request to the server side, and draws the calculation result. The technical invention focuses on using a web map server to calculate grid reachability, so as to achieve the effect of online public service facility reachability calculation and analysis.

[0052] Step 1) refer to the attached Figure 1 The overall process of the present application is as follows:

[0053] Refer to the attached Figure 2Fig. 1 is a schematic diagram of an overall architecture for implementing the present application. The server end connects different database resources, thereby providing different data sources for online public service facility accessibility calculation, including online collected road network data. The map server is responsible for performing specific calculation, which accepts requests and parameters sent by the client end and feeds back the results to the client end. Meanwhile, it also calls data from the database resources to participate in the online public service facility accessibility calculation. The client end needs to provide parameters for each accessibility calculation. Since the parameters can be dynamically adjusted and many basic operations are encapsulated in the present application, only the core parameter adjustment of the client end is needed. Only in this way can better interactive accessibility calculation be performed.

[0054] The present case selects ArcGIS server for online public service facility accessibility calculation. Therefore, the present application first encapsulates the improved grid cost accessibility calculation as an independent ArcToolbox in the ArcGIS software to achieve certain function encapsulation, thereby allowing the client end to dynamically adjust the core parameters to achieve the purpose of better online interactive public service facility accessibility analysis.

[0055] Step 1.1) encapsulate the custom grid cost accessibility calculation ArcToolbox tool by using the API for secondary development of ArcGIS. This step is based on the secondary development of ArcGIS 10.1-10.8 software version to encapsulate the custom grid cost accessibility calculation ArcToolbox tool.

[0056] Step 1.2) wherein the road cost grid data can be dynamically generated based on online map navigation data, and the road speed provided by the online map navigation is used in the road cost grid data calculation, which is mainly realized according to the following steps:

[0057] First, set a batch of OD point pair list for online map navigation query, each OD point pair including the geographic spatial coordinates of the starting point and the ending point; OD is the starting point to the ending point in traffic;

[0058] Second, obtain the road segment and the road segment speed of each OD point pair by using the travel API in the online map navigation interface. At the same time, the obtained road segment speed is assigned to the corresponding road segment.

[0059] Finally, convert the road segment obtained by online map navigation into road cost grid data, and the value of the road cost grid is formed according to the speed assigned to each road segment.

[0060] Step 1.3) add the closed road network and the priority road to the grid cost accessibility calculation through the rule 2 restrictive conditions;

[0061] Referring to the drawings Figure 3 The priority road passing rule refers to two roads intersecting in space, and some types of roads can actually pass through, such as a road passing through a railway underpass. Therefore, when the grid accessibility is calculated, it is ensured that the roads can pass through in the case of space intersection. This means that these roads are in the uppermost layer when the grid cost base map is calculated. In this example, the closed road system is mainly highways, subways, etc., which actually form a closed road blocking area, and in the running process, the vehicle cannot leave the road system at will. But the road above the closed road blocking area should be passable. Therefore, the uppermost layer is the non-closed road system in the city, which is also relative to the closed road system.

[0062] Step 2) Server-side accessibility calculation response function design.

[0063] Step 2.1) Deploy the ArcToolbox in step 1 to the server side and publish it as geoprogressingservice, and set it to asynchronous processing mode;

[0064] Step 2.2) The server side exposes the parameters of the geoprogressing service to the client; at least including the following parameters: different types of road networks, water systems, study area ranges, and source data for grid accessibility calculation.

[0065] Step 3) Client-side accessibility calculation parameter sending request function design.

[0066] Step 3.1) According to the parameters provided in step 2.2, the user sets the parameters;

[0067] Step 3.2) After the user sets the parameters, submit the calculation task to the server;

[0068] Referring to the drawings Figure 4 The accessibility analysis function panel schematic diagram includes different types of road networks, water systems, study area ranges, and source data for grid accessibility calculation.

[0069] Step 4) Client-side interactive adjustment of source function design.

[0070] Step 4.1) Convert the source data for grid accessibility calculation on the server side to GeoJson format data;

[0071] Step 4.2) The client uses the provided function to edit GeoJson format data to edit the source data for grid accessibility calculation, and saves the result as GeoJson format data;

[0072] Step 4.3) The client uses the edited and saved GeoJson format data as the source data for the new raster reachability calculation;

[0073] Step 4.4) The client resends the request, causing the server to recalculate the raster reachability.

[0074] Step 5) After receiving the results returned by the server, the client parses, draws, and statistically analyzes the results, and has the function of calculating and converting different travel modes.

[0075] Step 5.1) The client parses, draws, and analyzes the results returned by the server;

[0076] See appendix Figure 5 After the client sends the parameters for the request calculation and the server returns the results, the client completes the parsing and rendering of the reachability calculation results.

[0077] See appendix Figure 6 The client also completed the statistics of different accessibility results, such as the coverage area percentage for 0-5 minutes and 5-10 minutes.

[0078] Step 5.2) Further, according to the specified spatial range, conduct statistics on the residential land coverage of accessibility services;

[0079] Step 5.3) According to the conversion formula for walking and cycling time, statistically analyze the accessibility coverage under different modes of transportation. The default setting of this invention is: the conversion formula for walking and cycling time is that the bicycle speed is 5 times the walking speed.

[0080] See appendix Figure 7 This invention compares the differences in calculation results between closed and non-closed road systems.

[0081] See appendix Figure 8 This is a diagram showing the results of a public service facility accessibility calculation according to the present invention. (Attached) Figure 9 This displays the reachability results after source adjustment. The differences in calculation results after source adjustment are clearly visible. Because this invention provides online source adjustment functionality, the calculation has strong dynamic interactivity.

[0082] The present application also includes two innovative designs of the two aspects: the first is to take congestion information as a basic condition for grid accessibility calculation, so as to better simulate the influence of actual road travel conditions on accessibility. The second innovation is to simulate the travel accessibility distribution of different travel modes by speed equivalent conversion of different travel modes. This means that users only need to perform accessibility calculation of one traffic mode to obtain the results of different traffic modes. This has not been mentioned in the existing research results.

[0083] The above is a further detailed description of the present application in combination with specific preferred embodiments, and cannot be regarded as limiting the specific implementation of the present application to these descriptions. For ordinary skilled persons in the technical field to which the present application belongs, a number of simple deductions or substitutions can be made without departing from the concept of the present application, and all of them should be regarded as falling within the protection scope of the present application.

Claims

1. A method for analyzing the accessibility of public service facilities through online interactive means, characterized in that, Includes the following steps: Step 1: Improve raster cost accessibility calculation by encapsulating the improved calculation into a standalone ArcToolbox tool within ArcGIS software; details are as follows: Step 1.1: Dynamically generate road cost raster data based on online map navigation data, and use the road speed provided by online map navigation as the attribute value of the road cost raster data; Two restrictive conditions, namely closed road networks and priority roads, are added to the grid cost accessibility calculation; The method for dynamically generating road cost raster data based on online map navigation data includes the following steps: Step 1.1.1: Set up a batch list of OD point pairs for online map navigation queries. Each OD point pair includes the geospatial coordinates of the origin and destination. Step 1.1.2: Use the travel API in the online map navigation interface to obtain the queried road segments and speeds of each OD point pair; at the same time, assign the obtained road segment speeds to the corresponding road segments. Step 1.1.3: Convert the road segments obtained through online map navigation into road cost raster data. The value of the road cost raster is formed according to the speed assigned to each road segment. Step 1.2: Based on Step 1.1, utilize the ArcGIS secondary development API to encapsulate the ArcToolbox tool for raster cost accessibility calculation; Step 2: Deploy the ArcToolbox tool from Step 1 to the server and design the server-side reachability calculation response function; details are as follows: Step 2.1: Deploy the ArcToolbox tool from Step 1 to the server and publish it as a geoprocessing service, and set the geoprocessing service to asynchronous processing mode; Step 2.2: The server exposes the parameters of the geoprocessing service from Step 2.1 to the client; Step 3: Design the function to send requests for client reachability calculation parameters; details are as follows: Step 3.1: Encapsulate the parameters provided in Step 2.2 into a single request task package; Step 3.2: Submit the packaged results to the server as part of the computation task; Step 4: Design a client application to implement interactive source adjustment functionality and recalculate raster cost accessibility; details are as follows: Step 4.1: Convert the source data for raster cost reachability calculation pre-stored on the server into GeoJson format data; Step 4.2: The client uses the provided function to edit GeoJson format data to edit the source data for raster cost accessibility calculation and saves the result as GeoJson format data. Step 4.3: The client uses the edited and saved GeoJson format data as the source data for the new raster cost accessibility calculation; Step 4.4: The client resends the source data for the new raster cost reachability calculation obtained in Step 4.3, and the server also recalculates the raster cost reachability. Step 5: After receiving the new raster cost reachability calculation results from Step 4 returned by the server, the client parses, draws, and statistically analyzes the results, as follows: Step 5.1: The client parses, draws, and analyzes the new raster cost reachability calculation results returned by the server in Step 4.4; Step 5.2: Simultaneously, perform spatial overlay analysis between the new raster cost accessibility calculation results returned from Step 4.4 and residential land to determine the residential land coverage statistics of accessibility services.

2. The online interactive public service facility accessibility analysis method according to claim 1, characterized in that, The online interactive public service facility accessibility analysis method also includes: realizing the calculation and conversion function of different travel modes: according to the conversion formula of walking and cycling time, the accessibility coverage of different transportation modes is statistically analyzed.

3. The online interactive public service facility accessibility analysis method according to claim 1, characterized in that, Step 1 involves developing and encapsulating a custom ArcToolbox tool for calculating raster cost accessibility based on ArcGIS versions 10.1 to 10.

8.

4. The online interactive public service facility accessibility analysis method according to claim 1, characterized in that, In step 2.2, when the server exposes the parameters of a geoprocessing service to the client, the parameters include: different types of road networks, water systems, study area extent, and source data for raster cost accessibility calculation.

5. The online interactive public service facility accessibility analysis method according to claim 2, characterized in that, The formula for converting time between walking and cycling is that the speed of a bicycle is 5 times the walking speed.

Images