set_active
方法用于将指定grid的活动地层设置为指定阶段的地层。该函数可用于转换工作流程中的不同阶段。
Grid.set_active(stage: str, series: str, order: int)
stage
: str
series
: str
order
: int
# 导入Gempy并实例化Grid
import gempy as gp
geo_model = gp.create_model('Example')
# 在Grid对象上使用setDefaultFunctions方法
geo_model.grid.set_default_functions()
# 为Grid对象定义数据
gp.init_data(geo_model, [0, 10., 0, 10., 0, 10.], [50, 50, 50],
path_i='./data/input_data/simple_fault_model_orientations.csv',
path_o='./data/input_data/simple_fault_model_foliations.csv');
# 建模
gp.map_series_to_surfaces(geo_model, {"Fault_Series": 'fault',
"Strat_Series": ('rock2','rock1')},
remove_unused_series=True)
gp.set_interpolation_data(geo_model, compile_theano=True)
geo_model.surfaces.df
# 在第一个阶段第一系列的第二层将活动地层设置为'rock2'
geo_model.grid.set_active('first', 'Strat_Series', 2)
# 检查指定的活动地层是否正确
print(geo_model.surfaces.df)
# 输出:
# surface series order_surfaces color id
# 4 rock2 Strat_Series 2.0 #015482 5
# 5 rock1 Strat_Series 1.0 #9f0052 4
# 6 fault Fault_Series 1.0 #ffbe00 6
# 9 Basement Basement 1.0 #728f02 2
# 11 Rock2_Low Strat_Series 0.0 #7e1e9c 7