特許 7328650 位置情報に関連付けられたデータを利用してメッシュ統計を生成するためのシステム、方法、およびプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2023-08-08

(45)【発行日】2023-08-17

(54)【発明の名称】位置情報に関連付けられたデータを利用してメッシュ統計を生成するためのシステム、方法、およびプログラム

(51)【国際特許分類】

   G06F 17/18 20060101AFI20230809BHJP

   G06Q 30/0203 20230101ALI20230809BHJP

   G06Q 50/10 20120101ALI20230809BHJP

【ＦＩ】

G06F17/18 Z

G06Q30/0203

G06Q50/10

【請求項の数】 15

(21)【出願番号】P 2019068605

(22)【出願日】2019-03-29

(65)【公開番号】P2020166739

(43)【公開日】2020-10-08

【審査請求日】2022-03-23

【国等の委託研究の成果に係る記載事項】（出願人による申告）平成２７年度、国立研究開発法人科学技術振興機構、戦略的創造研究推進事業「グローバル・システムの持続可能性評価基盤に関する研究」委託研究、産業技術力強化法第１９条の適用を受ける特許出願

(73)【特許権者】

【識別番号】522107382

【氏名又は名称】佐藤 彰洋

(73)【特許権者】

【識別番号】513055610

【氏名又は名称】一般財団法人宇宙システム開発利用推進機構

(74)【代理人】

【識別番号】100078282

【弁理士】

【氏名又は名称】山本 秀策

(74)【代理人】

【識別番号】100107489

【弁理士】

【氏名又は名称】大塩 竹志

(74)【代理人】

【識別番号】100113413

【弁理士】

【氏名又は名称】森下 夏樹

(74)【代理人】

【識別番号】100181674

【弁理士】

【氏名又は名称】飯田 貴敏

(74)【代理人】

【識別番号】100181641

【弁理士】

【氏名又は名称】石川 大輔

(74)【代理人】

【識別番号】230113332

【弁護士】

【氏名又は名称】山本 健策

(72)【発明者】

【氏名】佐藤 彰洋

(72)【発明者】

【氏名】毛利 浩一郎

【審査官】坂庭 剛史

(56)【参考文献】

【文献】特開２００９－２９４４７９（ＪＰ，Ａ）

【文献】特開２０１８－２０５９５６（ＪＰ，Ａ）

【文献】国際公開第２００８／０８８０００（ＷＯ，Ａ１）

【文献】特開２００１－０９３０７７（ＪＰ，Ａ）

【文献】国際公開第２０１８／２０７８７８（ＷＯ，Ａ１）

【文献】特開平０９－１４７０９１（ＪＰ，Ａ）

【文献】特開平１１－３３９１７４（ＪＰ，Ａ）

【文献】特開２００３－０９１６２０（ＪＰ，Ａ）

【文献】特開２００３－３１６８０８（ＪＰ，Ａ）

【文献】特開２０１３－２２８７０１（ＪＰ，Ａ）

【文献】佐藤彰洋、西村正貴、椿 広計，世界メッシュコードに基づく各種データの可視化，ＦＩＴ２０１７ 第１６回情報科学技術フォーラム 講演論文集 第１分冊，日本，一般社団法人情報処理学会，2017年09月05日，pp.29-34

【文献】平戸崇博、醍醐市朗、松野泰也、足立芳寛，ＧＩＳデータを用いた鉄鋼蓄積量の推計，材料とプロセス［ＣＤ－ＲＯＭ］2009年第157回春期講演大会発表(2009年3月)第29回学生ポスターセッション要約集，Vol.22, No.2，日本，社団法人日本鉄鋼協会，2009年09月01日，p.1（PS-15），ISSN 1882-8922

【文献】寺木彰浩、阪田知彦，グリッド・サンプリングによる建物群配置評価に関する基礎的研究，第２０回ＧＩＳＡ学術研究発表大会 講演論文集［ＣＤ－ＲＯＭ］，日本，一般社団法人地理情報システム学会，2011年10月15日，Vol.20/2011，pp.1-4（F-3-3）

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ｇ０６Ｆ １７／１８

Ｇ０６Ｑ １０／００－９９／００

(57)【特許請求の範囲】

【請求項1】

位置情報に関連付けられた少なくとも１つのデータを利用してメッシュ統計を生成するためのシステムであって、前記位置情報に関連付けられたデータは、衛星画像に含まれる衛星画像データであり、前記システムは、
前記衛星画像を受信する受信手段と、
前記衛星画像データが属する少なくとも１つのメッシュを決定する決定手段であって、前記少なくとも１つのメッシュの各メッシュに少なくとも１つのデータが属しており、前記決定手段は、
前記衛星画像の四隅の座標を決定することと、
北西端のメッシュを特定することと、
前記北西端のメッシュの中央座標（latc，longc）を求めることと、
前記衛星画像を囲むメッシュの数Ｎ、Ｍを求めることと、
前記衛星画像を覆うメッシュｗ _ｉ の集合（１≦ｉ≦Ｎ×Ｍ）を求めることであって、ｗ _ｉ は、中心座標（lati，longi）を有するメッシュであり、中心座標（lati，longi）は、
（lati，longi）＝（latc＋（yy－１）×Δlat，longc＋（xx－１）×Δlong）
と表され、ここで、Δlatは、前記少なくとも１つのメッシュの緯度方向幅であり、Δlongは、前記少なくとも１つのメッシュの経度方向幅であり、１≦yy≦Ｍであり、１≦xx≦Ｎである、ことと
によって、前記少なくとも１つのメッシュを決定する、決定手段と、
前記少なくとも１つのメッシュのうちの１つのメッシュに属する少なくとも１つのデータに基づいて、前記１つのメッシュ内の所与の事項の統計量を示すメッシュデータを生成する生成手段と、
前記メッシュデータを格納する格納手段と
を備える、システム。

【請求項2】

前記生成手段は、前記メッシュに属する前記データの統計量を算出することによって、前記メッシュデータを生成するように構成されている、請求項１に記載のシステム。

【請求項3】

前記メッシュデータは、前記メッシュ内の所与の事項の最大値、最小値、中央値、およびパーセンタイル値を含む、請求項１または請求項２に記載のシステム。

【請求項4】

前記生成手段は、前記複数のメッシュに対して並列的に前記メッシュデータを生成するように構成されている、請求項１～３のいずれか一項に記載のシステム。

【請求項5】

前記受信手段は、位置情報に関連付けられた第１のデータと位置情報に関連付けられた第２のデータとを受信し、
前記生成手段は、
（１）前記メッシュに属する第１のデータに基づいて第１のメッシュデータを生成することと、
（２）前記メッシュに属する第２のデータに基づいて第２のメッシュデータを生成することと、
（３）前記第１のメッシュデータと前記第２のメッシュデータとに演算処理を施すことにより複合メッシュデータを生成することと
を行うように構成されている、請求項１～３のいずれか一項に記載のシステム。

【請求項6】

前記格納手段に格納された複数のメッシュデータの値を集計することにより、所望の領域のメッシュ統計を提供する集計手段
をさらに備える、請求項１～５のいずれか一項に記載のシステム。

【請求項7】

前記集計手段は、
提供すべきメッシュ統計の指示を受信することと、
前記指示に基づいて、前記格納手段に格納されたメッシュデータのうち集計すべきメッシュデータを選択することと、
前記選択されたメッシュデータの値を集計することと
を行うように構成されている、請求項６に記載のシステム。

【請求項8】

前記集計手段は、所望の領域の形状に対する、前記メッシュデータが関連付けられたメッシュの貢献度に基づいて、前記メッシュデータの値を集計することにより、前記所望の領域のメッシュ統計を提供するように構成されている、請求項６または請求項７に記載のシステム。

【請求項9】

前記集計手段は、前記所望の領域の形状と、前記メッシュデータが関連付けられたメッシュとの包含関係に基づいて、前記所望の領域の形状に対する前記メッシュの貢献度を決定する、請求項８に記載のシステム。

【請求項10】

前記集計手段は、第１の事項の統計量を示す第１のメッシュデータと、第２の事項の統計量を示す第２のメッシュデータとに演算処理を施すことによって、第３の事項の統計量を示す第３のメッシュデータに基づくメッシュ統計を提供するように構成されている、請求項７～９のいずれか一項に記載のシステム。

【請求項11】

前記衛星画像は、地表形状の標高を表す標高画像、夜間光強度を表す夜間光画像、雲密度を表す雲画像、輝度温度を表す輝度温度画像、人工構造物を特定する人工構造物特定画像のうちの少なくとも１つを含む、請求項１～１０のいずれか一項に記載のシステム。

【請求項12】

位置情報に関連付けられた少なくとも１つのデータを利用してメッシュ統計を生成するための方法であって、前記位置情報に関連付けられたデータは、衛星画像に含まれる衛星画像データであり、前記方法は、
前記衛星画像を受信することと、
前記衛星画像データが属する少なくとも１つのメッシュを決定することであって、前記少なくとも１つのメッシュの各メッシュに少なくとも１つのデータが属しており、前記決定することは、
前記衛星画像の四隅の座標を決定することと、
北西端のメッシュを特定することと、
前記北西端のメッシュの中央座標（latc，longc）を求めることと、
前記衛星画像を囲むメッシュの数Ｎ、Ｍを求めることと、
前記衛星画像を覆うメッシュｗ _ｉ の集合（１≦ｉ≦Ｎ×Ｍ）を求めることであって、ｗ _ｉ は、中心座標（lati，longi）を有するメッシュであり、中心座標（lati，longi）は、
（lati，longi）＝（latc＋（yy－１）×Δlat，longc＋（xx－１）×Δlong）
と表され、ここで、Δlatは、前記少なくとも１つのメッシュの緯度方向幅であり、Δlongは、前記少なくとも１つのメッシュの経度方向幅であり、１≦yy≦Ｍであり、１≦xx≦Ｎである、ことと
によって、前記少なくとも１つのメッシュを決定することを含む、ことと、
前記少なくとも１つのメッシュのうちの１つのメッシュに属する少なくとも１つのデータに基づいて、前記１つのメッシュ内の所与の事項の統計量を示すメッシュデータを生成することと、
前記メッシュデータを格納することと
を含む、方法。

【請求項13】

格納された複数のメッシュデータの値を集計することにより、所望の領域のメッシュ統計を提供することをさらに含む、請求項１２に記載の方法。

【請求項14】

位置情報に関連付けられた少なくとも１つのデータを利用してメッシュ統計を生成するためのプログラムであって、前記位置情報に関連付けられたデータは、衛星画像に含まれる衛星画像データであり、前記プログラムは、システムにおいて実行され、前記システムは、プロセッサとデータベース部とを備え、前記プログラムは、
前記衛星画像を受信することと、
前記衛星画像データが属する少なくとも１つのメッシュを決定することであって、前記少なくとも１つのメッシュの各メッシュに少なくとも１つのデータが属しており、前記決定することは、
前記衛星画像の四隅の座標を決定することと、
北西端のメッシュを特定することと、
前記北西端のメッシュの中央座標（latc，longc）を求めることと、
前記衛星画像を囲むメッシュの数Ｎ、Ｍを求めることと、
前記衛星画像を覆うメッシュｗ _ｉ の集合（１≦ｉ≦Ｎ×Ｍ）を求めることであって、ｗ _ｉ は、中心座標（lati，longi）を有するメッシュであり、中心座標（lati，longi）は、
（lati，longi）＝（latc＋（yy－１）×Δlat，longc＋（xx－１）×Δlong）
と表され、ここで、Δlatは、前記少なくとも１つのメッシュの緯度方向幅であり、Δlongは、前記少なくとも１つのメッシュの経度方向幅であり、１≦yy≦Ｍであり、１≦xx≦Ｎである、ことと
によって、前記少なくとも１つのメッシュを決定することを含む、ことと、
前記少なくとも１つのメッシュのうちの１つのメッシュに属する少なくとも１つのデータに基づいて、前記１つのメッシュ内の所与の事項の統計量を示すメッシュデータを生成することと、
前記メッシュデータを前記データベース部に格納することと
を含む処理を実行することを前記プロセッサに行わせる、プログラム。

【請求項15】

前記処理は、
格納された複数のメッシュデータの値を集計することにより、所望の領域のメッシュ統計を提供すること
をさらに含む、請求項１４に記載のプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、位置情報に関連付けられたデータを利用してメッシュ統計を生成するためのシステム、方法、およびプログラムに関する。

