def serialize_grid(repo):
"""
This function creates a one cell grid + a property on a partial ijk grid in the repo to be serialized
"""
earth_model = repo.createEarthModel("f2060ce0-fa3d-11e5-8620-0002a5d5c51b",
"Grid")
earth_model_interp = repo.createEarthModelInterpretation(
earth_model, "f5cd7520-fa3d-11e5-b65b-0002a5d5c51b", "Grid interp")
# ONE CELL GRID
single_cell_ijk_grid = repo.createIjkGridExplicitRepresentation(
earth_model_interp, "e69bfe00-fa3d-11e5-b5eb-0002a5d5c51b",
"One unfaulted sugar cube", 1, 1, 1)
python_points = [
0.0, 0.0, 300.0, 700.0, 0.0, 350.0, 0.0, 150.0, 300.0, 700.0, 150.0,
350.0, 0.0, 0.0, 500.0, 700.0, 0.0, 550.0, 0.0, 150.0, 500.0, 700.0,
150.0, 550.0
]
resqml_points = fesapi.DoubleArray(24)
for i, value in enumerate(python_points):
resqml_points.setitem(i, value)
single_cell_ijk_grid.setGeometryAsCoordinateLineNodes(
fesapi.resqml20__PillarShape_vertical,
fesapi.resqml20__KDirection_down, False, resqml_points)
# ONE CONTINUOUS PROPERTY ON A PARTIAL 3 cells I=1 J=1 K=3 GRID
partial_ijk_grid = repo.createPartialIjkGridRepresentation(
"bc92a216-aa17-4a1f-9253-8b3ab094bf84", "partial grid")
continuous_prop = repo.createContinuousProperty(
partial_ijk_grid, "827d1074-5b16-4c7d-80bb-a8f3e1413753",
"My length property", 1, fesapi.resqml22__IndexableElement_cells,
fesapi.resqml20__ResqmlUom_m,
fesapi.resqml20__ResqmlPropertyKind_length)
prop_values = [1.1, 2.2, 3.3]
resqml_values = fesapi.FloatArray(3)
for i, value in enumerate(prop_values):
resqml_values.setitem(i, value)
continuous_prop.pushBackFloatHdf5Array3dOfValues(resqml_values, 1, 1, 3,
1.1, 3.3)
def serialize_grid(repo):
"""
This function create some grids in the repo to be serialized
"""
earth_model = repo.createEarthModel("f2060ce0-fa3d-11e5-8620-0002a5d5c51b",
"Grid")
earth_model_interp = repo.createEarthModelInterpretation(
earth_model, "f5cd7520-fa3d-11e5-b65b-0002a5d5c51b", "Grid interp")
# ONE SUGAR
single_cell_ijk_grid = repo.createIjkGridExplicitRepresentation(
earth_model_interp, "e69bfe00-fa3d-11e5-b5eb-0002a5d5c51b",
"One unfaulted sugar cube", 1, 1, 1)
python_points = [
0.0, 0.0, 300.0, 700.0, 0.0, 350.0, 0.0, 150.0, 300.0, 700.0, 150.0,
350.0, 0.0, 0.0, 500.0, 700.0, 0.0, 550.0, 0.0, 150.0, 500.0, 700.0,
150.0, 550.0
]
resqml_points = fesapi.DoubleArray(24)
for i, value in enumerate(python_points):
resqml_points.setitem(i, value)
single_cell_ijk_grid.setGeometryAsCoordinateLineNodes(
fesapi.resqml20__PillarShape_vertical,
fesapi.resqml20__KDirection_down, False, resqml_points)
def show_ijk_grid(ijk_grid):
"""
This function shows some ijk grids which have been deserialized
"""
show_all_metadata(ijk_grid)
if ijk_grid.getGeometryKind(
) != fesapi.resqml2_AbstractIjkGridRepresentation.NO_GEOMETRY:
most_complex_pillar_geom = ijk_grid.getMostComplexPillarGeometry()
if most_complex_pillar_geom == fesapi.resqml20__PillarShape_vertical:
print("Most complex pillar geometry is vertical")
else:
if most_complex_pillar_geom == fesapi.resqml20__PillarShape_straight:
print("Most complex pillar geometry is straight")
else:
if most_complex_pillar_geom == fesapi.resqml20__PillarShape_curved:
print("Most complex pillar geometry is curved")
k_gap_count = ijk_grid.getKGapsCount()
print("K Gap Count = ", k_gap_count)
if ijk_grid.getGeometryKind(
) == fesapi.resqml2_AbstractIjkGridRepresentation.LATTICE:
print("This 3d grid has a lattice geometry.")
else:
if ijk_grid.getGeometryKind() == \
fesapi.resqml2_AbstractIjkGridRepresentation.PARAMETRIC:
print("This 3d grid has a parametric geometry.")
else:
if ijk_grid.getGeometryKind() == \
fesapi.resqml2_AbstractIjkGridRepresentation.EXPLICIT:
print("This 3d grid has an explicit geometry.")
else:
print("This 3d grid has an unknown geometry.")
patch_count = ijk_grid.getPatchCount()
print("Patch Count =", patch_count)
nb_xyz_points = ijk_grid.getXyzPointCountOfAllPatches()
print("XYZ points count :", nb_xyz_points)
for patch_index in range(patch_count):
nb_xyz_points_in_patch = ijk_grid.getXyzPointCountOfPatch(
patch_index)
xyz_points_in_patch = fesapi.DoubleArray(nb_xyz_points_in_patch *
3)
ijk_grid.getXyzPointsOfPatch(patch_index, xyz_points_in_patch)
for vertex_index in range(nb_xyz_points_in_patch):
x = xyz_points_in_patch.getitem(vertex_index * 3)
y = xyz_points_in_patch.getitem(vertex_index * 3 + 1)
z = xyz_points_in_patch.getitem(vertex_index * 3 + 2)
print("Vertex ", vertex_index, " : ", x, " ", y, " ", z)
xyz_points = fesapi.DoubleArray(nb_xyz_points * 3)
ijk_grid.getXyzPointsOfAllPatches(xyz_points)
ijk_grid.loadSplitInformation()
for cell_corner in range(8):
pt_index = ijk_grid.getXyzPointIndexFromCellCorner(
0, 0, 0, cell_corner)
print("Cell 0,0,0 corner ", cell_corner, " is at index ", pt_index)
print("Cell 0,0,0 corner ", cell_corner, " is ",
xyz_points.getitem(pt_index * 3), " ",
xyz_points.getitem(pt_index * 3 + 1), " ",
xyz_points.getitem(pt_index * 3 + 2))
ijk_grid.unloadSplitInformation()
else:
print("This 3d grid has no geometry.")