JP6972667B2 - Regional characteristic prediction method, regional characteristic prediction device and regional characteristic prediction program - Google Patents

Regional characteristic prediction method, regional characteristic prediction device and regional characteristic prediction program Download PDF

Info

Publication number
JP6972667B2
JP6972667B2 JP2017109187A JP2017109187A JP6972667B2 JP 6972667 B2 JP6972667 B2 JP 6972667B2 JP 2017109187 A JP2017109187 A JP 2017109187A JP 2017109187 A JP2017109187 A JP 2017109187A JP 6972667 B2 JP6972667 B2 JP 6972667B2
Authority
JP
Japan
Prior art keywords
weight matrix
spatial weight
map information
mesh
regional characteristic
Prior art date
Legal status (The legal status 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 status listed.)
Active
Application number
JP2017109187A
Other languages
Japanese (ja)
Other versions
JP2018205956A (en
Inventor
恭子 只木
均 小森谷
克仁 中澤
哲義 塩田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujitsu Ltd
Original Assignee
Fujitsu Ltd
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 Fujitsu Ltd filed Critical Fujitsu Ltd
Priority to JP2017109187A priority Critical patent/JP6972667B2/en
Publication of JP2018205956A publication Critical patent/JP2018205956A/en
Application granted granted Critical
Publication of JP6972667B2 publication Critical patent/JP6972667B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Instructional Devices (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Description

本発明は、地域特性予測方法、地域特性予測装置及び地域特性予測プログラムに関する。 The present invention relates to a regional characteristic prediction method, a regional characteristic prediction device, and a regional characteristic prediction program.

地域における人の流れや自然の変化など、地域における時間的な変化を予測する場合、対象とする地域が小さくなるほどその周辺の地域の影響を無視できなくなってくる。そのため、従来、計量統計において、空間重み行列を用いる方法が知られている。 When predicting temporal changes in a region, such as the flow of people and changes in nature, the smaller the target area, the more the influence of the surrounding area cannot be ignored. Therefore, conventionally, a method using a spatial weight matrix is known in metric statistics.

特開2012−141953号公報Japanese Unexamined Patent Publication No. 2012-141953

空間重み行列には、地域が隣接しているかどうかで決定するものや、ある距離範囲に入っているかどうかで決定するものがあるが、いずれの場合にも、ある地域に隣接していることによる影響は空間重み行列において等方的に表現される。 Some spatial weight matrices are determined by whether or not the regions are adjacent, and some are determined by whether or not they are within a certain distance range. In either case, it depends on whether or not they are adjacent to a certain region. The effect is expressed isotropically in the spatial weight matrix.

しかしながら、ある地域に隣接していても、隣接する地域それぞれが同程度の影響を受けるとは限らない。したがって、隣接する地域の影響が等方的に表現された空間重み行列を用いても、地域における時間的な変化を精度よく予測することができない可能性がある。 However, even if they are adjacent to a certain area, each adjacent area is not necessarily affected to the same extent. Therefore, even if a spatial weighting matrix that isotropically expresses the influence of adjacent regions is used, it may not be possible to accurately predict temporal changes in the regions.

本発明は、地域の特性を精度よく予測することが可能な地域特性予測方法、地域特性予測装置及び地域特性予測プログラムを提供することを目的とする。 An object of the present invention is to provide a regional characteristic prediction method, a regional characteristic prediction device, and a regional characteristic prediction program capable of accurately predicting regional characteristics.

一つの態様では、地域特性予測方法は、地図情報を格納する第1格納部を参照して、前記地図情報に含まれる細分化された地域間の隣接関係を示す空間重み行列を生成し、地図情報と対応付けて前記細分化された地域間の結び付きに影響を与える要素の情報を格納する第2格納部を参照して、前記空間重み行列を補正し、補正した前記空間重み行列と、地図情報と対応付けられた予測対象の特性に関するデータとに基づいて、前記予測対象の特性を予測する、処理をコンピュータが実行するものである。 In one embodiment, the regional characteristic prediction method refers to a first storage unit that stores map information to generate a spatial weighting matrix that indicates the adjacency between the subdivided regions contained in the map information, and maps. The spatial weight matrix is corrected, and the corrected spatial weight matrix and the map are referred to with reference to the second storage unit that stores the information of the elements that affect the connection between the subdivided regions in association with the information. The computer executes a process of predicting the characteristics of the prediction target based on the information and the data related to the characteristics of the prediction target associated with the information.

一つの側面として、地域の特性を精度よく予測することができる。 As one aspect, it is possible to accurately predict the characteristics of a region.

第1の実施形態に係る地域特性予測装置のハードウェア構成を概略的に示す図である。It is a figure which shows schematic the hardware composition of the area characteristic prediction apparatus which concerns on 1st Embodiment. 地域特性予測装置の機能ブロック図である。It is a functional block diagram of a regional characteristic prediction device. 図3(a)は、メッシュ地図DBに含まれるメッシュデータの一例を示す図であり、図3(b)は、メッシュの中心座標の一例を示す図である。FIG. 3A is a diagram showing an example of mesh data included in the mesh map DB, and FIG. 3B is a diagram showing an example of the center coordinates of the mesh. 図4(a)は、路線DBに含まれる路線重ね合わせメッシュデータを示す図であり、図4(b)は、路線DBに含まれる路線データを示す図である。FIG. 4A is a diagram showing route superimposed mesh data included in the route DB, and FIG. 4B is a diagram showing route data included in the route DB. 属性データの一例を示す図である。It is a figure which shows an example of attribute data. 地域特性予測装置の処理を示すフローチャートである。It is a flowchart which shows the processing of a region characteristic prediction apparatus. 図7(a)は、図6の説明に用いる9つのメッシュを示す図であり、図7(b)は、図7(a)のメッシュにおける空間重み行列を示す図である。7 (a) is a diagram showing nine meshes used in the explanation of FIG. 6, and FIG. 7 (b) is a diagram showing a spatial weight matrix in the mesh of FIG. 7 (a). 図8(a)は、図6の説明に用いる路線重ね合わせメッシュデータを示す図であり、図8(b)は、図8(a)を用いて補正した空間重み行列を示す図である。8 (a) is a diagram showing line superposition mesh data used in the explanation of FIG. 6, and FIG. 8 (b) is a diagram showing a spatial weight matrix corrected using FIG. 8 (a). 図9(a)は、比較例における予測結果を示すグラフであり、図9(b)は、第1の実施形態における予測結果を示すグラフである。FIG. 9A is a graph showing the prediction result in the comparative example, and FIG. 9B is a graph showing the prediction result in the first embodiment. 実測値を示すグラフである。It is a graph which shows the measured value. 図10と、図9(a)及び図9(b)とを比較した結果を示す表である。It is a table which shows the result of having compared FIG. 10 with FIG. 9 (a) and FIG. 9 (b). 第2の実施形態で用いる川が重ね合わせられたメッシュデータを示す図である。It is a figure which shows the mesh data which superposed the river used in the 2nd Embodiment. 図13(a)は、比較例における予測結果を示すグラフであり、図13(b)は、第2の実施形態における予測結果を示すグラフである。FIG. 13A is a graph showing the prediction result in the comparative example, and FIG. 13B is a graph showing the prediction result in the second embodiment. 第2の実施形態における実測値を示すグラフである。It is a graph which shows the measured value in 2nd Embodiment. 変形例を示す図である。It is a figure which shows the modification.

《第1の実施形態》
以下、地域特性予測装置の第1の実施形態について、図1〜図11に基づいて詳細に説明する。
<< First Embodiment >>
Hereinafter, the first embodiment of the regional characteristic prediction device will be described in detail with reference to FIGS. 1 to 11.

図1には、地域特性予測装置10のハードウェア構成が概略的に示されている。地域特性予測装置10は、例えばPC(Personal Computer)であり、図1に示すように、CPU(Central Processing Unit)90、ROM(Read Only Memory)92、RAM(Random Access Memory)94、記憶部(ここではHDD(Hard Disk Drive))96、ネットワークインタフェース97、表示部93、入力部95、及び可搬型記憶媒体用ドライブ99等を備えている。これら地域特性予測装置10の構成各部は、バス98に接続されている。地域特性予測装置10では、ROM92あるいはHDD96に格納されているプログラム(地域特性予測プログラムを含む)、或いは可搬型記憶媒体用ドライブ99が可搬型記憶媒体91から読み取ったプログラム(地域特性予測プログラムを含む)をCPU90が実行することにより、図2に示す、各部の機能が実現される。なお、図2においては、HDD96等に格納されている各種DB、データについても図示されている。 FIG. 1 schematically shows the hardware configuration of the regional characteristic prediction device 10. The regional characteristic prediction device 10 is, for example, a PC (Personal Computer), and as shown in FIG. 1, a CPU (Central Processing Unit) 90, a ROM (Read Only Memory) 92, a RAM (Random Access Memory) 94, and a storage unit ( Here, an HDD (Hard Disk Drive) 96, a network interface 97, a display unit 93, an input unit 95, a portable storage medium drive 99, and the like are provided. Each component of the regional characteristic prediction device 10 is connected to the bus 98. In the regional characteristic prediction device 10, a program stored in the ROM 92 or the HDD 96 (including the regional characteristic prediction program) or a program read from the portable storage medium 91 by the portable storage medium drive 99 (including the regional characteristic prediction program). ) Is executed by the CPU 90, so that the functions of each part shown in FIG. 2 are realized. Note that FIG. 2 also shows various DBs and data stored in the HDD 96 and the like.

図2は、地域特性予測装置10のCPU90がプログラムを実行することにより実現される機能を示す機能ブロック図である。図2に示すように、CPU90がプログラムを実行することで、生成部としての空間重み行列生成部12、補正部14、及び予測部16としての機能が実現されている。 FIG. 2 is a functional block diagram showing a function realized by the CPU 90 of the regional characteristic prediction device 10 executing a program. As shown in FIG. 2, by executing the program by the CPU 90, the functions as the spatial weight matrix generation unit 12, the correction unit 14, and the prediction unit 16 as the generation unit are realized.

空間重み行列生成部12は、第1格納部としてのメッシュ地図DB22に格納されているメッシュデータを参照して、空間重み行列を生成する。ここで、メッシュ地図DB22には、図3(a)に示すようなメッシュデータが格納されている。メッシュデータは、地図上に配置された、例えば縦横が1kmの矩形状のメッシュに関する情報であり、各メッシュの寸法や、各メッシュの中心座標(X,Y)を含む。なお、各メッシュにより地図に含まれる領域が複数の地域に細分化されている。空間重み行列生成部12は、メッシュデータから、図3(b)に示すような各メッシュの中心座標を取得する。そして、中心座標間の距離に基づいて、各メッシュの隣接関係を特定し、特定した隣接関係(隣接しているかどうか)に基づいて空間重み行列を生成する。 The spatial weight matrix generation unit 12 generates a spatial weight matrix with reference to the mesh data stored in the mesh map DB 22 as the first storage unit. Here, the mesh map DB 22 stores mesh data as shown in FIG. 3A. The mesh data is information about a rectangular mesh arranged on a map, for example, 1 km in length and width, and includes the dimensions of each mesh and the center coordinates (X, Y) of each mesh. The area included in the map is subdivided into a plurality of areas by each mesh. The spatial weight matrix generation unit 12 acquires the center coordinates of each mesh as shown in FIG. 3B from the mesh data. Then, the adjacency relationship of each mesh is specified based on the distance between the center coordinates, and the spatial weight matrix is generated based on the specified adjacency relationship (whether or not they are adjacent).

補正部14は、第2格納部としての路線DB24を参照して、空間重み行列を補正する。ここで、路線DB24には、図4(a)に示すような路線重ね合わせメッシュデータと、図4(b)に示すような路線データと、が含まれている。路線重ね合わせメッシュデータは、図4(a)に示すように、メッシュデータと路線の配置とを重ね合わせたデータである。なお、図4(a)では、2本の路線がメッシュデータに重ね合わされている。路線データは、図4(a)の路線重ね合わせメッシュデータから生成されるデータであり、図4(b)に示すように、「路線名」のフィールドと「メッシュID」のフィールドを有している。路線データでは、路線の名称と、路線が通過するメッシュの識別子(メッシュID)とが関連付けられている。 The correction unit 14 corrects the spatial weight matrix with reference to the line DB 24 as the second storage unit. Here, the route DB 24 includes route superimposed mesh data as shown in FIG. 4A and route data as shown in FIG. 4B. As shown in FIG. 4A, the line superimposition mesh data is data in which the mesh data and the line arrangement are superposed. In FIG. 4A, the two lines are superimposed on the mesh data. The route data is data generated from the route superimposition mesh data of FIG. 4 (a), and has a field of "route name" and a field of "mesh ID" as shown in FIG. 4 (b). There is. In the route data, the name of the route and the identifier (mesh ID) of the mesh through which the route passes are associated with each other.

補正部14では、2つの隣接するメッシュを同一の路線が通過している場合に、これら2つの隣接するメッシュに対応する空間重み行列の要素を補正する。例えば、2つのメッシュに対応する空間重み行列の要素をα倍(αは1よりも大きい値)する。このような補正を行うことで、路線の存在によるメッシュ間の結び付きの大きさ(強さ)を、空間重み行列において表現することができる。 The correction unit 14 corrects the elements of the spatial weight matrix corresponding to these two adjacent meshes when the same line passes through the two adjacent meshes. For example, the elements of the spatial weight matrix corresponding to the two meshes are multiplied by α (α is a value larger than 1). By making such a correction, the magnitude (strength) of the connection between the meshes due to the existence of the line can be expressed in the spatial weight matrix.

予測部16は、補正部14により補正された空間重み行列と、属性データ26とに基づいて、各地域の特性(例えば、各メッシュにおけるひとり世帯数)に関するデータの予測を行う。ここで、図5には、属性データ26の一例として、各年における各メッシュのひとり世帯数のデータが示されている。予測部16は、図5の属性データ26を用いて、例えば次の年の各メッシュにおけるひとり世帯数を予測する。なお、属性データ26は、予測部16が予測する特性(予測対象の特性)に関する情報を地図情報(各メッシュ)と対応付けたデータであるといえる。 The prediction unit 16 predicts data regarding the characteristics of each region (for example, the number of single households in each mesh) based on the spatial weight matrix corrected by the correction unit 14 and the attribute data 26. Here, FIG. 5 shows data on the number of single households in each mesh in each year as an example of the attribute data 26. The prediction unit 16 predicts, for example, the number of single households in each mesh of the next year by using the attribute data 26 of FIG. It can be said that the attribute data 26 is data in which information related to the characteristics (characteristics of the prediction target) predicted by the prediction unit 16 is associated with map information (each mesh).

(地域特性予測装置10の処理について)
次に、地域特性予測装置10の処理について、図6のフローチャートに沿って、その他図面を適宜参照しつつ詳細に説明する。なお、図6の処理が行われる前提として、メッシュ地図DB22及び路線DB24には、特性を予測したい地域に関する情報が格納されているものとする。また、属性データ26についても、入力部95を介して、ユーザにより入力済みか、メッシュ地図DB22内にdbfファイルとして含まれているものとする。なお、本処理の説明においては、説明の簡素化のため、メッシュデータが図7(a)のようなデータであり、路線重ね合わせメッシュデータが図8(a)のようなデータであるものとする。
(Regarding the processing of the regional characteristic prediction device 10)
Next, the processing of the regional characteristic prediction device 10 will be described in detail with reference to other drawings as appropriate along the flowchart of FIG. As a premise that the processing of FIG. 6 is performed, it is assumed that the mesh map DB 22 and the route DB 24 store information about the area for which the characteristics are to be predicted. Further, it is assumed that the attribute data 26 has already been input by the user via the input unit 95 or is included as a dbf file in the mesh map DB 22. In the description of this process, for the sake of simplification of the explanation, the mesh data is the data as shown in FIG. 7 (a), and the route superimposition mesh data is the data as shown in FIG. 8 (a). do.

図6の処理では、まずステップS12において、空間重み行列生成部12が、メッシュの中心座標を取得する。具体的には、図7(a)のメッシュデータから、各メッシュの識別子であるメッシュIDと、中心座標(X,Y)とを取得する(図3(b)参照)。 In the process of FIG. 6, first, in step S12, the spatial weight matrix generation unit 12 acquires the center coordinates of the mesh. Specifically, the mesh ID, which is an identifier of each mesh, and the center coordinates (X, Y) are acquired from the mesh data of FIG. 7 (a) (see FIG. 3 (b)).

次いで、ステップS14では、空間重み行列生成部12が、メッシュ間の隣接関係を決定する。本実施形態では、2つのメッシュの中心座標間の距離がメッシュの寸法と一致する場合に2つのメッシュが隣接していると判断する。例えば、図7(a)に示すように9つのメッシュ1〜9がある場合、メッシュ5と隣接するメッシュは、2,4,6,8となる。また、メッシュ1と隣接するメッシュは、2,4となる。なお、上記説明では、空間重み行列生成部12は、メッシュ間の隣接関係を2つのメッシュの中心座標間の距離に基づいて判断する場合について説明したが、これに限られるものではない。例えば、空間重み行列生成部12は、2つのメッシュが線で接しているか否かを判断し、線で接している場合に隣接すると判断してもよい。 Next, in step S14, the spatial weight matrix generation unit 12 determines the adjacency relationship between the meshes. In the present embodiment, when the distance between the center coordinates of the two meshes matches the dimensions of the meshes, it is determined that the two meshes are adjacent to each other. For example, when there are nine meshes 1 to 9 as shown in FIG. 7A, the meshes adjacent to the mesh 5 are 2, 4, 6, and 8. The number of meshes adjacent to mesh 1 is 2 and 4. In the above description, the space weight matrix generation unit 12 has described the case where the adjacency relationship between the meshes is determined based on the distance between the center coordinates of the two meshes, but the present invention is not limited to this. For example, the spatial weight matrix generation unit 12 may determine whether or not the two meshes are in contact with each other by a line, and if they are in contact with each other by a line, it may be determined that they are adjacent to each other.

次いで、ステップS16では、空間重み行列生成部12が、空間重み行列を生成する。具体的には、空間重み行列生成部12は、ステップS14において決定したメッシュ間の隣接関係を行列として表す。図7(a)の場合であれば、空間重み行列生成部12は、図7(b)に示すような空間重み行列を生成する。図7(b)の空間重み行列では、例えば、行「1」については、メッシュ1と隣接するメッシュの番号の列「2」、「4」にある行列要素の値を1とし、その他を0としている。また、例えば、行「2」については、メッシュ2と隣接するメッシュの番号の列「1」、「3」、「5」にある行列要素の値を1とし、その他を0としている。以下、同様にして、各行列要素の値が設定されている。 Next, in step S16, the spatial weight matrix generation unit 12 generates a spatial weight matrix. Specifically, the spatial weight matrix generation unit 12 represents the adjacency relationship between the meshes determined in step S14 as a matrix. In the case of FIG. 7 (a), the spatial weight matrix generation unit 12 generates the spatial weight matrix as shown in FIG. 7 (b). In the spatial weight matrix of FIG. 7B, for example, for row "1", the values of the matrix elements in the columns "2" and "4" of the mesh numbers adjacent to mesh 1 are set to 1, and the others are set to 0. It is supposed to be. Further, for example, for row "2", the values of the matrix elements in the columns "1", "3", and "5" of the mesh numbers adjacent to mesh 2 are set to 1, and the others are set to 0. Hereinafter, the values of each matrix element are set in the same manner.

なお、ステップS16の処理が終了すると、以降のステップS18〜S30において、補正部14により、空間重み行列を補正する処理が実行される。 When the process of step S16 is completed, the process of correcting the spatial weight matrix is executed by the correction unit 14 in the subsequent steps S18 to S30.

ステップS18に移行すると、補正部14は、未特定のメッシュを特定する。例えば、補正部14は、図7(a)のメッシュ1を特定したものとする。 When the process proceeds to step S18, the correction unit 14 identifies an unspecified mesh. For example, it is assumed that the correction unit 14 specifies the mesh 1 in FIG. 7 (a).

次いで、ステップS20では、補正部14が、特定したメッシュに路線が通っているか否かを判断する。この場合、補正部14は、図8(a)の路線重ね合わせメッシュデータを参照し、メッシュ1に路線が通っているか否かを判断する。図8(a)の例では特定されているメッシュ1に路線が通っていないため、このステップS20の判断は否定され、ステップS30に移行する。ステップS30に移行すると、補正部14は、全てのメッシュを特定済みか否かを判断する。このステップS30の判断が否定された場合には、ステップS18に戻る。ステップS18に戻ると、未特定のメッシュを特定する。ここでは、例えば、メッシュ2を特定したものとする。 Next, in step S20, the correction unit 14 determines whether or not the route passes through the specified mesh. In this case, the correction unit 14 refers to the line superimposition mesh data of FIG. 8A, and determines whether or not the line passes through the mesh 1. In the example of FIG. 8A, since the route does not pass through the specified mesh 1, the determination in step S20 is denied, and the process proceeds to step S30. When the process proceeds to step S30, the correction unit 14 determines whether or not all the meshes have been specified. If the determination in step S30 is denied, the process returns to step S18. Returning to step S18, the unspecified mesh is specified. Here, for example, it is assumed that the mesh 2 is specified.

次いで、ステップS20に移行すると、特定したメッシュ2に路線が通っているか否かを判断する。図8(a)の例では、メッシュ2に路線が通っているため、ステップS20の判断は肯定され、ステップS22に移行する。 Next, when the process proceeds to step S20, it is determined whether or not the route passes through the specified mesh 2. In the example of FIG. 8A, since the line passes through the mesh 2, the determination in step S20 is affirmed, and the process proceeds to step S22.

ステップS22に移行すると、補正部14は、未選択の隣接メッシュを選択する。ここでは、メッシュ2に隣接するメッシュ1,3,5のうちの1つ(例えばメッシュ1)を選択する。 When the process proceeds to step S22, the correction unit 14 selects an unselected adjacent mesh. Here, one of meshes 1, 3, and 5 adjacent to mesh 2 (for example, mesh 1) is selected.

次いで、ステップS24では、補正部14は、選択した隣接メッシュに同一の路線が通っているか否かを判断する。ここで、現時点で選択されている隣接メッシュ1にはステップS18で特定したメッシュ2と同一の路線が通っていないため、ステップS24の判断は否定されて、ステップS28に移行する。ステップS28では、補正部14が、全ての隣接メッシュを選択済みか否かを判断するが、ここでは、メッシュ1,3,5のうちのメッシュ1を選択したのみであるので、判断は否定され、ステップS22に戻る。 Next, in step S24, the correction unit 14 determines whether or not the same line passes through the selected adjacent mesh. Here, since the same route as the mesh 2 specified in step S18 does not pass through the adjacent mesh 1 currently selected, the determination in step S24 is denied and the process proceeds to step S28. In step S28, the correction unit 14 determines whether or not all adjacent meshes have been selected, but here, since only mesh 1 of meshes 1, 3 and 5 is selected, the determination is denied. , Return to step S22.

ステップS22に戻ると、補正部14は、未選択の隣接メッシュとして、例えば、メッシュ3を選択する。次いで、ステップS24では、補正部14は、選択した隣接メッシュ3にステップS18で特定したメッシュ2と同一の路線が通っているか否かを判断する。ここで、現時点で選択されている隣接メッシュ3にはメッシュ2と同一の路線が通っている。したがって、この場合にはステップS24の判断が肯定され、ステップS26に移行する。 Returning to step S22, the correction unit 14 selects, for example, mesh 3 as an unselected adjacent mesh. Next, in step S24, the correction unit 14 determines whether or not the same line as the mesh 2 specified in step S18 passes through the selected adjacent mesh 3. Here, the same line as the mesh 2 runs through the adjacent mesh 3 currently selected. Therefore, in this case, the determination in step S24 is affirmed, and the process proceeds to step S26.

ステップS26に移行した場合、補正部14は、対応する行列要素にα(α>1)を積算する。すなわち、補正部14は、図7(b)の上から2行目、左から3列目の行列要素「1」をα(=1×α)とする(図8(b)参照)。 When the process proceeds to step S26, the correction unit 14 integrates α (α> 1) into the corresponding matrix element. That is, the correction unit 14 sets α (= 1 × α) to the matrix element “1” in the second row from the top and the third column from the left in FIG. 7 (b) (see FIG. 8 (b)).

次いで、ステップS28では、補正部14が、全ての隣接メッシュを選択済みか否かを判断するが、ここでは、メッシュ1,3,5のうちのメッシュ1,3を選択したのみであるので、判断は否定され、ステップS22に戻る。 Next, in step S28, the correction unit 14 determines whether or not all adjacent meshes have been selected, but here, only meshes 1 and 3 out of meshes 1, 3 and 5 have been selected. The judgment is denied, and the process returns to step S22.

ステップS22に戻ると、補正部14は、未選択の隣接メッシュとして、例えば、メッシュ5を選択する。次いで、ステップS24では、補正部14は、選択した隣接メッシュ5にメッシュ2と同一の路線が通っているか否かを判断する。ここで、現時点で選択されている隣接メッシュ5にはメッシュ2と同一の路線が通っている。したがって、この場合にはステップS24の判断が肯定され、ステップS26に移行する。 Returning to step S22, the correction unit 14 selects, for example, mesh 5 as an unselected adjacent mesh. Next, in step S24, the correction unit 14 determines whether or not the same line as the mesh 2 passes through the selected adjacent mesh 5. Here, the same line as the mesh 2 runs through the adjacent mesh 5 currently selected. Therefore, in this case, the determination in step S24 is affirmed, and the process proceeds to step S26.

ステップS26に移行した場合、補正部14は、対応する行列要素にα(α>1)を積算する。すなわち、補正部14は、図7(b)の上から2行目、左から5列目の行列要素「1」をα(=1×α)とする(図8(b)参照)。 When the process proceeds to step S26, the correction unit 14 integrates α (α> 1) into the corresponding matrix element. That is, the correction unit 14 sets α (= 1 × α) to the matrix element “1” in the second row from the top and the fifth column from the left in FIG. 7 (b) (see FIG. 8 (b)).

その後はステップS28に移行するが、全ての隣接メッシュを選択済みであるので、ステップS28の判断は肯定され、ステップS30に移行する。ステップS30に移行すると、前述したように、補正部14は、全てのメッシュを特定済みか否かを判断するが、ここでの判断が否定されると、ステップS18に戻る。 After that, the process proceeds to step S28, but since all the adjacent meshes have been selected, the judgment of step S28 is affirmed, and the process proceeds to step S30. When the process proceeds to step S30, as described above, the correction unit 14 determines whether or not all the meshes have been specified, but if the determination here is denied, the process returns to step S18.

以降、補正部14は、メッシュ3〜9を順次特定し(S18)、特定したメッシュに隣接するメッシュを順次選択しつつ(S22)、特定したメッシュと選択した隣接メッシュとに同一の路線が通っているか否かを判断する(S24)。そして、同一の路線が通っていた場合には、対応する行列要素にαを積算すること(S26)で、空間重み行列を補正する。 After that, the correction unit 14 sequentially identifies the meshes 3 to 9 (S18), sequentially selects the meshes adjacent to the specified meshes (S22), and the same line passes through the specified mesh and the selected adjacent meshes. It is determined whether or not it is (S24). Then, when the same line passes, the spatial weight matrix is corrected by integrating α into the corresponding matrix elements (S26).

なお、図8(a)の例の場合、図8(b)に示すように空間重み行列が補正されるようになっている。 In the case of the example of FIG. 8A, the spatial weight matrix is corrected as shown in FIG. 8B.

上述したように全てのメッシュが特定済みとなり、ステップS30の判断が肯定されると、ステップS32に移行し、予測部16が、属性データ26と補正後の空間重み行列を用いて地域の特性を予測する。予測部16は、例えば、統計モデリングのためのプログラミング言語(例えばStan)を用いたベイズ推定を実施する機能を有し、属性データ26と補正後の空間重み行列に基づいて、地域の特性の予測を実行する。例えば、属性データ26が、過去における各メッシュでのひとり世帯数の推移を示すデータである場合には、予測部16は、属性データ26と補正後の空間重み行列とから、例えば各メッシュにおける翌年のひとり世帯数を予測する。この場合、2つの隣接するメッシュ間を同一の路線が通っている場合に、2つのメッシュに対応する空間重み行列の行列要素の値を大きく補正しているので、予測部16は、補正後の空間重み行列を用いることで、路線の存在によるメッシュ間の影響(例えば人の分布が路線に沿って広がる傾向にあること)を考慮して地域の特性を予測することができるようになっている。 As described above, when all the meshes have been specified and the judgment of step S30 is affirmed, the process proceeds to step S32, and the prediction unit 16 uses the attribute data 26 and the corrected spatial weight matrix to determine the characteristics of the region. Predict. The prediction unit 16 has a function of performing Bayesian estimation using, for example, a programming language for statistical modeling (for example, Stan), and predicts the characteristics of the region based on the attribute data 26 and the corrected spatial weight matrix. To execute. For example, when the attribute data 26 is data showing the transition of the number of single households in each mesh in the past, the prediction unit 16 uses the attribute data 26 and the corrected spatial weight matrix, for example, the next year in each mesh. Predict the number of single households. In this case, when the same line passes between two adjacent meshes, the value of the matrix element of the spatial weight matrix corresponding to the two meshes is greatly corrected, so that the prediction unit 16 has corrected the corrected value. By using the spatial weight matrix, it is possible to predict the characteristics of the region in consideration of the influence between the meshes due to the existence of the line (for example, the distribution of people tends to spread along the line). ..

以上のように、ステップS32までの処理が終了すると、図6の全処理が終了する。 As described above, when the processing up to step S32 is completed, all the processing in FIG. 6 is completed.

なお、図4(a)に示すように、分析対象の複数のメッシュの中に、複数の路線が通っている場合もある。このような場合には、着目する路線を異ならせて、図6のステップS18〜S30の処理を路線分だけ繰り返すこととすればよい。このとき、路線ごとにαを異ならせてもよい。例えば、在来線ほどαを大きくし、新幹線のような地域間に影響を及ぼす可能性の低い路線ほどαを小さくすることとしてもよい。 As shown in FIG. 4A, a plurality of routes may pass through the plurality of meshes to be analyzed. In such a case, the routes of interest may be different, and the processes of steps S18 to S30 in FIG. 6 may be repeated for the routes. At this time, α may be different for each line. For example, it may be possible to increase α for conventional lines and decrease α for lines such as the Shinkansen, which are less likely to affect areas.

図9(a)には、図3(a)のようにある自治体に設定したメッシュそれぞれにおけるひとり世帯数の平成7年(1995年)から平成17年(2005年)までのデータ(図5の属性データ26)と、空間重み行列(補正せず)を用いて、平成22年(2010年)の各メッシュにおけるひとり世帯数を予測した結果(平成22年のみ予測結果、その他は属性データ26に含まれる実測値)が示されている。また、図9(b)には、第1の実施形態のように補正した空間重み行列を用いて、図9(a)と同様に平成22年(2010年)の各メッシュにおけるひとり世帯数を予測した結果が示されている。さらに、図10には、各メッシュのひとり世帯数の実測値が示されている。 FIG. 9 (a) shows the data of the number of single households in each mesh set in a certain municipality as shown in FIG. 3 (a) from 1995 (1995) to 2005 (2005) (FIG. 5). Attribute data 26) and the result of predicting the number of single households in each mesh in 2010 (2010) using the spatial weight matrix (without correction) (prediction result only in 2010, others are attribute data 26) Measured values included) are shown. Further, in FIG. 9 (b), the number of single households in each mesh in 2010 (2010) is shown in the same manner as in FIG. 9 (a) by using the spatial weight matrix corrected as in the first embodiment. The predicted results are shown. Further, FIG. 10 shows an actually measured value of the number of single households in each mesh.

