KR101721114B1 - Method for Determining the Size of Grid for Clustering on Multi-Scale Web Map Services using Location-Based Point Data - Google Patents

Method for Determining the Size of Grid for Clustering on Multi-Scale Web Map Services using Location-Based Point Data Download PDF

Info

Publication number
KR101721114B1
KR101721114B1 KR1020160079891A KR20160079891A KR101721114B1 KR 101721114 B1 KR101721114 B1 KR 101721114B1 KR 1020160079891 A KR1020160079891 A KR 1020160079891A KR 20160079891 A KR20160079891 A KR 20160079891A KR 101721114 B1 KR101721114 B1 KR 101721114B1
Authority
KR
South Korea
Prior art keywords
point data
scale
grid
map
user terminal
Prior art date
Application number
KR1020160079891A
Other languages
Korean (ko)
Inventor
유기윤
이영민
양종현
Original Assignee
서울대학교산학협력단
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 서울대학교산학협력단 filed Critical 서울대학교산학협력단
Priority to KR1020160079891A priority Critical patent/KR101721114B1/en
Application granted granted Critical
Publication of KR101721114B1 publication Critical patent/KR101721114B1/en

Links

Images

Classifications

    • G06F17/30241
    • G06F17/30598
    • G06F17/30601
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/01Social networking

Abstract

The present invention relates to a method capable of determining a proper size of a grid while minimizing a modifiable areal unit problem (MAUP) when clustering a point data including positional information on a multi-scale web map. That is, the present invention considers quantitative aspects and qualitative aspects in clustering, and as a result, it is possible to minimize the influence of the MAUP. In order to consider quantitative aspects, a proper number of grids to be expressed for each scale may be calculated by using Topfers Radical Law, which is one of map generalization operators. In order to consider qualitative aspects, a size of grids, which can maintain the intrinsic distribution characteristic of the data even during the clustering process is performed for each scale, is determined by using Morans I, which is one of the global spatial relevance indexes. Thus, it is possible for users of a map service to obtain the same visual information about the distribution characteristic of the original data even when the scale of the map changes and the influence of the MAUP can be minimized.

Description

위치정보가 포함된 포인트 데이터를 다축척의 웹 지도상에 클러스터링하기 위하여 격자의 크기를 결정하는 방법{Method for Determining the Size of Grid for Clustering on Multi-Scale Web Map Services using Location-Based Point Data}[0001] The present invention relates to a method for determining a size of a grid for clustering point data including position information on a multi-scale web map,

본 발명은 위치정보가 포함된 포인트 데이터를 다축척의 웹 지도상에 축척별로 클러스터링하여 표현할 때, 공간단위 임의성의 문제(Modifiable Areal Unit Problem, MAUP)를 최소화하기 위한 방법에 대한 것으로, 더욱 상세하게는 클러스터링 시에 양적인 측면을 고려하기 위하여 지도 일반화 연산식인 Topfer's Radical Law와, 질적인 측면을 고려하여 공간적 자기상관성 측정 지표인 Moran's I를 이용하여 격자의 크기를 결정하는 방법에 관한 것이다.The present invention relates to a method for minimizing a Modifiable Areal Unit Problem (MAUP) when clustered point data including positional information is clustered on a multi-scale web map by scale. More specifically, Is a method for determining the size of a grid using Topfer's Radical Law, which is a map generalization equation, and Moran's I , which is a spatial autocorrelation measure index, in consideration of qualitative aspects in order to consider quantitative aspects in clustering.

GPS를 내장한 모바일 스마트 기기의 사용이 보편화되면서 트위터, 페이스북, 인스타그램 등과 같은 위치 기반 소셜 네트워크(Location-Based Social Network, LBSN) 서비스에 대한 사용이 활발히 이루어지고 있다.As the use of mobile smart devices with GPS is becoming commonplace, location-based social networks (LBSN) services such as Twitter, Facebook, Instagram are being actively used.

이러한 데이터를 다축척의 지도 서비스 상에서 효과적으로 시각화하기 위해서는 포인트들을 클러스터링하여 표현해야 할 필요성이 있다. LBSN 서비스를 통해 생성된 포인트 데이터는 오염원, 기상관측점 등과 같은 기존 포인트 데이터들에 비해 사람들이 많이 몰리는 도심지와 같은 특정 지역에 특히 집중되어 나타나며, 실시간으로 많은 양의 데이터가 생성된다는 특징이 있다. 따라서 해당 데이터를 지도상에 효율적으로 나타내기 위해서는 빠른 계산이 가능한 클러스터링 기법이 적용되어야 하고, 이러한 점에서 격자 기반의 클러스터링 기법이 적절하다. 이 기법은 대상 공간을 일정한 격자 구조로 구성하고, 모든 클러스터링 과정을 해당 격자 내에서 실행하는 방식으로, 객체 수에 대해 독립적이고, 미리 계산된 격자의 수에 의존하기 때문에 계산량이 적고 처리 속도가 빠르다는 장점이 있다.In order to effectively visualize such data on a multi-scale map service, there is a need to cluster and represent points. The point data generated through the LBSN service is especially concentrated in a specific area such as a downtown area where people are more crowded than conventional point data such as pollution sources and weather observation points, and a large amount of data is generated in real time. Therefore, in order to efficiently display the data on the map, a clustering technique capable of quick calculation must be applied, and a lattice-based clustering technique is appropriate in this respect. In this method, the target space is composed of a constant grid structure and all the clustering processes are executed in the corresponding grid. Since the method is independent of the number of objects and depends on the number of pre-calculated grids, the calculation amount is low and the processing speed is high .

더욱이 LBSN 서비스에서 생성된 데이터에 대해 격자 기반 클러스터링 기법을 적용하기 위해서는 지도의 축척별로 적절한 격자의 크기를 결정해야 하는데, 이에 대한 기준은 정해진 것이 없고, 데이터의 성격과 활용 목적에 따라 달라지기 때문에 통상 연구자의 주관이 개입될 수밖에 없다. 이때, 공간 단위를 설정하는 방법에 따라 분석 결과가 달라지는 MAUP가 발생한다. 새로운 공간 단위를 구축하는 문제는 지역을 구분하는 것과 밀접한 관련이 있는데, 이는 작은 스케일의 지역을 큰 스케일의 지역으로 합역하는 것을 의미한다. 그런데 여기서 합역하는 방식에 따라 동일한 지역을 다양한 방식으로 나누는 것이 가능하다. 즉, 합역 양식에 의해서 동일 지역을 다양한 스케일별로 생성할 수도 있고, 같은 스케일 내에서 서로 다른 구획을 가지는 지역을 생성할 수도 있는 것이다. 이러한 합역 과정에서 발생되는 문제를 MAUP라고 하며, 이는 면적 기반의 공간 단위를 다룰 때 필연적으로 발생하므로 격자 기반 클러스터링 기법을 사용할 때에는 반드시 MAUP의 영향을 고려해야 하는 실정이다.Furthermore, in order to apply the grid-based clustering technique to the data generated by the LBSN service, the appropriate grid size must be determined according to the scale of the map. Since there is no standard for this, The researcher 's subjectivity must be intervened. At this time, a MAUP occurs in which the analysis result is changed according to the method of setting the spatial unit. The problem of building a new unit of space is closely related to geographical separation, which means that small scale areas are merged into large scale areas. However, it is possible to divide the same area into various ways according to the way of merging. That is, it is possible to generate the same area for various scales by the joint style, or to create an area having different compartments within the same scale. The problem that arises in the process of merging is referred to as MAUP, which is inevitable when dealing with an area-based spatial unit. Therefore, the influence of MAUP must be considered when using the grid-based clustering method.

1998. 한국 GIS 학회지. 6(1). pp. 11-23.(점패턴분석을 이용한 수치지형도의 점사상 일반화). 유근배Journal of the Korean Geographical Society. 6 (1). pp. Point Mapping Generalization of Digital Topographies Using Point Pattern Analysis. Ryu

본 발명은 상기 실정을 감안하여, 위치정보를 포함하는 포인트 데이터를 격자 기반으로 클러스터링하여 표현하고자 할 때, 공간단위 임의성의 문제(Modifiable Areal Unit Problem, MAUP)를 최소화하면서 적절한 격자의 크기를 결정할 수 있는 방법을 제시하기 위한 것이 목적이다.SUMMARY OF THE INVENTION In view of the foregoing, it is an object of the present invention to provide a method and apparatus for determining a proper grid size while minimizing a MAID (Modified Arial Problem) The purpose of this is to present a method.

