numpy.meshgrid

numpy.meshgrid(*xi, **kwargs) [source]

Возвращает матрицы координат из векторов координат.

Создаёт N-мерные массивы координат для векторных вычислений N-мерных скалярных/векторных полей по N-мерным сеткам, заданным одномерными массивами координат x1, x2,…, xn.

Изменено в версии 1.9: Разрешены 1-мерные и 0-мерные случаи.

Параметры:

x1, x2,…, xn : array_like Одномерные массивы, представляющие координаты сетки. indexing : {‘xy’, ‘ij’}, optional Индексация декартовых координат (‘xy’, по умолчанию) или матричная индексация (‘ij’). Подробности см. в примечаниях. Новое в версии 1.7.0. sparse : bool, optional Если True, возвращается разреженная сетка для экономии памяти. По умолчанию False. Новое в версии 1.7.0. copy : bool, optional Если False, возвращаются представления исходных массивов для экономии памяти. По умолчанию True. Обратите внимание, что sparse=False, copy=False , вероятно, вернёт несмежные массивы. Кроме того, более чем один элемент широковещательного массива может ссылаться на одно местоположение в памяти. Если вам нужно записывать в массивы, сделайте копии предварительно. Новое в версии 1.7.0.

Возвращает:

X1, X2,…, XN : ndarray Для векторов x1, x2,…, ‘xn’ с длинами Ni=len(xi), возвращает массивы формы (N1, N2, N3,...Nn), если indexing=’ij’, или массивы формы (N2, N1, N3,...Nn), если indexing=’xy’. Элементы xi повторяются для заполнения матрицы вдоль первой размерности для x1, второй для x2 и так далее.

См. также

index_tricks.mgrid

Строит многомерную «сетку» с использованием индексной нотации.

index_tricks.ogrid

Строит открытую многомерную «сетку» с использованием индексной нотации.

Примечания

Эта функция поддерживает обе схемы индексации через ключевое слово indexing. Передача строки ‘ij’ возвращает сетку с матричной индексацией, а ‘xy’ — с декартовой индексацией. В двумерном случае с входами длины M и N выходы имеют форму (N, M) для индексации ‘xy’ и (M, N) для индексации ‘ij’. В трёхмерном случае с входами длины M, N и P выходы имеют форму (N, M, P) для индексации ‘xy’ и (M, N, P) для индексации ‘ij’. Разница проиллюстрирована следующим фрагментом кода:

xv, yv = np.meshgrid(x, y, sparse=False, indexing='ij')
for i in range(nx):
    for j in range(ny):
        # treat xv[i,j], yv[i,j]

xv, yv = np.meshgrid(x, y, sparse=False, indexing='xy')
for i in range(nx):
    for j in range(ny):
        # treat xv[j,i], yv[j,i]

В 1-мерном и 0-мерном случае ключи indexing и sparse не имеют эффекта.

Примеры

>>> nx, ny = (3, 2)
>>> x = np.linspace(0, 1, nx)
>>> y = np.linspace(0, 1, ny)
>>> xv, yv = np.meshgrid(x, y)
>>> xv
array([[0. , 0.5, 1. ],
       [0. , 0.5, 1. ]])
>>> yv
array([[0.,  0.,  0.],
       [1.,  1.,  1.]])
>>> xv, yv = np.meshgrid(x, y, sparse=True)  # make sparse output arrays
>>> xv
array([[0. ,  0.5,  1. ]])
>>> yv
array([[0.],
       [1.]])

meshgrid очень полезен для вычисления функций на сетке.

>>> import matplotlib.pyplot as plt
>>> x = np.arange(-5, 5, 0.1)
>>> y = np.arange(-5, 5, 0.1)
>>> xx, yy = np.meshgrid(x, y, sparse=True)
>>> z = np.sin(xx**2 + yy**2) / (xx**2 + yy**2)
>>> h = plt.contourf(x,y,z)
>>> plt.show()