図11には、図9(a)の予測値と図10の実測値との差(差:比較例)、及び図9(b)の予測値と図10の実測値との差(差:本願)を示している。図11の表から、差:比較例の絶対値の平均を求めると、78.8となり、差:本願の絶対値の平均を求めると27.0となった。この結果、第1の実施形態のように補正した空間重み行列を用いることで、ひとり世帯数を精度よく予測することができることが分かった。 11 shows the difference between the predicted value of FIG. 9A and the measured value of FIG. 10 (difference: comparative example), and the difference between the predicted value of FIG. 9B and the measured value of FIG. 10 (difference:). This application) is shown. From the table of FIG. 11, the difference: the average of the absolute values of the comparative examples was 78.8, and the difference: the average of the absolute values of the present application was 27.0. As a result, it was found that the number of single households can be accurately predicted by using the spatial weighting matrix corrected as in the first embodiment.

以上、詳細に説明したように、第1の実施形態によると、空間重み行列生成部12は、メッシュ地図DB22を参照して、細分化された地域(メッシュ)間の隣接関係を示す空間重み行列を生成し(S16)、補正部14は、路線DB24を参照して、空間重み行列を補正する(S18〜S30)。そして、予測部16は、補正した空間重み行列と、属性データとに基づいて、地域の特性を予測する。これにより、第1の実施形態では、メッシュ間の隣接関係のみならず、隣接するメッシュ間の結び付きに路線が与える影響(非対称性)を考慮して、属性データから地域の特性を予測することができる。したがって、地域の特性を精度よく予測することが可能である。 As described above in detail, according to the first embodiment, the spatial weight matrix generation unit 12 refers to the mesh map DB 22 and indicates the spatial weight matrix showing the adjacency relationship between the subdivided regions (mesh). (S16), and the correction unit 14 corrects the spatial weight matrix with reference to the line DB 24 (S18 to S30). Then, the prediction unit 16 predicts the characteristics of the region based on the corrected spatial weight matrix and the attribute data. As a result, in the first embodiment, it is possible to predict the characteristics of the region from the attribute data in consideration of not only the adjacency relationship between the meshes but also the influence (asymmetry) of the route on the connection between the adjacent meshes. can. Therefore, it is possible to accurately predict the characteristics of the region.

