Pythonを用いて等高線図を作成する方法について
xy座標とその点の高さをまとめた下のようなデータから、Pythonを用いて等高線図を作成するにはどのようにコードを作れば良いでしょうか?
(データのイメージ)
x y 高さ
0 0 2
1 1 5
2 2 3
3 3 2
0 1 1
1 0 4
2 0 3
xとyの範囲を決めてその範囲内で等間隔にメッシュを切り、格子点における高さを何らかの関数で与えてその高さの分布を等高線図として作成するというような方法については例がいくつか見つかるのですが、上記のようなことはどのようにして実現できるでしょうか。
アドバイス頂ければ幸いです。
No.1ベストアンサー
- 回答日時:
scatter,contour,surface等で等高線図を色々と作れます。
以下、scatterの場合のプログラム例です。色々と試せる様にコメント(#)にしている
ところがあります。グラフ2例添付しました。
以下のサンプルプログラムでどんな風になるか、試せると思います。
from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt
import numpy as np
#x = np.arange(-3, 3, 0.2)
#y = np.arange(-3, 3, 0.2)
#X, Y = np.meshgrid(x, y)
#Z = np.sin(X)+ np.cos(Y)
X=[0,1,2,3,0,1,2];
Y=[0,1,2,3,1,0,0];
Z=[2,5,3,2,1,4,3]
fig = plt.figure()
ax = Axes3D(fig)
#ax.plot_wireframe(X,Y,Z) # ワイヤーフレームプロット
#ax.plot_surface(X, Y, Z, rstride=1, cstride=1)# サーフェイス
#ax.plot3D(np.ravel(X),np.ravel(Y),np.ravel(Z))
#ax.contour3D(X,Y,Z)
#ax.contourf3D(X,Y,Z)
ax.scatter3D(np.ravel(X),np.ravel(Y),np.ravel(Z))
plt.xlabel('x'); plt.ylabel('y');
plt.show()
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
関連するカテゴリからQ&Aを探す
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・【大喜利】【投稿~11/12】 急に朝起こしてきた母親に言われた一言とは?
- ・好きな和訳タイトルを教えてください
- ・うちのカレーにはこれが入ってる!って食材ありますか?
- ・好きな「お肉」は?
- ・あなたは何にトキメキますか?
- ・おすすめのモーニング・朝食メニューを教えて!
- ・「覚え間違い」を教えてください!
- ・とっておきの手土産を教えて
- ・「平成」を感じるもの
- ・秘密基地、どこに作った?
- ・【お題】NEW演歌
- ・カンパ〜イ!←最初の1杯目、なに頼む?
- ・一回も披露したことのない豆知識
- ・これ何て呼びますか
- ・チョコミントアイス
- ・初めて自分の家と他人の家が違う、と意識した時
- ・「これはヤバかったな」という遅刻エピソード
- ・これ何て呼びますか Part2
- ・許せない心理テスト
- ・この人頭いいなと思ったエピソード
- ・牛、豚、鶏、どれか一つ食べられなくなるとしたら?
- ・あなたの習慣について教えてください!!
- ・ハマっている「お菓子」を教えて!
- ・高校三年生の合唱祭で何を歌いましたか?
- ・【大喜利】【投稿~11/1】 存在しそうで存在しないモノマネ芸人の名前を教えてください
- ・好きなおでんの具材ドラフト会議しましょう
- ・餃子を食べるとき、何をつけますか?
- ・あなたの「必」の書き順を教えてください
- ・ギリギリ行けるお一人様のライン
- ・10代と話して驚いたこと
- ・家の中でのこだわりスペースはどこですか?
- ・つい集めてしまうものはなんですか?
- ・自分のセンスや笑いの好みに影響を受けた作品を教えて
- ・【お題】引っかけ問題(締め切り10月27日(日)23時)
- ・大人になっても苦手な食べ物、ありますか?
- ・14歳の自分に衝撃の事実を告げてください
- ・架空の映画のネタバレレビュー
- ・「お昼の放送」の思い出
- ・昨日見た夢を教えて下さい
- ・ちょっと先の未来クイズ第4問
- ・【大喜利】【投稿~10/21(月)】買ったばかりの自転車を分解してひと言
- ・メモのコツを教えてください!
- ・CDの保有枚数を教えてください
- ・ホテルを選ぶとき、これだけは譲れない条件TOP3は?
- ・家・車以外で、人生で一番奮発した買い物
- ・人生最悪の忘れ物
- ・【コナン30周年】嘘でしょ!?と思った○○周年を教えて【ハルヒ20周年】
- ・あなたの習慣について教えてください!!
- ・都道府県穴埋めゲーム
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
プロダクションコードとは?
-
vba listviewにおけるtextのAli...
-
初級者の質問――time_tについて
-
C言語primeについて
-
COBOLのS9タイプからXタイプへ...
-
Arduino nano(アルドゥイーノ ...
-
指定したフォルダーを最前面表...
-
GitHubで拾ったプログラムが動...
-
プログラミングの問題と解答例...
-
DataGridViewで指定列の背景色...
-
無限ループ中にある任意のキー...
-
テキストの暗号化
-
VBAのプログラムが動いているか...
-
VBコンソールアプリ(.NET Frame...
-
2値化した画像の黒い部分のピク...
-
arduino mpu6050 計測できない
-
URLの最初の「h」を略す
-
C言語でのソケット通信のclose
-
XPathの関数「text()="value"」...
-
原子同士の結合をプログラムで...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
プロダクションコードとは?
-
指定したフォルダーを最前面表...
-
Arduino nano(アルドゥイーノ ...
-
vba listviewにおけるtextのAli...
-
C言語primeについて
-
世界一美しいソースコード
-
COBOLのS9タイプからXタイプへ...
-
マイコンからプログラムを読み...
-
無限ループ中にある任意のキー...
-
2値化した画像の黒い部分のピク...
-
初級者の質問――time_tについて
-
c言語のサイコロを100回振って...
-
Excel VBA じゃんけん
-
コンパイルエラーの意味
-
VBAにてTIFFから画像を取込みた...
-
arduino mpu6050 計測できない
-
C言語で1000までの素数を求める...
-
VB6 画面サイズ変更について
-
Labviewって・・・
-
signalハンドラに引数を渡したい
おすすめ情報