【背景技術】

【0002】

我が国の領域にメッシュを割り当て、統計調査で収集された既存の統計データを各メッシュに割り当てることにより生成される地域メッシュ統計が知られている。この地域メッシュ統計を、世界の領域に拡張する「世界メッシュ統計」を生成する試みがなされている（非特許文献１、非特許文献２を参照）。しかしながら、世界メッシュ統計を効率的に生成する手法は見出されていない。

【先行技術文献】

【非特許文献】

【0003】

【文献】佐藤彰洋、椿広計、「データ利活用基盤とグローバル観光サービスデザイン」、第9回横幹連合コンファレンス、2018.10.6-7、電気通信大学

【文献】世界メッシュ研究所、“メッシュ統計の利活用と世界メッシュ統計への拡張”、インターネット＜URL: http://www.fttsus.jp/worldgrids/ja/document-ja/

【発明の概要】

【発明が解決しようとする課題】

【0004】

本発明は、メッシュ統計を効率的に生成することが可能なシステム、方法、およびプログラムを提供することを目的とする。

【課題を解決するための手段】

【0005】

本発明の一実施形態において、位置情報に関連付けられた少なくとも１つのデータを利用してメッシュ統計を生成するためのシステムが提供され、前記システムは、前記少なくとも１つのデータを受信する受信手段と、前記少なくとも１つのデータが属する少なくとも１つのメッシュを決定する決定手段であって、前記少なくとも１つのメッシュの各メッシュに少なくとも１つのデータが属している、決定手段と、前記少なくとも１つのメッシュのうちの１つのメッシュに属する少なくとも１つのデータに基づいて、前記１つのメッシュ内の所与の事項の統計量を示すメッシュデータを生成することによって、メッシュデータを生成する生成手段と、前記メッシュデータを格納する格納手段とを備える。

【0006】

一実施形態において、前記生成手段は、前記メッシュに属する前記データの統計量を算出することによって、前記メッシュデータを生成するように構成されている。

【0007】

一実施形態において、前記メッシュデータは、前記メッシュ内の所与の事項の最大値、最小値、中央値を含む。

【0008】

一実施形態において、前記生成手段は、前記複数のメッシュに対して並列的に前記メッシュデータを生成するように構成されている。

【0009】

一実施形態において、前記受信手段は、位置情報に関連付けられた第１のデータと位置情報に関連付けられた第２のデータとを受信し、前記生成手段は、（１）前記メッシュに属する第１のデータに基づいて第１のメッシュデータを生成することと、（２）前記メッシュに属する第２のデータに基づいて第２のメッシュデータを生成することと、（３）前記第１のメッシュデータと前記第２のメッシュデータとに演算処理を施すことにより複合メッシュデータを生成することとを行うように構成されている。

【0010】

一実施形態において、前記格納手段に格納された複数のメッシュデータの値を集計することにより、所望の領域のメッシュ統計を提供する集計手段をさらに備える。

【0011】

一実施形態において、前記集計手段は、提供すべきメッシュ統計の指示を受信することと、前記指示に基づいて、前記格納手段に格納されたメッシュデータのうち集計すべきメッシュデータを選択することと、前記選択されたメッシュデータの値を集計することとを行うように構成されている。

【0012】

一実施形態において、前記集計手段は、所望の領域の形状に対する、前記メッシュデータが関連付けられたメッシュの貢献度に基づいて、前記メッシュデータの値を集計することにより、前記所望の領域のメッシュ統計を提供するように構成されている。

【0013】

一実施形態において、前記集計手段は、前記所望の領域の形状と、前記メッシュデータが関連付けられたメッシュとの包含関係に基づいて、前記所望の領域の形状に対する前記メッシュの貢献度を決定する。

【0014】

一実施形態において、前記集計手段は、第１の事項の統計量を示す第１のメッシュデータと、第２の事項の統計量を示す第２のメッシュデータとに演算処理を施すことによって、第３の事項の統計量を示す第３のメッシュデータに基づくメッシュ統計を提供するように構成されている。

【0015】

一実施形態において、前記位置情報に関連付けられたデータは、衛星センサから取得された衛星画像データである。

【0016】

一実施形態において、前記衛星画像は、地表形状の標高を表す標高画像、夜間光強度を表す夜間光画像、雲密度を表す雲画像、輝度温度を表す輝度温度画像、人工構造物を特定する人工構造物特定画像のうちの少なくとも１つを含む。

【0017】

本発明の一実施形態において、位置情報に関連付けられた少なくとも１つのデータを利用してメッシュ統計を生成するための方法が提供され、前記方法は、位置情報に関連付けられた少なくとも１つのデータを受信することと、複数のメッシュのうち、前記少なくとも１つのデータが属する少なくとも１つのメッシュを決定することであって、前記複数のメッシュは、位置情報に関連付けられており、前記少なくとも１つのメッシュの各メッシュに少なくとも１つのデータが属している、ことと、前記少なくとも１つのメッシュのうちの１つのメッシュに属する少なくとも１つのデータに基づいて、前記１つのメッシュ内の所与の事項の統計量を示すメッシュデータを生成することと、前記メッシュデータを格納することとを含む。

【0018】

本発明の一実施形態において、位置情報に関連付けられた少なくとも１つのデータを利用してメッシュ統計を生成するためのプログラムが提供され、前記プログラムは、システムにおいて実行され、前記システムは、プロセッサとデータベース部とを備え、前記プログラムは、位置情報に関連付けられた少なくとも１つのデータを受信することと、複数のメッシュのうち、前記少なくとも１つのデータが属する少なくとも１つのメッシュを決定することであって、前記複数のメッシュは、位置情報に関連付けられており、前記少なくとも１つのメッシュの各メッシュに少なくとも１つのデータが属している、ことと、前記少なくとも１つのメッシュのうちの１つのメッシュに属する少なくとも１つのデータに基づいて、前記１つのメッシュ内の所与の事項の統計量を示すメッシュデータを生成することと、前記メッシュデータを前記データベース部に格納することとを含む処理を実行することを前記プロセッサに行わせる。

【0019】

本発明の一実施形態において、少なくとも１つのメッシュデータを利用するためのシステムが提供され、前記システムは、少なくとも１つのメッシュデータを受信する受信手段と、前記受信された少なくとも１つのメッシュデータの値を集計することにより、所望の領域のメッシュ統計を提供する集計手段であって、前記集計手段は、前記所望の領域の形状に対する、前記少なくとも１つのメッシュデータが関連付けられたメッシュの貢献度に基づいて、前記少なくとも１つのメッシュデータの値を集計することにより、前記所望の領域のメッシュ統計を提供するように構成されている、集計手段とを備える。

【0020】

本発明の一実施形態において、少なくとも１つのメッシュデータを利用するための方法が提供され、前記方法は、少なくとも１つのメッシュデータを受信することと、前記受信された少なくとも１つのメッシュデータの値を集計することにより、所望の領域のメッシュ統計を提供することであって、前記提供することは、前記所望の領域の形状に対する、前記少なくとも１つのメッシュデータが関連付けられたメッシュの貢献度に基づいて、前記少なくとも１つのメッシュデータの値を集計することにより、前記所望の領域のメッシュ統計を提供することを含む、こととを含む。

【0021】

本発明の一実施形態において、複数のメッシュデータを利用するためのプログラムが提供され、前記プログラムは、システムにおいて実行され、前記システムは、プロセッサを備え、前記プログラムは、少なくとも１つのメッシュデータを受信することと、前記受信された少なくとも１つのメッシュデータの値を集計することにより、所望の領域のメッシュ統計を提供することであって、前記提供することは、前記所望の領域の形状に対する、前記少なくとも１つのメッシュデータが関連付けられたメッシュの貢献度に基づいて、前記少なくとも１つのメッシュデータの値を集計することにより、前記所望の領域のメッシュ統計を提供することを含む、こととを含む処理を実行することを前記プロセッサに行わせる。

【発明の効果】

【0022】

本発明によれば、メッシュ統計を効率的に生成することが可能なシステム、方法、およびプログラムを提供することが可能である。これにより、世界メッシュ統計を効率的に生成することが可能である。

【図面の簡単な説明】

【0023】

【図1】本発明の位置情報に関連付けられたデータを利用してメッシュ統計を生成するためのシステム１０００の構成の一例を示す図。

【図2】プロセッサ１２０の構成の一例を示す図。

【図3】メッシュの種類を説明する図。

【図4】プロセッサ１２０の代替実施形態であるプロセッサ１２０’の構成の一例を示す図。

【図5】メッシュデータを格納する第２のデータベース部４００のデータ構成の一例を示す図。

【図6】本発明の位置情報に関連付けられたデータを利用してメッシュ統計を生成するためのシステム１０００による処理６００の一例を示すフローチャート。

【図7】ステップＳ６０３において、生成手段１２２が、位置情報に関連付けられた複数のデータ（第１のデータおよび第２のデータ）に基づいて複合メッシュデータを生成する処理の一例を示すフローチャート。

【図8】決定手段１２１によって行われる処理８００の一例を示すフローチャート。

【図9】集計手段２２１によって行われる処理９００の一例を示すフローチャート。

【図10】ポリゴンＶとメッシュｗとの包含関係を説明する図。

【図11】集計手段２２１によって行われる処理１１００の一例を示すフローチャート。

【発明を実施するための形態】

【0024】

１．定義
本明細書において「メッシュ」とは、位置情報に関連付けられた領域のことをいう。位置情報は、例えば、緯度および経度である。メッシュは、複数のメッシュにより空間（または平面）を充填することができる外形形状を有する領域である。

【0025】

本明細書において「メッシュデータ」とは、１つのメッシュ内の所与の事項の統計量を示すデータのことを言う。所与の事項は、統計量として表すことが可能な任意の事項であり得、示強性量であってもよいし、示量性量であってもよい。

【0026】

本明細書において「メッシュ統計」とは、複数のメッシュデータに基づく統計量のことを言う。

【0027】

一例において、「メッシュ」が矩形領域であるとき、或る地域に複数のメッシュが適用されると、その地域が矩形の網の目で区分される。そのときの１つのメッシュ内の統計量を示すデータが「メッシュデータ」であり、複数のメッシュの複数のメッシュデータに基づく統計量が「メッシュ統計」である。例えば、１つのメッシュ内の人口の統計量（例えば、平均値、最大値、最小値、中央値等）が人口に関する「メッシュデータ」となり、或る地域内の複数のメッシュについて、人口に関するメッシュデータの値を集計した統計量が、その地域の人口に関する「メッシュ統計」となる。

【0028】

２．位置情報に関連付けられたデータを利用したメッシュ統計の生成
本発明の発明者は、メッシュ統計を効率的に生成する手法を見出した。それは、位置情報に関連付けられたデータを利用して、メッシュ統計を生成する手法である。位置情報に関連付けられたデータを利用してメッシュ統計を生成する手法は、既存の統計データを複数のメッシュのそれぞれに対応付けて割り当てる従来の手法とは対照的に、位置情報に関連付けられた複数のデータが属するそれぞれのメッシュを決定し、決定されたメッシュのそれぞれから統計データを生成することができる。すなわち、メッシュと関連付けられた統計データを生成することができるため、統計データとメッシュとを関連付ける手間が省かれ、メッシュ統計を効率的に生成することができるのである。

【0029】

位置情報に関連付けられたデータは、例えば、衛星画像データであり得る。衛星画像データは、衛星から取得された画像（衛星画像）に含まれるデータである。衛星画像データは、衛星センサから取得される。現在、世界中の衛星画像が利用可能であるため、世界中の衛星画像を利用することにより、世界メッシュ統計を効率的に生成することができる。

【0030】

衛星画像には、その目的に応じた情報がデータとして含まれており、衛星画像に含まれる情報に対応したメッシュ統計を生成することができる。

【0031】

例えば、衛星画像の一例である標高画像（デジタル標高モデル）は、各画素値がその場所の地表形状の標高を表す。この標高画像を利用することによって、標高についてのメッシュ統計を生成することができる。

【0032】

例えば、衛星画像の一例である夜間光強度を表す夜間光画像は、各画素値がその場所の夜間光強度を表す。この夜間光画像を利用することによって、夜間光強度についてのメッシュ統計を生成することができる。

【0033】

例えば、衛星画像の一例である雲密度を表す雲画像は、各画素値がその場所の雲密度を表す。この雲画像を利用することによって、雲密度についてのメッシュ統計を生成することができる。

