(*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.`}};
In[]:=
# test from WL to Python
from wolframclient.language import wl
data = wl.test
from wolframclient.language import wl
data = wl.test
In[]:=
data
Out[]=
{{120.41,23.9252,1.},{120.759,24.3829,1.},{121.452,24.9825,16.},{121.515,25.1059,13.},{121.513,25.0632,8.},{120.678,24.0996,5.},{121.202,25.0603,13.},{120.425,22.5641,9.},{120.338,22.5658,20.},{121.527,25.0624,15.},{120.333,22.6891,5.},{120.545,23.7119,6.},{121.53,25.0206,21.},{120.293,22.6749,2.},{121.205,24.9542,11.},{120.218,23.0483,1.},{120.247,23.4653,3.},{121.089,24.7406,11.},{120.617,24.1622,2.},{120.569,24.2256,0.},{120.641,24.152,3.},{121.366,25.0786,10.},{120.412,22.4795,6.},{121.459,25.013,19.},{121.579,25.05,14.},{120.288,22.6326,1.},{120.308,22.6054,3.},{120.685,23.913,12.},{120.488,22.6731,4.},{120.531,22.8836,4.},{120.82,24.5653,3.},{121.304,24.9957,15.},{121.76,25.1292,9.},{120.349,23.7575,2.},{121.449,25.1645,14.},{121.039,24.9001,5.},{121.481,25.069,12.},{120.972,24.8056,15.},{121.538,24.9772,25.},{121.433,25.038,15.},{120.346,23.5548,2.},{120.317,23.3056,3.},{120.328,22.7337,0.},{121.69,25.1797,10.},{121.508,25.0465,11.},{120.542,24.066,3.},{120.358,22.6274,8.},{120.561,22.5231,0.},{120.469,24.1317,1.},{120.306,22.7575,5.},{120.899,24.697,1.},{121.216,24.8639,5.},{120.742,24.2566,2.},{121.083,25.0355,12.},{120.203,23.7175,3.},{120.203,22.9846,2.}}
In[]:=
# test from WL to Python
from wolframclient.language import wl
data = wl.test
#Python Code
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(
data[:, 0],
data[:, 1],
data[:, 2],
variogram_model="spherical",
drift_terms=["regional_linear"]
)
z, ss = UK.execute("grid", gridx, gridy)
plt.imshow(z)
plt.show()
from wolframclient.language import wl
data = wl.test
#Python Code
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(
data[:, 0],
data[:, 1],
data[:, 2],
variogram_model="spherical",
drift_terms=["regional_linear"]
)
z, ss = UK.execute("grid", gridx, gridy)
plt.imshow(z)
plt.show()
Out[]=
Failure
