Takuya Miyashita

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風速データ(京大生存圏データベース)
  
  wspd.mon.mean.nc

解答例(参考) †

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