我有一個問題需要你的幫助。現在我有一個網格數據。它由多個不同的值組成。結構如下圖所示。 實際圖形事實上，我想要得到的圖形應該是由平滑的填充曲面組成的。如下所示。 我想要的圖形我提供了一個特定的數據集。請參考。 樣本數據集。這個數據集是“surfer”軟件定義的asc格式。您可以使用“pykrige”來閱讀。使用以下代碼進行渲染import sysimport matplotlib.pyplot as pltimport pykrige.kriging_tools as ktdef func(cmd): ascFile = cmd[1] ascFileContent = kt.read_asc_grid(ascFile) Z = ascFileContent[0] X = ascFileContent[1] Y = ascFileContent[2] print(Z) fig = plt.figure() ax = fig.add_subplot() contour = plt.contourf(X,Y,Z) #contour = plt.pcolormesh(X,Y,Z) plt.show()if __name__ == "__main__": cmd=["CONTOUR_KYZQ"] #gridFile cmd.append(r"D:\out.grd") func(cmd)我使用了各種方法，例如 pcolormesh、contourf。生成的圖形呈鋸齒狀。我怎么解決這個問題。請幫我。多謝。
1 回答
慕尼黑8549860
TA貢獻1818條經驗 獲得超11個贊
用于scipy.ndimage.zoom平滑圖像。這是我對數字 5 使用平滑，但您可以嘗試其他數字。
import matplotlib.pyplot as plt
import pykrige.kriging_tools as kt
import scipy.ndimage
ascFileContent = kt.read_asc_grid("out.grd")
Z = ascFileContent[0]
X = ascFileContent[1]
Y = ascFileContent[2]
fig = plt.figure()
plt.rcParams['figure.figsize'] = [20, 14]
new_z = scipy.ndimage.zoom(Z, 5)
contour = plt.contourf(np.linspace(X[0], X[-1], len(X) * 5),
np.linspace(Y[0], Y[-1], len(Y) * 5),
new_z)
plt.show()
如果您放大特定區域（如此代碼片段所示），您會注意到您的數據本質上是“方格的”。所以即使平滑了，也有一些區域看起來非常像網格（見下文）。
ascFileContent = kt.read_asc_grid("out.grd")
Z = ascFileContent[0]
X = ascFileContent[1]
Y = ascFileContent[2]
fig = plt.figure()
plt.rcParams['figure.figsize'] = [20, 14]
ZOOM_FACTOR = 5
new_z = scipy.ndimage.zoom(Z, ZOOM_FACTOR)
contour = plt.contourf(np.linspace(X[0], X[-1], len(X) * ZOOM_FACTOR),
np.linspace(Y[0], Y[-1], len(Y) * ZOOM_FACTOR),
new_z)
plt.xlim((106, 107))
plt.ylim((24, 26))
plt.show()
添加回答
舉報
0/150
提交
取消