【0034】

位置情報に関連付けられたデータは、例えば、ポイントデータであり得る。ポイントデータは、点の位置情報と、その点に関連付けられた情報とを有するデータである。近年インターネット上で大量の位置情報付きデータポイントデータが収集され公開されており、これらのデータを利用することにより、世界メッシュ統計を効率的に生成することができる。

【0035】

ポイントデータは、例えば、会社の所在地（位置情報）およびその従業員数を示す情報（その点に関連付けられた情報）を有するデータ、宿泊施設の所在地（位置情報）およびその収容人数または宿泊した人数を示す情報（その点に関連付けられた情報）を有するデータ等多岐にわたる。ポイントデータには、その目的に応じた情報がデータとして含まれており、ポイントデータに含まれる情報に対応したメッシュ統計を生成することができる。

【0036】

位置情報に関連付けられたデータは、例えば、ポリゴンデータであり得る。ポリゴンデータは、複数の点によって表される領域の位置情報と、その領域に関連付けられた情報とを有するデータである。

【0037】

ポリゴンデータは、例えば、或る地区の位置情報（領域の位置情報）およびその地区内に居住している人数を示す情報（その領域に関連付けられた情報）を有するデータ、浸水想定区域図に含まれるデータ（すなわち、単位区域の位置情報（領域の位置情報）および単位区域毎の想定浸水深を示す情報（その点に関連付けられた情報）を有するデータ）等多岐にわたる。ポリゴンデータには、その目的に応じた情報がデータとして含まれており、ポリゴンデータに含まれる情報に対応したメッシュ統計を生成することができる。

【0038】

位置情報に関連付けられたデータは、例えば、データ源から得られた生データであってもよいし、加工されたデータであってもよい。例えば、位置情報に関連付けられたデータが衛星画像データである場合、衛星画像データは、衛星のセンサで観測された源記録であってもよいし、源記録を緯度経度直交系に写像し地図上に重ねて利用できるようにオルソ加工を施したデータであってもよいし、源記録から特徴量を抽出したデータであってもよい。

【0039】

さらに、生成された複数のメッシュ統計に基づいて、所望のメッシュ統計を得ることができる。例えば、所望の領域（例えば、国、都道府県、市区町村等の行政界区画等）に含まれるメッシュの統計データの値を集計することによって、その領域のメッシュ統計を生成することができる。生成されたメッシュ統計は、国の機関、地方自治体、民間の会社等によって種々の用途のために利用され得る。

【0040】

例えば、位置情報に関連付けられたデータが衛星画像データである場合、衛星画像の一例である標高画像から生成された複数のメッシュ統計を所望の領域について集計することにより、標高についてのその領域のメッシュ統計（例えば、平均標高、標高分布、単位面積当たりの標高等）を生成することができる。このようなメッシュ統計は、例えば、国土計画における基本フレーム、グランドデザインに活用されることができる。

【0041】

例えば、衛星画像の一例である夜間光画像から生成された複数のメッシュ統計を所望の領域について集計することにより、夜間光強度についてのその領域のメッシュ統計（例えば、平均夜間光強度、夜間光強度分布、単位面積当たりの夜間光強度等）を生成することができる。このようなメッシュ統計は、人口把握、特に夜間人口把握のために活用されることができる。

【0042】

例えば、衛星画像の一例である雲画像から生成された複数のメッシュ統計を所望の領域について集計することにより、雲密度についてのその領域のメッシュ統計（例えば、平均雲密度、雲密度分布、単位面積当たりの雲密度等）を生成することができる。さらに、略リアルタイムの頻度で更新される雲画像を利用すると、略リアルタイムの頻度で更新可能な時系列のメッシュ統計も生成することができる。このようなメッシュ統計は、気象状況の把握、気象予報、ひいては災害予測のために活用されることができる。

【0043】

例えば、位置情報に関連付けられたデータがポイントデータである場合、会社の所在地およびその従業員数を示す情報を有する複数のデータを生成された複数のメッシュ統計を所望の領域について集計することにより、その領域内で働いている人数についてのメッシュ統計（例えば、平均人数、人数分布等）を生成することができる。このようなメッシュ統計は、例えば、労働環境調査等に活用されることができる。

【0044】

例えば、位置情報に関連付けられたデータがポイントデータである場合、宿泊施設の所在地およびその収容人数を示す情報を有する複数のデータから生成された複数のメッシュ統計を所望の領域について集計することにより、その領域内の宿泊収容人数についてのメッシュ統計（例えば、平均宿泊可能人数、宿泊可能人数分布等）を生成することができる。このようなメッシュ統計は、例えば、宿泊施設分布調査等に活用されることができる。

【0045】

例えば、位置情報に関連付けられたデータがポリゴンデータである場合、浸水想定区域図に含まれるデータから生成された複数のメッシュ統計を所望の領域について集計することにより、想定浸水深についてのその領域のメッシュ統計（例えば、平均浸水深、浸水深分布等）を生成することができる。このようなメッシュ統計は、例えば、危険区域の把握、災害予測のために活用されることができる。

【0046】

このように位置情報に関連付けられたデータを利用してメッシュ統計を生成することにより、既存の統計データが存在しない地域であっても、衛星画像から抽出可能な統計量に関するものにはなるが、メッシュ統計を生成することが可能になる。

【0047】

さらに、メッシュ統計は、１つのメッシュ統計と他のメッシュ統計とを組み合わせることができるため、既存のメッシュ統計から、新たな観点でのメッシュ統計を導出することができる。これにより、種々の観点のメッシュ統計を無限に創出することができるようになる。

【0048】

例えば、ユーザが、旅行客が自然派であるか都会派であるかのプレファレンスを知りたいとき、夜間光画像から生成された複数のメッシュ統計と、宿泊施設の位置情報と宿泊した人数を示す複数のポイントデータから生成された複数のメッシュ統計とを組み合わせることにより、明るい場所に宿泊した人数および暗い場所に宿泊した人数を示すメッシュ統計を導出することができる。これにより、明るい場所に宿泊した人数（例えば、街中のホテルに宿泊した都会派の旅行者の人数）と暗い場所に宿泊した人数（例えば、山中のロッジに宿泊した自然派の旅行者の人数）とを比較することができ、旅行者のプレファレンスを知ることができるようになる。これは、例えば、マーケティング戦略立案等のために利用され得る。

【0049】

例えば、ユーザが、旅行客が山派であるか海派であるかのプレファレンスを知りたいとき、地表の物理状態を示す土地被覆画像から生成された複数のメッシュ統計と、宿泊施設の位置情報と宿泊した人数を示す複数のポイントデータから生成された複数のメッシュ統計とを組み合わせることにより、周囲に森林が多い場所に宿泊した人数および海が近い場所に宿泊した人数を示すメッシュ統計を導出することができる。これにより、周囲に森林が多い場所に宿泊した人数（例えば、山中のロッジに宿泊した山派の旅行者の人数）と海が近い場所に宿泊した人数（例えば、海沿いの民宿に宿泊した海派の旅行者の人数）とを比較することができ、旅行者のプレファレンスを知ることができるようになる。これは、例えば、マーケティング戦略立案等のために利用され得る。

【0050】

例えば、ユーザが、標高が高い場所に住んでいる人数を知りたいとき、夜間光画像から生成された複数のメッシュ統計と、標高画像から生成された複数のメッシュ統計とを組み合わせることにより、標高毎の人口を示すメッシュ統計を導出することができる。これにより、例えば１０００ｍ以上の高地に住んでいる人数、１５００ｍ以上の高地に住んでいる人数等を知ることができるようになる。これは、例えば、居住環境調査等のために利用され得る。

【0051】

位置情報に関連付けられたデータを利用してメッシュ統計を生成することは、例えば、以下に説明する位置情報に関連付けられたデータを利用してメッシュ統計を生成するためのシステムによって実現され得る。

【0052】

以下、図面を参照しながら、本発明の実施の形態を説明する。

【0053】

３．位置情報に関連付けられたデータを利用してメッシュ統計を生成するためのシステムの構成
図１は、本発明の位置情報に関連付けられたデータを利用してメッシュ統計を生成するためのシステム１０００の構成の一例を示す。

【0054】

システム１０００は、第１のコンピュータ装置１００と、第２のコンピュータ装置２００と、第１のコンピュータ装置１００に接続された第１のデータベース部３００と、第１のコンピュータ装置１００および第２のコンピュータ装置２００に接続された第２のデータベース部４００を備える。第１のコンピュータ装置１００は、位置情報に関連付けられたデータから複数のメッシュデータを生成するための装置であり得、第２のコンピュータ装置２００は、複数のメッシュデータを集計することによりメッシュ統計を生成するための装置であり得る。

【0055】

第１のコンピュータ装置１００は、インターフェース部１１０と、プロセッサ１２０と、メモリ１３０とを備える。

【0056】

インターフェース部１１０は、第１のコンピュータ装置１００の外部とデータのやり取りを行う。第１のコンピュータ装置１００のプロセッサ１２０は、インターフェース部１１０を介して、第１のコンピュータ装置１００の外部からデータを受信することが可能であり、第１のコンピュータ装置１００の外部にデータを送信することが可能である。インターフェース部１１０は、任意の形式でデータのやり取りを行うことができる。

【0057】

インターフェース部１１０は、第１のコンピュータ装置１００の外部からデータを受信することが可能であるように構成されている受信手段を備え得る。受信手段が、どのような態様でデータを受信するかは問わない。受信手段は、例えば、第１のコンピュータ装置１００の外部からネットワークを介してデータを受信してもよいし、第１のコンピュータ装置１００に接続された記憶媒体（例えば、ＵＳＢメモリ、光ディスク等）または第１のデータベース部３００または第２のデータベース部４００からデータを受信してもよい。ネットワークを介してデータを受信する場合は、ネットワークの種類は問わない。受信手段は、例えば、Ｗｉ－ｆｉ等の無線ＬＡＮを利用してデータを受信してもよいし、インターネットを介してデータを受信してもよい。

【0058】

インターフェース部１１０の受信手段は、例えば、位置情報に関連付けられたデータを受信する。位置情報に関連付けられたデータは、例えば、第１のコンピュータ装置１００に接続された第１のデータベース部３００に格納され得、受信手段は、第１のデータベース３００から位置情報に関連付けられたデータを受信してもよいし、システム１０００の外部から位置情報に関連付けられたデータを受信してもよい。位置情報に関連付けられたデータは、例えば、衛星画像データであってもよいし、ポイントデータであってもよいし、ポリゴンデータであってもよい。例えば、衛星画像には衛星センサからの複数の衛星画像データが含まれているため、受信手段は、衛星画像を受信することを介して、複数の衛星画像データを受信するようにしてもよい。衛星画像は、例えば、地表形状の標高を表す標高画像、夜間光強度を表す夜間光画像、雲密度を表す雲画像、輝度温度を表す輝度温度画像、人工構造物を特定する人工構造物特定画像のうちの少なくとも１つを含むが、これらに限定されない。例えば、標高画像は、例えば、ＪＡＸＡ（宇宙航空研究開発機構）によって公開されており、このモデルは、陸域観測技術衛星「だいち（ＡＬＯＳ）」搭載のパンクロマチック立体視センサ（ＰＲＩＳＭ）で観測した衛星画像に基づいている。例えば、標高画像は、経済産業省および米国航空宇宙局（ＮＡＳＡ）によっても提供されており、このモデルは、資源探査用将来型センサＡＳＴＥＲ（Ａｄｖａｎｃｅｄ Ｓｐａｃｅｂｏｒｎｅ Ｔｈｅｒｍａｌ Ｅｍｉｓｓｉｏｎ ａｎｄ Ｒｅｆｌｅｃｔｉｏｎ Ｒａｄｉｏｍｅｔｅｒ）「ＡＳＴＥＲ」を用いて測定された数値地形データ（ＡＳＴＥＲ全球３次元地形データ）に基づいている。

【0059】

衛星画像は、例えば、各画素値がその場所の特性を表し得る。例えば、地表形状の標高を表す標高画像は、各画素値がその場所の地表形状の標高を表す。例えば、夜間光強度を表す夜間光画像は、各画素値がその場所の夜間光強度を表す。例えば、雲密度を表す雲画像は、各画素値がその場所の雲密度を表す。輝度温度を表す輝度温度画像は、各画素値がその場所の輝度温度を表す。例えば、人工構造物を特定する人工構造物特定画像は、各画素値がその場所の人工構造物の有無または種類を表す。

