この文書の現在のバージョンと選択したバージョンの差分を表示します。
両方とも前のリビジョン 前のリビジョン 次のリビジョン | 前のリビジョン | ||
py_scenalio [2018/11/12 22:11] ohno |
— (現在) | ||
---|---|---|---|
ライン 1: | ライン 1: | ||
- | ==== メッシュ温暖化シナリオデータ利用方法 ==== | ||
- | 温暖化シナリオデータは将来を予測したものではないので、将来の特定の年のデータを処理して得た結果には意味がありません。将来の数十年分のデータを処理し、現在条件の数十年分のデータを処理し、それらに統計処理を加えたうえで比較して初めて意味のある結果を得ることができます。場合によっては、数十年分の本当の現在気象値のデータを処理した結果も比較に加える必要があります。つまり、温暖化影響を評価するには気象データを大量に処理する技術が必要です。\\ | ||
- | このため、農研機構は、このデータを取得するエクセルシートは用意しません。メッシュ温暖化シナリオデータをPythonプログラムに読み込むための関数「GetSceData」を用意しています。この関数は、メッシュ農業気象データ利用ツール**AMD_tools3**に搭載されています。以下にその説明を記します。 | ||
- | 関数名:GetSceData | ||
- | 概要:気候予測シナリオデータを、気象データをデータ配信サーバーまたはローカルファイルから取得する関数。 | ||
- | 書式: | ||
- | GetSceData(element, timedomain, lalodomain, model, scenam, area=None, namuni=False, url='https://amd.rd.naro.go.jp/opendap/AMS') | ||
- | 引数(必須): | ||
- | element:気象要素記号で、'TMP_mea'などの文字列で与える | ||
- | timedomain:取得するデータの時間範囲で、['2008-05-05', '2008-05-05'] | ||
- | のような文字列の2要素リストで与える。特定の日のデータを | ||
- | 取得するときは、二カ所に同じ日付を与える。 | ||
- | lalodomain:取得するデータの緯度と経度の範囲で、 | ||
- | [36.0, 40.0, 130.0, 135.0] のように緯度,緯度,経度,経度の順で指定する。 | ||
- | 特定地点のデータを取得するときは、緯度と経度にそれぞれ同じ値を与える。 | ||
- | model:気候モデルの記号で、'MIROC5'などの文字列で与える | ||
- | scenam:排出シナリオ等の記号で、'RCP8.5'などの文字列で与える | ||
- | 引数(必要に応じ指定): | ||
- | namuni:True => 気象要素の正式名称と単位を取り出す。戻り値の数は2つ増えて6つになる。 | ||
- | False => 気象要素の正式名称を取り出さない。戻り値の数は4つ(気象値、時刻、緯度、経度)。 | ||
- | area:データを読み出すエリア(Area1~Area6)を指定する。省略した場合は自動的に選ばれる。 | ||
- | url:データファイルの場所を指定する。省略した場合はデータ配信サーバーに読みに行く。 | ||
- | ローカルにあるファイルを指定するときは、AreaN(N=1~6)の直上(通常は"・・・/AMS")を指定する。 | ||
- | 戻り値: | ||
- | 第1戻り値:指定した気象要素の三次元データ。[時刻、緯度、経度]の次元を持つ。 | ||
- | 第2戻り値:切り出した気象データの時刻の並び。Pythonの時刻オブジェクトの一次元配列である。 | ||
- | 第3戻り値:切り出した気象データの緯度の並び。実数の一次元配列である。 | ||
- | 第4戻り値:切り出した気象データの経度の並び。実数の一次元配列である。 | ||
- | 第5戻り値(namuni=Trueのときのみ):気象データの正式名称。文字列である。 | ||
- | 第6戻り値(namuni=Trueのときのみ):気象データの単位。文字列である。 | ||
- | | ||
- | 使用例:MIROC5モデルで予測したRCP8.5シナリオにおける、北緯35度、東経135度の地点の2020年~2030年の日最高気温を取得する場合。 | ||
- | import AMD_Tools3 as AMD | ||
- | model = 'MIROC5' | ||
- | scenario = 'RCP8.5' | ||
- | timedomain = ['2020-01-01', '2030-12-31'] | ||
- | lalodomain = [35.0, 35.0, 135.0, 135.0] | ||
- | Tm, tim, lat, lon = AMD.GetSceData('TMP_max', timedomain, lalodomain,model, scenario) | ||
- | |||
- | ===サンプルプログラム=== | ||
- | * {{ ::testscena.py |testscena.py}}:トップページの図を作成するサンプルプログラム。 | ||
- | |||
- | [[start#気候変化シナリオも搭載されています。|戻る]] |