(*testdata*)​​test={{120.409653`,23.925175`,1.`},{120.758833`,24.382942`,1.`},{121.451861`,24.982528`,16.`},{121.5145`,25.105917`,13.`},{121.513311`,25.0632`,8.`},{120.677689`,24.099611`,5.`},{121.201811`,25.060344`,13.`},{120.425311`,22.564136`,9.`},{120.337736`,22.565833`,20.`},{121.526528`,25.062361`,15.`},{120.332631`,22.689056`,5.`},{120.544994`,23.711853`,6.`},{121.529556`,25.020608`,21.`},{120.292917`,22.674861`,2.`},{121.204981`,24.954208`,11.`},{120.218333`,23.048333`,1.`},{120.24735`,23.465308`,3.`},{121.088903`,24.740644`,11.`},{120.616917`,24.162197`,2.`},{120.568794`,24.225628`,0.`},{120.641092`,24.151958`,3.`},{121.365703`,25.07857`,10.`},{120.41175`,22.4795`,6.`},{121.458667`,25.012972`,19.`},{121.578611`,25.05`,14.`},{120.288086`,22.632567`,1.`},{120.307564`,22.605386`,3.`},{120.685306`,23.913`,12.`},{120.488033`,22.673081`,4.`},{120.530542`,22.883583`,4.`},{120.8202`,24.565269`,3.`},{121.304417`,24.995667`,15.`},{121.760056`,25.129167`,9.`},{120.348742`,23.757547`,2.`},{121.449239`,25.1645`,14.`},{121.038653`,24.900142`,5.`},{121.481028`,25.06895`,12.`},{120.972075`,24.805619`,15.`},{121.537778`,24.977222`,25.`},{121.4325`,25.037972`,15.`},{120.345531`,23.554839`,2.`},{120.31725`,23.305633`,3.`},{120.328289`,22.733667`,0.`},{121.689881`,25.179667`,10.`},{121.507972`,25.046503`,11.`},{120.541519`,24.066`,3.`},{120.358083`,22.627392`,8.`},{120.561175`,22.523108`,0.`},{120.469061`,24.131672`,1.`},{120.305689`,22.757506`,5.`},{120.898572`,24.696969`,1.`},{121.21635`,24.863869`,5.`},{120.741711`,24.256586`,2.`},{121.082761`,25.035503`,12.`},{120.202842`,23.717533`,3.`},{120.202617`,22.984581`,2.`}};
(*PythonCODE*)​​
〉
from pykrige.uk import UniversalKriging
import numpy as np
import matplotlib.pyplot as plt

gridx = np.arange(120.203, 122.76, 0.001)
gridy = np.arange(22.4795, 25.1797, 0.001)​
​
UK = UniversalKriging(
test[:, 0],
test[:, 1],
test[:, 2],
variogram_model="spherical",
drift_terms=["regional_linear"]
)
​
z, ss = UK.execute("grid", gridx, gridy)
plt.imshow(z)
plt.show()