【0060】

ポイントデータは、点の位置情報と、その点に関連付けられたデータとを有する任意のデータである。ポリゴンデータは、複数の点によって表される領域の位置情報と、その領域に関連付けられた情報とを有する任意のデータである。

【0061】

インターフェース部１１０の受信手段が受信したデータは、後続の処理のために、プロセッサ１２０に渡される。

【0062】

インターフェース部１１０は、第１のコンピュータ装置１００の外部にデータを出力することが可能であるように構成されている出力手段を備え得る。出力手段が、どのような態様で第１のコンピュータ装置１００からデータを出力するかは問わない。例えば、出力手段が表示画面である場合、表示画面にデータを出力するようにしてもよい。あるいは、出力手段がデータ書き込み装置である場合、第１のコンピュータ装置１００に接続された記憶媒体または第２のデータベース部４００にデータを書き込むことによってデータを出力するようにしてもよい。あるいは、出力手段が送信器である場合、送信器がネットワークを介して第１のコンピュータ装置１００の外部にデータを送信することにより出力してもよい。この場合、ネットワークの種類は問わない。例えば、送信器は、インターネットを介してデータを送信してもよいし、ＬＡＮを介してデータを送信してもよい。例えば、出力手段は、データの出力先のハードウェアまたはソフトウェアによって取り扱い可能な形式に変換して、または、データの出力先のハードウェアまたはソフトウェアによって取り扱い可能な応答速度に調整してデータを出力するようにしてもよい。

【0063】

インターフェース部１１０の出力手段は、例えば、生成されたメッシュデータを出力する。出力されたメッシュデータは、例えば、第２のデータベース部４００に送信され、第２のデータベース部４００に格納され得る。

【0064】

プロセッサ１２０は、第１のコンピュータ装置１００全体の動作を制御する。プロセッサ１２０は、メモリ１３０に格納されているプログラムを読み出し、そのプログラムを実行する。これにより、第１のコンピュータ装置１００を所望のステップを実行する装置として機能させることが可能である。プロセッサ１２０は、単一のプロセッサによって実装されてもよいし、複数のプロセッサによって実装されてもよい。プロセッサ１２０によって処理されたデータは、出力のために、インターフェース部１１０に渡される。

【0065】

メモリ１３０には、第１のコンピュータ装置１００における処理を実行するためのプログラムやそのプログラムの実行に必要とされるデータ等が格納されている。メモリ１３０には、例えば、衛星画像を利用してメッシュ統計を生成するためのプログラム（例えば、後述する図６、図７、図８に示される処理を実現するプログラム）が格納されている。メモリ１３０には、任意の機能を実装するアプリケーションが格納されていてもよい。ここで、プログラムをどのようにしてメモリ１３０に格納するかは問わない。例えば、プログラムは、メモリ１３０にプリインストールされていてもよい。あるいは、プログラムは、ネットワークを経由してダウンロードされることによってメモリ１３０にインストールされるようにしてもよい。メモリ１３０は、任意の記憶手段によって実装され得る。

【0066】

第２のコンピュータ装置２００は、インターフェース部２１０と、プロセッサ２２０と、メモリ２３０とを備える。

【0067】

インターフェース部２１０は、第２のコンピュータ装置２００の外部とデータのやり取りを行う。第２のコンピュータ装置２００のプロセッサ２２０は、インターフェース部２１０を介して、第２のコンピュータ装置２００の外部からデータを受信することが可能であり、第２のコンピュータ装置２００の外部にデータを送信することが可能である。インターフェース部２１０は、任意の形式でデータのやり取りを行うことができる。

【0068】

インターフェース部２１０は、第２のコンピュータ装置２００の外部からデータを受信することが可能であるように構成されている受信手段を備え得る。受信手段が、どのような態様でデータを受信するかは問わない。受信手段は、例えば、第２のコンピュータ装置２００の外部からネットワークを介してデータを受信してもよいし、第２のコンピュータ装置２００に接続された記憶媒体（例えば、ＵＳＢメモリ、光ディスク等）または第２のデータベース部４００からデータを受信してもよい。ネットワークを介してデータを受信する場合は、ネットワークの種類は問わない。受信手段は、例えば、Ｗｉ－ｆｉ等の無線ＬＡＮを利用してデータを受信してもよいし、インターネットを介してデータを受信してもよい。

【0069】

インターフェース部２１０の受信手段は、例えば、複数のメッシュデータを受信する。複数のメッシュデータは、例えば、第２のデータベース部４００に格納され得、受信手段は、第２のデータベース部４００から複数のメッシュデータを受信してもよいし、システム１０００の外部から複数のメッシュデータを受信してもよい。

【0070】

受信手段は、例えば、システム１０００の外部から、提供すべきメッシュ統計の指示を受信することができる。提供すべきメッシュ統計の指示は、例えば、ユーザから送信され得、ユーザが所望するメッシュ統計の情報を含み得る。

【0071】

インターフェース部２１０の受信手段が受信したデータは、後続の処理のために、プロセッサ２２０に渡される。

【0072】

インターフェース部２１０は、第２のコンピュータ装置２００の外部にデータを出力することが可能であるように構成されている出力手段を備え得る。出力手段が、どのような態様で第２のコンピュータ装置２００からデータを出力するかは問わない。例えば、出力手段が表示画面である場合、表示画面にデータを出力するようにしてもよい。あるいは、出力手段がデータ書き込み装置である場合、第２のコンピュータ装置２００に接続された記憶媒体または第２のデータベース部４００にデータを書き込むことによってデータを出力するようにしてもよい。あるいは、出力手段が送信器である場合、送信器がネットワークを介して第２のコンピュータ装置２００の外部にデータを送信することにより出力してもよい。この場合、ネットワークの種類は問わない。例えば、送信器は、インターネットを介してデータを送信してもよいし、ＬＡＮを介してデータを送信してもよい。例えば、出力手段は、データの出力先のハードウェアまたはソフトウェアによって取り扱い可能な形式に変換して、または、データの出力先のハードウェアまたはソフトウェアによって取り扱い可能な応答速度に調整してデータを出力するようにしてもよい。インターフェース部２１０の出力手段は、例えば、生成されたメッシュ統計を出力する。出力されたメッシュ統計は、例えば、第２のデータベース部４００に送信され、第２のデータベース部４００に格納され得る。

【0073】

プロセッサ２２０は、第２のコンピュータ装置２００全体の動作を制御する。プロセッサ２２０は、メモリ２３０に格納されているプログラムを読み出し、そのプログラムを実行する。これにより、第２のコンピュータ装置２００を所望のステップを実行する装置として機能させることが可能である。プロセッサ２２０は、単一のプロセッサによって実装されてもよいし、複数のプロセッサによって実装されてもよい。プロセッサ２２０によって処理されたデータは、出力のために、インターフェース部２１０に渡される。

【0074】

メモリ２３０には、第２のコンピュータ装置２００における処理を実行するためのプログラムやそのプログラムの実行に必要とされるデータ等が格納されている。メモリ２３０には、例えば、衛星画像を利用してメッシュ統計を生成するためのプログラム（例えば、後述する図９に示される処理を実現するプログラム）が格納されている。メモリ２３０には、任意の機能を実装するアプリケーションが格納されていてもよい。ここで、プログラムをどのようにしてメモリ２３０に格納するかは問わない。例えば、プログラムは、メモリ２３０にプリインストールされていてもよい。あるいは、プログラムは、ネットワークを経由してダウンロードされることによってメモリ２３０にインストールされるようにしてもよい。メモリ２３０は、任意の記憶手段によって実装され得る。

【0075】

第１のコンピュータ装置１００に接続されている第１のデータベース部３００には、例えば、衛星画像が格納され得る。第１のデータベース部３００には、例えば、第２のコンピュータ装置２００によって出力されたデータ（例えば、生成されたメッシュデータ、生成されたメッシュ統計）が格納され得る。生成されたメッシュデータは、そのメッシュを識別する情報（例えば、地域メッシュ統計のために用いられている「地域メッシュコード」、もしくは、それを世界規模に拡張した「世界メッシュコード」、および／または、位置情報（例えば、メッシュの中央の座標、もしくは、メッシュの北西端の座標および南東端の座標）等）と関連付けられて格納され得る。「地域メッシュコード」は、日本全土に複数のメッシュを割り当て、各メッシュを識別できるようにするために付されたコードである。「世界メッシュコード」は、「地域メッシュコード」を世界規模に拡張したものであり、世界全土に複数のメッシュを割り当て、各メッシュを識別できるようにするために付されたコードである。緯度および経度の位置情報により、一意に「地域メッシュコード」または「世界メッシュコード」を特定することができる。

【0076】

第１のコンピュータ装置１００に接続されている第１のデータベース部３００には、例えば、位置情報に関連付けられたデータ（例えば、衛星画像データ、ポイントデータ、ポリゴンデータ）が格納され得る。

【0077】

第２のデータベース部４００には、例えば、第１のコンピュータ装置１００によって出力されたデータ（例えば、生成されたメッシュデータ）が格納され得る。生成されたメッシュデータは、そのメッシュを識別する情報（例えば、地域メッシュ統計のために用いられている「地域メッシュコード」、もしくは、それを世界規模に拡張した「世界メッシュコード」、および／または、位置情報（例えば、メッシュの中央の座標、もしくは、メッシュの北西端の座標および南東端の座標）等）と関連付けられて格納され得る。生成されたメッシュデータは、例えば、複数の統計量（最小値、最大値、中央値、平均値）を有することができる。第２のデータベース４００には、例えば、第２のコンピュータ装置２００によって出力されたデータ（例えば、生成されたメッシュ統計）が格納され得る。

【0078】

第２のデータベース部４００には、例えば、種々の領域の形状データが格納され得る。種々の領域は、例えば、国、都道府県、市区町村等の行政界区画等であるが、これらに限定されない。形状データは、例えば、種々の領域の形状を多角形で近似したデータであってもよいし、曲線で近似したデータであってもよい。

【0079】

図１に示される例では、第１のコンピュータ装置１００、第２のコンピュータ装置２００、第１のデータベース部３００、第２のデータベース部４００が、それぞれ別個の構成要素として説明されているが、本発明はこれに限定されない。第１のコンピュータ装置１００、第２のコンピュータ装置２００、第１のデータベース部３００、第２のデータベース部４００のうちの少なくとも２つを同一の構成要素内に設けることも可能である。例えば、第１のコンピュータ装置１００および第２のコンピュータ装置２００を同一の構成要素とすることができる。このとき、第１のコンピュータ装置１００および第２のコンピュータ装置２００内のインターフェース部１１０、２１０、プロセッサ１２０、２２０、メモリ部１３０、２３０は、それぞれ統合されてもよいし、別個であってもよい。例えば、第１のデータベース部３００を第１のコンピュータ装置１００の内部に設けることができる。このとき、第１のデータベース部３００は、メモリ１３０を実装する記憶手段と同一の記憶手段によって実装されてもよいし、メモリ１３０を実装する記憶手段とは別の記憶手段によって実装されてもよい。いずれにせよ、第１のデータベース部３００は、第１のコンピュータ装置１００のための格納部として構成される。例えば、第２のデータベース部４００を第２のコンピュータ装置２００の内部に設けることができる。このとき、第２のデータベース部４００は、メモリ２３０を実装する記憶手段と同一の記憶手段によって実装されてもよいし、メモリ２３０を実装する記憶手段とは別の記憶手段によって実装されてもよい。いずれにせよ、第２のデータベース部４００は、第２のコンピュータ装置２００のための格納部として構成される。第１のデータベース部３００および第２のデータベース部４００の構成は、特定のハードウェア構成に限定されない。例えば、第１のデータベース部３００および第２のデータベース部４００の各々は、単一のハードウェア部品で構成されてもよいし、複数のハードウェア部品で構成されてもよい。例えば、第１のデータベース部３００および／または第２のデータベース部４００は、第１のコンピュータ装置１００および／または第２のコンピュータ装置２００の外付けハードディスク装置として構成されてもよいし、ネットワークを介して接続されるクラウド上のストレージとして構成されてもよい。

【0080】