본 발명은 상기 목적을 달성하기 위하여, (a) 사용자단말기의 포인트데이터수집모듈이 위치 기반 소셜 네트워크(Location-Based Social Network, LBSN) 서비스에서 제공하는 오픈 API(Open Application Programming Interface)를 이용하여 수집을 원하는 지역의 위치정보를 포함하는 포인트 데이터를 수집하는 단계; (b) 상기 사용자단말기의 최근린분석모듈이 최근린분석(Nearest Neighbor Analysis, NNA)의 최근린지수(Nearest Neighbor Index, R)를 이용하여 수집된 포인트 데이터들의 본래 공간적 분포 특성을 파악하고 분석하는 단계; (c) 상기 사용자단말기의 지도일반화연산모듈이 Topfer's Radical Law를 이용하여 계수값(C)에 따른 복수의 클러스터(Cluster, 격자)의 개수를 축척별로 산출하는 단계; (d) 상기 사용자단말기의 제어모듈은 상기 (a) 단계에서 수집을 원하는 지역을 대상으로 수집한 위치정보가 포함된 포인트 데이터와, 해당 범위에 대하여 다양한 크기별로 격자데이터를 생성하여 그 생성된 격자데이터에 대하여 공간결합(Spatial Join) 연산을 수행하는 단계; (e) 상기 사용자단말기의 제어모듈은 상기 (d) 단계에서 어떠한 포인트 데이터와 결합되지 않은 격자를 제외한 전체 격자의 개수를 산출하는 단계; (f) 상기 사용자단말기의 제어모듈은 상기 (e) 단계에서 산출한 격자의 개수와, 상기 (c) 단계에서 도출한 축척별, C값별 복수의 격자 개수 중에서 동일한 값이 있는지를 판단하되, 동일한 값이 있는 경우에는 축척별, C값별 복수의 모든 격자 개수와 각각 동일한 상기 (e) 단계에서의 격자 개수를 구할 때까지 연산을 반복 수행하고, 동일한 값이 없을 경우에는 상기 (d) 단계로 돌아가 다른 크기의 격자데이터와 포인트 데이터에 대해 공간결합 연산과 상기 (e) 단계를 반복적으로 수행하는 단계; (g) 상기 사용자단말기의 격자크기분포특성분석모듈은 Moran's I를 이용하여 상기 (f) 단계에서 수행된 결과에 대하여 축척별, C값별로 분포 특성을 분석하는 단계; (h) 상기 사용자단말기의 제어모듈은 상기 (b) 단계에서 도출된 포인트 데이터의 본래 분포 특성과 축척별 계수값(C)에 따른 상기 (g) 단계의 분포 특성 분석 결과를 비교하는 단계, 및 (i) 상기 사용자단말기의 제어모듈은 상기 (h) 단계에서 분석한 포인트 데이터의 본래 분포 특성을 가장 잘 유지한 계수값(C)을 최종적으로 도출하는 단계를 포함하여 이루어진, 위치정보가 포함된 포인트 데이터를 다축척의 웹 지도상에 클러스터링하기 위하여 격자의 크기를 결정하는 방법을 제공한 것이 특징이다.(A) the point data collection module of the user terminal collects data using an open application programming interface (API) provided by a location-based social network (LBSN) service; Collecting point data including location information of a desired area; (b) the current lean analysis module of the user terminal identifies and analyzes the original spatial distribution characteristics of the point data collected using the Nearest Neighbor Index (RN) of Nearest Neighbor Analysis (NNA) step; (c) calculating a number of a plurality of clusters according to a scale value (C) according to a scale value using Topfer's Radical Law, by the map generalization operation module of the user terminal; (d) The control module of the user terminal generates point data including positional information collected for a region desired to be collected in the step (a), grid data for various sizes for the corresponding range, Performing a Spatial Join operation on the data; (e) the control module of the user terminal calculates the total number of grids excluding the grids not combined with any point data in the step (d); (f) The control module of the user terminal determines whether the number of grids calculated in the step (e) is the same among a plurality of grid numbers for each scale and C value derived in the step (c) If there is a value, the operation is repeated until the number of grids in the step (e) that is equal to the number of all the plurality of grids for each scale and C value is obtained. If there is no value, the process returns to the step (d) Repeatedly performing the spatial combining operation and the step (e) on the grid data and the point data of different sizes; (g) analyzing a distribution characteristic of the grid size distribution characteristic of the user terminal by using the Moran's I , analyzing distribution characteristics of the result performed in the step (f) for each scale and C value; (h) the control module of the user terminal compares the distribution characteristic analysis result of the step (g) according to the original distribution characteristic of the point data derived in the step (b) and the coefficient value (C) (i) the control module of the user terminal finally deriving a coefficient value (C) that best preserves the original distribution characteristic of the point data analyzed in the step (h) It is characterized by providing a method of determining the size of a grid in order to cluster point data on a multi-scale web map.

또한, 본 발명에서, (j) 상기 사용자단말기의 제어모듈은 상기 (d) 단계에서 (i) 단계까지 클러스터링에 필요한 모든 축척에 대하여 계수값(C)별로 반복하는 과정을 거쳐 최종적인 결과는 축척별로 도출된 C값에 따른 격자의 한 변의 길이(격자의 크기)를 얻을 수 있도록 연산하는 단계를 더 포함할 수 있다.In the present invention, (j) the control module of the user terminal repeats the steps (d) to (i) for all the scales required for clustering for each count value C, And calculating the length (the size of the lattice) of one side of the lattice according to the derived C value.

또한, 본 발명에서, 상기 (b) 단계에서, 최근린지수(R)는 다음의 수학식,Further, in the present invention, in the step (b), the latest Lindsey number (R)

Figure 112016061609799-pat00001
Figure 112016061609799-pat00001

(여기서,

Figure 112016061609799-pat00002
는 각각의 포인트 데이터로부터 가장 가까운 포인트 데이터까지의 거리의 평균이고,
Figure 112016061609799-pat00003
i포인트 데이터로부터 가장 가까운 포인트 데이터까지의 거리이며, n은 전체 포인트의 개수이고,
Figure 112016061609799-pat00004
는 CSR 하에서 기대되는 최근린 거리이며, A는 수집을 원하는 지역의 전체 면적이다.)을 이용하여 도출할 수 있다.(here,
Figure 112016061609799-pat00002
Is the average of the distance from each point data to the closest point data,
Figure 112016061609799-pat00003
Is the distance from the i point data to the closest point data, n is the number of all points,
Figure 112016061609799-pat00004
Is the recent lean distance expected under CSR, and A is the total area of the area you want to collect).

또한, 본 발명에서, 상기 (c) 단계에서, Topfer's Radical Law는 다음의 수학식,Also, in the present invention, in the step (c), Topfer's Radical Law may be expressed by the following equation,

Figure 112016061609799-pat00005
Figure 112016061609799-pat00005

(여기서,

Figure 112016061609799-pat00006
는 결과 지도의 포인트 수이고,
Figure 112016061609799-pat00007
는 입력 지도의 포인트 수이며,
Figure 112016061609799-pat00008
는 입력 지도의 축척계수이고,
Figure 112016061609799-pat00009
는 결과 지도의 축척계수이며, 그리고 C 는 심볼계수이다.)을 이용하여 도출할 수 있다.(here,
Figure 112016061609799-pat00006
Is the number of points in the result map,
Figure 112016061609799-pat00007
Is the number of points in the input map,
Figure 112016061609799-pat00008
Is the scale factor of the input map,
Figure 112016061609799-pat00009
Is the scale factor of the result map, and C is the symbol coefficient).

또한, 본 발명에서, 상기 (g) 단계에서, Moran's I는 다음의 수학식,Further, in the present invention, in the step (g), Moran's I is calculated by the following equation

Figure 112016061609799-pat00010
Figure 112016061609799-pat00010

(여기서, nij로 이루어진 모든 공간 단위의 개수이고,

Figure 112016061609799-pat00011
는 공간가중행렬이며,
Figure 112016061609799-pat00012
i번째 변수(i번째 격자의 Joint_Count 값)이고,
Figure 112016061609799-pat00013
x의 평균이다.)을 이용하여 도출할 수 있다.(Where n is the number of all spatial units made up of i and j ,
Figure 112016061609799-pat00011
Is a space weighting matrix,
Figure 112016061609799-pat00012
Is the i- th variable (the Joint_Count value of the i-th grid)
Figure 112016061609799-pat00013
Is the average of x ).

본 발명에 따르면, 위치정보가 포함된 포인트 데이터를 다축척의 지도상에 격자 기반으로 클러스터링하여 표현하고자 할 때, 적절한 격자의 크기를 결정하는 방법을 제시함으로써, 클러스터링 시에는 양적인 측면과 질적인 측면을 고려하며, 이를 통해 결과적으로 공간단위 임의성의 문제(MAUP)의 영향을 최소화할 수 있다. 즉, 양적인 측면을 고려하기 위해서는 지도 일반화 연산자 중 하나인 Topfer's Radical Law를 이용하여 각 축척별로 표출되어야 할 적절한 격자의 개수를 계산하고, 또한, 질적인 측면을 고려하기 위해서는 전역적 공간 연관성 지표 중 하나인 Moran's I를 이용하여 축척별로 클러스터링 과정이 수행되는 동안에도 본래 데이터가 가지는 분포 특성을 유지할 수 있는 격자의 크기를 도출함으로써, 지도 서비스를 이용하는 사용자는 지도의 축척이 변화할 때에도 본래 데이터의 분포 특성에 대한 동일한 시각적 정보를 얻을 수 있을 뿐 아니라 공간단위 임의성의 문제(MAUP)의 영향도 최소화할 수 있는 이점이 있다.According to the present invention, when a point data including position information is cluster-based on a multi-scale map, a method of determining an appropriate size of a grid is proposed. In the clustering, quantitative and qualitative aspects , Which in turn can minimize the effect of space unit randomness problem (MAUP). That is, in order to consider the quantitative aspect, it is necessary to calculate the number of appropriate grids to be displayed for each scale by using Topfer's Radical Law, one of the map generalization operators, and to calculate the qualitative aspect, The user of the map service can obtain the distribution characteristic of the original data even when the scale of the map is changed by deriving the size of the grid capable of maintaining the distribution characteristic of the original data even while the clustering process is performed by using the Moran's I , It is possible to obtain not only the same visual information but also the effect of the problem of spatial randomness (MAUP).

