set_solution_to_topography

set_solution_to_topography是Gempy中Solution类的一个方法。它将计算出的地质模型解决方案与地形进行合并，生成新的地形，用于表示地面形态。

输入参数

• set_solutions_to_towerpgraphy(self, lith_block=None, name='new_topography', **kwargs)：该方法有三个参数：
  • lith_block (Pandas dataframe): 包含地质模型解决方案的详细信息的Dataframe。
  • name (str): 新生成地形的名称。
  • **kwargs: 其他自定义参数，可选。

输出

该方法将返回包含新生成地形的gempy.SurfacePoints对象。

方法描述

该方法将计算出的地质模型解决方案从空间数据结构中提取出来，将其与地形数据集合并。因此，它输出包含地形数据和地质模型解决方案的新SurfacePoints对象。

注意事项

• 在使用此方法之前，必须先使用interpolate()方法来计算地质模型。
• 地质模型解决方案必须存储在gempy.Layer对象的lith_block属性中。

示例

import gempy as gp

# 计算地质模型
geo_model = gp.create_model('NewModel')
gp.init_data(geo_model, [0, 10, 0, 10, 0, 10], [50, 50, 50])
gp.add_surfaces(geo_model, ['rock1', 'rock2'])
gp.add_surface_points(geo_model, { 'X': [1, 5, 2, 8],
                                   'Y': [1, 2, 8, 7],
                                   'Z': [1, 7, 2, 5],
                                   'surface': ['rock1', 'rock1', 'rock2', 'rock2']})
gp.set_is_fault(geo_model, [0], rearrange=False)

gp.set_interpolation_data(geo_model, output=['geology'], compile_theano=True,
                          theano_optimizer='fast_compile', verbose=[])

sol = gp.compute_model(geo_model, compute_mesh=True)

# 将解决方案与地形合并
new_topography = sol.set_solution_to_topography(lith_block=sol._lith_block)