欠損値の扱い

2003年3月　西澤誠也 2004年3月　西澤誠也

欠損値を持ったデータを扱う際 NArrayMiss Class を使うと便利です。

NArraMissクラス
DCLでお絵書き
GPhys

NArrayMiss クラス

欠損値を持った NArray の配列を扱う場合
missing1.rb
require "narray" ndim = 10 rmiss = -999.0 x = NArray.sfloat(ndim).indgen x[3] = rmiss y = NArray.sfloat(ndim).indgen(10) y[5] = rmiss z = NArray.sfloat(ndim) for n in 0..ndim-1 if x[n]!=rmiss && y[n]!=rmiss then z[n] = x[n]+y[n] else z[n] = rmiss end end p z
もしくは
missing2.rb
require "narray" ndim = 10 rmiss =-999.0 x = NArray.sfloat(ndim).indgen x[3] = rmiss y = NArray.sfloat(ndim).indgen(10) y[5] = rmiss z = NArray.sfloat(ndim).fill(rmiss) mask = x.ne(rmiss)&y.ne(rmiss) z[mask] = x[mask]+y[mask] p z
となります。

ruby ではすべてはオブジェクトであり配列自身に欠損値情報をもたせることが出来ます。
NArrayMiss クラスは NArray の配列とマスクをあわせ持ったものです。 NArrayMiss クラスを利用すると
missing3.rb
require "narray_miss" ndim = 10 x = NArrayMiss.sfloat(ndim).indgen! x.invalidation(3) y = NArrayMiss.sfloat(ndim).indgen!(10) y.invalidation(5) z = x+y p z
のように、欠損のことを気にすること無く演算ができます。

四則演算以外にもいくつかのメソッドが用意されています。ドキュメント参照のこと
NArrayMiss#% NArrayMiss#** NArrayMiss#sum NArrayMiss#accum NArrayMiss#min NArrayMiss#max NArrayMiss#mean NArrayMiss#stddev
missing4.rb
require "narray_miss" ndim = 10 x = NArrayMiss.sfloat(ndim).indgen! x.invalidation(3) for method in ["sum","min","max","mean","stddev"] str = 'print "'+method+' = ",x.'+method+',"\n"' eval(str) end

NArrayからNArrayMissを作る際は以下のようにします。
NArrayMiss.to_nam(narray,mask)

DCLでお絵書き

DCLでお絵書きする際に欠損値処理を有効にするためには
DCL::gllset("lmiss",true)
とします。

簡単な例をあげます。
missing5.rb
require "numru/dcl" require "narray_miss" include NumRu rmiss = -9.99e10 nt = 50 nz = 50 tmin, tmax = 0.0, 5.0 zmin, zmax = 20.0,50.0 t = NArray.sfloat(nt+1, 1).indgen! * (tmax-tmin)/nt z = NArray.sfloat( 1,nz+1).indgen! * (zmax-zmin)/nz t[10..12,0] = rmiss z[0,5..10] = rmiss t = NArrayMiss.to_nam(t,t.ne(rmiss)) z = NArrayMiss.to_nam(z,z.ne(rmiss)) uz = NMMath::exp(-0.2*z)*(z**0.5) tz = -2.0*NMMath::exp(-0.1*z) u = uz*NMMath::sin(3.0*(tz+t)) DCL::gropn(1) DCL::gllset("lmiss",true) DCL::grfrm DCL::grswnd(tmin, tmax, zmin, zmax) DCL::uspfit DCL::grstrf DCL::usdaxs DCL::udcntr(u) DCL::grcls

NetCDF,GrADSファイルの解析

NetCDF,GrADSのデータを扱う場合は、 GPhysを使いましょう。内部でNArrayMissが使われます。くわしくは GPhys の説明を見てください。

return