도 1은 본 발명에서 위치 기반 소셜 네트워크 서비스를 통해 생성되는 포인트 데이터를 예시적으로 나타낸 것이다.
도 2는 본 발명에서 포인트 데이터의 클러스터링과 클러스터를 예시적으로 나타낸 것이다.
도 3은 본 발명에서 지도상에 마커들이 클러스터링되지 않은 축척별의 모습을 예시적으로 나타낸 것이다.
도 4는 본 발명에서 포인트 데이터를 다축척의 웹 지도상에서 효과적으로 시각화하기 위한 클러스터링 전과 후를 비교한 모습을 나타낸 것이다.
도 5는 본 발명에 따른 줌 레벨별 스케일을 예시적으로 나타낸 도표이다.
도 6은 본 발명에서 공간단위 임의성의 문제(MAUP)에서의 스케일 효과와 구획효과를 예시적으로 나타낸 것이다.
도 7은 본 발명에 따른 실시 예로, 위치정보가 포함된 포인트 데이터를 다축척의 웹 지도상에 클러스터링하기 위하여 격자의 크기를 결정하는 방법의 구현을 위한 시스템을 나타낸 구성도이다.
도 8은 본 발명에 따른 위치정보가 포함된 포인트 데이터를 다축척의 웹 지도상에 클러스터링하기 위하여 격자의 크기를 결정하는 방법을 나타낸 흐름도이다.
도 9는 본 발명에서 C값에 따른 줌 레벨별 복수의 포인트 개수의 계산을 예시적으로 나타낸 것이다.
도 10은 본 발명에서 공간결합 연산을 예시적으로 나타낸 것이다.
도 11은 본 발명에서 공간결합 연산을 통해 클러스터 개수와 클러스터로 표현될 격자의 개수의 일치여부를 예시적으로 나타낸 것이다.
도 12는 본 발명에서 공간분포를 예시적으로 나타낸 것이다.
도 13은 본 발명에서 공간분포에 대한 공간가중행렬을 나타낸 도표이다.
1 illustrates an exemplary point data generated through a location-based social network service in the present invention.
Figure 2 illustrates clustering and clustering of point data in the present invention.
3 illustrates an example of scale-by-scale appearance in which markers are not clustered on a map in the present invention.
FIG. 4 shows a comparison between before and after clustering for effectively visualizing point data on a multi-scale web map in the present invention.
FIG. 5 is a diagram exemplarily showing a scale for each zoom level according to the present invention. FIG.
Figure 6 is an illustration of the scale effect and the compartment effect in the problem of space unit randomness (MAUP) in the present invention.
7 is a diagram illustrating a system for implementing a method of determining a size of a grid to cluster point data including positional information on a multi-scale web map according to an embodiment of the present invention.
8 is a flowchart illustrating a method of determining a size of a grid to cluster point data including positional information according to the present invention on a multi-scale web map.
9 illustrates an exemplary calculation of a plurality of points per zoom level according to the C value in the present invention.
10 illustrates an exemplary spatial coupling operation in the present invention.
11 illustrates an example of the number of clusters and the number of clusters to be represented by clusters according to the spatial combining operation in the present invention.
12 shows an exemplary spatial distribution in the present invention.
13 is a diagram showing a spatial weighting matrix for spatial distribution in the present invention.

이하, 본 발명에 따른 위치정보가 포함된 포인트 데이터를 다축척의 웹 지도상에 클러스터링하기 위하여 격자의 크기를 결정하는 방법에 관하여 첨부된 도면을 참조하여 상세하게 설명한다.Hereinafter, a method for determining a size of a grid for clustering point data including positional information according to the present invention on a multi-scale web map will be described in detail with reference to the accompanying drawings.

먼저, 본 발명에서 사용되는 특정의 용어에 대하여 정의한다.First, specific terms used in the present invention will be defined.

도 1에서, 트위터, 인스타그램, 페이북 등의 위치 기반 소셜 네트워크(Location-Based Social Network, LBSN) 서비스를 통해 생성되는 데이터를 포인트 데이터(Point Data)라 한다. 더욱이 LBSN 서비스 사용자들은 지오태깅(Geotagging) 기능을 이용하여 자신이 작성하는 글과 함께 위치정보(x, y 좌표)를 포함시킬 수 있으므로 해당 서비스를 통해 생성된 데이터는 포인트의 속성을 가지는 공간 객체로 볼 수 있다. 이때, 위치정보는 사용자의 현재 위치이거나 또는 업로드한 사진을 찍은 위치 또는 사용자가 직접 지정한 위치가 될 수 있다. 또한, 포인트 데이터가 주로 활용되는 곳으로, 포인트 데이터는 다양한 분야에서 활용될 수 있는데, 대표적인 예로 오염원의 위치, 기상관측점, 인구 등이 있다. 더욱이, 2000년대 초반에 스마트폰 등의 스마트 기기가 개발되면서 트위터, 인스타그램, 페이스북 등의 LBSN 서비스를 통해 생성된 포인트 데이터가 새롭게 등장했으며, 이와 같은 LBSN 포인트 데이터는 기존의 포인트 데이터에 비해 사람들이 많이 몰리는 도심지와 같은 곳에 집중되어 나타나며, 실시간으로 많은 양의 데이터가 생성된다는 특징이 있다.In FIG. 1, data generated through a Location-Based Social Network (LBSN) service such as Twitter, Instagram, and PayBook is referred to as Point Data. Furthermore, LBSN service users can include location information (x, y coordinates) along with their own articles by using the geotagging function. Therefore, the data generated through the service is a spatial object having attributes of points can see. At this time, the location information may be the current location of the user, the location where the uploaded photograph was taken, or the location designated by the user. In addition, point data is mainly used, and point data can be used in various fields. Typical examples are location of pollution sources, weather observation points, and population. Furthermore, as smart devices such as smart phones were developed in the early 2000s, point data generated through LBSN services such as Twitter, Instagram, and Facebook newly appeared, and such LBSN point data is compared with existing point data It is concentrated in the same places as the downtown area where people are attracted to many people, and a large amount of data is generated in real time.

또한, 도 2에서, 클러스터링(Clustering, 군집화)의 본래 의미는 물리적 혹은 추상적 객체들을 비슷한 객체로 그룹화하는 것이며, 이때 같은 클러스터(Cluster, 군집) 내의 객체들과는 유사하고, 다른 클러스터의 객체들과는 상이한 특성을 지닌다. 본 발명에서는 여러 클러스터링 기법 중에서도 격자 기반 클러스터링(Grid-Based Clustering) 기법을 사용하는데, 이는 같은 격자에 포함되는 포인트들을 하나의 클러스터로 묶어서 표현하는 방법이다. 더욱이 클러스터링과 클러스터의 차이로서, 클러스터링은 군집화하는 과정 자체를 의미하며, 클러스터는 클러스터링의 결과로 도출된 결과물을 의미한다. 따라서 본 발명에서는 격자 기반 클러스터링 기법을 기반으로 하되, 격자 자체를 클러스터로 취급하므로, 클러스터는 곧 단위격자를 의미한다.2, the original meaning of clustering is to group physical or abstract objects into similar objects, which are similar to objects in the same cluster (cluster), and have different characteristics from those of other clusters I have. In the present invention, among various clustering techniques, a grid-based clustering technique is used, which expresses points included in the same grid by grouping them into one cluster. Moreover, as a difference between clustering and clustering, clustering means the process of clustering, and clustering means the result of clustering. Therefore, the present invention is based on a lattice-based clustering technique, and since the lattice itself is treated as a cluster, the cluster means a unit lattice.

또한, 도 3에서, 지도 분야에서의 클러스터링은 주로 시각적인 목적으로 사용된다. 일례로, 다음(DAUM) 아고라와 네티즌들이 구글 맵(Google Maps)의 위치 정보를 기반으로 하여 SNS와 인터넷 게시판을 통해 자발적으로 제작한 '폭우피해지도'는 2011년 기록적인 폭우가 발생했을 당시 정보 공유 플랫폼으로서 유용했지만, 지도상에 표시된 마커들이 클러스터링 되지 않은 채 몰려 있거나 겹쳐 있어 소축척 하에서 정확한 정보를 직관적으로 파악하기에는 어려움이 있었다.Also, in Fig. 3, clustering in the map field is mainly used for visual purposes. For example, the following (DAUM) agora and netizens voluntarily produced 'Sufi Damage Map' based on SNS and Internet bulletin board based on the location information of Google Maps, Although useful as a shared platform, markers on the map were packed or unclustered without clustering, making it difficult to intuitively grasp accurate information under small scale.

따라서 포인트 데이터를 다축척의 웹 지도상에서 효과적으로 시각화하기 위해서는 클러스터링 과정이 필수적이라고 할 수 있다. 이는 도 4에서, 세계지도를 중심으로 클러스터링 전과 후의 모습을 비교한 사진을 참조하면 알 수 있다.Therefore, clustering is essential to effectively visualizing point data on a multi-scale web map. This can be seen from FIG. 4 with reference to a photograph of the world map before and after clustering.

