GRDBLEND

grdblend − 部分的に重なり合う複数のグリッドファイルを 1 つの大きなグリ ッド上に混合する

grdblend blendfile −Ggrdfile −Ixinc[unit][=|+][/yinc[unit][=|+]] −Rwest/east/south/north[r] [ −Nnodata ] [ −Q ] [ −Zscale ] [ −V ] [ −W ] [ −fcolinfo ]

grdblend はグリッド化されたファイルの一覧を読み込んでパラメータを混合し て コサインテイパー加重を使って他のグリッドを混合することによりバイナリ グリッドファイルを作る。 grdblend は格子の一部にデータが無いときには 報 告 する。そのような拘束されない格子はユーザにより指定された値に設定され る[デフォルトは NaN ]。 1 つ以上の値のある格子は重み付き平均値に設定 さ れる。

blendfile

混合させるグリッドファイルごとに 1 行を持つアスキーファイル。そ れぞれの行にはグリッドファイル名、内部領域のための −R の設定 、 及 び相対加重 wr が、スペース又はタブで区切られて入っていなけれ ばならない。結合した重み付けスキームでは、このグリッドはその 領 域の外側では重み = ゼロが、内側の領域では重み = wr が、そして境 界領域におけるこれらの端成分の間では 2 次元のコサインテイパーさ れ た重みが与えられる。しかし、もし wr が与えられると、テイパー の意味が反対になる(すなわち、領域内部の重みがゼロになる)。ア ス キ ーファイルが与えられないときは grdblend は標準入力から読み込 む。

−G

grdfile はバイナリ出力グリッドファイルの名前である。

−I

x_inc [及び随意で y_inc ] は格子間隔である。分を示すには m を、秒を示すには c を付けること。もし代わりに e, k, i, n のいず れかの単位が付け加えられれば、それぞれ刻み幅はメートル、キロ メ ー トル、マイル、海里のいずれかと見なされ、その領域の中央の緯度 における経度の等価な度単位に変換される（変換方法は ELLIPSOID に 依 存 す る ）。もし /y_inc が与えられて 0 に設定されていれば、 /x_inc と同じ値に再設定され、そうでない場合は緯度の度に変換され る 。 もし = が付加されていたら、対応する最大の x (東) 又は y ( 北) が与えられた刻み幅にちょうど一致するように微調整される [ デ フォルトでは刻み幅は与えられた領域に合うように微調整される]。最 後に、刻み幅を与える代わりに、与えられた刻み幅に + を付加するこ と により、必要な格子点の数を指定することができる。刻み幅は格子 点の数と領域から再計算される。再計算される刻み幅は格子保存グ リ ッ ドか画素保存グリッドのどちらを選ぶかに依存する。詳細は付録 B を参照すること。

−R

xmin, xmax, ymin, ymax は対象とする領域を指定する。地理的領 域 に対しては、これらの境界は西、東、南、北に対応し、小数値を含 む度又は [+-]dd:mm[:ss.xxx][W|E|S|N] のフォーマットで指定できる 。東西南北の代わりに地図座標の左下と右上を使うときは r を付ける こと。 2 つの略号 −Rg −Rd は全球領域(それぞれ経度 0/360 又 は -180/+180 と緯度 -90/+90 )を表す。暦時間座標系に対しては相対時 間(選択された TIME_EPOCH に対して相対的で、選択された TIME_UNIT を単位とする。 t を −JX|x に付けること)又は [date]T[clock] 形式 の絶対時刻(T を −JX|x に付けること)のどちらかを与えることができ る。少なくとも date と clock のどちらかが存在していなければなら ない。 T は常に必要である。 date の文字列は [-]yyyy[-mm[-dd]] ( グレゴリオ暦)又は yyyy[-Www[-d]] (ISO 週暦)の形式でなければなら ない。 clock の文字列は hh:mm:ss[.xxx] の形式でなければならない 。 区切り文字を利用する場合にはその型と位置が指定どおりでなけれ ばならない(ただし入出力と描画フォーマットは変更可能である)。

−N

データ無し。入力グリッドが無い格子にこの値を設定する[デフォ ルトは NaN ]。

−Q

grdraster で使うのに適したヘッダの無いグリッドファイルを作 る。出力グリッドファイルはネイティブフォーマットである必要が あ る(すなわち netCDF ではない)。

−V

冗長モードを選択する。標準エラー出力に経過報告を送る[デフォ ルトでは「黙って」実行する]。

−W

混合せず、各格子に使われた重みだけを出力する。このオプシ ョ ン はただ 1 つだけグリッドが与えられたときだけ有効である[デフォ ルトでは混合する]。

−Z

出力値をファイルに書き出す前に scale でスケールする[ 1 ]。

−f

入出力列の特殊フォーマットを使う(時間又は地 理 デ ー タ) 。 i(nput) 又は o(utput) を指定すること [デフォルトでは入力と出力 の両方]。コンマ区切りで 1 列以上(又は列の範囲)を与えること。 各 列 又 は 列 の範囲の項目に対して T (絶対暦時刻)、 t (選択された TIME_EPOCH に対する相対時間)、 x (経度)、 y (緯度)、 f ( 小 数) の いずれかを付けること。 −f[i|o]g は −f[i|o]0x,1y (地理座標)の 省略形である。

4 つのグリッド化されたファイル piece_?.grd からグリッドファイルを作るた めに次のような blendfile を作る。

piece_1.grd -R<subregion_1> 1
piece_2.grd -R<subregion_2> 1
piece_3.grd -R<subregion_3> 1
piece_4.grd -R<subregion_4> 1

そして次のように実行する。

grdblend blend.job −Gblend.grd −R<full_region> −I<dx/dy> -V

GMT(l), grd2xyz(l), grdedit(l)