update_default_kriging是Gempy中的一个函数,用于更新默认反距离权重克里金插值的参数。
def update_default_kriging(self, range_a=None, range_b=None, sill=None, nugget=None, K=None)
无返回值
可以通过传入参数,更新默认的反距离权重克里金插值方法的相关参数,从而轻松调整插值方法的效果。
import gempy as gp
geo_model = gp.create_model('example_model')
gp.init_data(geo_model, [0, 10, 0, 10, 0, 10], [50, 50, 50], default_density=3300)
gp.map_stack_to_surfaces(geo_model, {"Fault_Series": 'fault', "Strat_Series": ('sandstone_2', 'shale', 'sandstone_1', 'limestone', 'basement')}, remove_unused_series=True)
gp.set_interpolator(geo_model)
gp.update_default_kriging(range_a=10, range_b=10, sill=1, nugget=0.4, K=3)
上述代码中,我们首先创建了一个名为"example_model"的地质模型,然后初始化了该模型的数据,设置了默认的密度为3300。接着我们将地层系列映射为表面,并删除未使用的系列。接下来我们设置了该模型的插值器,并调用了update_default_kriging函数来更新默认反距离权重克里金插值方法的相关参数。这里我们指定了range_a的值为10,range_b的值为10,sill的值为1,nugget的值为0.4,K的值为3。
如果您需要更改插值方法的参数,请确保已经在映射地层系列和设置插值器之后调用了此函数,否则将无法生效。