また、第1の実施形態では、路線が通るメッシュ間の結びつきが強くなるように、すなわち路線が通るメッシュ間の重みが大きくなるように、空間重み行列を補正するため、人の分布が路線沿いに広がる傾向にあること等を考慮して、地域の特性を精度よく予測することができる。 Further, in the first embodiment, the distribution of people is along the line in order to correct the spatial weight matrix so that the connection between the meshes through which the line passes is strong, that is, the weight between the meshes through which the line passes is large. It is possible to accurately predict the characteristics of the area in consideration of the fact that it tends to spread to the area.

なお、上記第1の実施形態では、メッシュ間を同一の路線が通るか否かに基づいて、行列要素を補正する場合について説明したが、これに限らず、その他の交通網(例えばバスの運行経路や幹線道路など)がメッシュ間を通るか否かに基づいて行列要素を補正することとしてもよい。また、複数の交通網に基づいて行列要素を補正してもよく、その場合には、交通網の種別ごとに補正値(α)を変更することとしてもよい。さらに、路線に沿って人の分布が広がる傾向は都心ほど高く、幹線道路に沿って人の分布が広がる傾向は郊外ほど高いなどの地域性があるため、このような地域性を考慮して補正値(α)を調整することとしてもよい。 In the first embodiment, the case where the matrix element is corrected based on whether or not the same route passes between the meshes has been described, but the present invention is not limited to this, and other transportation networks (for example, bus operation) are described. Matrix elements may be corrected based on whether a route (such as a route or a highway) passes between meshes. Further, the matrix element may be corrected based on a plurality of transportation networks, and in that case, the correction value (α) may be changed for each type of transportation network. Furthermore, the tendency for the distribution of people to spread along the route is higher in the city center, and the tendency for the distribution of people to spread along the trunk road is higher in the suburbs. The value (α) may be adjusted.