好ましい実施形態では、第１のデータベース部３００と、第２のデータベース部４００とは異なる構成要素で実装され得る。第１のデータベース部３００に格納される位置情報関連付けられたデータは、比較的に取り扱いに注意を要するセンシティブなデータであり得る。位置情報に関連付けられたデータは、例えば、個人情報を含む場合があり、個人情報が漏洩するリスクがある。また、位置情報に関連付けられたデータは、機密情報を含む場合があり、悪意を持って利用されるリスクがある。例えば、衛星画像のうち、光源または熱源を示す衛星画像は、光源または熱源により、人が住んでいる場所を示す場合もあれば、原子力発電所等の核施設を示す場合もあり、軍事的に利用されるリスクがある。このようなセンシティブなデータへのアクセスを第１のコンピュータ装置１００に制限し、第１のコンピュータ装置１００によって生成された、取り扱いが容易なメッシュデータを第２のデータベース部４００に格納して公衆に利用可能とすることにより、このようなリスクを排除することができる。さらに、メッシュデータであれば、個人が特定されるリスクは低く、セキュリティが強化されるという利点がある。また、衛星画像等のデータ量が大きいデータであっても、メッシュデータ化されることにより、扱いが容易になり、処理速度が向上する。

【0081】

図２は、プロセッサ１２０の構成の一例を示す。

【0082】

プロセッサ１２０は、少なくとも、決定手段１２１と、生成手段１２２とを備える。

【0083】

決定手段１２１は、インターフェース部１１０によって受信された、位置情報に関連付けられた複数のデータが、それぞれ、複数のメッシュのうちどのメッシュに属するかを決定するように構成されている。位置情報に関連付けられた複数のデータは、例えば、少なくとも１つの衛星画像に含まれる衛星画像データであり、決定手段１２１は、少なくとも１つの衛星画像に複数のメッシュを割り当てることにより、衛星画像に含まれる複数の衛星画像データが属するそれぞれのメッシュを決定する。決定手段１２１は、例えば、図８を参照して後述する処理８００によって、受信された衛星画像に複数のメッシュを割り当てることにより、衛星画像に含まれる複数の衛星画像データが属するそれぞれのメッシュを決定するができる。

【0084】

位置情報に関連付けられた複数のデータは、例えば、複数のポイントデータであり、決定手段１２１は、複数のポイントデータが属するそれぞれのメッシュを決定する。決定手段１２１は、例えば、ポイントデータが示す位置情報（緯度および経度）から、対応する「地域メッシュコード」またはそれを世界規模に拡張した「世界メッシュコード」を計算することによって、ポイントデータが属するメッシュを決定する。

【0085】

位置情報に関連付けられた複数のデータは、例えば、複数のポリゴンデータであり、決定手段１２１は、複数のポリゴンデータが属するそれぞれのメッシュを決定する。決定手段１２１は、例えば、ポリゴンデータを形成する複数の点のそれぞれが示す位置情報（緯度および経度）から、対応する「地域メッシュコード」またはそれを世界規模に拡張した「世界メッシュコード」を計算することによって、ポリゴンデータの各点が属するそれぞれのメッシュを決定し、それぞれのメッシュ全体を、ポリゴンデータが属するメッシュであると決定する。

【0086】

メッシュは、複数のメッシュにより空間（または平面）を充填することができる外形形状を有する領域であり、例えば、メッシュは、矩形領域、または、多角形領域であり得る。複数のメッシュにより空間を充填することができる外形形状は、例えば、正三角形、正方形、長方形、平行四辺形、正六角形を含む。

【0087】

一実施形態において、メッシュは、例えば、矩形領域であり、このとき、メッシュの各辺は、経線または緯線に略平行である。メッシュには、複数の種類のメッシュが存在し得る。例えば、図３に示されるように、大きさが異なる「１次メッシュ」、「２次メッシュ」、「３次メッシュ（または基準メッシュ）」、「２分の１メッシュ」、「４分の１メッシュ」、「８分の１メッシュ」が存在する。

【0088】

「１次メッシュ」は、例えば、緯度方向幅４０分および経度方向幅１度を有する矩形領域である。「２次メッシュ」は、１次メッシュを緯度方向および経度方向に８等分してできる矩形領域であり、例えば、緯度方向幅５分および経度方向幅７．５分を有する矩形領域である。このとき、「１次メッシュ」は、１６個の「２次メッシュ」を包含することになり、１６個の「２次メッシュ」のメッシュデータを集計したメッシュ統計が、対応する「１次メッシュ」のメッシュデータとなり得る。

【0089】

「３次メッシュ」は、２次メッシュを緯度方向および経度方向に１０等分してできる矩形領域であり、例えば、緯度方向幅３０秒および経度方向幅４５秒を有する矩形領域である。このとき、「２次メッシュ」は、１００個の「３次メッシュ」を包含することになり、１００個の「３次メッシュ」のメッシュデータを集計したメッシュ統計が、対応する「２次メッシュ」のメッシュデータとなり得る。

【0090】

「２分の１メッシュ」は、３次メッシュを緯度方向および経度方向に２等分してできる矩形領域であり、例えば、緯度方向幅１５秒および経度方向幅２２．５秒を有する矩形領域である。このとき、「３次メッシュ」は、４個の「２分の１メッシュ」を包含することになり、４個の「２分の１メッシュ」のメッシュデータを集計したメッシュ統計が、対応する「３次メッシュ」のメッシュデータとなり得る。

【0091】

「４分の１メッシュ」は、２分の１メッシュを緯度方向および経度方向に２等分してできる矩形領域であり、例えば、緯度方向幅７．５秒および経度方向幅１１．２５秒を有する矩形領域である。このとき、「２分の１メッシュ」は、４個の「４分の１メッシュ」を包含することになり、４個の「４分の１メッシュ」のメッシュデータを集計したメッシュ統計が、対応する「２分の１メッシュ」のメッシュデータとなり得る。

【0092】

「８分の１メッシュ」は、４分の１メッシュを緯度方向および経度方向に２等分してできる矩形領域であり、例えば、緯度方向幅３．７５秒および経度方向幅５．６２５秒を有する矩形領域である。このとき、「４分の１メッシュ」は、４個の「８分の１メッシュ」を包含することになり、４個の「８分の１メッシュ」のメッシュデータを集計したメッシュ統計が、対応する「４分の１メッシュ」のメッシュデータとなり得る。

【0093】

メッシュは、例えば、地域メッシュ統計のために用いられている「地域メッシュコード」、または、それを世界規模に拡張した「世界メッシュコード」によって識別され得る。あるいは、例えば、メッシュの位置情報（例えば、メッシュの中央の座標、もしくは、メッシュの北西端の座標および南東端の座標）によって識別され得る。

【0094】

以下では、メッシュが矩形領域であるとして、説明する。

【0095】

例えば、位置情報に関連付けられたデータが、衛星画像データである場合、決定手段１２１が衛星画像データが属するメッシュを決定する。このとき、決定手段１２１は、衛星画像に任意の数のメッシュを割り当てることにより、衛星画像データが属するメッシュを決定する。メッシュの数は、例えば、衛星画像の緯度および経度の範囲と、緯度方向幅Δlatと、経度方向幅Δlongとによって決定される。例えば、衛星画像が、最大緯度maxlat、最小緯度minlat、最大経度maxlong、最小経度minlongを有する場合、メッシュの数は、Ｎ×Ｍによって算出され、ここで、Ｍは緯度方向のメッシュの数であり、Ｎは経度方向のメッシュの数であり、ＮおよびＭは以下の式によって算出される。
Ｍ＝_└(maxlat-minlat)/(Δlat)_┘
Ｎ＝_└(maxlong-minlong)/(Δlong)_┘

【0096】

緯度方向幅Δlatおよび経度方向幅Δlongは、生成しようとするメッシュに応じて変動し得る。例えば、１次メッシュを割り当てる場合、Δlatは、４０分であり、Δlongは、１度であり得る。例えば、２次メッシュを割り当てる場合、Δlatは、５分であり、Δlongは、７．５分であり得る。例えば、３次メッシュを割り当てる場合、Δlatは、３０秒であり、Δlongは、４５秒であり得る。例えば、２分の１メッシュを割り当てる場合、Δlatは、１５秒であり、Δlongは、２２．５秒であり得る。例えば、４分の１メッシュを割り当てる場合、Δlatは、７．５秒であり、Δlongは、１１．２５秒であり得る。例えば、８分の１メッシュを割り当てる場合、Δlatは、３．７５秒であり、Δlongは、５．６２５秒であり得る。

【0097】

生成手段１２２は、複数のメッシュのうちの各メッシュに対して、そのメッシュに属するデータに基づいてメッシュデータを生成するように構成されている。生成手段１２２は、例えば、複数のメッシュのうちの各メッシュに対して、メッシュ内に含まれるデータの統計量（例えば、最小値、最大値、中央値、平均値）を算出することによって、メッシュデータを生成する。例えば、位置情報に関連付けられたデータが、衛星画像データである場合、生成手段１２２は、衛星画像の画素値の統計量（例えば、最小値、最大値、中央値、平均値）を算出することによって、メッシュデータを生成する。

【0098】

例えば、対象のメッシュ内にｍ個の画素が含まれている場合、画素値の集合Ｏは、
Ｏ＝｛ｏ_１，ｏ_２，・・・ｏ_ｍ｝
で表される。ここで、ｏ_ｉは、画素ｉの画素値である。

【0099】

生成手段１２２は、画素値の集合Ｏについて、統計量（例えば、最小値、最大値、中央値、平均値）を算出することができる。例えば、衛星画像にマスク等が存在し、メッシュ内に画素値が存在しない場合には、統計量を算出しないようにしてもよい。

【0100】

生成手段１２２は、例えば、複数のメッシュの各々に対して直列的にメッシュデータを生成するようにしてもよい。好ましい実施形態では、複数のメッシュの各々に対して並列的にメッシュデータを生成し得る。並列的にメッシュデータを生成することにより、大きなデータ量を有するデータであっても、計算量および計算時間が大きくならないという利点があるからである。特に、衛星画像は多数のデータ点が１ファイルに存在することにより、衛星画像に対する処理は計算量が膨大になり扱いづらいという問題があるが、メッシュ毎に並列的にメッシュデータを生成することにより、多数のデータ点が１つのファイルに入った高分解能の衛星画像であっても計算量および計算時間が大きくならない。

【0101】

生成手段１２２は、例えば、複合メッシュデータを生成することが可能なように構成されることもできる。ここで、複合メッシュデータとは、複数のメッシュデータに基づいて生成されたメッシュデータのことである。複合メッシュデータは、基づいたメッシュデータが示す統計量とは異なる新たな統計量を示すことができるため、基づいたメッシュデータとは別の観点からの統計量分析を可能にする。例えば、生成手段１２２は、位置情報に関連付けられた複数の第１のデータに基づいて、第１の事項に関する第１のメッシュデータを生成することと、位置情報に関連付けられた複数の第２のデータに基づいて、第２の事項に関する第２のメッシュデータを生成することと、第１のメッシュデータと第２のメッシュデータとに演算処理を施すことにより、複合メッシュデータを生成することとを行うができる。

【0102】

例えば、生成手段１２２は、第１の衛星画像データおよび第２の衛星画像データ、または、第１のポイントデータおよび第２のポイントデータ、または、第１のポリゴンデータおよび第２のポリゴンデータに基づいて、複合メッシュデータを生成するようにしてもよいし、衛星画像データ、ポイントデータ、ポリゴンデータのうちの少なくとも２つに基づいて、複合メッシュデータを生成するようにしてもよい。例えば、生成手段１２２は、位置情報に関連付けられたデータに基づいて、第１の事項に関する第１のメッシュデータを生成することと、位置情報に関連付けられたデータ以外のソースから第２の事項に関する第２のメッシュデータを取得することと、第１のメッシュデータと第２のメッシュデータとに演算処理を施すことにより、複合メッシュデータを生成することも可能である。

【0103】

演算処理は、例えば、四則演算であり得る。四則演算のうち、和および差は、第１のメッシュデータと第２のメッシュデータとが同一単位系である場合に行うことができる。四則演算のうち、積および商は、第１のメッシュデータと第２のメッシュデータとが同一単位系であっても、異なる単位系であっても行うことができ、積または商の演算結果の単位系は、もとの２つの単位系の積または商となる。

【0104】

