(57)【要約】【課題】複数の方位をそれぞれ向いた各面に対する着雪量が方位毎に異なるような状況に応じて、着雪量を的確に予測することができる、着雪予測装置を提供する。【解決手段】着雪予測装置２は、所定の時間間隔で区切られた各時間Ｔ（ｔ）及び所定の大きさで区切られた各領域Ｒ（ｒ）に対する複数の気象要素１１０を含む気象予報データ１００を取得するデータ取得部２２０と、データ取得部２２０により取得された気象予報データ１００に基づいて、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）における着雪情報２００（着雪速度Ｍｓ、着雪重量Ｍｗ又は着雪高さＭｈ）を、複数の方位Ｄ（ｄ）をそれぞれ向いた複数の着雪面Ｓ（ｄ）の方位毎に算定する着雪情報算定部２２１と、を備える。【選択図】      図１
2020067343.tif 000002
【特許請求の範囲】
【請求項１】
所定の時間間隔で区切られた各時間及び所定の大きさで区切られた各領域に対する複数の気象要素を含む気象予報データを取得するデータ取得部と、
前記データ取得部により取得された前記気象予報データに基づいて、前記各時間及び前記各領域における着雪情報を、複数の方位をそれぞれ向いた複数の着雪面の前記方位毎に算定する着雪情報算定部と、を備える、
ことを特徴とする着雪予測装置。
【請求項２】
前記データ取得部は、前記複数の気象要素として、少なくとも気温、降水量、風向及び風速を含む前記気象予報データを取得し、
前記着雪情報算定部は、前記データ取得部により取得された前記気象予報データと、前記複数の方位及び前記風向の間の関係により規定される方位毎の着雪率とに基づいて、前記各時間及び前記各領域における前記着雪情報を前記方位毎に算定する、
ことを特徴とする請求項１に記載の着雪予測装置。
【請求項３】
前記着雪情報算定部は、
前記気象予報データに基づいて、前記各時間及び前記各領域における単位時間当たりの衝突重量を計算する衝突重量計算部と、
前記複数の方位及び前記風向の間の関係に基づいて、前記各時間及び前記各領域における前記着雪率を前記方位毎に計算する着雪率計算部と、
前記衝突重量計算部により計算された前記衝突重量と、前記着雪率計算部により計算された前記方位毎の前記着雪率とに基づいて、前記各時間及び前記各領域における単位時間当たりの着雪重量を示す着雪速度を前記着雪情報として前記方位毎に計算する着雪速度計算部と、
前記着雪速度計算部により計算された前記方位毎の前記着雪速度を所定の単位時間分、前記方位毎に累積することにより、前記所定の単位時間後の前記各領域における着雪重量又は着雪高さを前記着雪情報として前記方位毎に計算する着雪情報累積部と、を備える、
ことを特徴とする請求項２に記載の着雪予測装置。
【請求項４】
前記着雪情報算定部は、
前記気象予報データに基づいて、前記各時間及び前記各領域において所定の落雪条件を満たすか否かを判定する落雪判定部と、
前記落雪判定部により前記所定の落雪条件を満たすと判定された前記時間及び前記領域における前記着雪情報をリセットする着雪情報リセット部と、をさらに備える、
ことを特徴とする請求項１乃至請求項３のいずれか１項に記載の着雪予測装置。
【請求項５】
前記着雪情報算定部は、
前記気象予報データに基づいて、前記各時間及び前記各領域において所定の落雪条件を満たすか否かを前記方位毎に判定する落雪判定部と、
前記落雪判定部により前記所定の落雪条件を満たすと判定された前記時間及び前記領域における前記方位の前記着雪情報をリセットする着雪情報リセット部と、をさらに備える、
ことを特徴とする請求項１乃至請求項３のいずれか１項に記載の着雪予測装置。
【請求項６】
前記落雪判定部は、前記気象要素が前記複数の方位にそれぞれ及ぼす落雪への影響を考慮することにより、前記所定の落雪条件を満たすか否かを前記方位毎に判定する、
ことを特徴とする請求項５に記載の着雪予測装置。
【請求項７】
前記着雪情報算定部により算定された前記方位毎の前記着雪情報を比較可能に表示画面に表示させる表示情報を生成する表示情報生成部をさらに備える、
ことを特徴とする請求項１乃至請求項６のいずれか１項に記載の着雪予測装置。
【請求項８】
前記表示画面は、前記着雪情報算定部により算定された前記方位毎の前記着雪情報を比較可能に表示する複数方位比較表示領域を備え、
前記複数方位比較表示領域は、前記複数方位比較表示領域の中心部を基準として前記方位の各々に対応する各位置に配置された複数の着雪情報表示部を備え、前記複数の着雪情報表示部は、前記各位置に対応する前記方位の着雪情報をそれぞれ表示する、
ことを特徴とする請求項７に記載の着雪予測装置。
【請求項９】
コンピュータを、請求項１乃至請求項８のいずれか１項に記載の着雪予測装置として機能させるための着雪予測プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、着雪予測装置及び着雪予測プログラムに関する。
【背景技術】
【０００２】
強風湿型着雪は、日本の太平洋側を通過する南岸低気圧が湿った雪と強風をもたらすことで、多量の湿雪が、道路設備（交通信号機、道路標識、高架道路の側壁等）、鉄道設備（鉄道信号機、鉄道標識、踏切、高架線路の側壁等）、電力設備（電柱、鉄塔等）、建築物（屋根、壁等）及び看板等の屋外構造物並びに樹木に着雪し、道路障害及び鉄道障害、電力設備、建築物及び看板の損壊並びに倒木等の被害を引き起こすものである。
【０００３】
このような強風湿型着雪による被害を受ける地域は、四国・九州から関東・東海、東北・北海道まで広範囲に及ぶため、着雪による被害を未然に防ぐことを目的として、各地域の着雪量を予測する予測手法の開発が求められている。
【０００４】
例えば、特許文献１には、架線の着雪予測を行う対象領域を所定の大きさで区切った各メッシュについて、着雪量に比例する着雪ポテンシャルを、架線の走行方位毎に求める着雪予測装置が開示されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００８−１５７７４５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
屋外構造物や樹木に着雪する場合、着雪の対象となる屋外構造物や樹木の面がそれぞれ向いている各方位に応じて、着雪量は異なるものであり、例えば、降雪時に北風が吹いているような場合、北向きの面に対する着雪量は、南向きの面に対する着雪量よりも多くなるものである。
【０００７】
特に、信号機や標識は、それらが設置された場所の状況に応じて信号機や標識の表示面が様々な方位を向くように設置されているので、信号機や標識の表示面が着雪により視認困難となって、交通障害や鉄道障害が発生する可能性を判断するためには、着雪の対象となる信号機や標識の表示面がそれぞれ向いている各方位に対して着雪量をそれぞれ予測することが必要になる。
【０００８】
しかし、特許文献１に開示された着雪予測装置は、架線を着雪の対象として、架線への着雪量に比例する着雪ポテンシャルを、架線の走行方位毎に予測するものである。したがって、特許文献１に開示された着雪予測装置は、例えば、南北方向に架かる架線に対する着雪ポテンシャルや、東西方向に架かる架線に対する着雪ポテンシャルをそれぞれ予測するものではあるが、北向きの面に対する着雪量や、南向きの面に対する着雪量をそれぞれ予測するものではない。そのため、特許文献１に開示された着雪予測装置は、複数の方位をそれぞれ向いた各面に対する着雪量が方位毎に異なるような状況に応じて、着雪量を的確に予測することができない、という問題点があった。
【０００９】
本発明は、このような事情に鑑みてなされたものであって、複数の方位をそれぞれ向いた各面に対する着雪量が方位毎に異なるような状況に応じて、着雪量を的確に予測することができる、着雪予測装置及び着雪予測プログラムを提供することを目的とする。
【課題を解決するための手段】
【００１０】
本発明は、上記課題を解決するものであって、本発明の一実施形態に係る着雪予測装置は、
所定の時間間隔で区切られた各時間及び所定の大きさで区切られた各領域に対する複数の気象要素を含む気象予報データを取得するデータ取得部と、
前記データ取得部により取得された前記気象予報データに基づいて、前記各時間及び前記各領域における着雪情報を、複数の方位をそれぞれ向いた複数の着雪面の前記方位毎に算定する着雪情報算定部と、を備える、ことを特徴とする。
【００１１】
また、本発明の一実施形態に係る着雪予測装置は、
前記データ取得部は、前記複数の気象要素として、少なくとも気温、降水量、風向及び風速を含む前記気象予報データを取得し、
前記着雪情報算定部は、前記データ取得部により取得された前記気象予報データと、前記複数の方位及び前記風向の間の関係により規定される方位毎の着雪率とに基づいて、前記各時間及び前記各領域における前記着雪情報を前記方位毎に算定する、ことを特徴とする。
【００１２】
また、本発明の一実施形態に係る着雪予測装置は、
前記着雪情報算定部は、
前記気象予報データに基づいて、前記各時間及び前記各領域における単位時間当たりの衝突重量を計算する衝突重量計算部と、
前記複数の方位及び前記風向の間の関係に基づいて、前記各時間及び前記各領域における前記着雪率を前記方位毎に計算する着雪率計算部と、
前記衝突重量計算部により計算された前記衝突重量と、前記着雪率計算部により計算された前記方位毎の前記着雪率とに基づいて、前記各時間及び前記各領域における単位時間当たりの着雪重量を示す着雪速度を前記着雪情報として前記方位毎に計算する着雪速度計算部と、
前記着雪速度計算部により計算された前記方位毎の前記着雪速度を所定の単位時間分、前記方位毎に累積することにより、前記所定の単位時間後の前記各領域における着雪重量又は着雪高さを前記着雪情報として前記方位毎に計算する着雪情報累積部と、を備える、
ことを特徴とする。
【００１３】
また、本発明の一実施形態に係る着雪予測装置は、
前記着雪情報算定部は、
前記気象予報データに基づいて、前記各時間及び前記各領域において所定の落雪条件を満たすか否かを判定する落雪判定部と、
前記落雪判定部により所定の落雪条件を満たすと判定された前記時間及び前記領域における前記着雪情報をリセットする着雪情報リセット部と、をさらに備える、ことを特徴とする。
【００１４】
また、本発明の一実施形態に係る着雪予測装置は、
前記着雪情報算定部は、
前記気象予報データに基づいて、前記各時間及び前記各領域において所定の落雪条件を満たすか否かを前記方位毎に判定する落雪判定部と、
前記落雪判定部により所定の落雪条件を満たすと判定された前記時間及び前記領域における前記方位の前記着雪情報をリセットする着雪情報リセット部と、をさらに備える、ことを特徴とする。
【００１５】
また、本発明の一実施形態に係る着雪予測装置は、
前記落雪判定部は、前記気象要素が前記複数の方位にそれぞれ及ぼす落雪への影響を考慮することにより、前記所定の落雪条件を満たすか否かを前記方位毎に判定する、ことを特徴とする。
【００１６】
また、本発明の一実施形態に係る着雪予測装置は、
前記着雪情報算定部により算定された前記方位毎の前記着雪情報を比較可能に表示画面に表示させる表示情報を生成する表示情報生成部をさらに備える、ことを特徴とする。
【００１７】
また、本発明の一実施形態に係る着雪予測装置は、
前記表示画面は、前記着雪情報算定部により算定された前記方位毎の前記着雪情報を比較可能に表示する複数方位比較表示領域を備え、
前記複数方位比較表示領域は、前記複数方位比較表示領域の中心部を基準として前記方位の各々に対応する各位置に配置された複数の着雪情報表示部を備え、前記複数の着雪情報表示部は、前記各位置に対応する前記方位の着雪情報をそれぞれ表示する、ことを特徴とする。
【００１８】
また、本発明の一実施形態に係る着雪予測プログラムは、
コンピュータを、上記着雪予測装置として機能させることを特徴とする。
【発明の効果】
【００１９】
本発明の一実施形態に係る着雪予測装置及び着雪予測システムによれば、着雪情報算定部が、気象予報データに基づいて、各時間及び各領域における着雪情報を、複数の方位をそれぞれ向いた複数の着雪面の方位毎に算定するので、各方位を向いた複数の着雪面に対する着雪量が方位毎に異なるような状況に応じて、着雪量を的確に予測することができる。
【図面の簡単な説明】
【００２０】
【図１】本発明の実施形態に係る着雪予測システム１の一例を示す全体構成図である。
【図２】本発明の実施形態に係る着雪予測装置２を示すブロック図である。
【図３】本発明の実施形態に係る着雪情報算定部２２１を示すブロック図である。
【図４】本発明の実施形態に係る着雪予測装置２の動作（着雪予測処理）を示すフローチャートである。
【図５】本発明の実施形態に係る着雪予測装置２の動作（ステップＳ２００の風向計算処理の詳細）を示すフローチャートである。
【図６】本発明の実施形態に係る着雪予測装置２の動作（ステップＳ３００の着雪情報算定処理の詳細）を示すフローチャートである。
【図７】本発明の実施形態に係る着雪予測装置２の動作（ステップＳ４００の方位別算定処理の詳細）を示すフローチャートである。
【図８】本発明の実施形態に係る着雪予測装置２の動作（ステップＳ５００の落雪判定処理の詳細）を示すフローチャートである。
【図９】本発明の実施形態に係る着雪予測装置２の動作（表示情報生成処理）及びユーザ端末１１の動作を示すフローチャートである。
【図１０】本発明の実施形態に係る着雪予測システム１における初期表示画面３Ａを示す図である。
【図１１】本発明の実施形態に係る着雪予測システム１における選択表示画面３Ｂを示す図である。
【図１２】本発明の実施形態に係る着雪予測システム１における時系列グラフ表示画面４を示す図である。
【図１３】本発明の実施形態に係る着雪予測システム１における第１の複数方位時系列グラフ表示画面５Ａを示す図である。
【図１４】本発明の実施形態に係る着雪予測システム１における第２の複数方位時系列グラフ表示画面５Ｂを示す図である。
【図１５】本発明の実施形態に係る着雪予測装置２の動作の変形例（ステップＳ５００の落雪判定処理の変形例）を示すフローチャートである。
【図１６】本発明の実施形態に係る着雪予測システム１における第３の複数方位時系列グラフ表示画面５Ｃを示す図である。
【発明を実施するための形態】
【００２１】
以下、本発明の一実施形態について添付図面を参照しつつ説明する。
【００２２】
図１は、本発明の実施形態に係る着雪予測システム１の一例を示す全体構成図である。
【００２３】
着雪予測システム１は、気象予報データ１００を提供する気象予報システム１０と、気象予報システム１０により提供された気象予報データ１００を取得し、気象予報データ１００に基づいて着雪情報２００を予測する着雪予測装置２と、各ユーザが使用し、着雪予測装置２により予測された着雪情報２００を各ユーザに出力するユーザ端末１１と、各装置間を相互に通信可能に接続するネットワーク１２と、を備える。なお、ネットワーク１２は、有線でも無線でもよく、また、インターネット回線でも専用回線でもよい。
【００２４】
気象予報システム１０は、気象庁、外国の気象機関、国内外の気象会社等により運営されるシステムである。気象予報システム１０は、複数の観測地点で収集された現在又は過去の観測データに基づいて気象状況の推移を予測することにより、気象予報データ１００を生成し、その気象予報データ１００を所定の提供時刻（例えば、１日に８回、３時間毎の時刻）に外部に提供する。
【００２５】
気象予報データ１００は、ＧＰＶ（Grid Point Value：格子点値）データと呼ばれるものであり、各種の数値予報モデルを用いて、所定の時間間隔（例えば、３０分、１時間、３時間、６時間）で区切られた各時間Ｔ（ｔ）及び所定の大きさ（例えば、２ｋｍメッシュ、５ｋｍメッシュ、２０ｋｍメッシュ）で区切られた各領域Ｒ（ｒ）に対する複数の気象要素１１０を予測したデータである。数値予報モデルとしては、例えば、全球数値用法モデル（ＧＳＭ：Global Spectral Model）、メソ数値予報モデル（ＭＳＭ：Meso Scale Model）、局地数値予報モデル（ＬＦＭ：Local Forecast Model）等が用いられる。
【００２６】
複数の気象要素１１０は、気象状況を定量的に表現するための各種の要素であり、例えば、気温［℃］、相対湿度［％］、降水量［ｍｍ／ｈ］、風ベクトル（東西成分及び南北成分の２要素）［ｍ／ｓ］、日射量［Ｗ／ｍ２］、気圧［ｈＰａ］及び雲量［％］等である。
【００２７】
なお、本実施形態では、気象予報データ１００は、メソ数値予報モデル（ＭＳＭ）を用いて予測したＧＰＶデータであるものとし、現時点から１時間間隔で区切られた３９時間先までの各時間Ｔ（ｔ）、及び、５ｋｍメッシュで区切られた各領域Ｒ（ｒ）に対する複数の気象要素１１０として、気温Ｔｅｍｐ、相対湿度Ｈｒ、降水量Ｐｒ、風ベクトルＷ（東西成分ｗｕ及び南北成分ｗｖ）、及び、日射量Ｓｒを含むデータであるものとして説明する。そのため、変数ｔは、１〜３９（＝Ｔｍａｘ）の間の整数値であり、変数ｒは、１〜各領域に通し番号を割り振った場合の最大値（＝Ｒｍａｘ）の間の整数値である。
【００２８】
ユーザ端末１１は、汎用のデスクトップ型コンピュータ等の情報処理装置や、ユーザが携帯可能なノートブック型コンピュータ、タブレット、スマートフォン、携帯電話、ウェアラブル端末等の携帯端末機器により構成されている。
【００２９】
ユーザ端末１１は、表示画面に対するユーザの入力操作を受け付けて、表示情報生成要求を着雪予測装置２に送信するとともに、その表示情報生成要求に対する応答として、着雪予測装置２から表示情報を受信し、その表示情報に基づいて、表示画面に気象予報データ１００や着雪情報２００を表示する。
【００３０】
着雪予測装置２は、気象予報データ１００に基づいて、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）における着雪情報２００を、複数の方位Ｄ（ｄ）をそれぞれ向いた複数の着雪面Ｓ（ｄ）の方位毎に予測する装置である。
【００３１】
着雪情報２００は、雪が着雪面Ｓ（ｄ）に着雪する際の着雪量を示す情報であり、例えば、単位時間当たりの着雪重量（単位面積当たりに着雪した雪の重量）を示す着雪速度Ｍｓ［ｋｇ／ｍ２／ｈ］、単位面積当たりに着雪した雪の重量を累積値として示す着雪重量Ｍｗ［ｋｇ／ｍ２］、又は、着雪した雪の高さを累積値として示す着雪高さＭｈ［ｃｍ］で表される情報である。
【００３２】
複数の方位Ｄ（ｄ）は、互いに異なる方位であり、全方位を、例えば、４等分、８等分、１６等分、３２等分等で均等に分けた場合には、４方位、８方位、１６方位、３２方位等の各方位を示すものである。また、複数の方位は、ユーザにより任意の方位が指定されたものでもよい。
【００３３】
着雪面Ｓ（ｄ）は、雪が着雪する対象物を想定した仮想的な面であり、対象物が有する鉛直面、又は、水平面に対して傾斜した傾斜面を想定したものである。雪が着雪する対象物は、例えば、道路設備（交通信号機、道路標識、高架道路の側壁等）、鉄道設備（鉄道信号機、鉄道標識、踏切、高架線路の側壁等）、電力設備（電柱、鉄塔等）、建築物（屋根、壁等）及び看板等の屋外構造物並びに樹木である。
【００３４】
なお、本実施形態では、複数の方位Ｄ（ｄ）は、８方位とし、８つの着雪面Ｓ（ｄ）は、８方位の各方位Ｄ（ｄ）をそれぞれ向いたものとして説明する。そのため、変数ｄは、１から８までの整数値であり、方位Ｄ（１）＝「北」、方位Ｄ（２）＝「北東」、方位Ｄ（３）＝「東」、方位Ｄ（４）＝「南東」、方位Ｄ（５）＝「南」、方位Ｄ（６）＝「南西」、方位Ｄ（７）＝「西」、方位Ｄ（８）＝「北西」である。したがって、方位Ｄ（１）＝「北」を向いた着雪面Ｓ（１）に対する着雪情報２００は、例えば、表示面が「北」又は「北」と近似できるような方向を向いて設置された実際の信号機や標識において、その表示面にどの程度の雪が着雪するのかを予測した情報である。
【００３５】
図２は、本発明の実施形態に係る着雪予測装置２を示すブロック図である。着雪予測装置２は、例えば、汎用のデスクトップ型コンピュータやサーバ等の情報処理装置により構成されている。
【００３６】
着雪予測装置２は、キーボード、タッチパネル等により構成される入力部２０と、ＨＤＤ、メモリ等により構成される記憶部２１と、ＣＰＵ等のプロセッサにより構成される制御部２２と、ネットワーク１２との通信インターフェースである通信部２３と、ディスプレイ、スピーカ等により構成される出力部２４と、を備える。なお、入力部２０や出力部２４は省略してもよく、その場合には、別の情報処理装置が、着雪予測装置２の入力部２０や出力部２４として機能するように構成すればよい。
【００３７】
記憶部２１には、気象予報データ１００を蓄積する気象予報データベース２１０と、着雪情報算定部２２１（詳細は後述する）により算定された着雪情報２００を蓄積する着雪情報データベース２１１と、着雪予測プログラム２１２と、が記憶されている。
【００３８】
制御部２２は、記憶部２１に記憶された着雪予測プログラム２１２を実行することにより、データ取得部２２０、着雪情報算定部２２１、及び、表示情報生成部２２２として機能し、着雪予測処理（詳細は図４〜８にて後述する）、及び、表示情報生成処理（詳細は図９〜１４にて後述する）を行う。
【００３９】
データ取得部２２０は、気象予報システム１０から通信部２３を介して気象予報データ１００を取得する。また、データ取得部２２０は、取得した気象予報データ１００を気象予報データベース２１０に登録する。なお、データ取得部２２０は、気象予報データ１００が記憶された、例えば、ＵＳＢメモリ等の外部記憶装置に接続されて、外部記憶装置から気象予報データ１００を取得してもよいし、着雪予測装置２の管理者により入力部２０を介して入力された気象予報データ１００を取得してもよい。
【００４０】
着雪情報算定部２２１は、データ取得部２２０により取得された気象予報データ１００に基づいて、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）における着雪情報２００として、着雪速度Ｍｓ（ｔ，ｒ，ｄ）、着雪重量Ｍｗ（ｔ，ｒ，ｄ）及び着雪高さＭｈ（ｔ，ｒ，ｄ）を、複数の方位Ｄ（ｄ）をそれぞれ向いた複数の着雪面Ｓ（ｄ）の方位毎に算定する。また、着雪情報算定部２２１は、算定した着雪情報２００を着雪情報データベース２１１に登録する。
【００４１】
本実施形態では、着雪情報算定部２２１は、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）における着雪情報２００として、着雪速度Ｍｓ（ｔ，ｒ，ｄ）、着雪重量Ｍｗ（ｔ，ｒ，ｄ）及び着雪高さＭｈ（ｔ，ｒ，ｄ）を、８方位Ｄ（ｄ）をそれぞれ向いた８つの着雪面Ｓ（ｄ）の方位毎、すなわち、８方位毎に算定するものとして説明する。
【００４２】
図３は、本発明の実施形態に係る着雪情報算定部２２１を示すブロック図である。
【００４３】
着雪情報算定部２２１は、気象予報データ１００に基づいて、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）における風向Ｗｄ（ｔ、ｒ）及び風速Ｗｓ（ｔ、ｒ）を計算する風向風速計算部２２１ａと、気象予報データ１００に基づいて、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）において所定の湿雪降雪条件を満たすか否かを判定する湿雪降雪判定部２２１ｂと、気象予報データ１００に基づいて、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）における単位時間当たりの衝突重量Ｐｎ（ｔ，ｒ）［ｋｇ／ｈｒ］を計算する衝突重量計算部２２１ｃと、気象予報データ１００に基づいて、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）における着雪密度Ｐｓ（ｔ，ｒ）［ｋｇ／ｍ３］を計算する着雪密度計算部２２１ｄと、複数の方位Ｄ（ｄ）及び風向Ｗｄ（ｔ、ｒ）の間の関係に基づいて、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）における着雪率α（ｔ，ｒ，ｄ）［％］を方位毎に計算する着雪率計算部２２１ｅと、衝突重量計算部２２１ｃにより計算された衝突重量Ｐｎ（ｔ，ｒ）と着雪率計算部２２１ｅにより計算された方位毎の着雪率α（ｔ，ｒ，ｄ）とに基づいて、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）における単位時間当たりの着雪重量を示す着雪速度Ｍｓ（ｔ，ｒ，ｄ）を着雪情報２００として方位毎に計算する着雪速度計算部２２１ｆと、着雪速度計算部２２１ｆにより計算された方位毎の着雪速度Ｍｓ（ｔ，ｒ，ｄ）を所定の単位時間分、方位毎に累積することにより、所定の単位時間後の各領域Ｒ（ｒ）における着雪重量Ｍｗ（ｔ，ｒ，ｄ）及び着雪高さＭｈ（ｔ，ｒ，ｄ）を着雪情報２００として方位毎に計算する着雪情報累積部２２１ｇと、を備える。
【００４４】
また、着雪情報算定部２２１は、気象予報データ１００に基づいて、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）において所定の落雪条件を満たすか否かを判定する落雪判定部２２１ｈと、落雪判定部２２１ｈにより所定の落雪条件を満たすと判定された時間Ｔ（ｔ）及び領域Ｒ（ｒ）における着雪情報２００をリセットする着雪情報リセット部２２１ｉと、をさらに備える。
【００４５】
風向風速計算部２２１ａは、気象予報データ１００に含まれる時間Ｔ（ｔ）及び領域Ｒ（ｒ）に対する風ベクトルＷ（東西成分ｗｕ及び南北成分ｗｖ）に基づいて、風向Ｗｄ（ｔ、ｒ）及び風速Ｗｓ（ｔ、ｒ）を計算する。
【００４６】
湿雪降雪判定部２２１ｂは、気象予報データ１００に含まれる時間Ｔ（ｔ）及び領域Ｒ（ｒ）に対する気温Ｔｅｍｐ、相対湿度Ｈｒ及び降水量Ｐｒに基づいて、図６に示すフローチャートのステップＳ３０１〜Ｓ３０３（詳細は後述する）により、所定の湿雪降雪条件を満たすか否かを判定する。湿雪降雪判定部２２１ｂが、所定の湿雪降雪条件を満たすと判定した場合には、衝突重量計算部２２１ｃ及び着雪密度計算部２２１ｄによる処理が実行され、所定の湿雪降雪条件を満たさないと判定した場合には、落雪判定部２２１ｈによる処理が実行される。
【００４７】
衝突重量計算部２２１ｃは、風向風速計算部２２１ａにより計算された風速Ｗｓ（ｔ、ｒ）と、気象予報データ１００に含まれる時間Ｔ（ｔ）及び領域Ｒ（ｒ）に対する降水量Ｐｒとに基づいて、衝突重量Ｐｎ（ｔ，ｒ）を計算する。
【００４８】
着雪密度計算部２２１ｄは、風向風速計算部２２１ａにより計算された風速Ｗｓ（ｔ、ｒ）に基づいて、着雪密度Ｐｓ（ｔ，ｒ）を計算する。
【００４９】
着雪率計算部２２１ｅは、複数の方位Ｄ（ｄ）と、気象予報データ１００に含まれる時間Ｔ（ｔ）及び領域Ｒ（ｒ）に対する風向Ｗｄ（ｔ、ｒ）との間の関係に基づいて、図７に示すフローチャートのステップＳ４２０〜Ｓ４２７（詳細は後述する）により、着雪率α（ｔ，ｒ，ｄ）を方位毎に計算する。
【００５０】
また、着雪率計算部２２１ｅは、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）における衝突率β（ｔ，ｒ，ｄ）を方位毎に計算する。なお、本実施形態では、衝突率β（ｔ，ｒ，ｄ）は、「６０」［％］の定数であるものとして説明する。
【００５１】
着雪速度計算部２２１ｆは、衝突重量計算部２２１ｃにより計算された衝突重量Ｐｎ（ｔ，ｒ）と、着雪率計算部２２１ｅにより計算された方位毎の着雪率α（ｔ，ｒ，ｄ）及び衝突率β（ｔ，ｒ，ｄ）とに基づいて、時間Ｔ（ｔ）及び領域Ｒ（ｒ）における着雪速度Ｍｓ（ｔ，ｒ，ｄ）を方位毎に計算する。
【００５２】
着雪情報累積部２２１ｇは、時間Ｔ（ｔ）及び領域Ｒ（ｒ）における着雪重量Ｍｗ（ｔ，ｒ，ｄ）を方位毎に計算する着雪重量累積部２２１ｇ１と、時間Ｔ（ｔ）及び領域Ｒ（ｒ）における着雪高さＭｈ（ｔ，ｒ，ｄ）を方位毎に計算する着雪高さ累積部２２１ｇ２とを備える。
【００５３】
落雪判定部２２１ｈは、気象予報データ１００に含まれる時間Ｔ（ｔ）及び領域Ｒ（ｒ）に対する気温Ｔｅｍｐ、降水量Ｐｒ及び日射量Ｓｒに基づいて、図８に示すフローチャートのステップＳ５０１〜Ｓ５２５（詳細は後述する）により、所定の落雪条件を満たすか否かを判定する。
【００５４】
着雪情報リセット部２２１ｉは、落雪判定部２２１ｈにより所定の落雪条件を満たすと判定された時間Ｔ（ｔ）及び領域Ｒ（ｒ）における着雪情報２００として、着雪重量Ｍｗ（ｔ，ｒ，ｄ）及び着雪高さＭｈ（ｔ，ｒ，ｄ）を「０」にリセットする。
【００５５】
なお、本実施形態では、着雪情報リセット部２２１ｉは、着雪重量Ｍｗ（ｔ，ｒ，ｄ）及び着雪高さＭｈ（ｔ，ｒ，ｄ）をリセットする際に、着雪面Ｓ（ｄ）に着雪した雪の全てが落雪することを想定して、「０」にリセットするものとして説明するが、着雪面Ｓ（ｄ）に着雪した雪の一部が落雪することを想定した場合には、一定の割合（例えば、５０％）又は一定の値（例えば、５ｋｇ／ｍ２及び１０ｃｍ）を減算することによりリセットしてもよいし、着雪面Ｓ（ｄ）に着雪した雪が時間の経過に伴って段階的に落雪することを想定した場合には、時間が経過する毎に、一定の割合（例えば、１０％）又は一定の値（例えば、１ｋｇ／ｍ２及び２ｃｍ）を減算することによりリセットしてもよい。
【００５６】
表示情報生成部２２２は、ユーザ端末１１により送信された表示情報生成要求に応じて、気象予報データベース２１０及び着雪情報データベース２１１を参照することにより、各種の気象情報として、気象予報データ１００に含まれる気象要素１１０と、着雪情報算定部２２１により算定された風速Ｗｓ（ｔ，ｒ）、風向Ｗｄ（ｔ，ｒ）及び着雪情報２００とを、ユーザ端末１１の表示画面に表示させる表示情報を生成し、その表示情報をユーザ端末１１に通信部２３を介して送信する。
【００５７】
表示情報生成部２２２は、ユーザ端末１１に表示させる表示画面として、例えば、気象情報を表示する際に表示の対象とする表示日時及び表示要素を選択する初期表示画面３Ａ（詳細は図１０にて後述する）、初期表示画面３Ａにて選択された日時の気象情報を地図上に重畳して表示する選択表示画面３Ｂ（詳細は図１１にて後述する）、特定地点における気象情報の時間的な推移を表示する時系列グラフ表示画面４（詳細は図１２にて後述する）、及び、方位毎の着雪情報２００を比較可能に表示する第１乃至第３の複数方位時系列グラフ表示画面５Ａ〜５Ｃ（詳細は図１３〜１４、１６にて後述する）等を表示させる表示情報を生成する。
【００５８】
（着雪予測処理について）
以下に、着雪予測装置２が着雪情報２００を予測する際の動作（着雪予測処理）について、図４〜８を参照しつつ説明する。
【００５９】
図４は、本発明の実施形態に係る着雪予測装置２の動作（着雪予測処理）を示すフローチャートである。着雪予測装置２は、気象予報システム１０により気象予報データ１００が外部に提供される所定の提供時刻になると、図４に示すフローチャートに従って、着雪予測処理を開始する。なお、着雪予測装置２は、所定の提供時刻とは異なる任意のタイミングで着雪予測処理を開始してもよいし、着雪予測装置２の管理者から指示を受け付けた場合に着雪予測処理を開始してもよい。
【００６０】
まず、データ取得部２２０が、気象予報システム１０から通信部２３を介して気象予報データ１００を取得する（ステップＳ１）。気象予報データ１００は、上述したように、メソ数値予報モデル（ＭＳＭ）を用いて予測したＧＰＶデータであるものとし、１時間間隔で区切られた３９時間先までの各時間Ｔ（ｔ）、及び、５ｋｍメッシュで区切られた各領域Ｒ（ｒ）に対する複数の気象要素１１０として、気温Ｔｅｍｐ、相対湿度Ｈｒ、降水量Ｐｒ、風ベクトルＷ（東西成分ｗｕ及び南北成分ｗｖ）、及び、日射量Ｓｒを含むデータであるものとして説明する。
【００６１】
次に、データ取得部２２０は、ステップＳ１にて取得した気象予報データ１００を気象予報データベース２１０に登録する（ステップＳ２）。
【００６２】
次に、着雪情報算定部２２１は、気象予報データ１００の各時間Ｔ（ｔ）に対して、「ｔ」を変数とする繰り返しループＡを実行する（ステップＳ１０、Ｓ１１）。ステップＳ１０における「ｔ＝１，Ｔｍａｘ，１」は、変数ｔの初期値、終値、インクリメント値をそれぞれ示すものであり、着雪情報算定部２２１は、ステップＳ１０とステップＳ１１との間に挟まれたステップＳ２０〜Ｓ２１の処理を、変数ｔの初期値１に対して１ずつインクリメントしながら変数ｔが終値Ｔｍａｘ（＝３９）になるまで繰り返し行うことで、当該処理を各時間Ｔ（ｔ）に対して行う。
【００６３】
次に、着雪情報算定部２２１は、気象予報データ１００の各領域Ｒ（ｒ）に対して、「ｒ」を変数とする繰り返しループＢを実行する（ステップＳ２０、Ｓ２１）。ステップＳ２０における「ｒ＝１，Ｒｍａｘ，１」は、変数ｒの初期値、終値、インクリメント値をそれぞれ示すものであり、着雪情報算定部２２１は、ステップＳ２０とステップＳ２１との間に挟まれたステップＳ１００〜Ｓ６００の処理を、変数ｒの初期値１に対して１ずつインクリメントしながら変数ｒが終値Ｒｍａｘになるまで繰り返し行うことで、当該処理を各領域Ｒ（ｒ）に対して行う。
【００６４】
次に、風向風速計算部２２１ａは、ステップＳ２にて気象予報データベース２１０に登録された気象予報データ１００を参照し、気象予報データ１００に含まれる時間Ｔ（ｔ）及び領域Ｒ（ｒ）に対する風ベクトルＷ（東西成分ｗｕ及び南北成分ｗｖ）に基づいて、下記の式（１）により風速Ｗｓ（ｔ、ｒ）を計算する（Ｓ１００）。
【数１】
2020067343.tif 000003
【００６５】
次に、風向風速計算部２２１ａは、ステップＳ２にて気象予報データベース２１０に登録された気象予報データ１００を参照し、図５に示す風向計算処理により、気象予報データ１００に含まれる時間Ｔ（ｔ）及び領域Ｒ（ｒ）に対する風ベクトルＷ（東西成分ｗｕ及び南北成分ｗｖ）に基づいて、風向Ｗｄ（ｔ、ｒ）を計算する（Ｓ２００）。
【００６６】
図５は、本発明の実施形態に係る着雪予測装置２の動作（ステップＳ２００の風向計算処理の詳細）を示すフローチャートである。
【００６７】
まず、風向風速計算部２２１ａは、気象予報データ１００に含まれる時間Ｔ（ｔ）及び領域Ｒ（ｒ）に対する風ベクトルＷの東西成分ｗｕの値と「０」とを比較し、その比較した結果に応じて場合分けを行う（ステップＳ２１０）。
【００６８】
ステップＳ２１０にて、風向風速計算部２２１ａは、東西成分ｗｕの値が「０」と等しいと判断した場合（ステップＳ２１０で「ｗｕ＝０」）、風ベクトルＷの南北成分ｗｖの値と「０」とを比較し、その比較した結果に応じてさらに場合分けを行う（ステップＳ２２０）。
【００６９】
そして、ステップＳ２２０にて、風向風速計算部２２１ａは、南北成分ｗｖの値が「０」と等しいと判断した場合（ステップＳ２２０で「ｗｖ＝０」）、風向Ｗｄ（ｔ、ｒ）に「０」（無風）を代入し（ステップＳ２２１）、南北成分ｗｖの値が「０」より小さいと判断した場合（ステップＳ２２０で「ｗｖ＜０」）、風向Ｗｄ（ｔ、ｒ）に「２７０」（東風）を代入し（ステップＳ２２２）、南北成分ｗｖの値が「０」より大きいと判断した場合（ステップＳ２２０で「ｗｖ＞０」）、風向Ｗｄ（ｔ、ｒ）に「９０」（西風）を代入し（ステップＳ２２３）。
【００７０】
ステップＳ２１０にて、風向風速計算部２２１ａは、東西成分ｗｕの値が「０」より小さいと判断した場合（ステップＳ２１０で「ｗｕ＜０」）、風ベクトルＷの南北成分ｗｖの値と「０」とを比較し、その比較した結果に応じてさらに場合分けを行う（ステップＳ２３０）。
【００７１】
そして、ステップＳ２３０にて、風向風速計算部２２１ａは、南北成分ｗｖの値が「０」と等しいと判断した場合（ステップＳ２３０で「ｗｖ＝０」）、風向Ｗｄ（ｔ、ｒ）に「０」（北風）を代入し（ステップＳ２３１）、南北成分ｗｖの値が「０」より小さいと判断した場合（ステップＳ２３０で「ｗｖ＜０」）、風向Ｗｄ（ｔ、ｒ）に「ｔａｎ−１（ｗｕ／ｗｖ）×１８０／π」の計算値を代入し（ステップＳ２３２）、南北成分ｗｖの値が「０」より大きいと判断した場合（ステップＳ２３０で「ｗｖ＞０」）、風向Ｗｄ（ｔ、ｒ）に「ｔａｎ−１（ｗｕ／ｗｖ）×１８０／π＋３６０」の計算値を代入する（ステップＳ２３３）。
【００７２】
ステップＳ２１０にて、風向風速計算部２２１ａは、東西成分ｗｕの値が「０」より大きいと判断した場合（ステップＳ２１０で「ｗｕ＞０」）、風ベクトルＷの南北成分ｗｖの値と「０」とを比較し、その比較した結果に応じてさらに場合分けを行う（ステップＳ２４０）。
【００７３】
そして、ステップＳ２４０にて、風向風速計算部２２１ａは、南北成分ｗｖの値が「０」と等しいと判断した場合（ステップＳ２４０で「ｗｖ＝０」）、風向Ｗｄ（ｔ、ｒ）に「１８０」（南風）を代入し（ステップＳ２４１）、南北成分ｗｖの値が「０」と等しくないと判断した場合（ステップＳ２４０で「ｗｖ≠０」）、風向Ｗｄ（ｔ、ｒ）に「ｔａｎ−１（ｗｕ／ｗｖ）×１８０／π＋１８０」の計算値を代入する（ステップＳ２４２）。
【００７４】
以上のようにして、風向風速計算部２２１ａは、風向Ｗｄ（ｔ、ｒ）を計算し、図５に示す風向計算処理を終了する。
【００７５】
次に、図４に戻り、着雪情報算定部２２１は、ステップＳ２にて気象予報データベース２１０に登録された気象予報データ１００を参照し、図６に示す着雪情報算定処理により、時間Ｔ（ｔ）及び領域Ｒ（ｒ）の気象予報データ１００に基づいて、時間Ｔ（ｔ）及び領域Ｒ（ｒ）における着雪速度Ｍｓ（ｔ，ｒ，ｄ）、着雪重量Ｍｗ（ｔ，ｒ，ｄ）及び着雪高さＭｈ（ｔ，ｒ，ｄ）を８方位毎に算定する（Ｓ３００）。
【００７６】
図６は、本発明の実施形態に係る着雪予測装置２の動作（ステップＳ３００の着雪情報算定処理の詳細）を示すフローチャートである。
【００７７】
まず、湿雪降雪判定部２２１ｂは、気象予報データ１００に基づいて、時間Ｔ（ｔ）及び領域Ｒ（ｒ）において所定の湿雪降雪条件を満たすか否かを判定する（ステップＳ３０１〜Ｓ３０３）。
【００７８】
具体的には、湿雪降雪判定部２２１ｂは、気象予報データ１００に含まれる時間Ｔ（ｔ）及び領域Ｒ（ｒ）に対する気温Ｔｅｍｐが、−０．５℃〜２．５℃の範囲内にあるか否かの判定（ステップＳ３０１）、気象予報データ１００に含まれる時間Ｔ（ｔ）及び領域Ｒ（ｒ）に対する相対湿度Ｈｒが、８０％以上であるか否かの判定（ステップＳ３０２）、気象予報データ１００に含まれる時間Ｔ（ｔ）及び領域Ｒ（ｒ）に対する降水量Ｐｒが、０より大きいか否かの判定（ステップＳ３０３）を行い、全ての判定を満たす場合には（ステップＳ３０１で「Ｙｅｓ」、ステップＳ３０２で「Ｙｅｓ」、ステップＳ３０３で「Ｙｅｓ」）、所定の湿雪降雪条件を満たすとものとして、次のステップＳ３１０に進む。
【００７９】
次に、衝突重量計算部２２１ｃは、ステップＳ１００にて風向風速計算部２２１ａにより計算された風速Ｗｓ（ｔ、ｒ）と、気象予報データ１００に含まれる時間Ｔ（ｔ）及び領域Ｒ（ｒ）に対する降水量Ｐｒとに基づいて、時間Ｔ（ｔ）及び領域Ｒ（ｒ）における単位時間当たりの衝突重量Ｐｎ（ｔ，ｒ）を、下記の式（２）により計算する（ステップＳ３１０）。
【数２】
2020067343.tif 000004
ただし、落下速度：Ｆｓ，水の密度：Ｐｗ＝１
【００８０】
次に、着雪密度計算部２２１ｄは、ステップＳ１００にて風向風速計算部２２１ａにより計算された風速Ｗｓ（ｔ、ｒ）に基づいて、時間Ｔ（ｔ）及び領域Ｒ（ｒ）における着雪密度Ｐｓ（ｔ，ｒ）を、下記の式（３）により計算する（ステップＳ３２０）。
【数３】
2020067343.tif 000005
【００８１】
そして、着雪密度計算部２２１ｄは、計算した着雪密度Ｐｓ（ｔ，ｒ）が、８５０を超えていると判定した場合には（ステップＳ３２１で「Ｙｅｓ」）、着雪密度Ｐｓ（ｔ，ｒ）に「８５０」を代入する（ステップＳ３２２）。
【００８２】
次に、着雪情報算定部２２１は、ステップＳ２にて気象予報データベース２１０に登録された気象予報データ１００を参照し、図７に示す方位別参照処理により、時間Ｔ（ｔ）及び領域Ｒ（ｒ）の気象予報データ１００と、複数の方位Ｄ（ｄ）及び風向Ｗｄ（ｔ、ｒ）の間の関係により規定される方位毎の着雪率α（ｔ，ｒ，ｄ）とに基づいて、時間Ｔ（ｔ）及び領域Ｒ（ｒ）における着雪速度Ｍｓ（ｔ，ｒ，ｄ）、着雪重量Ｍｗ（ｔ，ｒ，ｄ）及び着雪高さＭｈ（ｔ，ｒ，ｄ）を８方位毎に算定する（Ｓ４００）。
【００８３】
図７は、本発明の実施形態に係る着雪予測装置２の動作（ステップＳ４００の方位別算定処理の詳細）を示すフローチャートである。
【００８４】
まず、着雪情報算定部２２１は、８方位の各方位Ｄ（ｄ）に対して、「ｄ」を変数とする繰り返しループＣを実行する（ステップＳ４１０、Ｓ４１１）。ステップＳ４１０における「ｄ＝１，Ｄｍａｘ，１」は、変数ｄの初期値、終値、インクリメント値をそれぞれ示すものであり、Ｄｍａｘは、８方位の場合には、「８」である。すなわち、着雪情報算定部２２１は、ステップＳ４１０とステップＳ４１１との間に挟まれたステップＳ４２０〜Ｓ４５０の処理を、変数ｄの初期値１に対して１ずつインクリメントしながら変数ｄが終値Ｄｍａｘになるまで繰り返し行うことで、当該処理を８方位の各方位Ｄ（ｄ）に対して行う。
【００８５】
次に、着雪率計算部２２１ｅは、方位Ｄ（ｄ）を向いた着雪面Ｓ（ｄ）に対する風向Ｗｄ（ｔ、ｒ）の入射角Φｄｅｇ（ｄ）を、下記の式（４）により計算する（ステップＳ４２０）。すなわち、入射角Φｄｅｇ（ｄ）は、風向Ｗｄ（ｔ、ｒ）と、方位Ｄ（ｄ）との間の角度差を示すものである。
【数４】
2020067343.tif 000006
ただし、ｍｏｄ（ｘ，ｙ）は、ｘをｙで除算したときの剰余を示す。
【００８６】
そして、着雪率計算部２２１ｅは、計算した入射角Φｄｅｇ（ｄ）が、２７０度以上であると判定した場合には（ステップＳ４２１で「Ｙｅｓ」）、入射角Φｄｅｇ（ｄ）に「ａｂｓ（Φｄｅｇ（ｄ）−３６０」（ａｂｓは絶対値）の計算値を代入する（ステップＳ４２２）。
【００８７】
次に、着雪率計算部２２１ｅは、入射角Φｄｅｇ（ｄ）の単位を度数法から弧度法（ラジアン）に変換することにより、入射角Φｒａｄ（ｄ）を求める（ステップＳ４２３）。
【００８８】
次に、着雪率計算部２２１ｅは、入射角Φｒａｄ（ｄ）に基づいて、時間Ｔ（ｔ）及び領域Ｒ（ｒ）における方位Ｄ（ｄ）の着雪率α（ｔ，ｒ，ｄ）を、下記の式（５）により計算する（ステップＳ４２４）。すなわち、着雪率α（ｔ，ｒ，ｄ）は、風向Ｗｄ（ｔ、ｒ）と方位Ｄ（ｄ）との間の角度差が「０」の場合には、最大値の１００［％］と計算され、風向Ｗｄ（ｔ、ｒ）と方位Ｄ（ｄ）との間の角度差が大きくなるほど、小さくなるような値として計算される。なお、式（５）の係数ｋ１、ｋ２は、降雪の状況や、着雪面として想定する部材の材質（金属、樹脂、木材等）や傾斜面の傾斜角等に応じて適宜変更してもよい。
【数５】
2020067343.tif 000007
ただし、ｋ１＝０．０４０７，ｋ２＝３．９４４４
【００８９】
次に、着雪率計算部２２１ｅは、入射角Φｄｅｇ（ｄ）が、９０度より大きく、かつ、２７０度より小さいと判定した場合には（ステップＳ４２５で「Ｙｅｓ」）、着雪率α（ｔ，ｒ，ｄ）に「０」を代入する（ステップＳ４２６）。
【００９０】
次に、着雪率計算部２２１ｅは、衝突率β（ｔ，ｒ，ｄ）に６０［％］を示す「０．６」を代入する（ステップＳ４２７）。
【００９１】
次に、着雪速度計算部２２１ｆは、ステップＳ３１０にて衝突重量計算部２２１ｃにより計算された単位時間当たりの衝突重量Ｐｎ（ｔ，ｒ）と、ステップＳ４２４、Ｓ４２６にて着雪率計算部２２１ｅにより計算された時間Ｔ（ｔ）及び領域Ｒ（ｒ）における方位Ｄ（ｄ）の着雪率α（ｔ，ｒ，ｄ）と、ステップＳ４２７にて代入された衝突率β（ｔ，ｒ，ｄ）とに基づいて、下記の式（６）により、時間Ｔ（ｔ）及び領域Ｒ（ｒ）における方位Ｄ（ｄ）の着雪速度Ｍｓ（ｔ，ｒ，ｄ）を計算する（ステップＳ４３０）。
【数６】
2020067343.tif 000008
【００９２】
次に、着雪重量累積部２２１ｇ１は、時間Ｔ（ｔ−１）及び領域Ｒ（ｒ）における着雪重量Ｍｗ（ｔ−１，ｒ，ｄ）と、ステップＳ４３０にて着雪速度計算部２２１ｆにより計算された着雪速度Ｍｓ（ｔ，ｒ，ｄ）とに基づいて、下記の式（７）により、着雪高さＭｈ（ｔ−１，ｒ，ｄ）に着雪速度Ｍｓ（ｔ，ｒ，ｄ）を加算することにより、時間Ｔ（ｔ）及び領域Ｒ（ｒ）における方位Ｄ（ｄ）の着雪重量Ｍｗ（ｔ，ｒ，ｄ）を計算する（ステップＳ４４０）。
【数７】
2020067343.tif 000009
【００９３】
次に、着雪高さ累積部２２１ｇ２は、時間Ｔ（ｔ−１）及び領域Ｒ（ｒ）における着雪高さＭｈ（ｔ−１，ｒ，ｄ）と、ステップＳ４３０にて着雪速度計算部２２１ｆにより計算された着雪速度Ｍｓ（ｔ，ｒ，ｄ）と、ステップＳ３２０にて着雪密度計算部２２１ｄにより計算された着雪密度Ｐｓ（ｔ，ｒ）とに基づいて、下記の式（８）により、着雪高さＭｈ（ｔ−１，ｒ，ｄ）に、着雪速度Ｍｓ（ｔ，ｒ，ｄ）を着雪密度Ｐｓ（ｔ，ｒ）で除算した値を加算することにより、時間Ｔ（ｔ）及び領域Ｒ（ｒ）における方位Ｄ（ｄ）の着雪高さＭｈ（ｔ，ｒ，ｄ）を計算する（ステップＳ４５０）。
【数８】
2020067343.tif 000010
【００９４】
なお、上記の式（７）、（８）において、ｔ＝１の場合の着雪重量Ｍｗ（ｔ−１，ｒ，ｄ）及び着雪高さＭｈ（ｔ−１，ｒ，ｄ）は、現在時刻における着雪重量及び着雪高さに相当するため、実際の観測値を用いてもよいし、過去に行われた着雪予測処理の結果として、着雪情報データベース２１１に登録されている着雪重量及び着雪高さを用いてもよいし、それらのデータがない場合には「０」としてもよい。
【００９５】
以上のようにして、着雪情報算定部２２１は、繰り返しループＣによりステップＳ４２０〜Ｓ４５０をＤｍａｘ回（＝８回、すなわち、８方位分）繰り返すことにより、時間Ｔ（ｔ）及び領域Ｒ（ｒ）における着雪速度Ｍｓ（ｔ，ｒ，ｄ）、着雪重量Ｍｗ（ｔ，ｒ，ｄ）及び着雪高さＭｈ（ｔ，ｒ，ｄ）を８方位毎に算定し、図６〜７に示す着雪情報算定処理を終了する。
【００９６】
次に、図４に戻り、落雪判定部２２１ｈは、ステップＳ２にて気象予報データベース２１０に登録された気象予報データ１００を参照し、図８に示す落雪判定処理により、時間Ｔ（ｔ）及び領域Ｒ（ｒ）の気象予報データ１００に基づいて、時間Ｔ（ｔ）及び領域Ｒ（ｒ）において所定の落雪条件を満たすか否かを判定する（ステップＳ５００）。
【００９７】
図８は、本発明の実施形態に係る着雪予測装置２の動作（ステップＳ５００の落雪判定処理の詳細）を示すフローチャートである。
【００９８】
まず、落雪判定部２２１ｈは、気象予報データ１００に含まれる時間Ｔ（ｔ）が、昼の時間帯か夜の時間帯かを判定する（ステップＳ５０１）。
【００９９】
次に、落雪判定部２２１ｈは、昼の時間帯であると判定した場合には（ステップＳ５０１で「昼」）、気象予報データ１００に含まれる時間Ｔ（ｔ）及び領域Ｒ（ｒ）に対する気温Ｔｅｍｐが５℃以上で１時間以上継続しているか否かの判定（ステップＳ５１０）、気温Ｔｅｍｐが３℃以上で２時間以上継続しているか否かの判定（ステップＳ５１１）、気温Ｔｅｍｐが２℃以上で４時間以上継続しているか否かの判定（ステップＳ５１２）、気温Ｔｅｍｐが、１時間前の気温（時間Ｔ（ｔ−１）及び領域Ｒ（ｒ）に対する気温Ｔｅｍｐ）に比較して１．５℃以上上昇しているか否かの判定（ステップＳ５１３）、気象予報データ１００に含まれる時間Ｔ（ｔ）及び領域Ｒ（ｒ）に対する日射量Ｓｒが「０」でない時間（日射が有る時間）が１時間以上継続しているか否かの判定（ステップＳ５１４）、気象予報データ１００に含まれる時間Ｔ（ｔ）及び領域Ｒ（ｒ）に対する降水量Ｐｒが「０」である時間（無降水の時間）が６時間以上継続しているか否かの判定（ステップＳ５１５）を行い、いずかの判定を満たす場合には、所定の落雪条件を満たすとものとして、次のステップＳ５３０に進む。
【０１００】
一方、落雪判定部２２１ｈは、夜の時間帯であると判定した場合には（ステップＳ５０１で「夜」）、気象予報データ１００に含まれる時間Ｔ（ｔ）及び領域Ｒ（ｒ）に対する気温Ｔｅｍｐが５℃以上で２時間以上継続しているか否かの判定（ステップＳ５２０）、気温Ｔｅｍｐが３℃以上で４時間以上継続しているか否かの判定（ステップＳ５２１）、気温Ｔｅｍｐが２℃以上で８時間以上継続しているか否かの判定（ステップＳ５２２）、気温Ｔｅｍｐが、１時間前の気温（時間Ｔ（ｔ−１）及び領域Ｒ（ｒ）に対する気温Ｔｅｍｐ）に比較して１．５℃以上上昇しているか否かの判定（ステップＳ５２３）、気象予報データ１００に含まれる時間Ｔ（ｔ）及び領域Ｒ（ｒ）に対する日射量Ｓｒが「０」でない時間（日射が有る時間）が１時間以上継続しているか否かの判定（ステップＳ５２４）、気象予報データ１００に含まれる時間Ｔ（ｔ）及び領域Ｒ（ｒ）に対する降水量Ｐｒが「０」である時間（無降水の時間）が１２時間以上継続しているか否かの判定（ステップＳ５２５）を行い、いずかの判定を満たす場合には、所定の落雪条件を満たすとものとして、次のステップＳ５３０に進む。
【０１０１】
なお、上記のステップＳ５１０〜Ｓ５１２、５１４〜５１５、Ｓ５２０〜Ｓ５２２、５２４〜５２５は、特定の条件を満たす時間が所定の時間以上継続しているか否かの判定をそれぞれ行うものである。例えば、ステップＳ５１２における「気温Ｔｅｍｐが２℃以上」という特定条件を満たす時間が「４時間」という所定の時間以上継続しているか否かの判定では、落雪判定部２２１ｈは、気象予報データ１００に含まれるデータのうち、判定の対象としている時間の気温（時間Ｔ（ｔ）及び領域Ｒ（ｒ）に対する気温Ｔｅｍｐ）を参照するだけでなく、判定の対象としている時間よりも１時間前の気温（時間Ｔ（ｔ−１）及び領域Ｒ（ｒ）に対する気温Ｔｅｍｐ）と、判定の対象としている時間よりも２時間前の気温（時間Ｔ（ｔ−２）及び領域Ｒ（ｒ）に対する気温Ｔｅｍｐ）と、判定の対象としている時間よりも３時間前の気温（時間Ｔ（ｔ−３）及び領域Ｒ（ｒ）に対する気温Ｔｅｍｐ）についても参照することにより、３時間前の気温、２時間前の気温、１時間前の気温、及び、判定の対象としている時間の気温という「４時間」分の気温が、全て「２℃以上」であるか否かを判定することで、ステップＳ５１２の判定を即時的に行うものである。
【０１０２】
また、上記のステップＳ５１３、Ｓ５２３における、気温Ｔｅｍｐが１時間前の気温に比較して１．５℃以上上昇しているか否かの判定では、落雪判定部２２１ｈは、気象予報データ１００に含まれるデータのうち、判定の対象としている時間の気温（時間Ｔ（ｔ）及び領域Ｒ（ｒ）に対する気温Ｔｅｍｐ）を参照するだけでなく、判定の対象としている時間よりも１時間前の気温（時間Ｔ（ｔ−１）及び領域Ｒ（ｒ）に対する気温Ｔｅｍｐ）についても参照することにより、１時間前の気温と、判定の対象としている時間の気温とを比較し、判定の対象としている時間の気温が１時間前の気温よりも１．５℃以上高いか否かを判定することで、ステップＳ５１３、Ｓ５２３の判定を即時的に行うものである。
【０１０３】
そして、着雪情報リセット部２２１ｉは、落雪判定部２２１ｈにより所定の落雪条件を満たすと判定された時間Ｔ（ｔ）及び領域Ｒ（ｒ）における着雪情報２００として、着雪重量Ｍｗ（ｔ，ｒ，ｄ）及び着雪高さＭｈ（ｔ，ｒ，ｄ）を「０」にリセットする（ステップＳ５３０）。すなわち、着雪情報リセット部２２１ｉは、８方位の全ての方位に対する着雪重量Ｍｗ（ｔ，ｒ，ｄ）及び着雪高さＭｈ（ｔ，ｒ，ｄ）を「０」にリセットし、図８に示す落雪判定処理を終了する。
【０１０４】
次に、図４に戻り、着雪情報算定部２２１は、ステップＳ３００の着雪情報算定処理（特にステップＳ４００の方位別算定処理）及びステップ５００の落雪判定処理で算定した時間Ｔ（ｔ）及び領域Ｒ（ｒ）における８方位毎の着雪速度Ｍｓ（ｔ，ｒ，ｄ）、着雪重量Ｍｗ（ｔ，ｒ，ｄ）及び着雪高さＭｈ（ｔ，ｒ，ｄ）を、着雪情報２００として着雪情報データベース２１１に登録する（ステップＳ６００）。
【０１０５】
以上のようにして、着雪情報算定部２２１は、繰り返しループＢによりステップＳ１００〜Ｓ６００をＲｍａｘ回繰り返すとともに、繰り返しループＡによりステップＳ２０〜Ｓ２１をＴｍａｘ回繰り返すことにより、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）における着雪速度Ｍｓ（ｔ，ｒ，ｄ）、着雪重量Ｍｗ（ｔ，ｒ，ｄ）及び着雪高さＭｈ（ｔ，ｒ，ｄ）を８方位毎に算定し、その算定した８方位毎の着雪速度Ｍｓ（ｔ，ｒ，ｄ）、着雪重量Ｍｗ（ｔ，ｒ，ｄ）及び着雪高さＭｈ（ｔ，ｒ，ｄ）を着雪情報データベース２１１に登録し、図４〜８に示す着雪予測処理を終了する。
【０１０６】
（表示情報生成処理と、表示画面について）
以下に、着雪予測装置２が表示情報を生成する際の動作（表示情報生成処理）と、着雪予測装置２により生成された表示情報に基づく表示画面について、図９〜１４を参照しつつ説明する。
【０１０７】
図９は、本発明の実施形態に係る着雪予測装置２の動作（表示情報生成処理）及びユーザ端末１１の動作を示すフローチャートである。着雪予測装置２は、着雪予測処理により気象予報データベース２１０及び着雪情報データベース２１１に表示の対象となるデータが登録されている場合に、表示情報生成処理を開始可能であり、着雪予測処理とは異なるタイミングで表示情報生成処理を開始してもよいし、着雪予測処理と並行して表示情報生成処理を開始してもよい。
【０１０８】
まず、ユーザ端末１１は、例えば、着雪予測装置２にアクセスするためのＵＲＬを指定するユーザのアクセス操作を受け付けると（ステップＳ７００）、その指定されたＵＲＬにアクセスするとともに、初期選択肢を含む表示情報生成要求を着雪予測装置２に送信する（ステップＳ７０１）。
【０１０９】
そして、着雪予測装置２の表示情報生成部２２２は、ユーザ端末１１により送信された表示情報生成要求に応じて、気象予報データベース２１０及び着雪情報データベース２１１を参照することにより、ユーザ端末１１に初期表示画面３Ａを表示させる表示情報を生成し、その生成した表示情報をユーザ端末１１に送信する（ステップＳ８００）。
【０１１０】
次に、ユーザ端末１１は、表示情報生成要求に対する応答として、着雪予測装置２から表示情報を受信し、その表示情報に基づいて、図１０に示す初期表示画面３Ａを表示する（ステップＳ７０２）。
【０１１１】
図１０は、本発明の実施形態に係る着雪予測システム１における初期表示画面３Ａを示す図である。初期表示画面３Ａは、画面左上側に配置された日時選択領域３０と、画面左下側に配置された要素選択領域３１と、画面右側に配置されて、各種の気象情報として、気象要素１１０及び着雪情報２００等を地図上に重畳して表示する気象情報表示領域３２と、を備える。なお、日時選択領域３０、要素選択領域３１及び気象情報表示領域３２を表示する際の配置やサイズは適宜変更してもよい。
【０１１２】
日時選択領域３０は、気象情報表示領域３２に気象情報を表示する際に、表示の対象とする表示日を、例えば、カレンダー形式で選択する日付選択部３００と、表示の対象とする表示時間を、スライダーバーや、停止ボタン、再生ボタン及び一時停止ボタンで選択する時間選択部３０１と、を備える。
【０１１３】
要素選択領域３１は、気象情報表示領域３２に気象情報を表示する際に、表示の対象とする表示要素を、例えば、コンボボックスで選択する要素選択部３１０と、表示の対象とする表示方位を、例えば、ラジオボタンで選択する方位選択部３１１と、要素選択部３１０で選択された表示要素とともに風向Ｗｄ（ｔ、ｒ）を表示するか否かを、例えば、チェックボックスで指定する風向指定部３１２と、を備える。
【０１１４】
要素選択部３１０は、コンボボックスによる表示要素の選択肢として、気象予報データ１００に含まれる気象要素１１０である気温Ｔｅｍｐ、相対湿度Ｈｒ、降水量Ｐｒ及び日射量Ｓｒと、着雪情報算定部２２１により風ベクトルＷ（東西成分Ｗｕ及び南北成分Ｗｖ）から計算された風速Ｗｓ（ｔ，ｒ）と、着雪情報算定部２２１により算定された着雪情報２００としての着雪速度Ｍｓ（ｔ，ｒ，ｄ）、着雪重量Ｍｗ（ｔ，ｒ，ｄ）及び着雪高さＭｈ（ｔ，ｒ，ｄ）とが選択可能に構成されている。
【０１１５】
方位選択部３１１は、要素選択部３１０にて、着雪速度Ｍｓ（ｔ，ｒ，ｄ）、着雪重量Ｍｗ（ｔ，ｒ，ｄ）及び着雪高さＭｈ（ｔ，ｒ，ｄ）のいずれかが選択された場合に有効となり、ラジオボタンによる表示方位の選択肢として、８方位のいずれかの方位が選択可能に構成されている。なお、方位選択部３１１は、表示方位の選択肢として、８方位の平均値をさらに選択可能に構成してもよい。
【０１１６】
気象情報表示領域３２は、その気象情報表示領域３２の全体に地図を表示する地図表示部３２０と、地図表示部３２０に表示する地図の表示範囲や拡大率を指定する地図指定部３２１と、を備える。
【０１１７】
地図表示部３２０は、地図指定部３２１により指定された表示範囲や拡大率に基づいて地図を表示するとともに、その地図上に、要素選択部３１０で選択された表示要素及び方位選択部３１１で選択された表示方位に対応する気象情報を重畳して表示するように構成されている。図１０に示す初期表示画面３Ａは、表示情報生成要求に含まれる初期選択肢として、日時選択領域３０で現時点からの１時間後の「２０１８年９月２０日１３：００」が選択され、要素選択部３１０で「着雪重量Ｍｗ（ｔ，ｒ，ｄ）」が選択され、方位選択部３１１で「北」が選択されている場合に表示される表示画面であり、地図表示部３２０は、方位「北」に対する着雪重量Ｍｗ（ｔ，ｒ，ｄ）のコンター図を地図上に重畳して表示する。なお、コンター図は、カラースケールのコンター図でもよいし、グレースケールのコンター図でもよい。
【０１１８】
また、地図表示部３２０は、地図上の任意の特定地点をマウスでダブルクリックするユーザの地点特定操作を受け付けることで、その特定地点に対する気象情報を時系列グラフで表示するように構成されている。
【０１１９】
なお、地図表示部３２０に表示する地図には、例えば、都道府県名、市町村名、山、川、湖等の地名、行政区画の境界線、道路や鉄道の路線図等の補助情報が表示されていてもよく、地図に表示する補助情報をユーザにより選択可能に構成されていてもよい。また、地図指定部３２１は、地図をマウスでドラッグする操作を受け付けることで、地図の表示範囲を指定してもよいし、地図上でマウスホイールを回転する操作を受け付けることで、地図の拡大率を指定してもよい。
【０１２０】
図１０に示す初期表示画面３Ａにおいて、ユーザ端末１１は、日時選択領域３０及び要素選択領域３１に対するユーザの選択操作を受け付けると（ステップＳ７１０）、その選択操作により選択された表示日時及び表示要素を含む表示情報生成要求を着雪予測装置２に送信する（ステップＳ７１１）。そして、着雪予測装置２の表示情報生成部２２２は、その表示情報生成要求に応じて気象予報データベース２１０及び着雪情報データベース２１１を参照することにより、ユーザ端末１１に選択表示画面３Ｂを表示させる表示情報を生成し、その生成した表示情報をユーザ端末１１に送信すると（ステップＳ８１０）、ユーザ端末１１は、その表示情報に基づいて、図１１に示す選択表示画面３Ｂを表示する（ステップＳ７１２）。
【０１２１】
図１１は、本発明の実施形態に係る着雪予測システム１における選択表示画面３Ｂを示す図である。選択表示画面３Ｂは、図１０に示す初期表示画面３Ａと基本的な画面構成を共通にしたものであり、初期表示画面３Ａにて、日時選択領域３０で現時点からの１８間後の「２０１８年９月２１日６：００」が選択され、要素選択部３１０で「着雪高さＭｈ（ｔ，ｒ，ｄ）」が選択され、方位選択部３１１で「西」が選択された場合に表示される表示画面である。
【０１２２】
地図表示部３２０は、地図指定部３２１により指定された表示範囲や拡大率に基づいて地図を表示するとともに、その地図上に、方位「西」に対する着雪高さＭｈ（ｔ，ｒ，ｄ）のコンター図を重畳して表示する。
【０１２３】
図１１に示す選択表示画面３Ｂにおいて、ユーザ端末１１は、日時選択領域３０及び要素選択領域３１に対するユーザの選択操作をさらに受け付けると、初期表示画面３Ａから選択表示画面３Ｂに切り換えたときと同様にして、その選択操作に応じた表示日時及び表示要素を表示する選択表示画面に切り換える。
【０１２４】
また、図１１に示す選択表示画面３Ｂにおいて、ユーザ端末１１は、気象情報表示領域３２に表示された地図上の任意の特定地点をマウスでダブルクリックするユーザの地点特定操作を受け付けると（ステップＳ７２０）、その地点特定操作で受け付けた特定地点及び表示要素を含む表示情報生成要求を着雪予測装置２に送信する（ステップＳ７２１）。そして、着雪予測装置２の表示情報生成部２２２は、その表示情報生成要求に応じて気象予報データベース２１０及び着雪情報データベース２１１を参照することにより、ユーザ端末１１に時系列グラフ表示画面４を表示させる表示情報を生成し、その生成した表示情報をユーザ端末１１に送信すると（ステップＳ８２０）、ユーザ端末１１は、その表示情報に基づいて、図１２に示す時系列グラフ表示画面４を表示する（ステップＳ７２２）。
【０１２５】
図１２は、本発明の実施形態に係る着雪予測システム１における時系列グラフ表示画面４を示す図である。時系列グラフ表示画面４は、図１１に示す選択表示画面３Ｂが表示された状態、すなわち、要素選択部３１０で「着雪高さＭｈ（ｔ，ｒ，ｄ）」が選択され、方位選択部３１１で「西」が選択された状態にて、地図上の任意の特定地点がマウスでダブルクリックされた場合に表示される表示画面である。
【０１２６】
時系列グラフ表示画面４は、画面右側に配置されて、時間を横軸とし、着雪高さＭｈを縦軸としたグラフに対して、特定地点を含む領域Ｒ（ｒ）における方位「西」に対する着雪高さＭｈ（ｔ，ｒ，ｄ）の時間的な推移を、例えば、折れ線４００で表示する時系列グラフ領域４０と、画面左側に配置された要素選択領域４１と、を備える。なお、図１２に示す時系列グラフ領域４０は、横軸として、例えば、２０１８年９月２１日０：００~１６：００の範囲（期間）を表示するものであるが、全ての期間を表示するようにしてもよいし、表示の対象とする期間をスクロールバーにより変更可能にしてもよいし、表示の対象とする期間の開始時刻及び終了時刻の少なくとも一方を変更可能にしてもよい。
【０１２７】
要素選択領域４１は、時系列グラフ領域４０に折れ線４００を表示する際に、折れ線４００による表示の対象とする表示要素を、例えば、複数選択可能なリストボックスで選択する要素選択部４１０と、表示の対象とする表示方位を、例えば、複数選択可能なリストボックスで選択する方位選択部４１１と、を備える。
【０１２８】
要素選択部４１０は、リストボックスによる表示要素の選択肢として、図１０〜１１に示す要素選択部３１０と同様の表示要素の１つ又は複数が選択可能に構成されている。方位選択部４１１は、要素選択部４１０にて、着雪速度Ｍｓ（ｔ，ｒ，ｄ）、着雪重量Ｍｗ（ｔ，ｒ，ｄ）及び着雪高さＭｈ（ｔ，ｒ，ｄ）のいずれかが選択された場合に有効となり、コンボボックスによる表示方位の選択肢として、８方位の１つ又は複数が選択可能に構成されている。
【０１２９】
図１２に示す時系列グラフ表示画面４において、ユーザ端末１１は、例えば、方位選択部４１１にて８方位のうち複数の方位を選択するユーザの複数方位選択操作を受け付けると（ステップＳ７３０）、その複数方位選択操作により選択された複数の表示方位を含む表示情報生成要求を着雪予測装置２に送信する（ステップＳ７３１）。そして、着雪予測装置２の表示情報生成部２２２は、その表示情報生成要求に応じて気象予報データベース２１０及び着雪情報データベース２１１を参照することにより、ユーザ端末１１に第１又は第２の複数方位時系列グラフ表示画面５Ａ、５Ｂを表示させる表示情報を生成し、その生成した表示情報をユーザ端末１１に送信すると（ステップＳ８３０）、ユーザ端末１１は、その表示情報に基づいて、図１３に示す第１の複数方位時系列グラフ表示画面５Ａ、又は、図１４に示す第２の複数方位時系列グラフ表示画面５Ｂを表示する（ステップＳ７３２）。
【０１３０】
図１３は、本発明の実施形態に係る着雪予測システム１における第１の複数方位時系列グラフ表示画面５Ａを示す図である。第１の複数方位時系列グラフ表示画面５Ａは、要素選択部４１０で「着雪高さＭｈ（ｔ，ｒ，ｄ）」が選択され、方位選択部４１１にて「北」、［東］、「南」及び［西］の４方位が選択された場合に表示される表示画面である。
【０１３１】
第１の複数方位時系列グラフ表示画面５Ａは、画面右側に上下方向に並べて配置された４つの時系列グラフ領域５０Ａ〜５０Ｄと、画面左側に配置されて、図１２に示す時系列グラフ表示画面４の要素選択領域４１と同様の機能を有する要素選択領域５１とを備える。
【０１３２】
４つの時系列グラフ領域５０Ａ〜５０Ｄは、特定地点を含む領域Ｒ（ｒ）における「北」、［東］、「南」及び［西］の４方位毎の着雪高さＭｈ（ｔ，ｒ，ｄ）の時間的な推移を示す４つの折れ線５００Ａ〜５００Ｄをそれぞれ表示することにより、領域Ｒ（ｒ）における方位毎の着雪重量Ｍｗ（ｔ，ｒ，ｄ）の時間的な推移を比較可能に表示する。
【０１３３】
図１４は、本発明の実施形態に係る着雪予測システム１における第２の複数方位時系列グラフ表示画面５Ｂを示す図である。第２の複数方位時系列グラフ表示画面５Ｂは、要素選択部４１０で「着雪高さＭｈ（ｔ，ｒ，ｄ）」が選択され、方位選択部４１１にて８方位の全てが選択された場合に表示される表示画面である。
【０１３４】
第２の複数方位時系列グラフ表示画面５Ｂは、画面右側に配置された１つの時系列グラフ領域５０と、画面左側に配置されて、図１２に示す時系列グラフ表示画面４の要素選択領域４１と同様の機能を有する要素選択領域５１と、を備える。
【０１３５】
１つの時系列グラフ領域５０は、特定地点を含む領域Ｒ（ｒ）における８方位毎の着雪高さＭｈ（ｔ，ｒ，ｄ）の時間的な推移を示す８つの折れ線５００Ａ〜５００Ｈを表示することにより、領域Ｒ（ｒ）における方位毎の着雪高さＭｈ（ｔ，ｒ，ｄ）の時間的な推移を比較可能に表示する。
【０１３６】
また、第２の複数方位時系列グラフ表示画面５Ｂは、時系列グラフ領域５０における横軸の各時間Ｔ（ｔ）に対応する位置に配置された６つの複数方位比較表示領域５２Ａ〜５２Ｆを備える。なお、図１４に示す第２の複数方位時系列グラフ表示画面５Ｂでは、６つの複数方位比較表示領域５２Ａ〜５２Ｆを３時間間隔で配置したが、任意の時間間隔で配置してもよいし、１つの複数方位比較表示領域に対して、６つの複数方位比較表示領域５２Ａ〜５２Ｆを１秒間隔で時系列順に表示することでアニメーションのように表示してもよい。
【０１３７】
複数方位比較表示領域５２Ａ〜５２Ｆの各々は、複数方位比較表示領域５２Ａ〜５２Ｆの中心部５２１を基準として８方位の各々に対応する各位置に配置された８つの着雪情報表示部５２０Ａ〜５２０Ｈを備える。８つの着雪情報表示部５２０Ａ〜５２０Ｈは、着雪情報表示部５２０Ａ〜５２０Ｈの各位置に対応する８方位毎の着雪高さＭｈ（ｔ，ｒ，ｄ）をそれぞれ表示することにより、時間Ｔ（ｔ）及び領域Ｒ（ｒ）における方位毎の着雪高さＭｈ（ｔ，ｒ，ｄ）を比較可能に表示する。
【０１３８】
複数方位比較表示領域５２Ａ〜５２Ｆの各々は、例えば、正八角形とし、その正八角形の内側に同心状の正八角形を配置し、外側の正八角形の各頂点と内側の正八角形の各頂点とを直線でそれぞれ結ぶことにより、８つの着雪情報表示部５２０Ａ〜５２０Ｈが、中心部５２１からみて８方位の各々に対応する各位置に配置される。そして、８つの着雪情報表示部５２０Ａ〜５２０Ｈは、着雪情報表示部５２０Ａ〜５２０Ｈの各位置に対応する８方位毎の着雪高さＭｈ（ｔ，ｒ，ｄ）を、例えば、その着雪高さＭｈ（ｔ，ｒ，ｄ）の値が大きくなるほど、中心部５２１から放射方向に長くなるようなバーグラフ５２２により表示する。なお、複数方位比較表示領域５２Ａ〜５２Ｆは、正八角形に代えて円形でもよいし、バーグラフ５２２による表示に代えてレーダーチャートやポーラチャートによる表示でもよい。
【０１３９】
なお、図１３〜１４に示す第１及び第２の複数方位時系列グラフ表示画面５Ａ、５Ｂにおいて、要素選択部４１０にて表示要素として、例えば、気温Ｔｅｍｐがさらに選択された場合には、第１の複数方位時系列グラフ表示画面５Ａは、気温Ｔｅｍｐに対応する時系列グラフ領域をさらに備え、その時系列グラフ領域に、気温Ｔｅｍｐの時間的な推移を示す折れ線を表示するようにしてもよいし、第２の複数方位時系列グラフ表示画面５Ｂは、時系列グラフ領域５０に、気温Ｔｅｍｐの時間的な推移を示す折れ線をさらに表示するようにしてもよい。
【０１４０】
また、図１０に示す初期表示画面３Ａ又は図１１に示す選択表示画面３Ｂにおいて、地図表示部３２０は、地図上の任意の特定地点にてマウスで右クリックするユーザの操作を受け付けることで、その特定地点を含む領域Ｒ（ｒ）における方位毎の着雪情報２００を、図１４に示す複数方位比較表示領域５２Ａ〜５２Ｆを含むポップアップ画面で表示するように構成してもよい。
【０１４１】
また、図９に示すフローチャートでは、図１０に示す初期表示画面３Ａを表示した後、図１１〜１４の順番に表示画面を表示するものとして説明したが、表示画面を表示する順番は、図１１〜１４の順番に限られず、ユーザの操作に応じて表示画面を表示すればよい。
【０１４２】
（落雪判定処理の変形例について）
上記実施形態において、着雪情報算定部２２１は、落雪判定部２２１ｈと、着雪情報リセット部２２１ｉとを備えるものとして説明したが、落雪判定部２２１ｈ及び着雪情報リセット部２２１ｉの変形例として、落雪判定部２２１ｈは、気象予報データ１００に基づいて、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）において所定の落雪条件を満たすか否かを方位毎に判定し、着雪情報リセット部２２１ｉは、落雪判定部２２１ｈにより所定の落雪条件を満たすと判定された時間Ｔ（ｔ）及び領域Ｒ（ｒ）における特定の方位の着雪情報２００をリセットするようにしてもよい。
【０１４３】
具体的には、落雪判定部２２１ｈは、気象予報データベース２１０に登録された気象予報データ１００を参照し、図１５に示す落雪判定処理により、気象予報データ１００に含まれる各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）に対する気象要素１１０が複数の方位Ｄ（ｄ）にそれぞれ及ぼす落雪への影響を考慮することにより、所定の落雪条件を満たすか否かを方位毎に判定する。
【０１４４】
図１５は、本発明の実施形態に係る着雪予測装置２の動作の変形例（ステップＳ５００の落雪判定処理の変形例）を示すフローチャートである。図１５に示す落雪判定処理では、落雪判定部２２１ｈが、気象要素１１０が複数の方位Ｄ（ｄ）にそれぞれ及ぼす落雪への影響を考慮する方法の一例として、複数の方位Ｄ（ｄ）と、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）における太陽の方位角Ｓｄと、気象予報データ１００に含まれる各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）に対する日射量Ｓｒとの間の関係に基づいて、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）における所定の落雪条件を方位毎に設定し、気象予報データ１００に基づいて、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）においてその設定した所定の落雪条件を満たすか否かを方位毎に判定するものとして説明する。
【０１４５】
まず、落雪判定部２２１ｈは、気象予報データ１００に含まれる時間Ｔ（ｔ）が、昼の時間帯か夜の時間帯かを判定する（ステップＳ５４１）。
【０１４６】
次に、落雪判定部２２１ｈは、昼の時間帯であると判定した場合には（ステップＳ５４１で「昼」）、時間Ｔ（ｔ）及び領域Ｒ（ｒ）における太陽の方位角Ｓｄを特定する（ステップＳ５５０）。
【０１４７】
次に、落雪判定部２２１ｈは、８方位の各方位Ｄ（ｄ）に対して、「ｄ」を変数とする繰り返しループＤを実行する（ステップＳ５６０、Ｓ５６１）。ステップＳ５６０における「ｄ＝１，Ｄｍａｘ，１」は、変数ｄの初期値、終値、インクリメント値をそれぞれ示し、落雪判定部２２１ｈは、ステップＳ５６０とステップＳ５６１との間に挟まれたステップＳ５７０〜Ｓ５７２の処理を、変数ｄの初期値１に対して１ずつインクリメントしながら変数ｄが終値Ｄｍａｘになるまで繰り返し行うことで、当該処理を８方位の各方位Ｄ（ｄ）に対して行う。
【０１４８】
次に、落雪判定部２２１ｈは、方位Ｄ（ｄ）と、ステップＳ５５０で特定した太陽の方位角Ｓｄと、気象予報データ１００に含まれる時間Ｔ（ｔ）及び領域Ｒ（ｒ）に対する日射量Ｓｒとの間の関係に基づいて、所定の落雪条件を設定する（ステップＳ５７０）。
【０１４９】
すなわち、落雪判定部２２１ｈは、太陽の方位角Ｓｄと方位Ｄ（ｄ）との間の角度差が小さくなるほど（方位Ｄ（ｄ）を向いた着雪面Ｓ（ｄ）が太陽の方位角Ｓｄを向いているほど）、日射量Ｓｒが、方位Ｄ（ｄ）を向いた着雪面Ｓ（ｄ）に及ぼす落雪への影響は大きくなるように作用するため、図８に示すステップＳ５１０〜Ｓ５１５の落雪条件を緩和する方向（落雪条件を満たすと判定しやすくする方向）に変更することにより、所定の落雪条件を設定する。反対に、落雪判定部２２１ｈは、太陽の方位角Ｓｄと方位Ｄ（ｄ）との間の角度差が大きくなるほど、図８に示すステップＳ５１０〜Ｓ５１５の落雪条件をより厳しくする方向（落雪条件を満たすと判定しにくくする方向）に変更することにより、所定の落雪条件を設定する。
【０１５０】
例えば、ステップＳ５１２における気温Ｔｅｍｐが２℃以上で４時間以上継続しているか否かの判定では、４時間という継続時間を短くすることにより、落雪条件を緩和する方向に変更し、反対に、４時間という継続時間を長くすることにより、落雪条件を厳しくする方向に変更する。ステップＳ５１４における日射量Ｓｒが「０」でない時間が１時間以上継続しているか否かの判定では、１時間という継続時間を短くすることにより、落雪条件を緩和する方向に変更し、反対に、１時間という継続時間を長くすることにより、落雪条件を厳しくする方向に変更する。ステップＳ５１５における降水量Ｐｒが「０」である時間が６時間以上継続しているか否かの判定では、６時間という継続時間を短くすることにより、落雪条件を緩和する方向に変更し、反対に、６時間という継続時間を長くすることにより、落雪条件を厳しくする方向に変更する。
【０１５１】
そして、落雪判定部２２１ｈは、ステップＳ５７０で設定した所定の落雪条件を満たすか否かを判定し（ステップＳ５７１）、所定の落雪条件を満たすと判定した場合には（ステップＳ５７１で「Ｙｅｓ」）、着雪情報リセット部２２１ｉが、落雪判定部２２１ｈにより所定の落雪条件を満たすと判定された時間Ｔ（ｔ）及び領域Ｒ（ｒ）における方位Ｄ（ｄ）に対する着雪重量Ｍｗ（ｔ，ｒ，ｄ）及び着雪高さＭｈ（ｔ，ｒ，ｄ）を「０」にリセットする（ステップＳ５７２）。
【０１５２】
以上のようにして、落雪判定部２２１ｈ及び着雪情報リセット部２２１ｉは、繰り返しループＤによりステップＳ５７０〜Ｓ５７２をＤｍａｘ回（＝８回、すなわち、８方位分）繰り返すことにより、落雪判定部２２１ｈが、気象予報データ１００に含まれる時間Ｔ（ｔ）及び領域Ｒ（ｒ）に対する気象要素１１０が、複数の方位Ｄ（ｄ）にそれぞれ及ぼす落雪への影響を考慮することにより、所定の落雪条件を満たすか否かを方位毎に判定し、着雪情報リセット部２２１ｉが、落雪判定部２２１ｈにより所定の落雪条件を満たすと判定された時間Ｔ（ｔ）及び領域Ｒ（ｒ）における特定の方位に対する着雪情報２００をリセットする。
【０１５３】
このとき、表示情報生成部２２２は、ユーザ端末１１に表示させる表示画面として、第３の複数方位時系列グラフ表示画面５Ｃを表示させる表示情報を生成することにより、ユーザ端末１１は、その表示情報に基づいて、図１６に示す第３の複数方位時系列グラフ表示画面５Ｃを表示する。
【０１５４】
図１６は、本発明の実施形態に係る着雪予測システム１における第３の複数方位時系列グラフ表示画面５Ｃを示す図である。第３の複数方位時系列グラフ表示画面５Ｃは、図１４に示す第２の複数方位時系列グラフ表示画面５Ｂと基本的な画面構成を共通にしたものである。
【０１５５】
第３の複数方位時系列グラフ表示画面５Ｃは、落雪判定部２２１ｈにより、時刻「１２：００」に方位「南」及び「南西」に対して落雪条件を満たすと判定され、着雪情報リセット部２２１ｉにより、時刻「１２：００」における方位「南」及び「南西」の着雪情報２００がリセットされるととともに、落雪判定部２２１ｈにより、時刻「１５：００」に方位「西」に対して落雪条件を満たすと判定され、着雪情報リセット部２２１ｉにより、時刻「１５：００」における方位「西」の着雪情報２００がリセットされた場合に表示される表示画面である。
【０１５６】
一方、落雪判定部２２１ｈは、夜の時間帯であると判定した場合には（ステップＳ５４１で「夜」）、図８に示すステップＳ５２０〜Ｓ５２５と同様の落雪条件を、所定の落雪条件として設定する（ステップＳ５８０）。
【０１５７】
そして、落雪判定部２２１ｈは、ステップＳ５８０で設定した所定の落雪条件を満たすか否かを判定し（ステップＳ５８１）、所定の落雪条件を満たすと判定した場合には（ステップＳ５８１で「Ｙｅｓ」）、着雪情報リセット部２２１ｉが、落雪判定部２２１ｈにより所定の落雪条件を満たすと判定された時間Ｔ（ｔ）及び領域Ｒ（ｒ）における着雪重量Ｍｗ（ｔ，ｒ，ｄ）及び着雪高さＭｈ（ｔ，ｒ，ｄ）を「０」にリセットする（ステップＳ５８２）。すなわち、着雪情報リセット部２２１ｉは、８方位の全ての方位に対する着雪重量Ｍｗ（ｔ，ｒ，ｄ）及び着雪高さＭｈ（ｔ，ｒ，ｄ）を「０」にリセットし、図１５に示す落雪判定処理を終了する。
【０１５８】
以上のように、上記実施形態に係る着雪予測装置２において、着雪情報算定部２２１は、気象予報データ１００に基づいて、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）における着雪情報２００を、複数の方位Ｄ（ｄ）をそれぞれ向いた複数の着雪面Ｓ（ｄ）の方位毎に算定するので、各方位を向いた複数の着雪面に対する着雪量が方位毎に異なるような状況に応じて、着雪量を的確に予測することができる。
【０１５９】
また、着雪情報算定部２２１は、気象予報データ１００と、複数の方位Ｄ（ｄ）及び風向Ｗｄ（ｔ，ｒ）の間の関係により規定される方位毎の着雪率α（ｔ，ｒ，ｄ）とに基づいて、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）における着雪情報２００を方位毎に算定するので、風向が各方位を向いた複数の着雪面にそれぞれ及ぼす着雪への影響を方位毎の着雪率に反映することで、方位毎の着雪情報を的確に算定することができる。
【０１６０】
また、着雪情報算定部２２１は、気象予報データ１００に基づいて、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）における風向Ｗｄ（ｔ、ｒ）及び風速Ｗｓ（ｔ、ｒ）を計算する風向風速計算部２２１ａと、気象予報データ１００に基づいて、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）における単位時間当たりの衝突重量Ｐｎ（ｔ，ｒ）［ｋｇ／ｈｒ］を計算する衝突重量計算部２２１ｃと、複数の方位Ｄ（ｄ）及び風向Ｗｄ（ｔ、ｒ）の間の関係に基づいて、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）における着雪率α（ｔ，ｒ，ｄ）［％］を方位毎に計算する着雪率計算部２２１ｅと、衝突重量計算部２２１ｃにより計算された衝突重量Ｐｎ（ｔ，ｒ）と、着雪率計算部２２１ｅにより計算された方位毎の着雪率α（ｔ，ｒ，ｄ）とに基づいて、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）における単位時間当たりの着雪重量を示す着雪速度Ｍｓ（ｔ，ｒ，ｄ）を着雪情報２００として方位毎に計算する着雪速度計算部２２１ｆと、着雪速度計算部２２１ｆにより計算された方位毎の着雪速度Ｍｓ（ｔ，ｒ，ｄ）を所定の単位時間分、方位毎に累積することにより、所定の単位時間後の各領域Ｒ（ｒ）における着雪重量Ｍｗ（ｔ，ｒ，ｄ）及び着雪高さＭｈ（ｔ，ｒ，ｄ）を着雪情報２００として方位毎に計算する着雪情報累積部２２１ｇと、を備えるので、着雪情報として、着雪速度と、着雪重量又は着雪高さとを方位毎に算定することができる。
【０１６１】
また、着雪情報算定部２２１は、気象予報データ１００に基づいて、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）において所定の落雪条件を満たすか否かを判定する落雪判定部２２１ｈと、落雪判定部２２１ｈにより所定の落雪条件を満たすと判定された時間Ｔ（ｔ）及び領域Ｒ（ｒ）における着雪情報２００をリセットする着雪情報リセット部２２１ｉと、をさらに備えるので、着雪した雪が落雪することを考慮することで、実際の着雪量に即した着雪情報を算定することができる。
【０１６２】
また、着雪情報算定部２２１は、気象予報データ１００に基づいて、各時間Ｔ（ｔ）及び各領域Ｒ（ｒ）において所定の落雪条件を満たすか否かを方位毎に判定する落雪判定部２２１ｈと、落雪判定部２２１ｈにより所定の落雪条件を満たすと判定された時間Ｔ（ｔ）及び領域Ｒ（ｒ）における方位の着雪情報２００をリセットする着雪情報リセット部２２１ｉと、をさらに備えるので、着雪した雪が落雪することを方位毎に考慮することで、実際の着雪量に即した着雪情報を算定することができる。
【０１６３】
また、落雪判定部２２１ｈは、気象要素１１０が複数の方位Ｄ（ｄ）にそれぞれ及ぼす落雪への影響を考慮することにより、所定の落雪条件を満たすか否かを方位毎に判定するので、気象要素が各方位にそれぞれ及ぼす落雪への影響が異なることにより、各方位を向いた複数の着雪面に着雪した雪が、方位毎に異なる条件で落雪するような状況に応じて、気象要素が各方位にそれぞれ及ぼす落雪への影響を所定の落雪条件に反映することで、実際の着雪量に即した着雪情報を算定することができる。
【０１６４】
また、表示情報生成部２２２は、方位毎の着雪情報２００を比較可能に第１又は第２の複数方位時系列グラフ表示画面５Ａ、５Ｂに表示させる表示情報を生成するので、着雪量が方位毎にどの程度異なるのかを視覚的に把握することができる。
【０１６５】
また、第２の複数方位時系列グラフ表示画面５Ｂは、方位毎の着雪情報２００を比較可能に表示する複数方位比較表示領域５２Ａ〜５２Ｆを備え、複数方位比較表示領域５２Ａ〜５２Ｆは、複数方位比較表示領域５２Ａ〜５２Ｆの中心部５２１を基準として方位の各々に対応する各位置に配置された複数の着雪情報表示部５２０Ａ〜５２０Ｈを備え、複数の着雪情報表示部５２０Ａ〜５２０Ｈは、各位置に対応する方位の着雪情報２００をそれぞれ表示するので、各方位と、各方位に対する着雪情報とが、着雪情報を表示する位置で対応付けられることで、着雪量が方位毎にどの程度異なるのかを視覚的に把握することができる。
【０１６６】
（他の実施形態）
以上、本発明の実施形態について説明したが、本発明は上記実施形態に限定されるものではなく、本発明の技術的思想を逸脱しない範囲で適宜変更可能である。
【０１６７】
例えば、上記実施形態では、気象予報データ１００は、風速Ｗｓ及び風向Ｗｄを示す気象要素１１０として、風ベクトルＷ（東西成分ｗｕ及び南北成分ｗｖ）を含むものとして説明したが、気象予報データ１００が、気象要素１１０として、風速Ｗｓ及び風向Ｗｄを含むようにしてしてもよく、その場合には、風向風速計算部２２１ａを省略してもよい。
【０１６８】
また、上記実施形態では、着雪予測装置２は、単一の情報処理装置で構成されているものとして説明したが、複数の情報処理装置で構成されていてもよく、例えば、データ取得部２２０及び着雪情報算定部２２１を備えた計算用サーバと、気象予報データベース２１０及び着雪情報データベース２１１を備えたデータサーバと、表示情報生成部２２２を備えた表示用サーバとで構成されていてもよい。
【０１６９】
また、上記実施形態では、着雪予測プログラム２１２は、記憶部２１に記憶されたものとして説明したが、インストール可能な形式又は実行可能な形式のファイルでＣＤ−ＲＯＭ、ＤＶＤ等のコンピュータで読み取り可能な記録媒体に記録されて提供されてもよいし、ネットワーク１２に接続されたコンピュータ上に格納し、ネットワーク経由でダウンロードさせることにより提供されてもよい。
【符号の説明】
【０１７０】
１…着雪予測システム、２…着雪予測装置、
３Ａ…初期表示画面、３Ｂ…選択表示画面、４…時系列グラフ表示画面、
５Ａ…第１の複数方位時系列グラフ表示画面、
５Ｂ…第２の複数方位時系列グラフ表示画面、
５Ｃ…第３の複数方位時系列グラフ表示画面、
１０…気象予報システム、１１…ユーザ端末、１２…ネットワーク、
２０…入力部、２１…記憶部、２２…制御部、２３…通信部、２４…出力部、
３０…日時選択領域、３１…要素選択領域、３２…気象情報表示領域、
４０…時系列グラフ領域、４１…要素選択領域、
５０、５０Ａ〜５０Ｄ…時系列グラフ領域、５１…要素選択領域、
５２Ａ〜５２Ｆ…複数方位比較表示領域、
１００…気象予報データ、１１０…気象要素、
２００…着雪情報、２１０…気象予報データベース、
２１１…着雪情報データベース、２１２…着雪予測プログラム、
２２０…データ取得部、２２１…着雪情報算定部、２２１ａ…風向風速計算部、
２２１ｂ…湿雪降雪判定部、２２１ｃ…衝突重量計算部、２２１ｄ…着雪密度計算部、
２２１ｅ…着雪率計算部、２２１ｆ…着雪速度計算部、２２１ｇ…着雪情報累積部、
２２１ｇ１…着雪重量累積部、２２１ｇ２…着雪高さ累積部、２２１ｈ…落雪判定部、
２２１ｉ…着雪情報リセット部、２２２…表示情報生成部、
３００…日付選択部、３０１…時間選択部、３１０…要素選択部、
３１１…方位選択部、３１２…風向指定部、３２０…地図表示部、
３２１…地図指定部、４００…折れ線、４１０…要素選択部、４１１…方位選択部、
５００Ａ〜５００Ｈ…折れ線、５２０Ａ〜５２０Ｈ…着雪情報表示部、
５２１…中心部、５２２…バーグラフ
【図１】
2020067343.tif 000011
【図２】
2020067343.tif 000012
【図３】
2020067343.tif 000013
【図４】
2020067343.tif 000014
【図５】
2020067343.tif 000015
【図６】
2020067343.tif 000016
【図７】
2020067343.tif 000017
【図８】
2020067343.tif 000018
【図９】
2020067343.tif 000019
【図１０】
2020067343.tif 000020
【図１１】
2020067343.tif 000021
【図１２】
2020067343.tif 000022
【図１３】
2020067343.tif 000023
【図１４】
2020067343.tif 000024
【図１５】
2020067343.tif 000025
【図１６】
2020067343.tif 000026