なお、例えばメッシュを通る路線の長さを特定し、路線がわずかにしか通っていないメッシュの予測結果の精度が低いような場合には、予測結果の精度が上がるように、路線の長さに応じてαを調整してもよい。また、メッシュを通る路線の長さが所定の長さ未満である場合に補正を行わないようにしてもよい。 In addition, for example, if the length of the line passing through the mesh is specified and the accuracy of the prediction result of the mesh that the line passes only slightly is low, the length of the line is set so that the accuracy of the prediction result is improved. Α may be adjusted accordingly. Further, if the length of the line passing through the mesh is less than a predetermined length, the correction may not be performed.

《第2の実施形態》
次に、第2の実施形態について、図12〜図14に基づいて説明する。
<< Second Embodiment >>
Next, the second embodiment will be described with reference to FIGS. 12 to 14.

上記第1の実施形態では、同一の路線が通っている隣接するメッシュの結び付きが強くなるように、空間重み行列を補正する場合について説明したが、第2の実施形態では、河川による地域の分断(往来を阻害する要素)による影響を加味して、空間重み行列を補正する例について説明する。 In the first embodiment described above, a case where the spatial weight matrix is corrected so that the connection between adjacent meshes passing through the same line is strengthened has been described, but in the second embodiment, the division of the area by the river has been described. An example of correcting the spatial weight matrix in consideration of the influence of (elements that hinder traffic) will be described.

