明細書 :高解像度化装置、高解像度化方法、空間データ高解像度化方法、およびプログラム

発行国	日本国特許庁(JP)
公報種別	公開特許公報(A)
公開番号	特開2019-191936 (P2019-191936A)
公開日	令和元年10月31日(2019.10.31)
発明の名称または考案の名称	高解像度化装置、高解像度化方法、空間データ高解像度化方法、およびプログラム
国際特許分類	G06F  17/17        (2006.01) G06F  16/00        (2019.01)
ＦＩ	G06F 17/17 G06F 17/30 220Z
請求項の数または発明の数	7
出願形態	OL
全頁数	31
出願番号	特願2018-084055 (P2018-084055)
出願日	平成30年4月25日(2018.4.25)
発明者または考案者	【氏名】田中 佑典 【氏名】岩田 具治 【氏名】倉島 健 【氏名】戸田 浩之 【氏名】田中 利幸
出願人	【識別番号】000004226 【氏名又は名称】日本電信電話株式会社 【識別番号】504132272 【氏名又は名称】国立大学法人京都大学
個別代理人の代理人	【識別番号】110001519、【氏名又は名称】特許業務法人太陽国際特許事務所
審査請求	未請求
テーマコード	5B056
Ｆターム	5B056BB52 5B056HH03
要約	【課題】異なる粒度を持つ補助的な空間データを用いて、空間データを精度よく高解像度化することができるようにする。 【解決手段】空間補間部１５０が、値の種類を示す複数のドメインの各々について、当該ドメインの補助的な空間データを、任意の空間補間法を用いてターゲットとする粒度の当該補助的な空間データに補間し、パラメータ推定部１６０が、変換対象の空間データと、複数のドメインの各々についてのターゲットとする当該粒度の当該補助的な空間データの各々とに基づいて、当該変換対象の空間データをターゲットとする当該粒度の当該変換対象の空間データに変換するための回帰モデルのパラメータを推定する。 【選択図】図１
特許請求の範囲	【請求項1】 空間を分割するエリアの単位を示す粒度でのエリアごとに紐付く値からなる変換対象の空間データを、値の種類が異なる補助的な空間データを用いてターゲットとする前記粒度の前記変換対象の空間データに変換するための回帰モデルを学習する高解像度化装置であって、 前記値の種類を示す複数のドメインの各々について、前記ドメインの前記補助的な空間データを、任意の空間補間法を用いて前記ターゲットとする前記粒度の前記補助的な空間データに補間する空間補間部と、 前記変換対象の空間データと、前記複数のドメインの各々についての前記ターゲットとする前記粒度の前記補助的な空間データの各々とに基づいて、前記回帰モデルのパラメータを推定するパラメータ推定部と、 を含む高解像度化装置。 【請求項2】 前記空間補間部は、前記複数のドメインの各々について、前記空間補間法としてガウス過程回帰を用いて、前記ターゲットとなる前記粒度の前記補助的な空間データとして、前記ターゲットとする前記粒度でのエリアごとに前記ドメインの値の平均及び分散を求め、 前記パラメータ推定部は、前記変換対象の空間データと、前記空間補間部により前記複数のドメインの各々について求められた、前記ターゲットとする前記粒度でのエリア毎に前記ドメインの値の平均及び分散とに基づいて、前記回帰モデルのパラメータを推定する 請求項１記載の高解像度化装置。 【請求項3】 前記変換対象の空間データと、前記複数のドメインの各々についての前記ターゲットとする前記粒度の前記補助的な空間データの各々と、前記回帰モデルとに基づいて、前記ターゲットとする前記粒度の前記変換対象の空間データを算出する高解像度データ算出部 を更に含む請求項１又は２記載の高解像度化装置。 【請求項4】 空間を分割するエリアの単位を示す粒度でのエリアごとに紐付く値からなる変換対象の空間データを、値の種類が異なる補助的な空間データを用いてターゲットとする前記粒度の前記変換対象の空間データに変換するための回帰モデルを学習する高解像度化方法であって、 空間補間部が、前記値の種類を示す複数のドメインの各々について、前記ドメインの前記補助的な空間データを、任意の空間補間法を用いて前記ターゲットとする前記粒度の前記補助的な空間データに補間し、 パラメータ推定部が、前記変換対象の空間データと、前記複数のドメインの各々についての前記ターゲットとする前記粒度の前記補助的な空間データの各々とに基づいて、前記回帰モデルのパラメータを推定する 高解像度化方法。 【請求項5】 前記空間補間部は、前記複数のドメインの各々について、前記空間補間法としてガウス過程回帰を用いて、前記ターゲットとなる前記粒度の前記補助的な空間データとして、前記ターゲットとする前記粒度でのエリアごとに前記ドメインの値の平均及び分散を求め、 前記パラメータ推定部は、前記変換対象の空間データと、前記空間補間部により前記複数のドメインの各々について求められた、前記ターゲットとする前記粒度でのエリア毎に前記ドメインの値の平均及び分散とに基づいて、前記回帰モデルのパラメータを推定する 請求項４記載の高解像度化方法。 【請求項6】 高解像度データ算出部が、前記変換対象の空間データと、前記複数のドメインの各々についての前記ターゲットとする前記粒度の前記補助的な空間データの各々と、前記回帰モデルとに基づいて、前記ターゲットとする前記粒度の前記変換対象の空間データを算出する 請求項４又は５記載の高解像度化方法。 【請求項7】 コンピュータを、請求項１～３の何れか１項記載の高解像度化装置の各部として機能させるためのプログラム。
発明の詳細な説明	【技術分野】 【0001】 本発明は、高解像度化装置、高解像度化方法、およびプログラムに係り、特に、空間データを高解像度化するための高解像度化装置、高解像度化方法、およびプログラムに関する。 【背景技術】 【0002】 近年、政府や企業などの組織が、都市環境や事業の改善を目的として、様々な種類の空間データ（貧困度、大気汚染度、犯罪数、人口、交通量など）を収集し、公開している。空間データとは、エリア（緯度経度、住所、地域など）とそれに紐づく何らかの値とのペアで与えられるデータを指す。 【0003】 このような空間データは、収集コストが高く空間的に密なデータを収集することが難しいため、空間を分割する単位を広く取ったエリア、すなわち粒度のある程度粗いエリアにおいて、平均化されて提供されることが多い。 【0004】 しかし、より良い都市環境の改善のためには、高解像度の空間データが望ましい。例えば、貧困度が高いエリアや大気汚染度が高いエリアを詳細に絞りこむことで、粗い粒度の空間データよりも適切な政策の決定等の介入が可能となる。 【0005】 したがって、粗い粒度で収集された低解像度の空間データを、細かい粒度で表わされる高解像度の空間データに変換することが考えられる。 【0006】 例えば、ターゲットとする低解像度の空間データとは別に、他の種類の高解像度の空間データを用意し、それらを粒度の高解像度化のために補助的に用いて学習した回帰モデルを用いることによって、ターゲットとする空間データの細かい粒度における値を予測する技術が存在する（非特許文献１）。 【0007】 また、例えば、空間相関（位置が近い場合に、その位置に紐づくサンプル値も近い値をとるという性質）と前述の回帰モデルとの両方の効果を加味して、ターゲットとする空間データの細かい粒度における値を予測する技術が存在する（非特許文献２）。 【先行技術文献】 【0008】 【非特許文献1】C. Smith-Clarke, A. Mashhadi and L. Capra, “Poverty on the Cheap: Estimating Poverty Maps Using Aggregated Mobile Communication Networks”, Proceedings of the SIGCHI Conference on Human Factors in Computing Systems, 2014, pp. 511-520. 【非特許文献2】N. W. Park, “Spatial Downscaling of TRMM Precipitation Using Geostatistics and Fine Scale Environmental Variables”, Advances in Meteorology, volume 2013, article ID 237126, 2013. 【発明の概要】 【発明が解決しようとする課題】 【0009】 上記の従来技術（非特許文献１、２）は、ターゲットとする空間データが持つ空間相関を考慮しつつ、他の種類の解像度の高い空間データを補助的なデータとして用いて学習した回帰モデルに基づいて、ターゲットとする空間データの細かい粒度における値を予測することが可能である。 【0010】 ここで、従来技術では、“補助的な空間データの粒度（住所、地域などの集計単位）は互いに等しい”ことが仮定されている。 【0011】 しかし、この仮定は、より多くの種類の空間データを補助的に用いる場合には、必ずしも適切ではない。なぜなら、政府や企業などが収集する空間データは、各組織の方針に基づいて独自の粒度で集計されるからである。 【0012】 例えば、政府による国勢調査では、住所ごとの細かい粒度の人口密度データが得られるが、犯罪発生率データは警察の管轄エリア毎に得られる場合などが考えられる。 【0013】 このように、従来技術では、空間データ毎に粒度が異なる場合には適用困難である、という問題があった。 【0014】 本発明は、上記の点に鑑みてなされたものであり、異なる粒度を持つ補助的な空間データを用いて、空間データを精度よく高解像度化することができる高解像度化装置、高解像度化方法、およびプログラムを提供することを目的とする。 【課題を解決するための手段】 【0015】 本発明に係る高解像度化装置は、空間を分割するエリアの単位を示す粒度でのエリアごとに紐付く値からなる変換対象の空間データを、値の種類が異なる補助的な空間データを用いてターゲットとする前記粒度の前記変換対象の空間データに変換するための回帰モデルを学習する高解像度化装置であって、前記値の種類を示す複数のドメインの各々について、前記ドメインの前記補助的な空間データを、任意の空間補間法を用いて前記ターゲットとする前記粒度の前記補助的な空間データに補間する空間補間部と、前記変換対象の空間データと、前記複数のドメインの各々についての前記ターゲットとする前記粒度の前記補助的な空間データの各々とに基づいて、前記回帰モデルのパラメータを推定するパラメータ推定部と、を備えて構成される。 【0016】 また、本発明に係る高解像度化方法は、空間を分割するエリアの単位を示す粒度でのエリアごとに紐付く値からなる変換対象の空間データを、値の種類が異なる補助的な空間データを用いてターゲットとする前記粒度の前記変換対象の空間データに変換するための回帰モデルを学習する高解像度化方法であって、空間補間部が、前記値の種類を示す複数のドメインの各々について、前記ドメインの前記補助的な空間データを、任意の空間補間法を用いて前記ターゲットとする前記粒度の前記補助的な空間データに補間し、パラメータ推定部が、前記変換対象の空間データと、前記複数のドメインの各々についての前記ターゲットとする前記粒度の前記補助的な空間データの各々とに基づいて、前記回帰モデルのパラメータを推定する。 【0017】 本発明に係る高解像度化装置及び高解像度化方法によれば、空間補間部が、空間を分割するエリアの単位を示す粒度でのエリアごとに紐付く値の種類を示す複数のドメインの各々について、当該ドメインの補助的な空間データを、任意の空間補間法を用いてターゲットとする粒度の当該補助的な空間データに補間する。 【0018】 そして、パラメータ推定部が、当該変換対象の空間データと、複数のドメインの各々についての当該ターゲットとする当該粒度の当該補助的な空間データの各々とに基づいて、当該粒度の当該変換対象の空間データに変換するための回帰モデルのパラメータを推定する。 【0019】 このように、値の種類を示す複数のドメインの各々について、当該ドメインの補助的な空間データを、任意の空間補間法を用いてターゲットとする粒度の当該補助的な空間データに補間し、変換対象の空間データと、複数のドメインの各々についてのターゲットとする当該粒度の当該補助的な空間データの各々とに基づいて、当該変換対象の空間データをターゲットとする当該粒度の当該変換対象の空間データに変換するための回帰モデルのパラメータを推定することにより、異なる粒度を持つ補助的な空間データを用いて、空間データを精度よく高解像度化することができる。 【0020】 また、本発明に係る高解像度化装置の前記空間補間部は、前記複数のドメインの各々について、前記空間補間法としてガウス過程回帰を用いて、前記ターゲットとなる前記粒度の前記補助的な空間データとして、前記ターゲットとする前記粒度でのエリアごとに前記ドメインの値の平均及び分散を求め、前記パラメータ推定部は、前記変換対象の空間データと、前記空間補間部により前記複数のドメインの各々について求められた、前記ターゲットとする前記粒度でのエリア毎に前記ドメインの値の平均及び分散とに基づいて、前記回帰モデルのパラメータを推定することができる。 【0021】 また、本発明に係る高解像度化方法の前記空間補間部は、前記複数のドメインの各々について、前記空間補間法としてガウス過程回帰を用いて、前記ターゲットとなる前記粒度の前記補助的な空間データとして、前記ターゲットとする前記粒度でのエリアごとに前記ドメインの値の平均及び分散を求め、前記パラメータ推定部は、前記変換対象の空間データと、前記空間補間部により前記複数のドメインの各々について求められた、前記ターゲットとする前記粒度でのエリア毎に前記ドメインの値の平均及び分散とに基づいて、前記回帰モデルのパラメータを推定することができる。 【0022】 また、本発明に係る高解像度化装置は、前記変換対象の空間データと、前記複数のドメインの各々についての前記ターゲットとする前記粒度の前記補助的な空間データの各々と、前記回帰モデルとに基づいて、前記ターゲットとする前記粒度の前記変換対象の空間データを算出する高解像度データ算出部を備えることができる。 【0023】 また、本発明に係る高解像度化方法は、高解像度データ算出部が、前記変換対象の空間データと、前記複数のドメインの各々についての前記ターゲットとする前記粒度の前記補助的な空間データの各々と、前記回帰モデルとに基づいて、前記ターゲットとする前記粒度の前記変換対象の空間データを算出することができる。 【0024】 本発明に係るプログラムは、上記の高解像度化装置の各部として機能させるためのプログラムである。 【発明の効果】 【0025】 本発明の高解像度化装置、高解像度化方法、およびプログラムによれば、異なる粒度を持つ補助的な空間データを用いて、空間データを精度よく高解像度化することができる。 【図面の簡単な説明】 【0026】 【図1】本発明の実施の形態に係る高解像度化装置の構成を示すブロック図である。 【図2】本発明の実施の形態における粒度の例（上図）及び出力例（下図）である。 【図3】本発明の実施の形態に係る高解像度化装置の空間データ高解像度化処理ルーチンを示すフローチャートである。 【図4】本発明の実施の形態に係る高解像度化装置の検索処理ルーチンを示すフローチャートである。 【発明を実施するための形態】 【0027】 以下、本発明の実施の形態について図面を用いて説明する。 【0028】 ＜本発明の実施の形態に係る高解像度化装置の概要＞ まず、本発明の実施形態の概要について説明する。 【0029】 本実施形態では、様々な粒度を持つ様々なドメイン（種類）の空間データを補助的に用いて、変換対象とする粗い粒度の空間データとの間の統計的な関係を回帰モデルにより学習する。 【0030】 具体的には、クリギング（参考文献１）やガウス過程回帰（参考文献２）等の空間補間法を、様々な粒度を持つ補助的な空間データに適用することにより、これらの補助的な空間データの粒度を均一にする。そして、それらを入力として、変換対象の空間データとの統計的な関係を回帰モデルにより学習する。 ［参考文献１］T. Hengl, G. B. M. Heuvelink and A. Stein, “A Generic Framework for Spatial Prediction of Soil Variables based on Regression-Kriging”, Geoderma, volume 120, issues 11-2, 2004, pages 75-93. ［参考文献２］C. E. Rasmussen and C. K. I. Williams, Gaussian Processes for Machine Learning, MIT Press, 2006. 【0031】 このような構成により、粒度が互いに異なる補助的な空間データが与えられた場合においても、回帰モデルを学習可能である。 【0032】 例えば、住所ごとの人口密度データと、警察の管轄エリアごとの犯罪発生率データのように、集計化される単位が異なる場合においても、それらのデータを用いて回帰モデルを学習し、ターゲットとする粗い粒度の空間データの高解像度化に利用することができる。 【0033】 また、補助的な空間データに対して空間補間法を適用し、任意の点におけるサンプル値を予測することができるため、ターゲットとする粒度を任意に設定することができる。 【0034】 例えば、１００ｍ×１００ｍのグリッドで区切られたエリア集合等をターゲットとする粒度として設定することも可能である。これにより政府や企業等が要求する粒度に基づくエリア集合に対して、高解像度なデータを算出することができる。 【0035】 本実施形態では、様々な粒度を持つ補助的な空間データに対してガウス過程回帰を用いた空間補間を行う場合について説明する。ガウス過程回帰を用いることで、任意の点における予測値と分散（予測値に対する信頼度）を計算し、この予測値に対する信頼度を加味することで、効果的に回帰モデルを学習する。 【0036】 回帰モデルの学習では、補助的に用いる空間データのそれぞれに対する重要度を学習する。ここで、重要度とは、補助的な空間データのそれぞれが、ターゲットとする高解像度な空間データをどの程度説明できるかの度合いを表す。この重要度はターゲットとする空間データとの相関のみによって決まるわけではなく、補助的に用いる空間データの粒度の粗さにも依存する。 【0037】 例えば、変換対象とする空間データＡ、補助的な空間データＢ、Ｃがあったとし、ＡとＢと、及びＡとＣとの間の相関がそれぞれ全く同じだったとする。このとき、Ｂの空間粒度は細かく、Ｃの空間粒度は粗いとすると、高解像度な空間データをより説明できるのは細かい粒度の空間データＢであると考えられる。 【0038】 このように、補助的に用いるそれぞれの空間データの重要度を、それらの持つ空間粒度の違いを加味しつつ回帰モデルの推定が可能であるというメリットがある。 【0039】 したがって、補助的に用いる空間データの粒度の違いを加味しつつ効果的に回帰モデルを学習でき、学習済みの回帰モデルを用いることにより変換対象とする粗い粒度の空間データを高解像度に変換することができる。 【0040】 ＜本発明の実施の形態に係る高解像度化装置の構成＞ 図１を参照して、本発明の実施の形態に係る高解像度化装置１０の構成について説明する。図１は、本発明の実施の形態に係る高解像度化装置１０の構成を示すブロック図である。 【0041】 高解像度化装置１０は、ＣＰＵと、ＲＡＭと、後述する高解像度化処理ルーチン及び検索処理ルーチンを実行するためのプログラムを記憶したＲＯＭとを備えたコンピュータで構成され、機能的には次に示すように構成されている。 【0042】 図１に示すように、本実施形態に係る高解像度化装置１０は、補助的空間データ格納部１００と、ターゲット粒度エリア中心格納部１１０と、変換対象空間データ格納部１２０と、粒度エリア対応表格納部１３０と、操作部１４０と、空間補間部１５０と、パラメータ推定部１６０と、高解像度データ算出部１７０と、高解像度データ格納部１８０と、入力部１９０と、検索部２００と、出力部２１０とを備えて構成される。 【0043】 補助的空間データ格納部１００、ターゲット粒度エリア中心格納部１１０、変換対象空間データ格納部１２０、粒度エリア対応表格納部１３０、及び高解像度データ格納部１８０は、Ｗｅｂサーバや、データベースを具備するデータベースサーバ等であり、要求にしたがって、データを読み出し、該当のデータを要求元に送信する。 【0044】 補助的空間データ格納部１００は、変換対象の空間データの高解像度化のために補助的に利用され得る補助的な空間データを、値の種類を示すドメインごとに格納する。 【0045】 具体的には、入力空間全体を表す集合を とする。ここで、 はユークリッド空間の部分集合であり、例えば対象とする都市全体を表現するものである。 【0046】 また、ドメインを表す集合を とし、ドメイン の空間データが紐づくエリア集合を とする。ここで、 は の分割であり、 の要素であるエリアの和集合は に一致する。 【0047】 また、エリア におけるドメインｓの空間データの値を と表す。 【0048】 そして、ドメインsの空間データをまとめて 次元ベクトル とする。ここで、 は に含まれるエリアの数を表す。 【0049】 また、ドメインsのエリア における中心座標を とし、まとめて とする。 【0050】 すなわち、補助的空間データ格納部１００は、 を補助的な空間データとして格納している。 【0051】 ターゲット粒度エリア中心格納部１１０は、ターゲットとする粒度のエリア集合に対する各エリアｊの中心座標の集合 を、粒度ごとに格納する。 【0052】 ここで、粒度とは、空間を分割するエリアの単位である。例えば、図２上図に示すように、ターゲットとなり得る粒度であるターゲット粒度候補が、複数用意されており、ターゲット粒度候補Ａとターゲット粒度候補Ｂとでは、ターゲット粒度候補Ａの方が粗い粒度である。 【0053】 具体的には、ターゲットとする粒度のエリア集合を とする。細かいエリア における中心座標を とし、まとめて とする。 【0054】 ここで、 は、 に含まれるエリアの数を表す。 【0055】 変換対象空間データ格納部１２０は、変換対象とする粗い粒度の空間データ を格納する。 【0056】 ここでは、変換対象の空間データの粗い粒度のエリア集合を とする。また、粗いエリア におけるターゲットの空間データの値を と表し、まとめて 次元ベクトル とする。 【0057】 ここで、 は、 に含まれるエリアの数を表す。 【0058】 粒度エリア対応表格納部１３０は、ターゲットとする粒度の各エリアが、粗い粒度のどのエリアに含まれるかを表す対応表 を格納する。 【0059】 ここで、 は、粗い粒度のエリアｉに含まれるターゲットとする粒度のエリアｊの集合を表す。 【0060】 操作部１４０は、補助的空間データ格納部１００、ターゲット粒度エリア中心格納部１１０、変換対象空間データ格納部１２０、及び粒度エリア対応表格納部１３０が格納するデータに対するユーザからの各種操作を受け付ける。 【0061】 各種操作とは、例えば、格納された情報を登録、修正、削除する操作等である。操作部１４０の入力手段は、キーボードやマウス、メニュー画面、タッチパネルによるもの等、なんでもよい。また、操作部１４０は、マウス等の入力手段のデバイスドライバや、メニュー画面の制御ソフトウェアで実現される。 【0062】 空間補間部１５０は、複数のドメインの各々について、当該ドメインの補助的な空間データを、任意の空間補間法を用いてターゲットとする粒度の当該補助的な空間データに補間する。 【0063】 本実施形態では、複数のドメインの各々について、空間補間法としてガウス過程回帰を用いて、ターゲットとなる粒度の補助的な空間データとして、ターゲットとする粒度でのエリアごとに当該ドメインの値の平均及び分散を求める場合について説明する。 【0064】 具体的には、空間補間部１５０は、以下に説明するように、補助的空間データ格納部１００に格納されたドメインｓの補助的な空間データ を、様々な粒度を持つ空間データに補間するためのドメインｓにおける空間データに関するハイパーパラメータを求めた後、粒度毎に、各ドメインｓにおける連続表現（予測分布）として、ターゲット粒度エリア中心格納部１１０に格納された各エリア中心での平均（予測値）及び分散を推定する。 【0065】 まず、 をある座標とし、 をドメインｓにおける補助的な空間データに対する潜在的な関数とする。 はガウス過程に従うと仮定し、平均は０、相関関数は下記式（１）とする。 【0066】 【数1】 【0067】 ここで、 はある点 の周りへの点への相関の範囲を決めるスケールパラメータであり、 は相関の大きさ（ｍａｇｎｉｔｕｄｅ）を決める分散パラメータである。式（１）はｓｑｕａｒｅｄ－ｅｘｐｏｎｅｎｔｉａｌ ｋｅｒｎｅｌと呼ばれ、空間座標におけるあるデータの類似度を測るためにも最もよく用いられる相関関数の一つである。 【0068】 ドメインｓのエリア集合における各エリアの中心座標 を表す訓練データが与えられたとしたとき、 の結合分布は多次元ガウス分布で表すことができ、下記式（２）と書ける。 【0069】 【数2】 【0070】 ここで、 は 行列であり、各要素は である。 【0071】 次に、ドメインｓの空間データ は、ガウスノイズが加えられて得られるものと仮定する。このとき、ドメインｓの空間データ は、条件付き確率である下記式（３）に従う。 【0072】 【数3】 【0073】 ここで、 は、単位行列であり、 はノイズに対する分散パラメータである。 【0074】 ここで、上記式（２）及び式（３）の共役性により、 を解析的に積分消去することができる。このため、ドメインｓの空間データ の周辺尤度が、下記式（４）により定義される。 【0075】 【数4】 【0076】 そして、空間補間部１５０は、式（４）について、周辺尤度最大化に基づいて、ハイパーパラメータである を推定する。式（４）の両辺の対数をとった周辺化対数尤度関数は、下記式（５）で表すことができる。 【0077】 【数5】 【0078】 すなわち、空間補間部１５０は、式（５）を最大にするようなハイパーパラメータ を求める。ここで、最適化手法は何を用いてもよいが、例えば、ＢＦＧＳ法（参考文献３）を用いることができる。 ［参考文献３］D.C. Liu and J. Nocedal, “On the Limited Memory BFGS Method for Large Scale Optimization”, Mathematical Programming, vol. 45, 1989, pp. 503-528. 【0079】 次に、空間補間部１５０は、推定したハイパーパラメータ を用いて、粒度毎に、ターゲット粒度エリア中心格納部１１０に格納された各エリアｊの中心座標 での予測分布を計算する。 を、あるターゲット粒度のテストデータとしたとき、各テストサンプル における予測値 の結合予測分布は、下記式（６）となる。 【0080】 【数6】 【0081】 ここで、 は、平均（予測平均）であり、下記式（７）により計算することができ、 は分散（相関行列）であり、下記式（８）により計算することができる。 【0082】 【数7】 【0083】 【数8】 【0084】 ここで、 は、 行列であり、 の各要素は、訓練データ の各訓練サンプル と、テストデータ の各テストサンプル との相関関数（式（１））の値である。ただし、推定したハイパーパラメータ を用いる。 【0085】 また、 は、 行列であり、 の各要素は、テストデータ の各テストサンプル 間の相関関数（式（１））の値である。ただし、推定したハイパーパラメータ を用いる。 【0086】 このように、空間補間部１５０は、上記の処理を全てのドメインｓについて、ハイパーパラメータ を推定し、ターゲットとする粒度毎に、平均 及び分散 を求めることにより、ドメインｓの補助的な空間データを、ターゲットとする粒度の当該補助的な空間データに補間する。 【0087】 そして、空間補間部１５０は、求めた全てのドメインｓにおける平均 及び分散 を、パラメータ推定部１６０及び高解像度データ算出部１７０に渡す。 【0088】 パラメータ推定部１６０は、変換対象の空間データと、複数のドメインの各々についてのターゲットとする粒度でのエリア毎にドメインの値の平均及び分散とに基づいて、ターゲットとする粒度の変換対象の空間データに変換するための回帰モデルのパラメータを推定する。 【0089】 具体的には、パラメータ推定部１６０は、まず、 を変換対象とする空間データに対する潜在的な関数とし、 はガウス過程に従うものとし、平均関数を下記式（９）とする。 【0090】 【数9】 【0091】 ここで、 はドメインｓの空間データに対する重みパラメータであり、 はバイアスパラメータである。また、相関関数は上記式（１）と同様のものを用いる。ただし、相関関数のハイパーパラメータを とする。 【0092】 次に、パラメータ推定部１６０は、ターゲットとする粒度におけるエリア集合の各エリアｊの中心座標 に対する変換対象の空間データのサンプル値を とする。 【0093】 ここで、 は、変換対象の空間データを、ターゲットとする粒度の変換対象の空間データに変換した空間データであり、以下高解像度データ と呼ぶ。 【0094】 各ドメインｓにおける補助的な空間データの予測値が与えられたとしたとき、 の条件付き確率は、下記式（１０）と書ける。 【0095】 【数10】 【0096】 ここで、 であり、 は の相関行列である。また、 を数値１を並べた列ベクトルであるとして、 とした。 【0097】 変換対象の空間データ は、 を構成する高解像度データ を平均化して表されると仮定すると、 は下記式（１１）の条件付き分布から生成されるとみなすことができる。 【0098】 【数11】 【0099】 ここで、 はノイズに対する分散パラメータであり、 は の平均化行列である。 【0100】 また、パラメータ推定部１６０は、粒度エリア対応表格納部１３０に格納されている を用いて、 の要素 は、 を満たす場合は であり、満たさない場合は０とする。 【0101】 次に、パラメータ推定部１６０は、モデルパラメータである重みパラメータ 及び空間相関を制御するパラメータ を推定する。 【0102】 具体的には、変換対象空間データ格納部１２０に格納されている変換対象の空間データ の周辺化尤度は下記式（１２）で表すことができる。 【0103】 【数12】 【0104】 ここで、 とした。また、 は の共分散行列であり、各要素 は下記式（１３）で表すことができる。 【0105】 【数13】 【0106】 ここで、 はクロネッカーのデルタ関数であり、Ａ＝Ｂのとき１の値を、それ以外のとき０の値をとる関数である。また、上記式（１３）の２行目における第２項は、 の回帰残差に対する分散を表す。ここに、 が加えられていることに注目する。各ドメインｓにおける補助的な空間データを用いて、空間補間部１５０により求められた平均 に対する分散 が の回帰残差に対する分散に影響する。これにより、補助的な空間データを用いて推定した予測の信頼度、すなわち、空間補間の信頼度を加味しつつ、重みパラメータ を推定することができる。 【0107】 そして、パラメータ推定部１６０は、上記式（１２）を最大化するように、重みパラメータ 及び空間相関を制御するパラメータ を推定する。 【0108】 ここで、式（１２）の両辺に対数をとった周辺化対数尤度関数は、下記式（１４）と書くことができる。 【0109】 【数14】 【0110】 そこで、パラメータ推定部１６０は、上記式（１４）を最大にするような、重みパラメータ 及び空間相関を制御するパラメータ を推定する。最適化手法は何を用いてもよいが、例えば、ＢＦＧＳ法（参考文献３）を用いることができる。 【0111】 そして、パラメータ推定部１６０は、推定した重みパラメータ 及び空間相関を制御するパラメータ を、高解像度データ算出部１７０に渡す。 【0112】 高解像度データ算出部１７０は、変換対象の空間データと、複数のドメインの各々についてのターゲットとする粒度の補助的な空間データの各々と、回帰モデルとに基づいて、ターゲットとする粒度の変換対象の空間データを算出する。 【0113】 具体的には、高解像度データ算出部１７０は、変換対象空間データ格納部１２０に格納された変換対象とする空間データ 、空間補間部１５０により求められた各ドメインにおける補助的な空間データの平均 及び分散 、及びターゲット粒度エリア中心格納部１１０に格納されたターゲットとする粒度のエリア集合に対する各エリアｊの中心座標の集合 と、粒度エリア対応表格納部１３０に格納された対応表 と、パラメータ推定部１６０で推定された重みパラメータ 及び空間相関を制御するパラメータ とを用いて、ターゲットとする粒度の変換対象の空間データである高解像度データ の予測分布 を求めることにより、高解像度データ を算出する。 【0114】 予測分布 は、下記式（１５）である。 【0115】 【数15】 【0116】 ここで、 は平均であり、下記式（１６）で表すことができる。 【0117】 【数16】 【0118】 また、 は相関行列であり、下記式（１７）で表すことができる。 【0119】 【数17】 【0120】 そして、高解像度データ算出部１７０は、算出したターゲットとする粒度の変換対象の空間データである高解像度データ の平均 及び相関行列 を、高解像度データ格納部１８０に格納する。 【0121】 高解像度データ格納部１８０は、高解像度データ算出部１７０により算出されたターゲットとする粒度の変換対象の空間データである高解像度データ の平均 及び相関行列 を、格納する。 【0122】 そして、高解像度化装置１０は、空間補間部１５０、パラメータ推定部１６０、及び高解像度データ算出部１７０の処理を、ターゲット粒度候補の全てについて行う。 【0123】 入力部１９０は、ターゲット粒度候補のＩＤであるターゲット粒度ＩＤの入力を受け付ける。 【0124】 そして、入力部１９０は、受け付けたターゲット粒度ＩＤを、検索部２００に渡す。 【0125】 なお、入力部１９０にターゲット粒度ＩＤを入力する手段は、キーボードやマウス、メニュー画面、タッチパネルによるもの等、なんでもよい。また、入力部１９０は、マウス等の入力手段のデバイスドライバや、メニュー画面の制御ソフトウェアで実現され得る。 【0126】 検索部２００は、入力部１９０により受け付けたターゲット粒度ＩＤにより指定された粒度のエリア集合に対して、変換対象の空間データを高解像度化した高解像度データ を検索する。 【0127】 具体的には、図２上図に示すように、入力部１９０が受け付けたターゲット粒度ＩＤと同一のＩＤであるターゲット粒度候補の変換対象の空間データである高解像度データ の平均 を、高解像度データ格納部１８０から検索する。 【0128】 そして、検索部２００は、検索結果である を出力部２１０に渡す。 【0129】 出力部２１０は、検索部２００から取得した を高解像度データとして出力する。 【0130】 例えば、図２下図に示すように、粗い粒度の変換対象の空間データを、ターゲット粒度候補Ａの持つ粒度に高解像度データを出力する。 【0131】 また、上記ではディスプレイへの表示を例にしたが、出力部２１０は、プリンタへの印字、音出力、外部装置への送信等の出力をする構成としてもよい。また、出力部２１０は、ディスプレイやスピーカ等の出力デバイスを含む構成としてもよい。また、出力部２１０は、出力デバイスのドライバソフトまたは、出力デバイスのドライバソフトと出力デバイス等で実現され得る。 【0132】 ＜本発明の実施の形態に係る高解像度化装置の作用＞ 図３は、本発明の実施の形態に係る高解像度化処理ルーチンを示すフローチャートである。 【0133】 高解像度化装置１０において、高解像度化処理を行う命令が実行されると、図３に示す高解像度化処理ルーチンが実行される。 【0134】 まず、ステップＳ１００において、空間補間部１５０は、複数のドメインの各々について、当該ドメインの補助的な空間データを様々な粒度の当該補助的な空間データに補間するためのハイパーパラメータを推定する。 ステップＳ１１０において、高解像度化装置１０が、１番目のターゲット粒度ＩＤを選択する。 【0135】 ステップＳ１２０において、空間補間部１５０は、複数のドメインの各々について、当該ドメインの補助的な空間データについて上記ステップＳ１００で推定されたハイパーパラメータを用いて、当該ドメインの補助的な空間データを、任意の空間補間法を用いてターゲットとする粒度の当該補助的な空間データに補間する。 【0136】 ステップＳ１３０において、パラメータ推定部１６０は、変換対象の空間データと、複数のドメインの各々についてのターゲットとする粒度でのエリア毎にドメインの値の平均及び分散とに基づいて、ターゲットとする粒度の変換対象の空間データに変換するための回帰モデルのパラメータを推定する。 【0137】 ステップＳ１４０において、高解像度データ算出部１７０は、変換対象の空間データと、複数のドメインの各々についてのターゲットとする粒度の補助的な空間データの各々と、回帰モデルとに基づいて、ターゲットとする粒度の変換対象の空間データを算出する。 【0138】 ステップＳ１５０において、高解像度データ格納部１８０は、上記ステップＳ１４０により算出されたターゲットとする粒度の変換対象の空間データである高解像度データの平均及び相関行列を、格納する。 【0139】 ステップＳ１６０において、高解像度化装置１０は、全てのターゲット粒度ＩＤについて、処理を行ったか否かを判定する。 【0140】 全てのターゲット粒度ＩＤについて処理を行っていない場合（ステップＳ１６０のＮＯ）、ステップＳ１７０において、高解像度化装置１０は、次のターゲット粒度ＩＤを選択し、上記ステップＳ１２０～Ｓ１６０の処理を繰り返す。 【0141】 一方、全てのターゲット粒度ＩＤについて処理を行っている場合（ステップＳ１６０のＹＥＳ）、処理を終了する。 【0142】 また、図４は、本発明の実施の形態に係る検索処理ルーチンを示すフローチャートである。 【0143】 ステップＳ２００において、ターゲットとする粒度のＩＤであるターゲット粒度ＩＤの入力を受け付ける。 【0144】 ステップＳ２１０において、検索部２００は、上記ステップＳ２００により受け付けたターゲット粒度ＩＤにより指定された粒度のエリア集合に対して、変換対象の空間データを高解像度化した高解像度データを検索する。 【0145】 ステップＳ２２０において、出力部２１０は、検索部２００から取得した高解像度データを出力する。 【0146】 以上説明したように、本発明の実施形態に係る高解像度化装置によれば、値の種類を示す複数のドメインの各々について、当該ドメインの補助的な空間データを、任意の空間補間法を用いてターゲットとする粒度の当該補助的な空間データに補間し、変換対象の空間データと、複数のドメインの各々についてのターゲットとする当該粒度の当該補助的な空間データの各々とに基づいて、当該変換対象の空間データをターゲットとする当該粒度の当該変換対象の空間データに変換するための回帰モデルのパラメータを推定するため、異なる粒度を持つ補助的な空間データを用いて、空間データを精度よく高解像度化することができる。 【0147】 なお、本発明は、上述した実施の形態に限定されるものではなく、この発明の要旨を逸脱しない範囲内で様々な変形や応用が可能である。 【0148】 例えば、上記実施形態では、予め全ての粒度について高解像度データ の平均 及び相関行列 を、高解像度データ格納部１８０に格納しておく構成としたが、ターゲット粒度ＩＤが入力される度に学習済みモデルで高解像度化して出力する構成としてもよい。 【0149】 また、入力部１９０が、変換対象となる空間データ 及びターゲット粒度ＩＤを受け付け、これを変換対象空間データ格納部１２０に格納し、入力されたターゲット粒度ＩＤの粒度をターゲットとする粒度として、変換対象となる空間データを高解像度化し、これを出力する構成としてもよい。 【0150】 また、本願明細書中において、プログラムが予めインストールされている実施形態として説明したが、当該プログラムを、コンピュータ読み取り可能な記録媒体に格納して提供することも可能である。 【符号の説明】 【0151】 １０ 高解像度化装置 １００ 補助的空間データ格納部 １１０ ターゲット粒度エリア中心格納部 １２０ 変換対象空間データ格納部 １３０ 粒度エリア対応表格納部 １４０ 操作部 １５０ 空間補間部 １６０ パラメータ推定部 １７０ 高解像度データ算出部 １８０ 高解像度データ格納部 １９０ 入力部 ２００ 検索部 ２１０ 出力部
図面	【図1】 【図2】 【図3】 【図4】