Я не нашел ответа на этот вопрос: у меня есть сетка, определенная в текстовом файле с четырьмя столбцами: (долгота, широта, глубина, скольжение).Каждая строка является точкой сетки.
Я могу сгенерировать точечную диаграмму этих точек, используя следующий простой код:
# Main imports:
import numpy as np
from pylab import *
from mpl_toolkits.mplot3d import Axes3D
# Read the grid:
points = np.loadtxt("grid.txt")
# Retrieve parameters from the grid:
lon = points[:,0]
lat = points[:,1]
depth = points[:,2]
slip = points[:,3]
# 3-D plot of the model:
fig = figure(1)
ax = fig.add_subplot(111, projection='3d')
p = ax.scatter(lon, lat, depth, c=slip, vmin=0, vmax=max(slip), s=30, edgecolor='none', marker='o')
fig.colorbar(p)
title("Published finite fault in 3-D")
ax.set_xlabel("Longitude [degrees]")
ax.set_ylabel("Latitude [degrees]")
ax.set_zlabel("Depth [km]")
ax.invert_zaxis()
jet()
grid()
show()
И я получаю следующий рисунок:
Что я хочу сделать, так это уметь интерполировать эти точки, чтобы создать «непрерывную» сетку поверхности и отобразить ее на двухмерных и трехмерных графиках.Поэтому как-то я должен учитывать все (долг, широта, глубина, скольжение) в интерполяции.Буду признателен за ваши предложения.Заранее спасибо!