昔作ったCSVで3Dデータ

今はどこにデータがあるのか分からないのでrequestsライブラリを使ってデータを得る事が出来ないのですが、昔保存した上高地の地形データのCSVを表示するプログラムです。

from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt
import numpy as np
import csv

f = open("tikei.csv")
reader = csv.reader(f)
data = []
tmp = []
for row in reader:
    data.append(row)
    
data = np.array(data)
x = range(len(data[0]))
y = range(len(data))
x = np.array(x, dtype=float)
y = np.array(y, dtype=float)

X, Y = np.meshgrid(x, y)
Z = np.array(data, dtype=float)

fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
ax.plot_wireframe(X, Y, Z, color="#0000FF", rstride=1, cstride=1, linewidth=0.1)
ax.set_zlim(0, 3776)
plt.show()

とこんな感じでCSVの行を横、列を縦、値を高さとすることで地形データを表示するプログラムでした。
ちなみに3行追加すると等高線が作れます。

plt.contourf(Z, cmap="rainbow")
plt.colorbar()
plt.show()

という昔作ったプログラムでした。

ちなみに機械学習で決定境界を作るプログラムはこれを応用して作っています。