例えば、第１のメッシュデータが衛星画像データから生成された人工物占有面積に関するデータであり、第２のメッシュデータが衛星画像データから生成された森林面積に関するデータである場合、第１のメッシュデータと第２のメッシュデータとの和を取ることによって、人工物と森林とによって占有されている面積に関する複合メッシュデータを得ることができる。例えば、複数のメッシュのうちの１つのメッシュについて、第１の衛星画像の画素値の平均値を算出することによりそのメッシュ内の人工物占有面積が約５０，０００ｍ^２であると算出され、第２の衛星画像の画素値の平均値を算出することによりそのメッシュ内の森林面積が約３０，０００ｍ^２であると算出された場合、人工物と森林とによって占有されている面積に関するそのメッシュの複合メッシュデータは、約８０，０００ｍ^２となる。

【0105】

例えば、第１のメッシュデータが衛星画像データから生成された森林面積に関するデータであり、第２のメッシュデータがポリゴンデータから生成された針葉樹林に関するデータである場合、第１のメッシュデータと第２のメッシュデータとの差を取ることによって、針葉樹林ではない森林の面積に関する複合メッシュデータを得ることができる。例えば、複数のメッシュのうちの１つのメッシュについて、第１の衛星画像の画素値の平均値を算出することによりそのメッシュ内の森林面積が約３０，０００ｍ^２であると算出され、ポリゴンデータのデータの平均値を算出することによりそのメッシュ内の針葉樹林面積が約１０，０００ｍ^２であると算出された場合、針葉樹林ではない森林の面積に関するそのメッシュの複合メッシュデータは、約２０，０００ｍ^２となる。

【0106】

例えば、第１のメッシュデータがポイントデータから生成された宿泊施設の稼働率に関するデータであり、第２のメッシュデータがポイントデータから生成された宿泊施設の収容人数に関するデータである場合、第１のメッシュデータと第２のメッシュデータとの積を取ることによって、宿泊者数に関する複合メッシュデータを得ることができる。例えば、複数のメッシュのうちの１つのメッシュについて、第１のポイントデータのデータの平均値を算出することによりそのメッシュ内の宿泊施設の稼働率が約３５％であると算出され、第２のポイントデータのデータの平均値を算出することによりそのメッシュ内の宿泊施設の収容人数が約１０，０００人であると算出された場合、宿泊者数に関するそのメッシュの複合メッシュデータは、約３，５００人となる。

【0107】

例えば、第１のメッシュデータがポリゴンデータから生成された第２次産業の売上額に関するデータであり、第２のメッシュデータがポリゴンデータから生成された第２次産業就業者数に関するデータである場合、第１のメッシュデータと第２のメッシュデータとの商を取ることによって、就業者１人当たりの第２次産業売上額に関する複合メッシュデータを得ることができる。例えば、複数のメッシュのうちの１つのメッシュについて、第１のポリゴンデータのデータの平均値を算出することによりそのメッシュ内の第２次産業の売上額が約８億万円であると算出され、第２のポリゴンデータのデータの平均値を算出することによりそのメッシュ内の第２次産業就業者数が約１００人であると算出された場合、就業者１人当たりの第２次産業売上額に関するそのメッシュの複合メッシュデータは、約８００万円となる。

【0108】

生成手段１２２によって生成されたメッシュデータ（または複合メッシュデータ）は、インターフェース部１１０に渡され、インターフェース部１１０を介して装置１００の外部に出力される。出力されたデータは、例えば、第１のデータベース部３００に格納され得る。

【0109】

図４は、プロセッサ２２０の構成の一例を示す。プロセッサ２２０は、少なくとも、集計手段２２１を備える。

【0110】

集計手段２２１は、インターフェース部２１０によって受信された複数のメッシュデータの値を集計することにより、所望の領域のメッシュ統計を生成するように構成されている。インターフェース部２１０によって受信される複数のメッシュデータは、第１のコンピュータ装置１００によって生成され第２のデータベース部４００に格納された複数のメッシュデータであってもよいし、システム１０００の外部から送信された複数のメッシュデータであってもよい。例えば、集計手段２２１は、所望の領域に対する、複数のメッシュデータが関連付けられたそれぞれのメッシュの包含関係に基づいて、複数のメッシュデータの値を集計することができる。所望の領域は、例えば、国、都道府県、市区町村等の行政界区画であり得る。このような領域は、複雑な外形形状を有し得るが、集計手段２２１は、この領域を例えばポリゴン（多角形）で近似し、ポリゴンとメッシュとの包含関係から、複数のメッシュデータの値を集計し得る。所望の領域のポリゴンとメッシュとの包含関係に基づいて複数のメッシュデータの値を集計する処理は、図９を参照して後述する処理９００に説明される。

【0111】

集計手段２２１は、例えば、インターフェース部２１０によって受信された、提供すべきメッシュ統計の指示を受信するようにしてもよい。集計手段２２１は、受信した指示に基づいて、メッシュ統計を生成するように、複数のメッシュデータを集計することができる。集計手段２２１は、例えば、受信された指示に基づいて、第２のデータベース部４００に格納されている複数のメッシュデータを検索し、集計すべき複数のメッシュデータを選択し、選択された複数のメッシュデータの値を集計することができる。受信手段２２１は、例えば、受信された指示に基づいて、受信された複数のメッシュデータから集計すべき複数のメッシュデータを選択し、選択された複数のメッシュデータの値を集計することができる。これにより、集計手段２２１は、ユーザによって要求されたメッシュ統計を生成することができる。

【0112】

集計手段２２１が複数のメッシュデータの値を集計することによって生成されたメッシュ統計は、インターフェース部１１０に渡され、インターフェース部１１０を介して装置１００の外部に出力される。出力されたデータは、例えば、第１のデータベース部３００に格納され得る。

【0113】

図５は、複数のメッシュデータを格納する第２のデータベース部４００のデータ構成の一例を示す。

【0114】

第２のデータベース部４００は、例えば、同一の事項についてのメッシュデータをまとめて格納し得る。第２のデータベース部４００は、例えば、第１の事項（例えば、標高）についての第１のメッシュデータを格納する第１のデータベース５０１と、第２の事項（例えば、雲密度）についての第２のメッシュデータを格納する第２のデータベース５０２と、第３の事項（例えば、夜間光強度）についての第３のメッシュデータを格納する第３のデータベース５０３を備え得る。例えば、第１のメッシュデータは、標高画像から生成されたメッシュデータであり、第２のメッシュデータは、雲画像から生成されたメッシュデータであり、第３のメッシュデータは、夜間光画像から生成されたメッシュデータであり得る。

【0115】

メッシュデータは、例えば、最大値、最小値、中央値のうちの少なくとも１つを含み得、好ましくは、最大値、最小値、および中央値を含み得る。最大値、最小値、および中央値を含むことにより、生成されるメッシュ統計の幅が広がるからである。メッシュデータは、例えば、最大値、最小値、および中央値に加えて、平均値も含むようにしてもよい。メッシュデータは、例えば、最大値、最小値、および中央値に加えて、他の四分位点（例えば２５パーセンタイル（第一四分位点）、７５パーセンタイル（第三四分位点））を含むようにしてもよい。

【0116】

メッシュデータは、そのメッシュを識別する情報（例えば、地域メッシュコードもしくは世界メッシュコード、および／または、位置情報等）と関連付けられて格納され得る。本例では、メッシュの北西端の座標（緯度，経度）および南東端の座標（緯度，経度）によって表される位置情報と関連付けられて格納されている。座標の値は、正の値が北緯または東経を表し、負の値が南緯または西経を表し得る。

【0117】

第１のデータベース５０１、第２のデータベース５０２、および第３のデータベース５０３は、それぞれ同一のサイズのメッシュについてのメッシュデータを格納しているが、異なるサイズのメッシュについてのメッシュデータを格納するデータベースを備えるようにしてもよい。

【0118】

上述した図２、図４に示される例では、プロセッサ１２０、２２０の各構成要素が同一のプロセッサ１２０またはプロセッサ２２０内に設けられているが、本発明はこれに限定されない。プロセッサ１２０またはプロセッサ２２０の各構成要素が、複数のプロセッサ部に分散される構成も本発明の範囲内である。このとき、複数のプロセッサ部は、同一のハードウェア部品内に位置してもよいし、近傍または遠隔の別個のハードウェア部品内に位置してもよい。

【0119】

上述した図１に示される例では、第１のコンピュータ装置１００の各構成要素が第１のコンピュータ装置１００内に設けられており、第２のコンピュータ装置２００の各構成要素が第２のコンピュータ装置２００内に設けられているが、本発明はこれに限定されない。第１のコンピュータ装置１００の各構成要素のいずれかが第１のコンピュータ装置１００の外部に設けられることも可能であり、第２のコンピュータ装置２００の各構成要素のいずれかが第２のコンピュータ装置２００の外部に設けられることも可能である。例えば、プロセッサ１２０、メモリ１３０のそれぞれが別々のハードウェア部品で構成されている場合には、各ハードウェア部品が任意のネットワークを介して接続されてもよい。このとき、ネットワークの種類は問わない。各ハードウェア部品は、例えば、ＬＡＮを介して接続されてもよいし、無線接続されてもよいし、有線接続されてもよい。第１のコンピュータ装置１００および第２のコンピュータ装置２００は、特定のハードウェア構成には限定されない。例えば、プロセッサ１２０、２２０をデジタル回路ではなくアナログ回路によって構成することも本発明の範囲内である。第１のコンピュータ装置１００および第２のコンピュータ装置２００の構成は、その機能を実現できる限りにおいて上述したものに限定されない。

【0120】

４．位置情報に関連付けられたデータを利用してメッシュ統計を生成するためのシステムによる処理
図６は、本発明の位置情報に関連付けられたデータを利用してメッシュ統計を生成するためのシステム１０００による処理６００の一例を示す。

【0121】

ステップＳ６０１では、第１のコンピュータ装置１００のインターフェース部１１０が、位置情報に関連付けられた少なくとも１つのデータを受信する。位置情報に関連付けられたデータは、例えば、衛星画像データ、ポイントデータ、ポリゴンデータであり得る。衛星画像データを受信する場合、衛星画像を受信することにより、衛星画像に含まれる衛星画像データ（例えば、画素値）を抽出するようにしてもよい。

【0122】

衛星画像は、例えば、地表形状の標高を表す標高画像、夜間光強度を表す夜間光画像、雲密度を表す雲画像、輝度温度を表す輝度温度画像、人工構造物を特定する人工構造物特定画像のうちの少なくとも１つを含むが、これらに限定されない。

【0123】

受信された位置情報に関連付けられたデータは、インターフェース部１１０からプロセッサ１２０に渡され、第１のコンピュータ装置１００のプロセッサ１２０はこれを受信する。

【0124】

ステップＳ６０２では、プロセッサ１２０が、複数のメッシュのうち、受信された位置情報に関連付けられたデータが属するメッシュを決定する。これは、例えば、上述した決定手段１２１によって行われる。決定手段１２１は、例えば、後述する図８に示される処理８００によって、受信された衛星画像に複数のメッシュを割り当てることにより、衛星画像データが属するメッシュを決定するができる。決定手段１２１は、例えば、ポイントデータの位置情報から地域メッシュコードまたは世界メッシュコードを算出することによって、受信されたポイントデータが属するメッシュを決定することができる。決定手段１２１は、例えば、ポリゴンデータを形成する複数の点のそれぞれが示す位置情報から地域メッシュコードまたは世界メッシュコードを算出することによって、受信されたポリゴンデータが属するメッシュを決定することができる。ステップＳ６０１で位置情報に関連付けられた複数のデータが受信された場合には、複数のデータのそれぞれについて、属するメッシュを決定することができる。

【0125】

ステップＳ６０３では、プロセッサ１２０が、少なくとも１つのデータに基づいて、そのデータが属するメッシュのメッシュデータを生成する。これは、例えば、上述した生成手段１２２によって行われる。

【0126】

例えば、生成手段１２２は、複数のメッシュのうちの第１のメッシュについて、第１のメッシュに含まれる位置情報に関連付けられたデータの集合Ｏ_１を取得し、Ｏ_１の統計量（例えば、最小値、最大値、中央値、平均値）を算出することによって、第１のメッシュのメッシュデータを生成する。その後、または、これと並列的に、生成手段１２２は、複数のメッシュのうちの第２のメッシュについて、第２のメッシュに含まれる位置情報に関連付けられたデータの集合Ｏ_２を取得し、Ｏ_２の統計量（例えば、最小値、最大値、中央値、平均値）を算出することによって、第２のメッシュのメッシュデータを生成する。その後、または、これと並列的に、生成手段１２２は、複数のメッシュのうちの第ｎのメッシュについて、第ｎのメッシュに含まれる位置情報に関連付けられたデータの集合Ｏ_ｎを取得し、Ｏ_ｎの統計量（例えば、最小値、最大値、中央値、平均値）を算出することによって、第ｎのメッシュのメッシュデータを生成する。このように、生成手段１２２は、複数のメッシュの全てに対して、直列的または並列的にメッシュデータを生成することができる。複数のメッシュデータを並列的に生成することにより、大きなデータ量を有するデータ、特に、多数のデータ点が１つのファイルに入った高分解能の衛星画像であっても計算時間が大きくならない。本明細書において、「並列的」は、「同時に」の意味に限定されず、少なくとも部分的に同時に行われる場合も含むものとする。

