使用 Julia 中的 GeoStats 包进行 3D 克里金插值的断言错误

问题描述

我正在尝试使用 GeoStats.jl 包构建克里金曲面插值示例。 我提供了 x 和 y 坐标作为输入，z 坐标作为值来训练函数。然后我尝试预测其他输入集的 z 值。 目标是预测拓扑插值，示例代码如下：

using GeoStats,DataFrames
using Variography,Plots

#create your interpolator object
OK = OrdinaryKriging(GaussianVariogram())

# fit it to the data:

x_train = [1.0,1.0,10.0,2.5,7.5] # x coordinates train
y_train = [1.0,5.0,5.0] # y coordinates train
z = [2.0,2.0,4.5,4.5] # z coordinates train

table_train = DataFrame(:x=>x_train,:y=>y_train)

? = georef(table_train,[x y]')

scatter(?,z,label="train points")

krig = KrigingEstimators.fit(OK,?,z) # fit function

# predictions:

x_test = [1.5,1.5,11.0,3.5,5.5] # x coordinates test
y_test = [1.4,10.5,10.6,4.0,6.0,8.0] # y coordinates test

table_test = DataFrame(:xt=>x_test,:yt=>y_test)

test = georef(table,[xt yt]')

result = []
variance =[]

for i in test
    μ,σ² = KrigingEstimators.predict(krig,[i])
    push!(result,μ)
    push!(variance,σ²)
end
df_krig_vario = DataFrame(:predict=>result,:real=>y_test,:variance=>variance)
println(first(df_krig_vario,5))
mean_var = sum(variance)/length(variance)
println("")
println("mean variance is $mean_var")
test = reshape(test,length(test),1)
plot!(test,y_test,label="actual")
plot!(test,result,label="predict",legend=:bottomright,title="Gaussian Variogram")

但是，当我运行此代码时，它会引发以下错误：

AssertionError: invalid variable name

Stacktrace:
 [1] macro expansion at C:\Users\User\.julia\packages\GeoStatsBase\aI9iy\src\plotrecipes\data.jl:14 [inlined]
 [2] apply_recipe(::AbstractDict{Symbol,Any},::AbstractData,::Any) at C:\Users\User\.julia\packages\RecipesBase\92zOw\src\RecipesBase.jl:282
 [3] _process_userrecipes!(::Any,::Any,::Any) at C:\Users\User\.julia\packages\RecipesPipeline\VEk89\src\user_recipe.jl:36
 [4] recipe_pipeline!(::Any,::Any) at C:\Users\User\.julia\packages\RecipesPipeline\VEk89\src\RecipesPipeline.jl:70
 [5] _plot!(::Plots.Plot,::Any) at C:\Users\User\.julia\packages\Plots\KN96c\src\plot.jl:172
 [6] #plot#129 at C:\Users\User\.julia\packages\Plots\KN96c\src\plot.jl:58 [inlined]
 [7] scatter(::Any,::Vararg{Any,N} where N; kw::Any) at C:\Users\User\.julia\packages\RecipesBase\92zOw\src\RecipesBase.jl:403
 [8] top-level scope at In[58]:18
 [9] include_string(::Function,::Module,::String,::String) at .\loading.jl:1091

我可以知道如何解决这个问题吗？

谢谢！

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）