다음으로, 도 5에서, 지도의 줌 레벨에 따른(축척별) 스케일을 표시하였는데, 이 모든 줌 레벨에 대하여 클러스터링 작업을 실시할 필요는 없다. 즉 예를 들어, 서울을 대상으로 클러스터링 작업을 수행한다고 한다면, 서울의 영역을 기준으로 봤을 때 줌 레벨 0부터 4까지는 너무 소축척이기 때문에 클러스터링은 무의미하다고 할 수 있고, 줌 레벨 17부터 19까지는 상대적으로 대축척에 해당하기 때문에 클러스터링으로 표현하는 것이 아닌 본래 데이터의 형태, 즉 개별 포인트의 형태로 표현하는 것이 더욱 적합하다고 할 수도 있다. 이렇게 되면 최종적으로는 줌 레벨 5부터 16까지가 클러스터링 대상이 된다. 따라서 클러스터링을 실시할 줌 레벨의 범위는 실험의 대상과 목적에 따라 달라지며, 작업자의 판단에 따를 수도 있다.Next, in FIG. 5, a scale (by scale) according to the zoom level of the map is displayed, but it is not necessary to perform the clustering operation for all the zoom levels. For example, if clustering is performed for Seoul, clustering is meaningless because zoom levels 0 to 4 are too small when viewed from the Seoul area, and zoom levels 17 to 19 are relative It may be said that it is more appropriate to express it in the form of data originally, that is, in the form of individual points, rather than being represented by clustering. Finally, zoom levels 5 through 16 are clustering targets. Therefore, the range of the zoom level to be clustering depends on the object and purpose of the experiment, and may be determined by the operator.

또한, 트위터, 인스타그램, 페이스북 등 위치정보를 포함하는 LBSN 데이터들은 해당 데이터가 참조하는 장소에 대한 좌표 정보를 포함하고 있기 때문에 이를 공간상에 표현할 때 포인트의 형태를 취하게 된다. 어떠한 공간적 범위, 예를 들어, 서울시를 대상으로 수집된 포인트 데이터를 특정 공간 단위, 즉, 행정동 또는 격자를 기준으로 분석하기 위해서는 포인트 속성의 원데이터를 폴리곤 속성의 데이터에 결합(공간결합)하여 사용해야 한다. 이때, 연구자는 분석에 사용할 공간 단위를 선택하거나 기존의 공간 단위를 조작하여 새로운 공간 단위를 만들어내야 하는데, 이에 대한 명확한 기준은 정해진 것이 없기 때문에 연구자의 판단에 따를 수밖에 없다. In addition, since LBSN data including position information such as tweeter, Instagram, and Facebook includes coordinate information on a place referred to by the corresponding data, it takes the form of a point when it is expressed in space. In order to analyze point data collected in a certain spatial range, for example, Seoul, based on a specific spatial unit, that is, an administrative district or a grid, the original data of the point attribute should be combined with the data of the polygon attribute do. At this time, the researcher must select a space unit to be used for analysis or create a new space unit by manipulating an existing space unit. Since there is no definite criterion for this, the researcher must judge it.

새로운 공간 단위를 구축하는 문제는 지역구분(Regionalization)과 밀접한 관련성을 가지는데, 이는 실질적으로 작은 스케일의 지역을 큰 스케일의 지역으로 합역(Spatial Aggregation)하는 것을 의미한다. 그런데 합역양식(Aggregation Scheme)에 따라 동일한 지역을 다양한 방식으로 나눌 수 있다. 즉, 합역양식에 의해 다양한 스케일의 지역을 생성할 수도 있고, 같은 스케일에서 서로 다른 구획을 가지는 지역을 생성할 수도 있는 것이다. 이러한 합역 과정에서 발생되는 문제점을 공간단위 임의성의 문제(Modifiable Areal Unit Problem, MAUP)라고 한다.The problem of establishing a new spatial unit is closely related to the regionalization, which means that the actual small scale area is spatially aggregated into a large scale area. However, the same area can be divided in various ways according to the aggregation scheme. In other words, it is possible to create regions of various scales according to the joint style, or to create regions having different divisions on the same scale. The problem that arises in the process of this integration is called the Modifiable Areal Unit Problem (MAUP).

이와 같이 MAUP은 분석에 사용되는 공간 단위를 어떻게 설정하느냐에 따라 분석 결과에 영향을 끼친다. 공간 단위는 상이한 스케일이나 구획을 만들기 위해서 합역되거나 형태가 변형될 수 있다. 즉, MAUP는 스케일 효과(Scale Effect)와 구획 효과(Zoning Effect)의 두 가지 관점을 갖는다.Thus, the MAUP affects the analysis results depending on how the spatial unit used for the analysis is set. Spatial units can be merged or modified to create different scales or segments. That is, MAUP has two viewpoints: a scale effect and a zoning effect.

더욱이 도 6에서, 스케일 효과는 서로 다른 스케일의 공간 단위를 대상으로 분석을 실시했을 때 스케일에 따라 결과가 달라지는 것을 의미하며, 구획 효과는 동일한 스케일에서 공간 단위의 재그룹화를 통해 발생되는 서로 다른 결과를 말한다.Further, in FIG. 6, the scale effect means that the results are different according to the scale when the analysis is performed on the spatial units of different scales, and the division effect is different results .

도 7은 본 발명에서 위치정보가 포함된 포인트 데이터를 다축척의 웹 지도상에 클러스터링하기 위하여 격자의 크기를 결정하는 방법을 구현하기 위한 시스템이다. 사용자단말기(10)는 테스크톱이나 노트북 등의 PC나 태블릿이나 스마트폰 등이 포함된다. 사용자단말기(10)에는 클러스터링 애플리케이션이 설치된다. 클러스터링 애플리케이션에는 제어모듈(11), 포인트데이터수집모듈(12), 최근린분석모듈(13), 지도일반화연산모듈(14) 및 격자크기분포특성분석모듈(15)을 포함한다. 그리고 사용자단말기(10)는 트위터, 인스타그램, 페이스북 등의 위치 기반 소셜 네트워크(LBSN)(20) 서비스에 통신망을 통해 접속된다. 클러스터링 애플리케이션은 상용프로그램인 ArcMap 또는 프로그램 언어를 통해 직접 알고리즘을 코딩한 프로그램에 의하여 이루어진다.7 is a system for implementing a method of determining a size of a grid in order to cluster point data including positional information on a multi-scale web map in the present invention. The user terminal 10 includes a PC, a tablet, a smart phone, and the like, such as a desk top or a notebook computer. The user terminal 10 is installed with a clustering application. The clustering application includes a control module 11, a point data collection module 12, a current lean analysis module 13, a map generalization calculation module 14 and a grid size distribution characteristic analysis module 15. The user terminal 10 is connected to a location based social network (LBSN) 20 service such as Twitter, Instagram, Facebook, etc. through a communication network. The clustering application is accomplished by a commercial program, ArcMap, or a program that directly codes the algorithm through a programming language.

다음으로, 본 발명에 따른 위치정보가 포함된 포인트 데이터를 다축척의 웹 지도상에 클러스터링하기 위하여 격자의 크기를 결정하는 방법에 관하여 도 8의 흐름도를 참조하여 설명한다.Next, a method for determining the size of a grid to cluster point data including positional information according to the present invention on a multi-scale web map will be described with reference to the flowchart of FIG.

