Python演習問題:NetCDFデータの読み込み & 動画作成

インポート・エイリアス設定

import numpy as np
import matplotlib.pyplot as plt
import matplotlib.animation as animation
import datetime
import netCDF4

問題

  1. NetCDFファイル,wspd.mon.mean.ncを読み込みなさい.
    ncdata = netCDF4.Dataset(filename, 'r', format='NETCDF4')
  2. ncファイルから変数lon, lat, time, wspdを取得しなさい.
    var_name = ncdata.variables['var_name'][:]
  3. 変数timeをpythonのdatetime型に変換しなさい.
    (補足1)変数timeの値は1800年1月1日0時0分0秒からの経過時間(単位:hour)である.
    (補足2)1800年1月1日から1900年1月1日までの日数は36524日である.
    datetime.datetime
    datetime.timedelta
  4. 変数wspdの最初の時刻wspd[0,:,:]のコンター図を描きなさい.
    fig = plt.figure()
    ax = fig.add_subplot(111)
    ax.pcolor(X, Y, C)
  5. 複数のステップを繰り返し描き,動画を作成しなさい.
    ani = animation.FuncAnimation(fig, function, fargs=(XX,XX,...), frames=n) # 方法①
    ani = animation.ArtistAnimation(fig, XX) # 方法②
    ani.save(filename)

データ

  • NCEP再解析データの月平均高度10m風速データ(京大生存圏データベース)

解答例(参考)

1〜24ステップの2年間分の動画を作成してみた.

wspd.gif

Front page   Edit Diff Attach Copy Rename Reload   New List of pages Search Recent changes   Help   RSS of recent changes
Last-modified: 2018-07-11 (Wed) 05:25:45 (314d)