【0127】

例えば、位置情報に関連付けられたデータ内にマスク等が存在し、メッシュ内にデータが存在しない場合には、統計量を算出しないようにしてもよい。このとき、そのメッシュについてのメッシュデータを生成する処理を終了してもよい。

【0128】

ステップＳ６０４では、プロセッサ１２０が、インターフェース部１１０を介して第２のデータベース部４００にメッシュデータを送信し、メッシュデータを第２のデータベース部４００に格納する。プロセッサ１２０は、例えば、ステップＳ６０３でメッシュデータを生成する度にステップＳ６０４でメッシュデータを第２のデータベース部４００に格納してもよいし、例えば、ステップＳ６０３で複数のメッシュのうちのいくつかについてメッシュデータを生成した後にステップＳ６０４でメッシュデータを第２のデータベース部４００に格納してもよいし、例えば、ステップＳ６０３で複数のメッシュの全てについてメッシュデータを生成した後にステップＳ６０４でメッシュデータを第２のデータベース部４００に格納してもよい。

【0129】

ステップＳ６０４では、メッシュデータに加えて、メッシュを識別する情報（例えば、地域メッシュコードもしくは世界メッシュコード、および／または、位置情報等）が、第２のデータベース部４００に送信されて、メッシュデータに関連付けられて第２のデータベース部４００に格納され得る。

【0130】

ステップＳ６０１でインターフェース部１１０が、位置情報に関連付けられた第１のデータと位置情報に関連付けられた第２のデータとを含む複数のデータを受信した場合、ステップＳ６０２では、決定手段１２１は、複数のメッシュのうち、第１のデータおよび第２のデータのそれぞれが属するメッシュを決定する。次いで、ステップＳ６０３では、生成手段１２２は、第１のデータおよび第２のデータに基づいて複合メッシュデータを生成することができる。なお、ステップＳ６０２で考慮される複数のメッシュは、第１のデータおよび第２のデータのそれぞれにおいて、同じ位置で区画される同じサイズのメッシュである必要がある。メッシュが区画される位置またはサイズがずれている場合には、それぞれのメッシュデータに演算処理を行うことができず、複合メッシュデータを生成することができなくなるからである。

【0131】

図７は、ステップＳ６０３において、生成手段１２２が、複数のデータ（位置情報に関連付けられた第１のデータおよび位置情報に関連付けられた第２のデータ）に基づいて複合メッシュデータを生成する処理の一例を示す。

【0132】

ステップＳ７０１では、生成手段１２２が、メッシュに属する位置情報に関連付けられた第１のデータに基づいて第１のメッシュデータを生成する。これは上述した処理と同様の処理であり得る。例えば、生成手段１２２は、複数のメッシュのうちの第１のメッシュについて、第１のメッシュに含まれる位置情報に関連付けられた第１のデータの集合Ｏ_１１を取得し、Ｏ_１１の統計量（例えば、最小値、最大値、中央値、平均値）を算出することによって、第１のメッシュの第１のメッシュデータを生成する。その後、または、これと並列的に、生成手段１２２は、複数のメッシュのうちの第２のメッシュについて、第２のメッシュに含まれる位置情報に関連付けられた第１のデータの集合Ｏ_１２を取得し、Ｏ_１２の統計量（例えば、最小値、最大値、中央値、平均値）を算出することによって、第２のメッシュの第１のメッシュデータを生成する。その後、または、これと並列的に、生成手段１２２は、複数のメッシュのうちの第ｎのメッシュについて、第ｎのメッシュに含まれる位置情報に関連付けられた第１のデータの集合Ｏ_１ｎを取得し、Ｏ_１ｎの統計量（例えば、最小値、最大値、中央値、平均値）を算出することによって、第ｎのメッシュの第１のメッシュデータを生成する。

【0133】

ステップＳ７０２では、生成手段１２２が、メッシュに属する位置情報に関連付けられた第２のデータに基づいて第２のメッシュデータを生成する。これも上述した処理と同様の処理であり得る。例えば、生成手段１２２は、複数のメッシュのうちの第１のメッシュについて、第１のメッシュに含まれる位置情報に関連付けられた第２のデータの集合Ｏ_２１を取得し、Ｏ_２１の統計量（例えば、最小値、最大値、中央値、平均値）を算出することによって、第１のメッシュの第２のメッシュデータを生成する。その後、または、これと並列的に、生成手段１２２は、複数のメッシュのうちの第２のメッシュについて、第２のメッシュに含まれる位置情報に関連付けられた第２のデータの集合Ｏ_２２を取得し、Ｏ_２２の統計量（例えば、最小値、最大値、中央値、平均値）を算出することによって、第２のメッシュの第２のメッシュデータを生成する。その後、または、これと並列的に、生成手段１２２は、複数のメッシュのうちの第ｎのメッシュについて、第ｎのメッシュに含まれる位置情報に関連付けられた第２のデータの集合Ｏ_２ｎを取得し、Ｏ_２ｎの統計量（例えば、最小値、最大値、中央値、平均値）を算出することによって、第ｎのメッシュの第２のメッシュデータを生成する。

【0134】

第１のメッシュデータと第２のメッシュデータとが生成されると、ステップＳ７０３において、生成手段１２２が、第１のメッシュデータと第２のメッシュデータとに演算処理を施すことにより、複合メッシュデータを生成する。演算処理は、例えば、四則演算であり得る。これにより、第１のメッシュデータおよび第２のメッシュデータとは異なる事項についての統計量を示す複合メッシュデータを生成することができる。これは、第１のメッシュデータおよび第２のメッシュデータとは別の観点からの統計量分析を可能にする。

【0135】

上述した例ではステップＳ７０１の後にステップＳ７０２が行われるように説明したが、本発明は、これらの順序に限定されない。ステップＳ７０２の後にステップＳ７０１が行われてもよいし、ステップＳ７０１とステップＳ７０２が並列的に行われてもよい。

【0136】

ステップＳ７０３は、例えば、ステップＳ７０１およびステップＳ７０２で第１のメッシュデータおよび第２のメッシュデータが生成される度に行われるようにしてもよいし、例えば、ステップＳ７０１およびステップＳ７０２で複数のメッシュのうちのいくつかについて第１のメッシュデータおよび第２のメッシュデータを生成した後に行われるようにしてもよいし、例えば、ステップＳ７０１およびステップＳ７０２で複数のメッシュの全てについて第１のメッシュデータおよび第２のメッシュデータを生成した後に行われるようにしてもよい。

【0137】

処理６００によって生成され、第１のデータベース部３００に格納された複数のメッシュデータは、メッシュ統計として利用されることができる。例えば、略リアルタイムで更新される衛星画像をソースに用いることにより、略リアルタイムの頻度で複数のメッシュデータを生成することも可能である。これにより、より最新の統計情報を需要者に提供することが可能になる。

【0138】

例えば、格納された複数のメッシュデータは、それらの値を、所望の領域に対する、メッシュデータが関連付けられたメッシュの貢献度に基づいて、集計手段２２１によって集計されることができ、これにより、所望の領域のメッシュ統計が生成される。貢献度は、例えば、所望の領域の形状と、メッシュデータが関連付けられたメッシュとの包含関係に基づいて決定されてもよいし、別のメッシュデータに基づいて決定されてもよい。例えば、集計手段２２１は、図９を参照して後述する処理９００によって、格納された複数のメッシュデータの値を集計することにより、所望の領域のメッシュ統計を生成することができる。

【0139】

図８は、決定手段１２１によって行われる処理８００の一例を示す。処理８００は、受信された１つの衛星画像に複数のメッシュを割り当てることにより、衛星画像データが属するメッシュを決定するための処理である。

【0140】

決定手段１２１が衛星画像を受信すると、ステップＳ８０１では、決定手段１２１が、受信された衛星画像の四隅の座標（緯度，経度）を求める。決定手段１２１は、例えば、衛星画像に含まれる位置情報を参照することにより、四隅の座標（緯度，経度）を求めることができる。決定手段１２１は、例えば、衛星画像に写っている、座標が既知の特徴との相対的な位置関係から、四隅の座標（緯度，経度）を求めることができる。

【0141】

例えば、衛星画像が、最小緯度minlat、最大緯度maxlat、最小経度minlong、最大経度maxlongを有する場合、四隅の座標は、（maxlat，minlong）、（minlat，minlong）、（maxlat，maxlong）、（minlat，maxlong）と表され得る。

【0142】

ステップＳ８０２では、決定手段１２１が、北西端のメッシュｗ_ＮＷを特定する。メッシュｗ_ＮＷは、座標（maxlat，minlong）を含むメッシュである。

【0143】

ステップＳ８０３では、決定手段１２１が、メッシュｗ_ＮＷの中央の座標（緯度，経度）を求める。メッシュｗ_ＮＷが最小緯度lat0、最大緯度lat1、最小経度long0、最大経度long1を有する場合、メッシュｗ_ＮＷの中央の座標（latc，longc）は、
（latc，longc）＝（（lat0＋lat1）／２，（long0＋long1）／２）
によって求められる。

【0144】

ここで、メッシュの緯度方向幅Δlatおよび経度方向幅Δlongとすると、
最小緯度lat0＝maxlat－Δlat
最大緯度lat1＝maxlat
最小経度long0＝minlong
最大経度long1＝minlong＋Δlong
として求められ得る。緯度方向幅Δlatおよび経度方向幅Δlongは、生成しようとするメッシュに応じて変動し得る。例えば、１次メッシュを割り当てる場合、Δlatは、４０分であり、Δlongは、１度であり得る。例えば、２次メッシュを割り当てる場合、Δlatは、５分であり、Δlongは、７．５分であり得る。例えば、３次メッシュを割り当てる場合、Δlatは、３０秒であり、Δlongは、４５秒であり得る。例えば、２分の１メッシュを割り当てる場合、Δlatは、１５秒であり、Δlongは、２２．５秒であり得る。例えば、４分の１メッシュを割り当てる場合、Δlatは、７．５秒であり、Δlongは、１１．２５秒であり得る。例えば、８分の１メッシュを割り当てる場合、Δlatは、３．７５秒であり、Δlongは、５．６２５秒であり得る。

【0145】

ステップＳ８０４では、決定手段１２１が、衛星画像を囲むメッシュの数Ｎ、Ｍを求める。ここで、Ｍは緯度方向のメッシュの数であり、Ｎは経度方向のメッシュの数であり、ＮおよびＭは以下の式によって算出される。
Ｍ＝_└(maxlat-minlat)/(Δlat)_┘
Ｎ＝_└(maxlong-minlong)/(Δlong)_┘

【0146】

ステップＳ８０５では、決定手段１２１が、衛星画像を覆うメッシュの集合Ｗ＝ｗ_１，ｗ_２，．．．，ｗ_ＮＭ、すなわち、メッシュｗ_ｉの集合（１≦ｉ≦Ｎ×Ｍ）を求める。

【0147】

ｗ_ｉは、中心座標（lati，longi）を有するメッシュであり、中心座標（lati，longi）は、
（lati，longi）＝（latc＋（yy－１）×Δlat，longc＋（xx－１）×Δlong）
と表され、ここで、１≦yy≦Ｍであり、１≦xx≦Ｎである。

【0148】

上述した処理によって割り当てられたメッシュの集合Ｗのそれぞれのメッシュは、中央の座標ならびに緯度方向幅Δlatおよび経度方向幅Δlongが分かっているので、メッシュの北西端の座標および南東端の座標を算出することが可能である。

【0149】

なお、処理８００は、衛星画像にメッシュを割り当てる処理の一例であり、他の処理によっても、衛星画像にメッシュを割り当てることができることが理解される。

【0150】

