Rotated Block Model

This example demonstrates how to create a rotated block model.

import tempfile

from pathlib import Path
from parq_blockmodel import ParquetBlockModel
import pyvista as pv

Create a parquet block model

We will create a block model with a shape of 3x3x3 blocks, each block having a size of 1.0 in all dimensions. We will apply rotation.

corner = (0.0, 0.0, 0.0)
block_size = (1.0, 1.0, 1.0)
shape = (2, 2, 2)

temp_dir: Path = Path(tempfile.gettempdir()) / "block_model_example"

temp_dir.mkdir(parents=True, exist_ok=True)
# create a temporary file path for the block model

pbm: ParquetBlockModel = ParquetBlockModel.create_demo_block_model(filename=temp_dir / "demo_block_model.parquet",
                                                                   block_size=block_size, corner=corner, shape=shape,
                                                                   axis_azimuth=30, axis_dip=0, axis_plunge=0,
                                                                   )

pbm

ParquetBlockModel(name=demo_block_model, path=/tmp/block_model_example/demo_block_model.pbm.parquet)

Check the block model

print("Block Model Path:", pbm.blockmodel_path)
print("Name:", pbm.name)
print("Axis Angles:", pbm.geometry.axis_angles)
print("Data Shape:", pbm.data.shape)
print("Data Head:\n", pbm.data.head())
print("Model Attributes:", pbm.attributes)

Block Model Path: /tmp/block_model_example/demo_block_model.pbm.parquet
Name: demo_block_model
Axis Angles: (29.999999999999996, 0.0, 0.0)
Data Shape: (8, 4)
Data Head:
                         index_c  index_f  depth depth_category
x         y        z
 0.183013 0.683013 0.5        0        0    1.5           deep
                   1.5        1        4    0.5        shallow
-0.316987 1.549038 0.5        2        2    1.5           deep
                   1.5        3        6    0.5        shallow
 1.049038 1.183013 0.5        4        1    1.5           deep
Model Attributes: ['index_c', 'index_f', 'depth', 'depth_category']

Visualise

p: pv.Plotter = pbm.plot(scalar='depth', grid_type="image")
p.show()