8

我正在尝试将一些数据放入 imshow() 图中。我的问题是数据不是作为 MxN 数组,而是作为 3xN 数组(x 和 y 坐标和值)。这些点不是按规则网格排列的,而是位于 [xmin,xmax,ymin 和 ymax]=[-pi/2,pi/2,0,3.5] 内。

In [117]: shape(data)
Out[117]: (3L, 102906L)

如何从该数据中获得漂亮的图像图?非常感谢您的帮助。

顺便说一句,数据表示杆表面上的温度值作为轴向和方位角位置的函数,想想 cfd 网格。

4

1 回答 1

12

我建议使用 griddata-method 进行插值。一个样本是:

import numpy as np
from matplotlib.mlab import griddata
import matplotlib.pyplot as plt

xs0 = np.random.random((1000)) * np.pi - np.pi/2
ys0 = np.random.random((1000)) * 3.5
zs0 = np.random.random((1000))

N = 30j
extent = (-np.pi/2,np.pi/2,0,3.5)

xs,ys = np.mgrid[extent[0]:extent[1]:N, extent[2]:extent[3]:N]

resampled = griddata(xs0, ys0, zs0, xs, ys)

plt.imshow(resampled.T, extent=extent)
plt.plot(xs0, ys0, "r.")
plt.plot(xs, ys, "b.")
plt.title("imshow for irregularly spaced data using griddata")
plt.show()

我想从您的 3*X 阵列到三个 X 阵列的过渡是显而易见的。

结果是:

样本

红点表示数据的“原始”位置,蓝点表示现在规则间隔的数据。

griddata 返回一个掩码数组。无法评估插值的所有点都被屏蔽,然后绘制为白色区域。

HTH,托尔斯滕

于 2013-01-03T14:11:57.463 回答