우선, 사용자단말기(10)의 포인트데이터수집모듈(12)이 위치 기반 소셜 네트워크(LBSN)(20)에서 제공하는 오픈 API(Open Application Programming Interface)를 이용하여 수집을 원하는 지역의 위치정보를 포함하는 포인트 데이터를 수집한다(S1). 그리고 사용자단말기(10)의 최근린분석모듈(13)이 최근린분석(Nearest Neighbor Analysis, NNA)의 최근린지수(Nearest Neighbor Index, R)를 이용하여 수집된 포인트 데이터들의 본래 공간적 분포 특성을 파악하고 분석한다(S2). 이 분석된 결과는 다음의 (S8) 단계에서 이용된다. 최근린분석은 포인트 데이터의 공간적 분포 특성을 파악하기 위한 분석법으로, 각각의 포인트 데이터로부터 가장 가까운 포인트 데이터까지의 거리의 평균값(

Figure 112016061609799-pat00014
)과 완전한 랜덤 분포(Complete Spatial Randomness, CSR) 하에서 기대되는 최근린 거리(
Figure 112016061609799-pat00015
)를 비교한다. 이때, 기대되는 최근린 거리(
Figure 112016061609799-pat00016
)에 대한 관측된 평균 최근린 거리(
Figure 112016061609799-pat00017
)의 비율을 R이라고 한다. 이러한 최근린지수(R)는 다음의 수학식 1과 같다.First, the point data collection module 12 of the user terminal 10 uses the open application programming interface (API) provided by the location-based social network (LBSN) The point data is collected (S1). Then, the recent lean analysis module 13 of the user terminal 10 grasps the inherent spatial distribution characteristics of the collected point data using the Nearest Neighbor Index (RN) of the nearest neighbor analysis (NNA) (S2). The analyzed result is used in the next step (S8). The recent lean analysis is an analysis method for grasping the spatial distribution characteristic of point data. The average value of the distance from each point data to the closest point data (
Figure 112016061609799-pat00014
) And the expected lean distance under Complete Spatial Randomness (CSR)
Figure 112016061609799-pat00015
). At this time, the expected lean distance (
Figure 112016061609799-pat00016
The observed mean recent lean distance (
Figure 112016061609799-pat00017
) Is referred to as R. This latest Lindis number (R) is given by the following equation (1).

Figure 112016061609799-pat00018
Figure 112016061609799-pat00018

여기서,

Figure 112016061609799-pat00019
는 각각의 포인트 데이터로부터 가장 가까운 포인트 데이터까지의 거리의 평균이고,
Figure 112016061609799-pat00020
i포인트 데이터로부터 가장 가까운 포인트 데이터까지의 거리이며, n은 전체 포인트의 개수이고,
Figure 112016061609799-pat00021
는 CSR 하에서 기대되는 최근린 거리이며, A는 수집을 원하는 지역의 전체 면적이다.here,
Figure 112016061609799-pat00019
Is the average of the distance from each point data to the closest point data,
Figure 112016061609799-pat00020
Is the distance from the i point data to the closest point data, n is the number of all points,
Figure 112016061609799-pat00021
Is the nearest lean distance expected under CSR, and A is the total area of the area you want to collect.

R은 각 포인트들의 실제 간격이 랜덤한 분포 패턴에서의 간격에 비해 어느 정도의 차이가 있는지를 나타내는 값으로, R이 1인 경우는 완전한 랜덤 분포 유형, R이 1보다 큰 경우는 분산된 분포 유형, 그리고 R이 1보다 작은 경우는 클러스터된 분포 유형을 보인다고 해석한다.R is a value indicating how much the actual spacing of points is different from the spacing in the random distribution pattern. If R is 1, it is a complete random distribution type. If R is greater than 1, , And if R is less than 1, the clustered distribution type is interpreted.

다음으로, 사용자단말기(10)의 지도일반화연산모듈(14)이 Topfer's Radical Law를 이용하여 계수값(C)에 따른 복수의 클러스터(Cluster, 격자)의 개수를 축척별로 산출한다(S3). 즉 지도일반화연산모듈(14)은 지도 일반화 연산식인 Topfer's radical law의 심볼계수(C) 값을 0.1 ~ 1.9까지 0.1씩 19단계로 증가시키며, 각각의 경우에 대해 축척별로 도출되어야 할 적절 포인트의 개수(

Figure 112016061609799-pat00022
)를 계산한다. Topfer's Radical Law는 다음의 수학식 2와 같다.Next, the map generalization operation module 14 of the user terminal 10 calculates the number of clusters according to the count value C by scale using Topfer's Radical Law (S3). That is, the map generalization operation module 14 increases the symbol coefficient ( C ) value of Topfer's radical law, which is a map generalization operation expression, from 0.1 to 1.9 in increments of 0.1 in 19 steps. In each case, the number of appropriate points (
Figure 112016061609799-pat00022
). Topfer's Radical Law is expressed by Equation 2 below.

Figure 112016061609799-pat00023
Figure 112016061609799-pat00023

여기서,

Figure 112016061609799-pat00024
는 결과 지도의 포인트 수이고,
Figure 112016061609799-pat00025
는 입력 지도의 포인트 수이며,
Figure 112016061609799-pat00026
는 입력 지도의 축척계수이고,
Figure 112016061609799-pat00027
는 결과 지도의 축척계수이며, 그리고 C 는 심볼계수이다. 바람직하게는 입력 지도의 포인트 수는 (S1) 단계에서의 수집된 포인트 데이터 수이고, 입력 지도의 축척은 사용자가 선택한 줌 레벨의 범위 중 가장 대축척이며, 그리고 입력 지도는 클러스터링을 하지 않기 때문에 결과 지도의 축척계수는 최대 19개가 될 수 있다. 즉 지도의 축척계수는 구글 맵(Google Maps)에서 줌 레벨별로 제공하는 축척값을 이용할 수 있으며, 이는 도 5에 도시되었다. 이를 이용하여 클러스터링이 필요한 축척에 대한 범위를 설정한다. 더욱이 C값의 범위와 지도의 축척은 사용자의 선택에 따라 조절이 가능하다. 다시 말해,
Figure 112016061609799-pat00028
Figure 112016061609799-pat00029
값은 도 5의 스케일을 이용한다.here,
Figure 112016061609799-pat00024
Is the number of points in the result map,
Figure 112016061609799-pat00025
Is the number of points in the input map,
Figure 112016061609799-pat00026
Is the scale factor of the input map,
Figure 112016061609799-pat00027
Is the scale factor of the result map, and C is the symbol coefficient. Preferably, the number of points in the input map is the number of collected point data in step (S1), the scale of the input map is the largest of the range of the zoom level selected by the user, and since the input map does not clustering, Can be a maximum of 19 scale factors. That is, the scale factor of the map can use a scale value provided for each zoom level in Google Maps, which is shown in FIG. Using this, we set the scope for the scale that needs clustering. Furthermore, the range of C values and the scale of the map can be adjusted according to the user's choice. In other words,
Figure 112016061609799-pat00028
Wow
Figure 112016061609799-pat00029
The values use the scale of FIG.

도 9는 C값에 따른 줌 레벨별 복수의 포인트 개수의 계산을 예시적으로 나타낸 것이다. 예를 들어, 줌 레벨 19에서 보인 포인트의 개수가 1,000개일 때, 줌 레벨 15에서 보여야 할 적절한 포인트의 수는 상기 수학식 2를 이용하여 다음과 같이 계산할 수 있다(C = 1일 때).Figure 9 is an exemplary illustration of the calculation of a plurality of points per zoom level according to the value of C; For example, when the number of points seen at zoom level 19 is 1,000, the number of appropriate points to be seen at zoom level 15 can be calculated using Equation (2) as follows (when C = 1).

Figure 112016061609799-pat00030
Figure 112016061609799-pat00030

즉, 줌 레벨 19에서 보인 포인트의 개수가 1,000개일 때, 줌 레벨 15에서 보여야 할 적절한 포인트의 수 250개가 된다. 이때, 도출된 포인트의 개수(

Figure 112016061609799-pat00031
)를 축척별(줌 레벨별)로 도출되어야 할 적절한 클러스터(격자)의 개수(n)로 간주한다.That is, when the number of points seen at zoom level 19 is 1,000, the appropriate number of points to show at zoom level 15 is 250. At this time, the number of derived points (
Figure 112016061609799-pat00031
) As the number ( n ) of appropriate clusters (grids) to be derived by scale (per zoom level).

본 발명에서는 포인트 데이터가 줌 레벨별로 클러스터링되어 표현될 때도 본래의 분포 특성을 유지하는 것을 중요시하기 때문에 줌 레벨별로 계산된 적절 포인트의 개수를 클러스터(격자)의 개수(n)로 간주한 것이다. 결국, 본 발명에서 구하고자 하는 것은, 상기 (S3) 단계에서 도출된 적절 클러스터(격자)의 개수(n)를 도출할 수 있는 클러스터(격자)의 크기, 즉, 한 변의 길이, l를 결정하는 것이다.In the present invention, it is important to maintain the original distribution characteristics even when the point data is clustered according to the zoom level. Therefore, the number of appropriate points calculated for each zoom level is regarded as the number ( n ) of clusters (grids). In other words, what is sought in the present invention is to determine the size of a cluster (grid), i.e., the length l of one side, from which the number n of appropriate clusters derived from the step S3 is derived will be.

따라서 개수를 알면 역으로 그 크기도 알 수 있지 않느냐고 할 수 있지만, 이는 단순한 치환의 문제가 아니다. 왜냐하면, 지도상에 모든 클러스터(격자)가 표출되는 것이 아니기 때문이다. 개중에는 클러스터로써 표현될 필요가 없는 격자가 있을 수 있기 때문이다. 이에 대해서는 (S4) 단계와 (S5) 단계에서 각각 설명한다.Therefore, if we know the number, we can say that we can not know its size, but this is not a matter of simple substitution. This is because not all clusters (grids) are displayed on the map. In some cases, there may be a grid that does not need to be represented as a cluster. This will be described in step S4 and step S5, respectively.

사용자단말기(10)의 제어모듈(11)은 상기 (S1) 단계에서 수집을 원하는 지역을 대상으로 수집한 위치정보가 포함된 포인트 데이터와, 해당 범위에 대하여 다양한 크기별로 격자데이터를 생성하여 그 생성된 격자데이터에 대하여 공간결합(Spatial Join) 연산을 수행한다(S4). 이때, 다양한 크기라는 것은 축척별로 생성한다는 의미가 아니라, 축척을 떠나서 랜덤으로 최대한 다양한 크기의 격자를 생성한다는 것을 의미한다. 이러한 격자데이터는 (S3) 단계에서의 포인트에 대한 격자와는 다른 것이다.The control module 11 of the user terminal 10 generates point data including positional information collected for a region desired to be collected in the step S1 and lattice data for various sizes for the corresponding range, Spatial Join operation is performed on the lattice data (S4). At this time, various sizes do not mean that they are generated according to the scale, but mean that they generate grids of various sizes as random as possible out of scale. This grid data is different from the grid for the point in step (S3).