第2の実施形態では、補正部14は、隣接する2つのメッシュの間に川が流れている場合に、その2つのメッシュに対応する行列要素の値をβ倍(β<1)するように補正する。例えば、図12に示すような川が重ね合わせられたメッシュデータにおいて、メッシュ5に隣接するメッシュは2,4,6,8となるが、メッシュ5とメッシュ2,6との間に川が存在し、メッシュ間を分断している。したがって、メッシュ5、2及びメッシュ5,6に対応する行列要素の値をβ倍する。また、その他のメッシュについても同様に補正を行う。そして、予測部16は、補正後の空間重み行列を用いて、地域の特性(例えば完全失業者数)を予測する。 In the second embodiment, the correction unit 14 causes the value of the matrix element corresponding to the two meshes to be multiplied by β (β <1) when the river flows between the two adjacent meshes. to correct. For example, in the mesh data in which rivers are superimposed as shown in FIG. 12, the number of meshes adjacent to mesh 5 is 2, 4, 6, 8, but there is a river between mesh 5 and mesh 2, 6. However, the mesh is divided. Therefore, the values of the matrix elements corresponding to meshes 5, 2 and meshes 5, 6 are multiplied by β. The same correction is performed for other meshes. Then, the prediction unit 16 predicts the characteristics of the region (for example, the number of unemployed people) using the corrected spatial weight matrix.