図９は、集計手段２２１によって行われる処理９００の一例を示す。処理９００は、所望の領域とメッシュとの包含関係に基づいて、第２のデータベース部４００に格納された複数のメッシュデータの値を集計することにより、所望の領域のメッシュ統計を提供するための処理である。

【0151】

集計手段２２１は、第２のデータベース部４００からインターフェース部１１０を介して所望の領域の形状データを受信しているものとする。また、集計手段２２１は、所望の領域の形状データによって定義されるポリゴンＶを覆い得るｍ個のメッシュｗ_ｉに対応するｍ個のメッシュデータの値ｘ_ｉおよびそれらのメッシュを識別する情報も受信しているものとする。ここで、１≦ｉ≦ｍであり、ｍは、２以上の整数である。

【0152】

処理９００は、ステップＳ９０１で開始し、ｉ＝１が定義される。すなわち、ｍ個のメッシュのうち、メッシュｗ_ｉ＝１についての処理が行われる。

【0153】

ステップＳ９０２では、集計手段２２１が、形状データによって定義されるポリゴンＶと、メッシュｗ_ｉとの包含関係を決定する。ポリゴンＶと、メッシュｗ_ｉとの包含関係は、例えば、図１０に説明される関係である。

【0154】

図１０は、ポリゴンＶとメッシュｗとの包含関係を説明する図である。

【0155】

（１）は、ポリゴンＶがメッシュｗに完全に内包されている場合を示している。この場合、Ｓ（ｗ）＞Ｓ（Ｖ）＝Ｓ（Ｖ∩ｗ）＞０となる。ここで、Ｓは面積を表す。すなわち、Ｓ（ｗ）は、メッシュｗの面積であり、Ｓ（Ｖ）は、ポリゴンＶの面積であり、Ｓ（Ｖ∩ｗ）は、メッシュｗとポリゴンＶとの共通部分の面積である。

【0156】

（２）は、ポリゴンＶがメッシュｗを完全に外包している場合を示している。この場合、Ｓ（Ｖ）≧Ｓ（ｗ）＝Ｓ（Ｖ∩ｗ）＞０となる。

【0157】

（３）は、ポリゴンＶがメッシュｗと共通部分を有している場合を示している。この場合、Ｓ（Ｖ），Ｓ（ｗ）＞Ｓ（Ｖ∩ｗ）＞０となる。

【0158】

（４）は、ポリゴンＶがメッシュｗと共通部分を有していない場合を示している。この場合、Ｓ（Ｖ∩ｗ）＝０となる。

【0159】

再び図９を参照して、ステップＳ９０３では、集計手段２２１が、ステップＳ９０２で決定された包含関係に基づいて、ポリゴンＶのメッシュｗ_ｉに対する貢献度ρ（Ｖ，ｗ_ｉ）を算出する。貢献度ρ（Ｖ，ｗ_ｉ）は、ポリゴンＶがメッシュｗ_ｉに対してどのくらいの割合を占めているかの程度を示す値であり、例えば、以下の式によって定義される。
ρ（Ｖ，ｗ_ｉ）＝Ｓ（Ｖ∩ｗ）／Ｓ（ｗ）

【0160】

すなわち、図１０を参照して説明した包含関係において、（２）ポリゴンＶがメッシュｗを完全に外包している場合、ρ（Ｖ，ｗ_ｉ）＝１となり、（４）ポリゴンＶがメッシュｗと共通部分を有していない場合、ρ（Ｖ，ｗ_ｉ）＝０となり、（１）ポリゴンＶがメッシュｗに完全に内包されている場合、および、（３）ポリゴンＶがメッシュｗと共通部分を有している場合、１＞ρ（Ｖ，ｗ_ｉ）＞０となる。

【0161】

貢献度が算出されると、ステップＳ９０４では、集計手段２２１が、算出された貢献度ρ（Ｖ，ｗ_ｉ）をメッシュデータの値ｘ_ｉに乗じる。

【0162】

ステップＳ９０５では、ｉがインクリメントされる。

【0163】

ステップＳ９０６では、集計手段２２１が、ｉとメッシュの数ｍとを比較し、ｉがｍよりも大きい場合、ステップＳ９０７に進む。ｉがｍ以下の場合、ステップＳ９０２に戻り、ｉがｍよりも大きくなるまで、すなわち、ｍ個のメッシュ全てに対してρ（Ｖ，ｗ_ｉ）×ｘ_ｉを算出するまで、ステップＳ９０２～ステップＳ９０６を繰り返す。

【0164】

ステップＳ９０７では、集計手段２２１が、ｍ個のメッシュ全てに対して算出されたρ（Ｖ，ｗ_ｉ）×ｘ_ｉを総和する。得られた値が、所望の領域のメッシュ統計となり得る。

【0165】

このようにして、所望の領域とメッシュとの包含関係に基づいて、第２のデータベース部４００に格納された複数のメッシュデータの値を集計することにより、所望の領域のメッシュ統計を提供することで、例えば、統計データが存在しない領域（例えば、国、都道府県、市区町村等の行政界区画等）であっても、衛星画像から抽出可能な統計量に関するものにはなるが、メッシュ統計を提供することができる。さらに、略リアルタイムで更新される衛星画像をソースに用いることにより、略リアルタイムの頻度で更新可能なメッシュ統計を提供することも可能である。これにより、より最新の統計情報を需要者に提供することが可能になる。

【0166】

上述した例では、第２のデータベース部４００に格納された複数のメッシュデータの値を集計することを説明したが、システム１０００の外部から受信された複数のメッシュデータの値を集計することも本発明の範囲内である。

【0167】

上述した例では、ｍ個のメッシュ全てに対して算出されたρ（Ｖ，ｗ_ｉ）×ｘ_ｉを総和（Σρ（Ｖ，ｗ_ｉ）×ｘ_ｉ）することによって、所望の領域のメッシュ統計を提供したが、例えば、貢献度ρ（Ｖ，ｗ_ｉ）が０ではないメッシュのメッシュデータｘ_ｉを総和（Σ_{（ρ（Ｖ，ｗｉ）≠０）}ｘ_ｉ）することによってメッシュ統計を提供するようにしてもよいし、貢献度ρ（Ｖ，ｗ_ｉ）が１であるメッシュのメッシュデータｘ_ｉを総和（Σ_{（ρ（Ｖ，ｗｉ）＝１）}ｘ_ｉ）することによってメッシュ統計を提供するようにしてもよい。Σρ（Ｖ，ｗ_ｉ）×ｘ_ｉが、真の統計値により近似した値を与え得るが、Σ_{（ρ（Ｖ，ｗｉ）≠０）}ｘ_ｉまたはΣ_{（ρ（Ｖ，ｗｉ）＝１）}ｘ_ｉもまた、真の統計値にある程度近似した値を与え得る。

【0168】

上述した例では、所望の領域とメッシュとの包含関係に基づいて貢献度を算出したが、貢献度は、包含関係に基づくものに限定されない。貢献度は、例えば、他のメッシュデータに基づいて算出されてもよい。例えば、貢献度は、衛星画像データから生成された人工物占有面積に関するメッシュデータに基づいて算出されてもよく、この貢献度は、人工物が密集している区域であるかどうかの指標となる。例えば、人口についてのメッシュデータの値を所望の領域について再集計する場合にこの貢献度を利用することにより、人工物が密集している区域ほど人口に対する貢献度が大きく、人工物が点在している区域ほど人口に対する貢献度が小さいものとして、人口についてのメッシュデータの値を集計することができる。これにより、結果として得られるメッシュ統計の精度がより高いものとなる。

【0169】

図１１は、集計手段２２１によって行われる処理１１００の一例を示す。処理１１００は、ユーザから受けた要求に基づいて、所望のメッシュ統計を提供するための処理である。

【0170】

第２のコンピュータ装置２００のインターフェース部２１０が、提供すべきメッシュ統計の指示を受信すると、ステップＳ１１０１では、集計手段２２１が、インターフェース部２１０から、提供すべきメッシュ統計の指示を受信する。提供すべきメッシュ統計の指示は、例えば、ユーザから送信された指示であり得る。

【0171】

ステップＳ１１０２では、集計手段２２１が、提供すべきメッシュ統計の指示に基づいて、メッシュデータのうち集計すべきメッシュデータを選択する。集計手段２２１は、例えば、第２のデータベース部４００に格納されている複数のメッシュデータを検索することにより、集計すべきメッシュデータを選択してもよいし、インターフェース部２１０を介して既に受信されている複数のメッシュデータから、集計すべきメッシュデータを選択してもよい。集計手段２２１は、例えば、提供すべきメッシュ統計のために必要な事項を示すメッシュデータを選択するようにしてもよいし、提供すべきメッシュ統計のために必要な事項の必要な統計量（例えば、最大値、最小値、中央値、平均値等）を示すメッシュデータを選択するようにしてもよい。

【0172】

ステップＳ１１０３では、集計手段２２１は、ステップＳ１１０２で選択されたメッシュデータの値を集計する。集計手段２２１は、例えば、ステップＳ１１０２で選択された第１の事項の統計量を示す第１のメッシュデータと、ステップＳ１１０２で選択された第２の事項の統計量を示す第２のメッシュデータとに演算処理を施す。これにより、提供すべきメッシュ統計が作成される。演算処理は、例えば、四則演算であり得る。

【0173】

例えば、旅行客が自然派であるか都会派であるかのプレファレンスを知りたいというユーザの要望から、明るい場所に宿泊した人数および暗い場所に宿泊した人数を示すメッシュ統計を提供すべき場合、ステップＳ１１０２で、集計手段２２１は、夜間光画像から生成された夜間光強度の平均値を示すメッシュデータと、宿泊施設のポイントデータから生成された宿泊人数の平均値または中央値を示すメッシュデータとを選択し得る。ステップＳ１１０３では、集計手段２２１は、宿泊人数の平均値または中央値を示すメッシュデータの値を夜間光強度の平均値を示すメッシュデータの値で除算することにより、夜間光強度当たりの宿泊人数を示すメッシュ統計を提供することができる。

【0174】

例えば、旅行客が山派であるか海派であるかのプレファレンスを知りたいというユーザの要望から、周囲に森林が多い場所に宿泊した人数および海が近い場所に宿泊した人数を示すメッシュ統計を提供すべき場合、ステップＳ１１０２で、集計手段２２１は、地表の物理状態を示す土地被覆画像から生成された土地利用状態を示すメッシュデータと、宿泊施設のポイントデータから生成された宿泊人数の平均値を示すメッシュデータとを選択し得る。ステップＳ１１０３では、集計手段２２１は、土地利用状態を示すメッシュデータの値に宿泊人数の平均値を示すメッシュデータの値を加算することにより、宿泊人数の平均値に土地利用状態のラベルを付したメッシュ統計を提供することができる。

【0175】

例えば、標高が高い場所に住んでいる人数を知りたいというユーザの要望から、標高毎の人口を示すメッシュ統計を提供すべき場合、ステップＳ１１０２で、集計手段２２１は、夜間光画像から生成された夜間光強度の平均値を示すメッシュデータと、標高画像から生成された平均標高を示すメッシュデータとを選択し得る。ステップＳ１１０３では、集計手段２２１は、夜間光強度の平均値を示すメッシュデータの値に平均標高を示すメッシュデータの値を加算することにより、夜間光強度に平均標高のラベルを付したメッシュ統計を提供することができる。

【0176】

図６、図７、図８、図９、図１１を参照して上述した例では、図６、図７、図８、図９、図１１に示される各ステップの処理は、プロセッサ１２０または２２０とメモリ１３０または２３０に格納されたプログラムとによって実現することが説明されたが、本発明はこれに限定されない。図６、図７、図８、図９、図１１に示される各ステップの処理のうちの少なくとも１つは、制御回路などのハードウェア構成によって実現されてもよい。

【0177】

本発明は、上述した実施形態に限定されるものではない。本発明は、特許請求の範囲によってのみその範囲が解釈されるべきであることが理解される。当業者は、本発明の具体的な好ましい実施形態の記載から、本発明の記載および技術常識に基づいて等価な範囲を実施することができることが理解される。

【産業上の利用可能性】

【0178】

本発明は、メッシュ統計を効率的に生成することが可能なシステム、方法、およびプログラム等を提供するものとして有用である。

【符号の説明】

【0179】

１００ コンピュータ装置
１１０ インターフェース部
１２０ プロセッサ
１３０メモリ
２００ データベース部
１０００ システム

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】