tương tự với scipy.interpolate.griddata?

Question

Tôi muốn nội suy một đám mây điểm 3D đã cho:

Tôi đã xem xét scipy.interpolate.griddata và kết quả chính xác là những gì tôi cần, nhưng như tôi hiểu, tôi cần nhập "griddata" như x = [[0,0,0],[1,1,1],[2,2,2]].

Nhưng đám mây điểm 3D của tôi không có giao diện lưới này - Giá trị x, y không hoạt động như lưới - dù sao chỉ có một giá trị z duy nhất cho mỗi giá trị x, y. *

Vì vậy, có một giải pháp thay thế cho scipy.interpolate.griddata cho đám mây không-trong-lưới-điểm của tôi?

* chỉnh sửa: "không nhìn lưới" có nghĩa là đầu vào của tôi trông như thế này:

x = [0,4,17] 
y = [-7,25,116] 
z = [50,112,47] 

Answer 1

Đây là một chức năng tôi sử dụng cho loại công cụ này:

from numpy import linspace, meshgrid 

def grid(x, y, z, resX=100, resY=100): 
    "Convert 3 column data to matplotlib grid" 
    xi = linspace(min(x), max(x), resX) 
    yi = linspace(min(y), max(y), resY) 
    Z = griddata(x, y, z, xi, yi) 
    X, Y = meshgrid(xi, yi) 
    return X, Y, Z 

Sau đó, sử dụng nó như này:

X, Y, Z = grid(x, y, z) 

Answer 2

Scipy có tài liệu với ví dụ cụ thể về cách sử dụng scipy.interpo late.griddata và họ giải thích chính xác những gì bạn đang yêu cầu. Nhìn đây: http://docs.scipy.org/doc/scipy/reference/generated/scipy.interpolate.griddata.html

Nói tóm lại, bạn làm điều này có được "dữ liệu lưới":

grid_x, grid_y = np.mgrid[0:1:100j, 0:1:200j] 

này sẽ làm cho một mạng lưới 100x200 trong đó khoảng 0-1 ở cả x và y hướng .

grid_x, grid_y = np.mgrid[-10:10:51j, 0:2:20j] 

Điều này sẽ tạo ra lưới 51x20 nằm trong khoảng từ -10 đến 10 theo hướng x và từ 0 đến 2 theo hướng y.

Bây giờ bạn phải sửa đầu vào cho scipy.interpolate.griddata.