図13(a)には、ある自治体に設定したメッシュそれぞれにおける完全失業者数(平成7年から平成17年)のデータと、空間重み行列(補正せず)を用いて平成22年の完全失業者数を予測した結果(比較例)が示されている。また、図13(b)には、図13(a)の場合と同様のデータと、補正した空間重み行列とを用いて、平成22年の完全失業者数を予測した結果(第2の実施形態)が示されている。さらに、図14には、平成7年から平成22年までの完全失業者数の推移(実測値)が示されている。 In FIG. 13 (a), the data of the number of unemployed people (1995 to 2005) in each mesh set in a certain municipality and the space weight matrix (without correction) are used to show the total unemployment in 2010. The result of predicting the number of people (comparative example) is shown. Further, in FIG. 13 (b), the result of predicting the number of unemployed people in 2010 using the same data as in the case of FIG. 13 (a) and the corrected spatial weight matrix (second implementation). Form) is shown. Further, FIG. 14 shows the transition (actual measurement value) of the number of unemployed people from 1995 to 2010.

これらの結果から、第1の実施形態と同様、比較例の予測値と実測値との差の絶対値の平均値を求めると、10.0となり、第2の実施形態の予測値と実測値との差の絶対値の平均値を求めると、9.6となった。すなわち、第2の実施形態のように空間重み行列を補正することで、予測精度に4%の向上が見られた。 From these results, as in the first embodiment, the average value of the absolute value of the difference between the predicted value and the measured value of the comparative example is 10.0, which is 10.0, which is the predicted value and the measured value of the second embodiment. The average value of the absolute values of the difference from and was 9.6. That is, by correcting the spatial weight matrix as in the second embodiment, the prediction accuracy was improved by 4%.

以上、説明したように、第2の実施形態によると、メッシュ間の結び付きに影響を与える要素として、メッシュ間の往来を阻害する要素(例えば川)の情報を用い、メッシュ間に川がある場合にメッシュ間の結び付きを弱めるように行列要素を補正するので、メッシュ間の往来を阻害する要素を考慮して、地域の特性を精度よく予測することができる。 As described above, according to the second embodiment, information on an element (for example, a river) that hinders traffic between meshes is used as an element that affects the connection between meshes, and there is a river between meshes. Since the matrix elements are corrected so as to weaken the connection between the meshes, it is possible to accurately predict the characteristics of the region in consideration of the factors that hinder the traffic between the meshes.

なお、上記実施形態では、地域の特性として、完全失業者数を予測する場合について説明したが、これに限らず、例えば、地域の特性として、植生や生物の生息域などを予測することとしてもよい。 In the above embodiment, the case of predicting the number of unemployed people as a characteristic of the region has been described, but the present invention is not limited to this, and for example, as a characteristic of the region, it is also possible to predict the habitat of vegetation and organisms. good.

なお、上記実施形態では、メッシュ間の往来を阻害する要素が川である場合を例にとり説明したが、これに限らず、メッシュ間の往来を阻害する要素は、例えばダムや飛行場、基地などであってもよい。 In the above embodiment, the case where the element that hinders the traffic between meshes is a river has been described as an example, but the present invention is not limited to this, and the element that hinders the traffic between meshes is, for example, a dam, an airfield, a base, or the like. There may be.

なお、上記第1、第2の実施形態では、予測部16がベイズ推定により地域の特性を予測する場合について説明したが、これに限らず、その他の方法により地域の特性を予測することとしてもよい。 In the first and second embodiments, the case where the prediction unit 16 predicts the characteristics of the region by Bayesian estimation has been described, but the present invention is not limited to this, and the characteristics of the region may be predicted by other methods. good.

なお、上記第1、第2の実施形態を組み合わせてもよい。すなわち、メッシュ間を交通網がつなぐ場合には、行列要素の値が大きくなるように補正し、メッシュ間に川などの往来を阻害する要素がある場合には、行列要素の値が小さくなるように補正してもよい。 The first and second embodiments may be combined. That is, when the traffic network connects the meshes, the value of the matrix element is corrected to be large, and when there is an element such as a river between the meshes that obstructs the traffic, the value of the matrix element is reduced. It may be corrected to.

なお、上記第1、第2の実施形態では、地域特性予測装置10が、図2の各機能を有する場合について説明したが、これに限られるものではない。例えば、図2の各機能を、図15に示すような、インターネットなどのネットワーク80に接続されたクラウドサーバ110が有していてもよい。この場合、利用者端末70から入力されたデータをクラウドサーバ110が受信し、クラウドサーバ110において図6の処理を実行することとしてもよい。そして、クラウドサーバ110から利用者端末70に処理結果(予測結果)を送信し、利用者端末70にて受信した処理結果を利用する。なお、クラウドサーバ110は、国内、海外のいずれに設置されてもよい。 In the first and second embodiments, the case where the regional characteristic prediction device 10 has each function of FIG. 2 has been described, but the present invention is not limited to this. For example, each function of FIG. 2 may be possessed by a cloud server 110 connected to a network 80 such as the Internet as shown in FIG. In this case, the cloud server 110 may receive the data input from the user terminal 70 and execute the process of FIG. 6 on the cloud server 110. Then, the processing result (prediction result) is transmitted from the cloud server 110 to the user terminal 70, and the processing result received by the user terminal 70 is used. The cloud server 110 may be installed in Japan or overseas.

なお、上記の処理機能は、コンピュータによって実現することができる。その場合、処理装置が有すべき機能の処理内容を記述したプログラムが提供される。そのプログラムをコンピュータで実行することにより、上記処理機能がコンピュータ上で実現される。処理内容を記述したプログラムは、コンピュータで読み取り可能な記録媒体(ただし、搬送波は除く)に記録しておくことができる。 The above processing function can be realized by a computer. In that case, a program that describes the processing content of the function that the processing device should have is provided. By executing the program on a computer, the above processing function is realized on the computer. The program describing the processing content can be recorded on a computer-readable recording medium (however, the carrier wave is excluded).

プログラムを流通させる場合には、例えば、そのプログラムが記録されたDVD(Digital Versatile Disc)、CD−ROM(Compact Disc Read Only Memory)などの可搬型記録媒体の形態で販売される。また、プログラムをサーバコンピュータの記憶装置に格納しておき、ネットワークを介して、サーバコンピュータから他のコンピュータにそのプログラムを転送することもできる。 When a program is distributed, it is sold in the form of a portable recording medium such as a DVD (Digital Versatile Disc) or a CD-ROM (Compact Disc Read Only Memory) on which the program is recorded. It is also possible to store the program in the storage device of the server computer and transfer the program from the server computer to another computer via the network.

プログラムを実行するコンピュータは、例えば、可搬型記録媒体に記録されたプログラムもしくはサーバコンピュータから転送されたプログラムを、自己の記憶装置に格納する。そして、コンピュータは、自己の記憶装置からプログラムを読み取り、プログラムに従った処理を実行する。なお、コンピュータは、可搬型記録媒体から直接プログラムを読み取り、そのプログラムに従った処理を実行することもできる。また、コンピュータは、サーバコンピュータからプログラムが転送されるごとに、逐次、受け取ったプログラムに従った処理を実行することもできる。 The computer that executes the program stores, for example, the program recorded on the portable recording medium or the program transferred from the server computer in its own storage device. Then, the computer reads the program from its own storage device and executes the processing according to the program. The computer can also read the program directly from the portable recording medium and execute the processing according to the program. Further, the computer can also sequentially execute the processing according to the received program each time the program is transferred from the server computer.

上述した実施形態は本発明の好適な実施の例である。但し、これに限定されるものではなく、本発明の要旨を逸脱しない範囲内において種々変形実施可能である。 The embodiments described above are examples of preferred embodiments of the present invention. However, the present invention is not limited to this, and various modifications can be made without departing from the gist of the present invention.

