Takuya Miyashita

Cartesian grid から Geographic grid への変換 †

状況 †

外部 function で生成した .grd や .nc などのファイルはデフォルトで Cartesian grid になっている場合があり，それを Geographic grid (緯度経度) 座標に変換したい．
緯度経度なのに Cartesian のままだと projection 方法によっては思うように図が描けなかったり面倒になったりする．
Cartesian か Geographic かは，grdinfo の表示から確認できる．

$ gmt grdinfo hogehoge.grd
...
hogehoge_cut_geo.grd: Remark: 
hogehoge.grd: Gridline node registration used [Cartesian grid]
hogehoge_cut_geo.grd: Grid file format: cf = GMT netCDF format (32-bit float, deprecated)
...

解決方法 †

共通オプションの -f を使う．
公式doc -> Common Options - GMT
grdconvert でも grdcut でも grdsample でも使い方は同じで，-fog をつけるだけ．
下記の hogehoge

gmt grdcut hogehoge.grd -R130/142/30/36 -Ghogehoge_cut_geo.grd -fog -V

オプションの -fog の o は output ファイルの指定を指す．
生成されたファイルを grdinfo で確認すると，

$ gmt grdinfo nankaiTrough_geo.grd 
...
hogehoge_cut_geo.grd: Remark: 
hogehoge_cut_geo.grd: Gridline node registration used [Geographic grid]
hogehoge_cut_geo.grd: Grid file format: cf = GMT netCDF format (32-bit float, deprecated)
...

のように Geographic grid に変わっているはず．