(57)【要約】  【課題】  メッシュ気象値の推定適用範囲を拡張してメッシュ気象値の推定精度を向上させる。  【解決手段】  固定気象観測地点の気象データで推定されるメッシュ気象値の適用範囲外の地点に気象センサ装置を設置し、この設置された気象センサ装置の観測データと、該気象センサ装置の設置地点に近い固定気象観測地点における既存の気象データとから、上記気象センサの設置地点における過去の気象データを推定し、これら気象センサの設置地点で推定される過去の気象データと固定気象観測地点の既存の気象データとを用いて地形因子解析により上記気象センサの設置地点におけるメッシュ気象値を推定することにより、気象センサ設置地点の予測メッシュ図が作成され、それに基づいて、固定気象観測地点の気象データのみで推定されるメッシュ気象値の適用範囲外の地域であってもきめ細かな精度の高い予測情報が得られる。  【選択図】      図１
2009257868.tif 000002
【特許請求の範囲】
【請求項１】
複数の固定気象観測地点で観測された気象データを用いて地形因子解析によりメッシュ気象値を推定する気象予測システムにおいて、
前記固定気象観測地点の気象データで推定されるメッシュ気象値の適用範囲外の地点に気象センサ装置を設置し、この設置された気象センサ装置の観測データと、該気象センサ装置の設置地点に近い前記固定気象観測地点における既存の気象データとから、上記気象センサの設置地点における過去の気象データを推定し、これら気象センサの設置地点で推定される過去の気象データと前記固定気象観測地点の既存の気象データとを用いて地形因子解析により上記気象センサの設置地点におけるメッシュ気象値を推定したことを特徴とする気象予測システム。
【請求項２】
前記気象センサ装置の観測データと、該気象センサ装置の設置地点に近い前記固定気象観測地点における既存の気象データとで相関式を作成し、上記気象センサ装置の観測期間の経過と連動して該気象センサ装置の作成を繰り返す請求項１記載の気象予測システム。
【請求項３】
前記固定気象観測地点の気象データで推定されるメッシュ気象値の適用範囲外の地点で異なる地域に気象センサ装置を複数台設置した請求項１または２記載の気象予測システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、例えば農業圃場の気象解析などに使用される気象予測システムに関する。
詳しくは、複数の固定気象観測地点で観測された気象データを用いて地形因子解析によりメッシュ気象値を推定する気象予測システムに関する。
【背景技術】
【０００２】
従来、この種の気象予測システムとしては、圃場の地理的情報を取得すると共に、この圃場内の複数箇所に分散配置したセンサ装置で測定した気象データ及び土壌データをリアルタイムに収集し、上記地理的情報を用いた地形因子解析で、農作業支援情報を予測するための一定区画毎のメッシュ気象値を推定することにより、その圃場の標高、起伏量、傾斜度とその向きなどの地形条件の影響を考慮した予測メッシュ図が作成され、それに基づいて地形が異なる圃場であってもきめ細かな精度の高い予測情報（作付け開始予測、生育・収量予測、地中温度予測、水温予測、土壌水分予測、病害虫発生予測）が提供されるようにした農作業決定支援用予測システムがある（例えば、特許文献１参照）。
【０００３】
【特許文献１】特開２００５−０８５０５９号公報（第６−９頁、図１−７）
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかし乍ら、このような従来の気象予測システムでは、メッシュ気象値の推定適用範囲が地形因子解析に用いた固定気象観測地点の地域範囲に制限されるため、その適用範囲外でメッシュ気象値を得ようとする場合には、気象メッシュ推定式を外挿することになり、推定精度が著しく低下するという問題があった。
このメッシュ気象値の推定精度が劣ると、それに基づいて得られる予測情報と実態との間に誤差が生じ、それにより、例えば農業圃場の気象解析に使用した場合には、病害虫発生予測や成熟期予測などの農業生産活動において種々の判断を求められる場面で、誤った判断を誘引する虞があった。
【０００５】
本発明のうち第一の発明は、メッシュ気象値の推定適用範囲を拡張してメッシュ気象値の推定精度を向上させることを目的としたものである。
第二の発明は、第一の発明の目的に加えて、気象センサの設置地点における過去の気象データの推定確率を向上させて更なるメッシュ気象値の精度向上を図ることを目的としたものである。
第三の発明は、第一の発明または第二の発明の目的に加えて、短い期間でメッシュ気象値の推定適用範囲を更に拡張してメッシュ気象値の推定精度を向上させることを目的としたものである。
【課題を解決するための手段】
【０００６】
前述した目的を達成するために、本発明のうち第一の発明は、固定気象観測地点の気象データで推定されるメッシュ気象値の適用範囲外の地点に気象センサ装置を設置し、この設置された気象センサ装置の観測データと、該気象センサ装置の設置地点に近い固定気象観測地点における既存の気象データとから、上記気象センサの設置地点における過去の気象データを推定し、これら気象センサの設置地点で推定される過去の気象データと固定気象観測地点の既存の気象データとを用いて地形因子解析により上記気象センサの設置地点におけるメッシュ気象値を推定したことを特徴とするものである。
第二の発明は、第一の発明の構成に、前記気象センサ装置の観測データと、該気象センサ装置の設置地点に近い前記固定気象観測地点における既存の気象データとで相関式を作成し、上記気象センサ装置の観測期間の経過と連動して該気象センサ装置の作成を繰り返す構成を加えたことを特徴とする。
第三の発明は、第一の発明または第二の発明の構成に、前記固定気象観測地点の気象データで推定されるメッシュ気象値の適用範囲外の地点で異なる地域に気象センサ装置を複数台設置した構成を加えたことを特徴とする。
【発明の効果】
【０００７】
本発明のうち第一の発明は、固定気象観測地点の気象データで推定されるメッシュ気象値の適用範囲外の地点に気象センサ装置を設置し、この設置された気象センサ装置の観測データと、該気象センサ装置の設置地点に近い固定気象観測地点における既存の気象データとから、上記気象センサの設置地点における過去の気象データを推定し、これら気象センサの設置地点で推定される過去の気象データと固定気象観測地点の既存の気象データとを用いて地形因子解析により上記気象センサの設置地点におけるメッシュ気象値を推定することにより、気象センサ設置地点の予測メッシュ図が作成され、それに基づいて、固定気象観測地点の気象データのみで推定されるメッシュ気象値の適用範囲外の地域であってもきめ細かな精度の高い予測情報が得られる。
従って、メッシュ気象値の推定適用範囲を拡張してメッシュ気象値の推定精度を向上させることができる。
その結果、固定気象観測地点の気象データで推定されるメッシュ気象値の適用範囲外の地域でメッシュ気象値を得ようとする場合には気象メッシュ推定式を外挿する従来のものに比べ、例えば農業圃場の気象解析に使用した場合には、病害虫発生や成熟期予測などの予測が的確になり、農業生産活動において種々の判断を求められる場面で、勘や経験と関係なく正確に判断できて、農業圃場の管理と高度な生産技術の獲得が図れる。
【０００８】
第二の発明は、第一の発明の効果に加えて、気象センサ装置の観測データと、該気象センサ装置の設置地点に近い前記固定気象観測地点における既存の気象データとで相関式を作成し、上記気象センサ装置の観測期間の経過と連動して該気象センサ装置の作成を繰り返すことにより、気象センサ装置の観測データと、その近くの固定気象観測地点に蓄積された過去の気象データとの相関係数が、上記気象センサ装置の観測期間が経過する度に改善される。
従って、気象センサの設置地点における過去の気象データの推定確率を向上させて更なるメッシュ気象値の精度向上を図ることができる。
【０００９】
第三の発明は、第一の発明または第二の発明の効果に加えて、固定気象観測地点の気象データで推定されるメッシュ気象値の適用範囲外の地点で異なる地域に気象センサ装置を複数台設置することにより、同時期に複数の異なる地域の予測メッシュ図が作成され、それらに基づいてきめ細かな精度の高い予測情報が得られる。
従って、短い期間でメッシュ気象値の推定適用範囲を更に拡張してメッシュ気象値の推定精度を向上させることができる。
【発明を実施するための最良の形態】
【００１０】
本発明の気象予測システムの実施形態は、対象地域内の固定気象観測地点に設置される既存の気象観測装置と、その対象地域外に設置される気象センサ装置と、これら既存の気象観測装置が設置される固定気象観測地点及び該気象センサ装置が設置される地点の地理的情報と、上記既存の気象観測装置で過去に観測された気象データ及び上記気象センサ装置の観測データを用いて地形因子解析によりメッシュ気象値を推定するためのソフトウエアからなる。
【００１１】
上記対象地域内の固定気象観測地点とは、設置場所を変えず長い年月に渡って観測可能な観測点を言い、そこに設置される既存の気象観測装置としては、例えば特開２００５−０８５０５９号公報に記載される農作業決定支援用予測システムのように、対象地域内に複数のセンサ装置を分散配置し、これらセンサ装置から所定期間に亘って観測した過去の気象データを使用することが好ましい。
それ以外には、上記対象地域内の固定気象観測地点に設置される既存の気象観測装置として、気象庁の地域気象観測システム（アメダス）などの観測所を利用し、これらの観測所から得られる過去の気象データを使用することも可能である。
【００１２】
この気象データの種類としては、気温、相対湿度、日射量、降水量、風向、風速、気圧などが含まれる。
また、特開２００５−０８５０５９号公報と同様に、対象地域内に分散配置されるセンサ装置で観測された過去の気象データを使用すれば、必要に応じて、各観測地点の地中温度や含水量や土壌成分の含有量などの土壌データを使用することも可能である。
【００１３】
上記対象地域外に設置される気象センサ装置としては、小型・軽量でしかもソーラー電源を有する無線通信方式のものが最適である。
更に、この気象センサ装置として、特開２００５−０８５０５９号公報と同様に例えば特開２００４−３０１５１１号公報に記載されるような無線中継機能を持ち、該気象センサ装置から無線送信される気象データを、既存のセンサ装置を順次経由して端末装置まで無線転送すれば、最も離れた気象センサ装置から端末装置までの間隔が無線到達距離以上になっても無線の送受信が正常に動作し、それにより、障害物の有無に関係なくどこからでも広い測定領域に亘って気象データを簡便かつ低コストで収集できるという利点がある。
また、上記無線通信方式としては、例えば小電力無線を使うタイプ、携帯電話を使うタイプやＰＨＳを使うが考えられ、これら無線送信の他に、部分的又は全体的にシリアルケーブルなどを介して端末装置と接続することも可能である。
【００１４】
上記対象地域内の固定気象観測地点及び上記対象地域外の気象センサ装置が設置される地点の地理的情報としては、国土地理院発行の数値地図ファイルなどの地理データを使用し、この地理データとして、最も細かな50ｍメッシュの標高データファイルを使用することが好ましい。
【００１５】
更に、上記地形因子解析によりメッシュ気象値を推定するためのソフトウエアとしては、地形要素計算ソフトウエアと、重回帰分析計算ソフトウエアと、気象メッシュ計算・作図ソフトウエアを有する。
【００１６】
この地形要素計算ソフトウエアは、上記固定気象観測地点及び気象センサ装置の設置地点の地理データとして得た標高データ（ｘ１）から、これらの地点における起伏量（ｘ２）、最大傾斜度（ｘ３）、その傾斜方向（ｘ４）、開放度などの地形の特徴を表現する要素を計算する。
【００１７】
上記重回帰分析計算ソフトウエアは、気象データを目的変数（ｙ）とし、上記地理データの地形要素を説明変数（ｘ１，ｘ２，ｘ３，ｘ４）とする重回帰式（ｙ＝ａ１・ｘ１＋ａ２・ｘ２＋ａ３・ｘ３＋ａ４・ｘ４＋ａ５）を作成し、その式の適合度検定を統計的手法で計算する。
なお、重回帰式中のａ１，ａ２，ａ３，ａ４，ａ５は係数である。
【００１８】
上記気象メッシュ計算・作図ソフトウエアは、上記重回帰分析計算ソフトウエアにより得られたメッシュ推定式を用いて、対象地域全体を分割したメッシュ毎に気象値を推定するソフトウエアであり、その気象要素としては気温、湿度、降水量、日射量、日照時間などの日別値である。
更に必要に応じて、この気象メッシュ計算・作図ソフトウエアにより計算した結果を、階級ごとに色分けし地域分布を作図することが好ましい。
【００１９】
そして、斯かる気象予測システムを用いてメッシュ推定式の適用地形範囲を拡張する方法を図１に示すフローチャートに沿って工程順に説明する。
先ず、既存の気象観測装置やアメダス観測所などが設置される固定気象観測地点の地形範囲以外の地点に、気象センサ装置を臨時的に設置し、気象データを収集する（ステップ１）。
【００２０】
次に、上記気象センサ装置から得られた観測データを用いて、該気象センサ装置の設置地点に最も近い上記固定気象観測地点における既存の気象データとの相関式を作成し（ステップ２）、この気象センサ装置が設置された臨時観測点の気象データを過去に遡って推定する（ステップ３）。
【００２１】
これに続いて、上記気象センサ装置が設置された臨時観測点で推定される過去の気象データと、上記固定気象観測地点の既存の気象データとを用いて地形因子解析を行い、それに加えて、新たに重回帰式を作成する。
それにより得られる重回帰式の適用範囲は、上記気象センサ装置が設置された臨時観測点の地形分だけ拡張される（ステップ４）。
【００２２】
最後に、上記固定気象観測地点における既存の気象データとの相関式で得られた上記気象センサ装置の臨時観測点の推定値を用いて、これら固定気象観測地点及び臨時観測点の地形範囲におけるメッシュ気象値を計算する（ステップ５）。
以下、本発明の一実施例を図面に基づいて説明する。
【実施例１】
【００２３】
この実施例は、図１〜図１１に示す如く、前記対象地域内の固定気象観測地点に設置される既存の気象観測装置として、例えば特開２００５−０８５０５９号公報に記載される農作業決定支援用予測システムのように、対象地域内に分散配置された複数のセンサ装置を使用し、これらセンサ装置が所定期間に亘って観測した過去の気象データを使用すると共に、対象地域外に一台の気象センサ装置が設置される場合を示すものであり、その具体的な手順を図１に示すフローチャートに基づいて説明する。
【００２４】
この特開２００５−０８５０５９号公報に記載される農作業決定支援用予測システムでは、対象地域の地理的情報を取得するステップ１１と、メッシュ区画を地形の複雑さの程度によって任意に選択して決定するステップ１２と、メッシュ区画毎の地形データを計算するステップ１３と、上記メッシュ区画毎に分散配置した複数のセンサ装置からのデータを収集するステップ１４と、各計測地点毎に月別・階級別の平均値を計算するステップ１５と、上述したデータを組み合わせて計測地点のデータによる地形因子解析を行うステップ１６と、メッシュ気象値を計算するステップ１７とからなる。
【００２５】
ここで、帯広市（2006年）の実例によって説明すれば、前記メッシュ推定式の適用メッシュ範囲（適用地域）を、平均標高＝20ｍ〜370ｍ、起伏量＝０ｍ〜22ｍ、最大傾斜度＝0.0〜0.22、最大傾斜の方向＝１：NNE〜８：N（時計回り）の地形範囲とする。
【００２６】
この対象地域内には、図２に示す如く、固定気象観測地点として７地点（Ａ，Ｂ，Ｃ，Ｄ，Ｅ，Ｆ，Ｇ）が設置されており、その地形範囲は、図３に示す如く、平均標高ｘ１＝79ｍ〜234ｍ、起伏量ｘ２＝１ｍ〜７ｍ、最大傾斜度ｘ３＝0.02〜0.06、最大傾斜の方向ｘ４＝１（NE）〜８（N）である。
上記の７地点（Ａ，Ｂ，Ｃ，Ｄ，Ｅ，Ｆ，Ｇ）は1996年に設置され、2006年までに約11年間の気象データが蓄積されている。
【００２７】
これら７地点の気象データを用いて、特開２００５−０８５０５９号公報に記載されるの手順により求めたメッシュ推定式の適用地形範囲は、これら固定気象観測地点の地形範囲内に制限される。
何故なら、上記地形範囲外の地形メッシュに対しては、気象メッシュ推定式を外挿で適用することになり、推定精度が著しく低下するからである。
そこで、その適用地形範囲外の地形に対して、高精度なメッシュ推定式を求める簡便法について以下に手順の詳細を説明する。
【００２８】
上記ステップ１７の後工程には、メッシュ推定式の適用地形範囲を拡張するかどうかのステップ１８があり、拡張しない場合にはそのまま終了するが、拡張する場合には、上述したステップ１〜５に進む。
【００２９】
先ず、地形条件がメッシュ推定式の適用地形範囲外の地点、即ち推定精度が低下する地点に前記気象センサ装置を設置し、この設置位置の気象データを実測して該気象データの収集を行う（ステップ１）。
今、この気象データの設置位置を図２に示すＪ地点とする。
Ｊ地点の地形条件は、図３に示す如く、平均標高＝278ｍ、起伏量＝16.0ｍ、最大傾斜度＝0.10、最大傾斜の方向＝６（NW）であり、このＪ地点と最も近い固定気象観測地点は図２に示す如く、Ｅ地点になる。
データ収集期間は、メッシュ推定式を適用する期間とし、例えば2006年５月１日から10月31日としている。
【００３０】
このＪ地点に設置される気象センサ装置として、特開２００５−０８５０５９号公報や特開２００４−３０１５１１号公報に記載されるセンサ装置と同様に小型・軽量なものを使用すれば、運搬が容易であるため、簡単に設置位置の移動ができるという利点がある。。
【００３１】
次に、上記気象センサ装置の設置地点（Ｊ地点）と、このＪ地点と最も近傍に存在する固定気象観測地点（Ｅ地点）との間で、相関式を作成する。
この相関式は、気象メッシュの計算に必要な気象要素毎に作成する。
また、相関式は期間を通して用いることができるように、通常の場合は一つの式で良いが、相関係数の劣る場合は月毎に作成する工夫が必要である。
【００３２】
そして、このように作成した相関式を用いて上記気象センサ装置の設置地点（Ｊ地点）の気象値を過去に遡って推定する。
固定気象観測地点の観測開始が1996年であるから、1996年から2005年までのＪ地点の過去データを作成することができる。
実際は上記気象センサ装置で観測をしていないが、気象データは存在する。
言い換えると、この気象センサ装置の設置地点（Ｊ地点）は、ヴァーチャルな気象観測地点といえる。
本発明の特許性は、短期間の気象観測を行って、その結果を近傍の観測地点の蓄積された実測データから過去に遡って推定し、メッシュ推定式の推定精度を向上させる考え方にある。
【００３３】
最後に、既存の固定気象観測地点（Ａ〜Ｇ地点）に新設した気象センサ装置の設置地点（Ｊ地点）を加えた合計８地点のデータ（1996年〜2006年までの11年分）で、再度地形因子解析によりメッシュ推定式を作成する。
【００３４】
ここで得られたメッシュ推定式の適用地形範囲は、平均標高＝79ｍ〜278ｍ、起伏量＝1ｍ〜16.0ｍ、最大傾斜度＝0.02〜0.10、最大傾斜の方向＝1（NE）〜8（N）となり、図３に示すメッシュ推定式の適用地形範囲が拡張される。
実例におけるメッシュ推定式の適用地形範囲の拡張前のメッシュ図を図４（ａ）に示し、拡張後のメッシュ図を図４（ｂ）に示しており、気象センサ装置の設置地点（Ｊ地点）が含まれる地域全体にメッシュ図が拡張されている。
【００３５】
このような方法でメッシュ推定式を改良した実例を気象要素別に図５〜図１１に示す。
これらは図２に示す如く、上記気象センサ装置の設置地点（Ｊ地点）と同じ適用地形範囲内のＵ地点に、精度検証用の気象センサ装置を設置し、その実測値と比較した結果を示している。
【００３６】
メッシュ推定式の修正前後における精度比較として、修正前の日平均気温を図５（ａ）に示し、修正後の日平均気温を図５（ｂ）に示している。
ここで、誤差の標準偏差を比較すると、修正前の推定誤差は1.5 ℃であったが、修正後には0.3 ℃まで小さくなり、メッシュ推定式の推定精度が大きく向上したことが解る。
【００３７】
同様に、修正前の日最高気温を図６（ａ）に示し、修正後の日最高気温を図６（ｂ）に示す。
修正前の日最低気温を図７（ａ）に示し、修正後の日最低気温を図７（ｂ）に示す。
修正前の日合計降水量を図８（ａ）に示し、修正後の日合計降水量を図８（ｂ）に示す。
修正前の日平均湿度を図９（ａ）に示し、修正後の日平均湿度を図９（ｂ）に示す。
修正前の日合計日射量を図１０（ａ）に示し、修正後の日合計日射量を図１０（ｂ）に示す。
修正前の日合計日照時間を図１１（ａ）に示し、修正後の日合計日照時間を図１１（ｂ）に示す。
これらについても、誤差の標準偏差を比較すると、修正前の推定誤差に比べて修正後には縮小され、それによりメッシュ推定式の推定精度が大きく向上したことが解る。
【００３８】
また、上記メッシュ推定式の適用地形範囲を更に拡張する場合は、翌年（2007年）に上記気象センサ装置を推定式の適用地形範囲外の場所に移設して観測する。
移設する場所の地形条件は、図３に示す如く、平均標高＝20ｍ〜79ｍか又は278ｍ〜370ｍ、起伏量＝０ｍ〜1ｍか又は16.0ｍ〜22ｍ、最大傾斜度＝0〜0.02か又は0.10〜0.22、最大傾斜の方向＝1（NE）〜８（N）である。
ここで、この移設ポイントを図２に示すＫ地点とし、その地形条件は平均標高＝350ｍ、起伏量＝19.0ｍ、最大傾斜度＝0.15、最大傾斜の方向＝４（S）とする。
【００３９】
このＫ地点の観測終了後、上述した手順により移設地点の過去データを遡って推定する。
この時、ヴァーチャル気象観測地点は、前年の観測ポイントを含めて２地点（Ｊ地点とＫ地点）となり、これらＪ地点とＫ地点の過去の気象データが揃ったら、再度、メッシュ推定式を作成し直す。
Ａ〜Ｇの７つの固定気象観測地点に、Ｊ地点とＫ地点の２つの移動観測地点を加えて、計９地点の1996年〜2007年までの12年分の過去データによりメッシュ推定式を作成する。
ここで得られたメッシュ推定式の適用地形範囲は、平均標高＝79ｍ〜350ｍ、起伏量＝1ｍ〜19.0ｍ、最大傾斜度＝0.02〜0.15、最大傾斜の方向＝１（NE）〜８（N）となり前年より更に拡張される。
【００４０】
それ以降、同様な方法を毎年繰り返すことによって、対象地域の全域をカバーする高精度なメッシュ推定式を作成することができる。
【００４１】
尚、前示実施例では、対象地域内の固定気象観測地点に設置される既存の気象観測装置として、特開２００５−０８５０５９号公報のように対象地域内に複数のセンサ装置を分散配置し、これらセンサ装置が所定期間に亘って観測した過去の気象データを使用したが、これに限定されず、例えば気象庁の地域気象観測システム（アメダス）などの観測所を利用し、これらの観測所から得られる過去の気象データを使用しても良い。
この場合も上述した実施例と同様な作用効果が得られる。
【００４２】
更に、対象地域外に一台の気象センサ装置が設置される場合を示したが、これに限定されず、対象地域外の異なる地形の場所に複数台の気象センサ装置を設置しても良い。
この場合には、同時期に複数の異なる地域のメッシュ図が作成され、更に短い年月で適用地形範囲を拡張することができ、高精度なメッシュ推定式の作成が可能となる。
【図面の簡単な説明】
【００４３】
【図１】本発明の一実施例を示す気象予測システムのフローチャートである。
【図２】気象観測地点の配置図である。
【図３】メッシュ推定式の適用地形範囲を示す表である。
【図４】メッシュ推定式の適用地形範囲のメッシュ図であり、（ａ）が拡張前を示し、（ｂ）が拡張後を示している。
【図５】メッシュ推定式の修正前後における精度比較したグラフであり、（ａ）が修正前の日平均気温を示し、（ｂ）が修正後の日平均気温を示している。
【図６】メッシュ推定式の修正前後における精度比較したグラフであり、（ａ）が修正前の日最高気温を示し、（ｂ）が修正後の日最高気温を示している。
【図７】メッシュ推定式の修正前後における精度比較したグラフであり、（ａ）が修正前の日最低気温を示し、（ｂ）が修正後の日最低気温を示している。
【図８】メッシュ推定式の修正前後における精度比較したグラフであり、（ａ）が修正前の日合計降水量を示し、（ｂ）が修正後の日合計降水量を示している。
【図９】メッシュ推定式の修正前後における精度比較したグラフであり、（ａ）が修正前の日平均湿度を示し、（ｂ）が修正後の日平均湿度を示している。
【図１０】メッシュ推定式の修正前後における精度比較したグラフであり、（ａ）が修正前の日合計日射量を示し、（ｂ）が修正後の日合計日射量を示している。
【図１１】メッシュ推定式の修正前後における精度比較したグラフであり、（ａ）が修正前の日合計日照時間を示し、（ｂ）が修正後の日合計日照時間を示している。
【図１】
2009257868.tif 000003
【図２】
2009257868.tif 000004
【図３】
2009257868.tif 000005
【図４】
2009257868.tif 000006
【図５】
2009257868.tif 000007
【図６】
2009257868.tif 000008
【図７】
2009257868.tif 000009
【図８】
2009257868.tif 000010
【図９】
2009257868.tif 000011
【図１０】
2009257868.tif 000012
【図１１】
2009257868.tif 000013