ROOT/(x,y,z)形式の3次元データを扱う

(x,y,z)形式の3次元データを扱う†

TGraph2D†

• ROOT; TGraph2D Class Reference

TGraph2D *g = new TGraph2D("input.dat", "%*s %lg %lg %lg"); // 読み込みファイル指定。%*sで読み飛ばす行。
g->Draw("TRI"); // グラフの描画
double z = g->Interpolate(22,1.2); // ドロネー補間により、任意のx,y上のzを返す
int n = g->GetN(); // 読み込んだデータ点の数
Double_t *x = g->GetX(); // 読み込んだデータ点を変数 *x に代入
printf("%f\n",x[3]); // 読み込んだデータ点の3番目のxを返す

等高線プロットを作る†

• 簡単にやるには

  g->Draw("CONT5");

  で出来る。デフォルトだと20本の等高線が引かれる。それぞれの等高線は、内部的にはTGraphになっている。
  • 03.Histograms.pdf
• 自力で1本づつ引くには

  ((TGraph*)((TList*)g->GetContourList(z))->At(0))->Draw("C");

  とする。
  • 但し、GetContourListを使う前に
• バグ?
  • z軸のlogスケールに設定されたキャンバスでDrawするとzの値がlog10を取った値に書き換えられるっぽい挙動をする。要検証。

Topへ戻る