なお、以上の第1、第2の実施形態の説明に関して、更に以下の付記を開示する。
(付記1) 地図情報を格納する第1格納部を参照して、前記地図情報に含まれる細分化された地域間の隣接関係を示す空間重み行列を生成し、
地図情報と対応付けて前記細分化された地域間の結び付きに影響を与える要素の情報を格納する第2格納部を参照して、前記空間重み行列を補正し、
補正した前記空間重み行列と、地図情報と対応付けられた予測対象の特性に関するデータとに基づいて、前記予測対象の特性を予測する、
処理をコンピュータが実行することを特徴とする地域特性予測方法。
(付記2) 前記細分化された地域間の結び付きに影響を与える要素の情報は、交通網に関する情報を含むことを特徴とする付記1に記載の地域特性予測方法。
(付記3) 前記補正する処理では、前記交通網が通る地域間の重みが大きくなるように前記空間重み行列を補正する、ことを特徴とする付記2に記載の地域特性予測方法。
(付記4) 前記細分化された地域間の結び付きに影響を与える要素の情報は、前記地域間の往来を阻害する要素の情報を含むことを特徴とする付記1〜3のいずれかに記載の地域特性予測方法。
(付記5) 前記補正する処理では、往来が阻害されている地域間の重みが小さくなるように前記空間重み行列を補正する、ことを特徴とする付記4に記載の地域特性予測方法。
(付記6) 地図情報を格納する第1格納部を参照して、前記地図情報に含まれる細分化された地域間の隣接関係を示す空間重み行列を生成する生成部と、
地図情報と対応付けて前記細分化された地域間の結び付きに影響を与える要素の情報を格納する第2格納部を参照して、前記空間重み行列を補正する補正部と、
補正した前記空間重み行列と、地図情報と対応付けられた予測対象の特性に関するデータとに基づいて、前記予測対象の特性を予測する予測部と、
を備える地域特性予測装置。
(付記7) 前記細分化された地域間の結び付きに影響を与える要素の情報は、交通網に関する情報を含むことを特徴とする付記6に記載の地域特性予測装置。
(付記8) 前記補正部は、前記交通網が通る地域間の重みが大きくなるように前記空間重み行列を補正する、ことを特徴とする付記7に記載の地域特性予測装置。
(付記9) 前記細分化された地域間の結び付きに影響を与える要素の情報は、前記地域間の往来を阻害する要素の情報を含むことを特徴とする付記6〜8のいずれかに記載の地域特性予測装置。
(付記10) 前記補正部は、往来が阻害されている地域間の重みが小さくなるように前記空間重み行列を補正する、ことを特徴とする付記9に記載の地域特性予測装置。
(付記11) 地図情報を格納する第1格納部を参照して、前記地図情報に含まれる細分化された地域間の隣接関係を示す空間重み行列を生成し、
地図情報と対応付けて前記細分化された地域間の結び付きに影響を与える要素の情報を格納する第2格納部を参照して、前記空間重み行列を補正し、
補正した前記空間重み行列と、地図情報と対応付けられた予測対象の特性に関するデータとに基づいて、前記予測対象の特性を予測する、
処理をコンピュータに実行させるための地域特性予測プログラム。
Regarding the above description of the first and second embodiments, the following additional notes will be further disclosed.
(Appendix 1) With reference to the first storage unit for storing the map information, a spatial weight matrix showing the adjacency relationship between the subdivided areas included in the map information is generated.
The spatial weight matrix is corrected by referring to the second storage unit that stores the information of the elements that affect the connection between the subdivided regions in association with the map information.
The characteristics of the prediction target are predicted based on the corrected spatial weight matrix and the data related to the characteristics of the prediction target associated with the map information.
A regional characteristic prediction method characterized by a computer performing processing.
(Appendix 2) The regional characteristic prediction method according to Appendix 1, wherein the information of the elements affecting the connection between the subdivided regions includes information on the transportation network.
(Supplementary Note 3) The regional characteristic prediction method according to Supplementary note 2, wherein in the correction process, the spatial weight matrix is corrected so that the weight between the areas through which the transportation network passes becomes large.
(Appendix 4) The description in any of the appendices 1 to 3, wherein the information of the elements that influence the connection between the subdivided regions includes the information of the elements that hinder the traffic between the regions. Regional characteristic prediction method.
(Supplementary Note 5) The regional characteristic prediction method according to Supplementary note 4, wherein in the correction process, the spatial weight matrix is corrected so that the weights between regions where traffic is hindered become small.
(Appendix 6) With reference to the first storage unit for storing map information, a generation unit for generating a spatial weight matrix indicating an adjacency relationship between subdivided regions included in the map information, and a generation unit.
A correction unit that corrects the spatial weight matrix with reference to a second storage unit that stores information on elements that affect the connection between the subdivided regions in association with the map information.
A prediction unit that predicts the characteristics of the prediction target based on the corrected spatial weight matrix and data related to the characteristics of the prediction target associated with the map information.
Regional characteristic prediction device equipped with.
(Appendix 7) The regional characteristic prediction device according to Appendix 6, wherein the information of the elements affecting the connection between the subdivided regions includes information on the transportation network.
(Supplementary Note 8) The regional characteristic prediction device according to Supplementary note 7, wherein the correction unit corrects the spatial weight matrix so that the weight between areas through which the transportation network passes becomes large.
(Appendix 9) The description in any of the appendices 6 to 8, wherein the information of the elements that influence the connection between the subdivided regions includes the information of the elements that hinder the traffic between the regions. Regional characteristic predictor.
(Supplementary Note 10) The regional characteristic prediction device according to Supplementary note 9, wherein the correction unit corrects the spatial weight matrix so that the weights between regions where traffic is hindered become small.
(Appendix 11) With reference to the first storage unit for storing the map information, a spatial weight matrix showing the adjacency relationship between the subdivided areas included in the map information is generated.
The spatial weight matrix is corrected by referring to the second storage unit that stores the information of the elements that affect the connection between the subdivided regions in association with the map information.
The characteristics of the prediction target are predicted based on the corrected spatial weight matrix and the data related to the characteristics of the prediction target associated with the map information.
A regional characteristic prediction program that allows a computer to perform processing.

10 地域特性予測装置
12 空間重み行列生成部(生成部)
14 補正部
16 予測部
22 メッシュ地図DB(第1格納部)
24 路線DB(第2格納部)
26 属性データ(予想対象の特性に関するデータ)
10 Regional characteristic prediction device 12 Spatial weight matrix generation unit (generation unit)
14 Correction unit 16 Prediction unit 22 Mesh map DB (1st storage unit)
24 Line DB (2nd storage)
26 Attribute data (data related to the characteristics of the forecast target)

Claims (8)