도 10에서, 포인트 데이터와 격자데이터를 겹친 포인트 + 격자의 속성 테이블을 보면 'Join_Count'라는 것이 생긴 것을 볼 수 있는데, 이는 공간결합 연산을 사용하는 경우에 필연적으로 발생되는 것으로, 이는 각 격자와 결합된 포인트 데이터의 개수를 의미한다.In FIG. 10, it can be seen that 'Join_Count' is generated when the attribute table of the point + lattice that overlaps the point data and the lattice data is generated. This is inevitably generated when the spatial join operation is used. And the number of point data.

또한, 사용자단말기(10)의 제어모듈(11)은 상기 (S4) 단계에서 어떠한 포인트 데이터와 결합되지 않은 격자를 제외한 전체 격자의 개수를 산출한다(S5). 즉, (S5) 단계에서는 상기 (S4) 단계에서 어떠한 포인트 데이터와도 결합되지 않은 격자를 제외한 전체 격자의 개수를 산출하는데, 그 의미는 도 10에서, 포인트 + 격자의 속성 테이블에서 Join_Count의 값이 0인 격자의 개수를 산출한다는 것이다. 그 이유는 어떠한 포인트와도 결합되지 않은 격자는 클러스터로 표현될 필요가 없어서 (S3) 단계에서 도출된 적절한 클러스터의 개수에 포함되지 않기 때문이다.In addition, the control module 11 of the user terminal 10 calculates the total number of grids excluding the grids not combined with any point data at step S4 (S5). That is, in step (S5), the total number of grids excluding the grid not combined with any point data is calculated in step S4, meaning that the value of Join_Count in the point + grid attribute table is The number of grid points of 0 is calculated. The reason is that a grid not combined with any point does not need to be represented by a cluster and is not included in the number of appropriate clusters derived in the step S3.