地図情報を格納する第1格納部を参照して、前記地図情報に含まれる細分化された地域間の隣接関係を示す空間重み行列を生成し、
地図情報と対応付けて前記細分化された地域間の結び付きに影響を与える要素の情報を格納する第2格納部を参照して、前記空間重み行列を補正し、
補正した前記空間重み行列と、地図情報と対応付けられた予測対象の特性に関するデータとに基づいて、前記予測対象の特性を予測する、
処理をコンピュータが実行することを特徴とする地域特性予測方法。
With reference to the first storage unit for storing the map information, a spatial weight matrix showing the adjacency relationship between the subdivided regions included in the map information is generated.
The spatial weight matrix is corrected by referring to the second storage unit that stores the information of the elements that affect the connection between the subdivided regions in association with the map information.
The characteristics of the prediction target are predicted based on the corrected spatial weight matrix and the data related to the characteristics of the prediction target associated with the map information.
A regional characteristic prediction method characterized by a computer performing processing.
前記細分化された地域間の結び付きに影響を与える要素の情報は、交通網に関する情報を含むことを特徴とする請求項1に記載の地域特性予測方法。 The regional characteristic prediction method according to claim 1, wherein the information on the factors that influence the connection between the subdivided regions includes information on the transportation network. 前記補正する処理では、前記交通網が通る地域間の重みが大きくなるように前記空間重み行列を補正する、ことを特徴とする請求項2に記載の地域特性予測方法。 The regional characteristic prediction method according to claim 2, wherein in the correction process, the spatial weight matrix is corrected so that the weight between the areas through which the transportation network passes becomes large. 前記細分化された地域間の結び付きに影響を与える要素の情報は、前記地域間の往来を阻害する要素の情報を含むことを特徴とする請求項1〜3のいずれか一項に記載の地域特性予測方法。 The region according to any one of claims 1 to 3, wherein the information of the element that influences the connection between the subdivided regions includes the information of the element that hinders the traffic between the regions. Characteristic prediction method. 前記補正する処理では、往来が阻害されている地域間の重みが小さくなるように前記空間重み行列を補正する、ことを特徴とする請求項4に記載の地域特性予測方法。 The regional characteristic prediction method according to claim 4, wherein in the correction process, the spatial weight matrix is corrected so that the weight between regions where traffic is hindered becomes small. 前記予測する処理では、補正した前記空間重み行列と、地図情報と対応付けられた各地域の所定の特性に関するデータとに基づいて、予測対象の地域の前記所定の特性に関するデータを予測する、ことを特徴とする請求項1〜5のいずれか一項に記載の地域特性予測方法。In the prediction process, the data regarding the predetermined characteristics of the region to be predicted is predicted based on the corrected spatial weight matrix and the data regarding the predetermined characteristics of each region associated with the map information. The regional characteristic prediction method according to any one of claims 1 to 5, wherein the method is characterized by. 地図情報を格納する第1格納部を参照して、前記地図情報に含まれる細分化された地域間の隣接関係を示す空間重み行列を生成する生成部と、
地図情報と対応付けて前記細分化された地域間の結び付きに影響を与える要素の情報を格納する第2格納部を参照して、前記空間重み行列を補正する補正部と、
補正した前記空間重み行列と、地図情報と対応付けられた予測対象の特性に関するデータとに基づいて、前記予測対象の特性を予測する予測部と、
を備える地域特性予測装置。
With reference to the first storage unit that stores the map information, a generation unit that generates a spatial weight matrix indicating the adjacency relationship between the subdivided regions included in the map information, and a generation unit.
A correction unit that corrects the spatial weight matrix with reference to a second storage unit that stores information on elements that affect the connection between the subdivided regions in association with the map information.
A prediction unit that predicts the characteristics of the prediction target based on the corrected spatial weight matrix and data related to the characteristics of the prediction target associated with the map information.
Regional characteristic prediction device equipped with.
地図情報を格納する第1格納部を参照して、前記地図情報に含まれる細分化された地域間の隣接関係を示す空間重み行列を生成し、
地図情報と対応付けて前記細分化された地域間の結び付きに影響を与える要素の情報を格納する第2格納部を参照して、前記空間重み行列を補正し、
補正した前記空間重み行列と、地図情報と対応付けられた予測対象の特性に関するデータとに基づいて、前記予測対象の特性を予測する、
処理をコンピュータに実行させるための地域特性予測プログラム。
With reference to the first storage unit for storing the map information, a spatial weight matrix showing the adjacency relationship between the subdivided regions included in the map information is generated.
The spatial weight matrix is corrected by referring to the second storage unit that stores the information of the elements that affect the connection between the subdivided regions in association with the map information.
The characteristics of the prediction target are predicted based on the corrected spatial weight matrix and the data related to the characteristics of the prediction target associated with the map information.
A regional characteristic prediction program that allows a computer to perform processing.
JP2017109187A 2017-06-01 2017-06-01 Regional characteristic prediction method, regional characteristic prediction device and regional characteristic prediction program Active JP6972667B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2017109187A JP6972667B2 (en) 2017-06-01 2017-06-01 Regional characteristic prediction method, regional characteristic prediction device and regional characteristic prediction program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2017109187A JP6972667B2 (en) 2017-06-01 2017-06-01 Regional characteristic prediction method, regional characteristic prediction device and regional characteristic prediction program

Publications (2)

Publication Number Publication Date
JP2018205956A JP2018205956A (en) 2018-12-27
JP6972667B2 true JP6972667B2 (en) 2021-11-24

Family

ID=64957150

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017109187A Active JP6972667B2 (en) 2017-06-01 2017-06-01 Regional characteristic prediction method, regional characteristic prediction device and regional characteristic prediction program

Country Status (1)

Country Link
JP (1) JP6972667B2 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109738604B (en) * 2019-01-14 2021-04-27 北京工商大学 Water area characteristic model establishing method based on space weight matrix
JP7328650B2 (en) 2019-03-29 2023-08-17 彰洋 佐藤 Systems, methods, and programs for generating mesh statistics utilizing data associated with location information
CN111914052B (en) * 2020-08-04 2022-11-29 中国人民解放军战略支援部队信息工程大学 Method and device for determining regional spatial incidence relation

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000322402A (en) * 1999-05-11 2000-11-24 Hitachi Ltd Stream of people analyzing method
JP4552099B2 (en) * 2001-03-09 2010-09-29 清水建設株式会社 Hospital user number prediction system
JP2003296540A (en) * 2002-03-29 2003-10-17 Kumagai Gumi Co Ltd Diagnostic area analyzing system
JP5024959B2 (en) * 2008-05-30 2012-09-12 学校法人早稲田大学 Optimal route search device, optimal route search device and program
JP2011158985A (en) * 2010-01-29 2011-08-18 Ohbayashi Corp Connectiveness evaluation method for supply facility
JP2014187624A (en) * 2013-03-25 2014-10-02 National Institute Of Information & Communication Technology Repair method and repair program of network
US20160063516A1 (en) * 2014-08-29 2016-03-03 The Nielsen Company (Us), Llc Methods and apparatus to estimate commercial characteristics based on geospatial data
JP6757561B2 (en) * 2015-09-17 2020-09-23 株式会社パスコ Trade area forecasting device and trade area forecasting program

Also Published As

Publication number Publication date
JP2018205956A (en) 2018-12-27

Similar Documents

Publication Publication Date Title
JP6972667B2 (en) Regional characteristic prediction method, regional characteristic prediction device and regional characteristic prediction program
Moon et al. Future changes in monsoon duration and precipitation using CMIP6
Chen et al. Global land projection based on plant functional types with a 1-km resolution under socio-climatic scenarios
CN106649331B (en) Business circle identification method and equipment
JP5871193B2 (en) Monitoring data analysis apparatus, monitoring data analysis method, and monitoring data analysis program
US8849738B2 (en) Predicting a user behavior number of a word
CN109522382B (en) Spatial data gridding statistical method and device
JP5871192B2 (en) Monitoring data analysis apparatus, monitoring data analysis method, and monitoring data analysis program
JP5421949B2 (en) Traffic volume prediction device, traffic volume prediction method and program
CN106202098A (en) The method and apparatus of click location in record and reduction webpage
Peng et al. Climate change multi-model projections in CMIP6 scenarios in Central Hokkaido, Japan
Han et al. Conservation prioritization with machine learning predictions for the black-necked crane Grus nigricollis, a flagship species on the Tibetan Plateau for 2070
Dubrovsky et al. Developing a reduced-form ensemble of climate change scenarios for Europe and its application to selected impact indicators
Spickett et al. Adaptation strategies for health impacts of climate change in Western Australia: Application of a Health Impact Assessment framework
Rich et al. Are North American bird species' geographic ranges mainly determined by climate?
CN110058939A (en) System Expansion method, device and equipment
JP2007108809A (en) Time-series prediction system, time-series prediction method, and time-series prediction program
Noël et al. Peak refreezing in the Greenland firn layer under future warming scenarios
CN112468984B (en) Method for selecting address of power wireless private network base station and related equipment
McNicol et al. Protected areas reduce deforestation and degradation and enhance woody growth across African woodlands
Lu et al. Predicting long-term urban growth in Beijing (China) with new factors and constraints of environmental change under integrated stochastic and fuzzy uncertainties
CN110991867A (en) Marketing resource adjustment method and device based on gross benefit rate influence factor
CN113449217A (en) Method and equipment for migration track, thermodynamic diagram generation and mesh point determination
CN115422318A (en) Business data analysis method and device, storage medium and computer equipment
CN106023131A (en) Cartographic information analysis method and device

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200310

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20210129

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210209

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210412

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20211005

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20211018

R150 Certificate of patent or registration of utility model

Ref document number: 6972667

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150