그리고 사용자단말기(10)의 제어모듈(11)은 상기 (S5) 단계에서 Join_Count가 0인 격자를 제외하고 산출한 격자의 개수와, 상기 (S3) 단계에서 도출한 축척별, C값별 복수의 격자 개수 중에서 동일한 값이 있는지를 판단하되, 동일한 값이 있는 경우에는 축척별, C값별 복수(줌 레벨 19인 경우에는 클러스터링을 하지 않으므로 최대 19X19개가 될 수 있음)의 모든 격자 개수와 각각 동일한 상기 (S5) 단계에서의 격자 개수를 구할 때까지 연산을 반복 수행하고, 동일한 값이 없을 경우에는 상기 (S4) 단계로 돌아가 다른 크기의 격자데이터와 포인트 데이터에 대해 공간결합 연산과 상기 (S5) 단계를 반복적으로 수행한다(S6). 이로써 축척별, C값별 모든 격자 개수에 대한 격자데이터의 크기가 정해지게 된다.The control module 11 of the user terminal 10 calculates the number of grids excluding the grid with Join_Count equal to 0 and the number of grids according to the scale and C value derived in the step (S3) If there is the same value, it is judged whether or not there is the same value. If there is the same value, it is judged whether or not the same value exists (S5 If there is no identical value, the process returns to step S4 to repeat the spatial combining operation on the grid data and the point data of different sizes and the step (S5) repeatedly until the number of grids in the step (S6). Thus, the size of the grid data for all the grid numbers for each scale and C value is determined.

한편, 도 11에서, 상기 (S3) 단계 내지 (S6) 단계까지를 보다 이해하기 쉽게 설명한 것이다. 즉, 상기 (S3) 단계에서, C값에 따른 줌 레벨(축척)별 클러스터(격자)의 개수(n)를 산출한다. 예컨대, C값이 19개이고, 줌 레벨이 총 10단계라고 한다면 총 190개의 n가 산출된다. 그리고 (S4) 단계 및 (S5) 단계에서, 포인터 데이터와 격자데이터를 공간결합으로 연산하였을 때에 격자의 개수가 9개이면 공간결합 결과로 Joint_Count가 0인 격자의 개수는 2개이고, 클러스터로 표현된 격자의 개수는 7개가 된다. 따라서 단계(S6)에서 클러스터의 개수(n)과 클러스터로 표현될 격자의 개수가 일치하면 (S7) 단계를 수행하고, 일치하지 않으면 (S4) 단계와 (S5) 단계를 반복 수행한다.On the other hand, in FIG. 11, the steps (S3) to (S6) are explained more easily. That is, in step (S3), the number (n) of clusters (grids) for each zoom level (scale) according to the C value is calculated. For example, if the C value is 19 and the zoom level is 10 in total, a total of 190 n is calculated. If the number of grids is nine when the pointer data and the lattice data are computed by spatial combination at step S4 and step S5, the number of grids whose Joint_Count is 0 as a result of spatial combination is 2, The number of grids is seven. Therefore, if the number n of clusters coincides with the number of grids to be represented by clusters in step S6, step S7 is performed. Otherwise, steps S4 and S5 are repeated.

또한, 상기 (S3) 단계에서 축척별, C값별 복수의 모든 격자 개수와 각각 동일한 상기 (S5) 단계에서의 격자 개수가 구해지면, 사용자단말기(10)의 격자크기분포특성분석모듈(15)은 Moran's I 을 이용하여 상기 (S6) 단계에서 수행된 결과에 대하여 축척별, C값별로 분포 특성을 분석한다(S7). 즉 사용자단말기(10)는 격자크기분포특성모듈(15)에서 Moran's I 을 이용하여 상기 (S6) 단계에서 수행된 결과에 대해 축척별, C값별로 분포 특성을 분석한다. 따라서 Moran's I 는 다음의 수학식 3과 같이 계산할 수 있다.If the number of grids in step S5 is equal to the total number of grids for each scale and C value in step S3, the grid size distribution characteristic analysis module 15 of the user terminal 10 calculates Moran's I is used to analyze the distribution characteristics of the results performed in step S6 for each scale and C value (S7). That is, the user terminal 10 analyzes the distribution characteristics of each of the scales and the C values for the result of the step (S6) using Moran's I in the lattice size distribution characteristic module 15. [ Therefore, Moran's I can be calculated as shown in Equation 3 below.

Figure 112016061609799-pat00032
Figure 112016061609799-pat00032

여기서, nij로 이루어진 모든 공간 단위의 개수이고,

Figure 112016061609799-pat00033
는 공간가중행렬이며,
Figure 112016061609799-pat00034
i번째 변수(i번째 격자의 Joint_Count 값)이고,
Figure 112016061609799-pat00035
x의 평균이다.Where n is the number of all spatial units made up of i and j ,
Figure 112016061609799-pat00033
Is a space weighting matrix,
Figure 112016061609799-pat00034
Is the i- th variable (the Joint_Count value of the i-th grid)
Figure 112016061609799-pat00035
Is the average of x .

이때, 상기

Figure 112016061609799-pat00036
는 도 13의 공간가중행렬로, 도 12의 공간분포의 예시에서, 이웃하면 1 아니면 0의 값을 가지며, 자기 자신은 이웃으로 취급하지 않기 때문에 대각선의 값은 0이 된다. Moran's I 는 공간적 자기상관성을 측정하기 위한 통계량으로, 거의 -1과 1 사이의 값을 가지며, -1에 가까울수록 분산된 분포, 0에 가까울수록 랜덤한 분포, 1에 가까울수록 클러스터된 분포를 보인다고 해석한다.At this time,
Figure 112016061609799-pat00036
Is a space weighting matrix of FIG. 13, in the example of the spatial distribution in FIG. 12, has a value of 1 or 0 in the neighborhood, and the value of the diagonal line is 0 because it does not treat itself as a neighbor. Moran's I is a statistic for measuring spatial autocorrelation. It has a value between -1 and 1, and the closer to -1, the more scattered, the closer to 0, the closer to 1, the more clustered Interpret.

다음으로, 사용자단말기(10)의 제어모듈(11)은 상기 (S2) 단계에서 도출된 포인트 데이터의 본래 분포 특성과 축척별, 계수값(C)에 따른 상기 (S7) 단계의 분포 특성 분석 결과를 비교한다(S8). 그리고 사용자단말기(10)의 제어모듈(11)은 상기 (S8) 단계에서 분석한 포인트 데이터의 본래 분포 특성을 가장 잘 유지한 계수값(C)을 최종적으로 도출한다(S9). 따라서 (S2) 단계에서 NNA의 R값을 이용하여 분석한 포인트 데이터의 본래 분포 특성과 (S7) 단계에서 Moran's I 를 통한 분석 결과를 비교하는 것이다. 예를 들어, R값이 1보다 작다면, 본래 포인트 데이터의 분포는 클러스터되어 있다는 의미이므로 (S7) 단계에서 계산한 Moran's I 의 값이 19개의 C값 중에서 1에 가까운 C값을 축척별로 최종 결정하게 되는 것이다. 그러므로 최종적으로 얻어지는 결과는 결정된 C값에 따른 축척별 적절한 클러스터(격자)의 한 변의 길이(l)이다. 즉, 사용자단말기(10)의 제어모듈(11)은 상기 (S4) 단계에서 (S9) 단계까지 클러스터링에 필요한 모든 축척에 대하여 계수값(C)별로 반복하는 과정을 거쳐 최종적인 결과는 축척별로 도출된 C값에 따른 적절한 격자((S6)단계에서 얻어진 (S4)단계의 격자데이터)의 한 변의 길이, 달리 말해 격자의 크기이다. 이로써 MAUP의 영향을 최소화할 수 있다.Next, the control module 11 of the user terminal 10 analyzes the distribution characteristics of the step S7 according to the original distribution characteristics of the point data derived in the step S2, (S8). Then, the control module 11 of the user terminal 10 finally derives a coefficient value C that best maintains the original distribution characteristic of the point data analyzed in the step S8 (S9). Therefore, the original distribution characteristic of the point data analyzed using the R value of the NNA at step S2 is compared with the analysis result through Moran's I at step S7. For example, if the R value is smaller than 1, since the distribution of the point data is originally meant to be clustered, the value of Moran's I calculated in the step (S7) It will be done. Therefore, the final result is the length ( l ) of one side of the appropriate cluster (grid) for each scale according to the determined C value. That is, the control module 11 of the user terminal 10 repeats the steps of (S4) to (S9) for all the scales required for clustering for each coefficient value C, (The lattice data of step (S4) obtained in the step (S6)) according to the value of C, that is, the length of one side of the lattice, in other words, the size of the lattice. This minimizes the effect of the MAUP.

이상의 설명에서 본 발명은 특정의 실시 예와 관련하여 도시 및 설명하였지만, 청구범위에 의해 나타난 발명의 사상 및 영역으로부터 벗어나지 않는 한도 내에서 다양한 개조 및 변화가 가능하다는 것을 이 기술분야에서 통상의 지식을 가진 자라면 누구나 쉽게 알 수 있을 것이다.While the invention has been shown and described with respect to the specific embodiments thereof, it will be understood by those skilled in the art that various changes and modifications may be made without departing from the spirit and scope of the invention as defined by the appended claims. Anyone with it will know easily.

10: 사용자단말기 11: 제어모듈 12: 포인트데이터수집모듈 13: 최근린분석모듈 14: 지도일반화연산모듈 15: 격자크기분포특성분석모듈 20: 위치 기반 소셜 네트워크(LBSN)10: user terminal 11: control module 12: point data collection module 13: current lean analysis module 14: map generalization calculation module 15: lattice size distribution characteristic analysis module 20: location based social network (LBSN)

Claims (5)

(a) 사용자단말기의 포인트데이터수집모듈이 위치 기반 소셜 네트워크(Location-Based Social Network, LBSN) 서비스에서 제공하는 오픈 API(Open Application Programming Interface)를 이용하여 수집을 원하는 지역의 위치정보를 포함하는 포인트 데이터를 수집하는 단계;
(b) 상기 사용자단말기의 최근린분석모듈이 최근린분석(Nearest Neighbor Analysis, NNA)의 최근린지수(Nearest Neighbor Index, R)를 이용하여 수집된 포인트 데이터들의 본래 공간적 분포 특성을 파악하고 분석하는 단계;
(c) 상기 사용자단말기의 지도일반화연산모듈이 Topfer's Radical Law를 이용하여 계수값(C)에 따른 복수의 클러스터(Cluster, 격자)의 개수를 축척별로 산출하는 단계;
(d) 상기 사용자단말기의 제어모듈은 상기 (a) 단계에서 수집을 원하는 지역을 대상으로 수집한 위치정보가 포함된 포인트 데이터와, 해당 범위에 대하여 다양한 크기별로 격자데이터를 생성하여 그 생성된 격자데이터에 대하여 공간결합(Spatial Join) 연산을 수행하는 단계;
(e) 상기 사용자단말기의 제어모듈은 상기 (d) 단계에서 어떠한 포인트 데이터와 결합되지 않은 격자를 제외한 전체 격자의 개수를 산출하는 단계;
(f) 상기 사용자단말기의 제어모듈은 상기 (e) 단계에서 산출한 격자의 개수와, 상기 (c) 단계에서 도출한 축척별, C값별 복수의 격자 개수 중에서 동일한 값이 있는지를 판단하되, 동일한 값이 있는 경우에는 축척별, C값별 복수의 모든 격자 개수와 각각 동일한 상기 (e) 단계에서의 격자 개수를 구할 때까지 연산을 반복 수행하고, 동일한 값이 없을 경우에는 상기 (d) 단계로 돌아가 다른 크기의 격자데이터와 포인트 데이터에 대해 공간결합 연산과 상기 (e) 단계를 반복적으로 수행하는 단계;
(g) 상기 사용자단말기의 격자크기분포특성분석모듈은 Moran's I를 이용하여 상기 (f) 단계에서 수행된 결과에 대하여 축척별, C값별로 분포 특성을 분석하는 단계;
(h) 상기 사용자단말기의 제어모듈은 상기 (b) 단계에서 도출된 포인트 데이터의 본래 분포 특성과 축척별 계수값(C)에 따른 상기 (g) 단계의 분포 특성 분석 결과를 비교하는 단계, 및
(i) 상기 사용자단말기의 제어모듈은 상기 (h) 단계에서 분석한 포인트 데이터의 본래 분포 특성을 가장 잘 유지한 계수값(C)을 최종적으로 도출하는 단계를 포함하여 이루어진, 위치정보가 포함된 포인트 데이터를 다축척의 웹 지도상에 클러스터링하기 위하여 격자의 크기를 결정하는 방법.
(a) The point data collection module of the user terminal uses the open application programming interface (API) provided by the location-based social network (LBSN) Collecting data;
(b) the current lean analysis module of the user terminal identifies and analyzes the original spatial distribution characteristics of the point data collected using the Nearest Neighbor Index (RN) of Nearest Neighbor Analysis (NNA) step;
(c) calculating a number of a plurality of clusters according to a scale value (C) according to a scale value using Topfer's Radical Law, by the map generalization operation module of the user terminal;
(d) The control module of the user terminal generates point data including positional information collected for a region desired to be collected in the step (a), grid data for various sizes for the corresponding range, Performing a Spatial Join operation on the data;
(e) the control module of the user terminal calculates the total number of grids excluding the grids not combined with any point data in the step (d);
(f) The control module of the user terminal determines whether the number of grids calculated in the step (e) is the same among a plurality of grid numbers for each scale and C value derived in the step (c) If there is a value, the operation is repeated until the number of grids in the step (e) that is equal to the number of all the plurality of grids for each scale and C value is obtained. If there is no value, the process returns to the step (d) Repeatedly performing the spatial combining operation and the step (e) on the grid data and the point data of different sizes;
(g) analyzing a distribution characteristic of the grid size distribution characteristic of the user terminal by using the Moran's I , analyzing distribution characteristics of the result performed in the step (f) for each scale and C value;
(h) the control module of the user terminal compares the distribution characteristic analysis result of the step (g) according to the original distribution characteristic of the point data derived in the step (b) and the coefficient value (C)
(i) the control module of the user terminal finally deriving a coefficient value (C) that best preserves the original distribution characteristic of the point data analyzed in the step (h) A method for determining the size of a grid to cluster point data on a multi-scale web map.
제1항에 있어서, (j) 상기 사용자단말기의 제어모듈은 상기 (d) 단계에서 (i) 단계까지 클러스터링에 필요한 모든 축척에 대하여 계수값(C)별로 반복하는 과정을 거쳐 최종적인 결과는 축척별로 도출된 C값에 따른 격자의 한 변의 길이(격자의 크기)를 얻을 수 있도록 연산하는 단계를 더 포함하는, 위치정보가 포함된 포인트 데이터를 다축척의 웹 지도상에 클러스터링하기 위하여 격자의 크기를 결정하는 방법.
2. The method of claim 1, wherein the control module of the user terminal repeats the process of (d) to (i) And calculating a length of one side of the lattice (size of the lattice) corresponding to the derived C values, wherein the step of calculating the size of the lattice in order to cluster the point data including the position information on the multi- / RTI >
제1항에 있어서, 상기 (b) 단계에서, 최근린지수(R)는 다음의 수학식,
Figure 112016061609799-pat00037

(여기서,
Figure 112016061609799-pat00038
는 각각의 포인트 데이터로부터 가장 가까운 포인트 데이터까지의 거리의 평균이고,
Figure 112016061609799-pat00039
i포인트 데이터로부터 가장 가까운 포인트 데이터까지의 거리이며, n은 전체 포인트의 개수이고,
Figure 112016061609799-pat00040
는 CSR 하에서 기대되는 최근린 거리이며, A는 수집을 원하는 지역의 전체 면적이다.)을 이용하여 도출하는, 위치정보가 포함된 포인트 데이터를 다축척의 웹 지도상에 클러스터링하기 위하여 격자의 크기를 결정하는 방법.
2. The method according to claim 1, wherein, in step (b), the latest Lindsey number (R)
Figure 112016061609799-pat00037

(here,
Figure 112016061609799-pat00038
Is the average of the distance from each point data to the closest point data,
Figure 112016061609799-pat00039
Is the distance from the i point data to the closest point data, n is the number of all points,
Figure 112016061609799-pat00040
Is the latest lean distance expected under the CSR, and A is the total area of the area desired to be collected). In order to cluster the point data including the position information on the multi-scale web map, How to determine.
제1항에 있어서, 상기 (c) 단계에서, Topfer's Radical Law는 다음의 수학식,
Figure 112016061609799-pat00041

(여기서,
Figure 112016061609799-pat00042
는 결과 지도의 포인트 수이고,
Figure 112016061609799-pat00043
는 입력 지도의 포인트 수이며,
Figure 112016061609799-pat00044
는 입력 지도의 축척계수이고,
Figure 112016061609799-pat00045
는 결과 지도의 축척계수이며, 그리고 C 는 심볼계수이다.)을 이용하여 도출하는, 위치정보가 포함된 포인트 데이터를 다축척의 웹 지도상에 클러스터링하기 위하여 격자의 크기를 결정하는 방법.
The method of claim 1, wherein, in step (c), Topfer's Radical Law is expressed by the following equation:
Figure 112016061609799-pat00041

(here,
Figure 112016061609799-pat00042
Is the number of points in the result map,
Figure 112016061609799-pat00043
Is the number of points in the input map,
Figure 112016061609799-pat00044
Is the scale factor of the input map,
Figure 112016061609799-pat00045
Is the scale factor of the result map, and C is the symbol factor). A method of determining the size of a grid to cluster point data including position information on a multi-scale web map.
제1항에 있어서, 상기 (g) 단계에서, Moran's I는 다음의 수학식,
Figure 112016061609799-pat00046

(여기서, nij로 이루어진 모든 공간 단위의 개수이고,
Figure 112016061609799-pat00047
는 공간가중행렬이며,
Figure 112016061609799-pat00048
i번째 변수(i번째 격자의 Joint_Count 값)이고,
Figure 112016061609799-pat00049
x의 평균이다.)을 이용하여 도출하는, 위치정보가 포함된 포인트 데이터를 다축척의 웹 지도상에 클러스터링하기 위하여 격자의 크기를 결정하는 방법.
2. The method of claim 1, wherein in step (g), Moran's I is calculated by the following equation:
Figure 112016061609799-pat00046

(Where n is the number of all spatial units made up of i and j ,
Figure 112016061609799-pat00047
Is a space weighting matrix,
Figure 112016061609799-pat00048
Is the i- th variable (the Joint_Count value of the i-th grid)
Figure 112016061609799-pat00049
Is the average of x .) A method for determining the size of a grid to cluster point data containing position information on a multi-scale web map.
KR1020160079891A 2016-06-27 2016-06-27 Method for Determining the Size of Grid for Clustering on Multi-Scale Web Map Services using Location-Based Point Data KR101721114B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020160079891A KR101721114B1 (en) 2016-06-27 2016-06-27 Method for Determining the Size of Grid for Clustering on Multi-Scale Web Map Services using Location-Based Point Data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020160079891A KR101721114B1 (en) 2016-06-27 2016-06-27 Method for Determining the Size of Grid for Clustering on Multi-Scale Web Map Services using Location-Based Point Data

Publications (1)

Publication Number Publication Date
KR101721114B1 true KR101721114B1 (en) 2017-03-30

Family

ID=58503047

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020160079891A KR101721114B1 (en) 2016-06-27 2016-06-27 Method for Determining the Size of Grid for Clustering on Multi-Scale Web Map Services using Location-Based Point Data

Country Status (1)

Country Link
KR (1) KR101721114B1 (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108615059A (en) * 2018-05-10 2018-10-02 中国人民解放军战略支援部队信息工程大学 A kind of lake automatically selecting method and device based on Dynamic Multiscale cluster
KR20190048790A (en) * 2017-10-31 2019-05-09 한국과학기술원 System and method for distributed processing spatial data
KR20190072105A (en) * 2017-12-15 2019-06-25 주식회사 상헌 Apparatus for providing 3d terrain and facility data using online map open api and method thereof
CN110399441A (en) * 2019-07-04 2019-11-01 佳都新太科技股份有限公司 A kind of massive point data aggregate rendering method, device, equipment and storage medium
CN112200355A (en) * 2020-09-30 2021-01-08 杭州鲁尔物联科技有限公司 Landslide prediction method, device, equipment and storage medium
CN113343518A (en) * 2021-05-21 2021-09-03 中国测绘科学研究院 Landscape space area integrity evaluation method and system based on adjacency relation analysis
US11893647B2 (en) 2017-04-27 2024-02-06 Snap Inc. Location-based virtual avatars

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013129580A1 (en) * 2012-02-28 2013-09-06 公立大学法人大阪府立大学 Approximate nearest neighbor search device, approximate nearest neighbor search method, and program
KR20140055663A (en) * 2012-11-01 2014-05-09 네이버비즈니스플랫폼 주식회사 Apparatus, method and computer readable recording medium of clustering data for displaying the information based on a lacation
US8983138B2 (en) * 2010-07-15 2015-03-17 Olympus Corporation Image processing device, information storage device, and image processing method
KR20160017016A (en) * 2016-01-27 2016-02-15 네이버 주식회사 Apparatus, method and computer readable recording medium of clustering data for displaying the information based on a lacation

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8983138B2 (en) * 2010-07-15 2015-03-17 Olympus Corporation Image processing device, information storage device, and image processing method
WO2013129580A1 (en) * 2012-02-28 2013-09-06 公立大学法人大阪府立大学 Approximate nearest neighbor search device, approximate nearest neighbor search method, and program
KR20140055663A (en) * 2012-11-01 2014-05-09 네이버비즈니스플랫폼 주식회사 Apparatus, method and computer readable recording medium of clustering data for displaying the information based on a lacation
KR20160017016A (en) * 2016-01-27 2016-02-15 네이버 주식회사 Apparatus, method and computer readable recording medium of clustering data for displaying the information based on a lacation

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
1998. 한국 GIS 학회지. 6(1). pp. 11-23.(점패턴분석을 이용한 수치지형도의 점사상 일반화). 유근배

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11893647B2 (en) 2017-04-27 2024-02-06 Snap Inc. Location-based virtual avatars
KR20190048790A (en) * 2017-10-31 2019-05-09 한국과학기술원 System and method for distributed processing spatial data
KR101980699B1 (en) 2017-10-31 2019-05-22 한국과학기술원 System and method for distributed processing spatial data
KR20190072105A (en) * 2017-12-15 2019-06-25 주식회사 상헌 Apparatus for providing 3d terrain and facility data using online map open api and method thereof
KR102046112B1 (en) * 2017-12-15 2019-11-18 주식회사 상헌 Apparatus for providing 3d terrain and facility data using online map open api and method thereof
CN108615059A (en) * 2018-05-10 2018-10-02 中国人民解放军战略支援部队信息工程大学 A kind of lake automatically selecting method and device based on Dynamic Multiscale cluster
CN108615059B (en) * 2018-05-10 2020-08-07 中国人民解放军战略支援部队信息工程大学 Lake automatic selection method and device based on dynamic multi-scale clustering
CN110399441A (en) * 2019-07-04 2019-11-01 佳都新太科技股份有限公司 A kind of massive point data aggregate rendering method, device, equipment and storage medium
CN112200355A (en) * 2020-09-30 2021-01-08 杭州鲁尔物联科技有限公司 Landslide prediction method, device, equipment and storage medium
CN112200355B (en) * 2020-09-30 2024-01-05 杭州鲁尔物联科技有限公司 Landslide prediction method, landslide prediction device, landslide prediction equipment and storage medium
CN113343518A (en) * 2021-05-21 2021-09-03 中国测绘科学研究院 Landscape space area integrity evaluation method and system based on adjacency relation analysis
CN113343518B (en) * 2021-05-21 2023-12-01 中国测绘科学研究院 Landscape space region integrity evaluation method and system based on adjacency analysis

Similar Documents

Publication Publication Date Title
KR101721114B1 (en) Method for Determining the Size of Grid for Clustering on Multi-Scale Web Map Services using Location-Based Point Data
US11042572B2 (en) System and method for spatial clustering using multiple-resolution grids
CN109033011B (en) Method and device for calculating track frequency, storage medium and electronic equipment
CN110019568B (en) Spatial clustering-based addressing method and device, computer equipment and storage medium
CN111858820B (en) Land property identification method and device, electronic equipment and storage medium
CN109241846B (en) Method and device for estimating space-time change of remote sensing image and storage medium
CN109376761A (en) The method for digging and device of a kind of address mark and its longitude and latitude
CN110019622A (en) A kind of real time mass map point polymerization and device
CN103888541A (en) Method and system for discovering cells fused with topology potential and spectral clustering
CN111475746B (en) Point-of-interest mining method, device, computer equipment and storage medium
Peters Quadtree-and octree-based approach for point data selection in 2D or 3D
CN109059949B (en) Shortest path calculation method and device
CN111125490A (en) Population activity number extraction method, device and medium based on Baidu thermodynamic diagram
Wang et al. A fast candidate viewpoints filtering algorithm for multiple viewshed site planning
CN114254950A (en) Telecommunication resource data processing method and device, electronic equipment and storage medium
CN109684425A (en) Declaration form station address matching process and device
WO2016037346A1 (en) Measuring and diagnosing noise in urban environment
CN112464040B (en) Graph structure recognition, visual display and display operation method and device
CN110674134B (en) Geographic information data storage method, query method and device
Beilschmidt et al. A linear-time algorithm for the aggregation and visualization of big spatial point data
KR101516631B1 (en) Visible neighbor query for spatial skyline in mobile environment
US20230306602A1 (en) Information generation method and apparatus, electronic device, and computer readable medium
WO2018004083A1 (en) Apparatus and method for estimating floating population using map search history
CN114610825A (en) Method and device for confirming associated grid set, electronic equipment and storage medium
CN114357102A (en) Road network data generation method and device

Legal Events

Date Code Title Description
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20200302

Year of fee payment: 4