{
"cells": [
{
"cell_type": "markdown",
"id": "early-wonder",
"metadata": {},
"source": [
"# Isosurfaces\n",
"\n",
"![](\"./images/iso_surface_app.png\")
\n",
"\n",
"With this application, users can create 3D isosurfaces around data values stored on any type of object.\n",
"\n",
"The application leverages the [Scikit-Image.measure.marching_cube](https://scikit-image.org/docs/dev/api/skimage.measure.html?highlight=marching_cubes#marching-cubes) routine.\n",
"\n",
"\n",
" \n",
"New user? Visit the [Getting Started](../installation.rst) page."
]
},
{
"cell_type": "markdown",
"id": "handy-detroit",
"metadata": {},
"source": [
"## Application\n",
"\n",
"The following sections provide details on the different parameters controlling the application. Interactive widgets shown below are for demonstration purposes only."
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "liberal-jacob",
"metadata": {
"scrolled": false
},
"outputs": [
{
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "a61a7f4df93f4baca35ce19dc7013560",
"version_major": 2,
"version_minor": 0
},
"text/plain": [
"HBox(children=(VBox(children=(VBox(children=(Label(value='Workspace', style=DescriptionStyle(description_width…"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"from geoapps.iso_surfaces.application import IsoSurface\n",
"\n",
"app = IsoSurface(geoh5=r\"../../../assets/FlinFlon.geoh5\")\n",
"app()"
]
},
{
"cell_type": "markdown",
"id": "suffering-platform",
"metadata": {},
"source": [
"## Project Selection\n",
"\n",
"Select and connect to an existing **geoh5** or **ui.json** project file containing data. "
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "internal-pixel",
"metadata": {},
"outputs": [
{
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "d66ceb63672a421985558c382a28ddcd",
"version_major": 2,
"version_minor": 0
},
"text/plain": [
"VBox(children=(Label(value='Workspace', style=DescriptionStyle(description_width='initial')), HBox(children=(F…"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"app.project_panel"
]
},
{
"cell_type": "markdown",
"id": "stretch-company",
"metadata": {},
"source": [
"See the [Project Panel](base_application.ipynb#Project-Panel) page for more details."
]
},
{
"cell_type": "markdown",
"id": "instructional-peace",
"metadata": {},
"source": [
"## Object/Data Selection\n",
"\n",
"List of objects and associated data pulled from the target workspace. The data fields are transfered to the newly created surface."
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "athletic-archives",
"metadata": {},
"outputs": [
{
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "910cfab0f5054010af4993befb57ef8b",
"version_major": 2,
"version_minor": 0
},
"text/plain": [
"VBox(children=(Dropdown(description='Object:', index=6, options=(['', None], ['Labels/Label3', UUID('39b2d390-…"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"app.data_panel"
]
},
{
"cell_type": "markdown",
"id": "according-smell",
"metadata": {},
"source": [
"## Iso-values\n",
"\n",
"Set the values to create a surface around. Define a list of fixed contours with the `Fixed contours` parameter. Define a range of contours with the `Contour min`, `Contour max`, and `Contour spacing` parameters."
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "organic-brunswick",
"metadata": {},
"outputs": [
{
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "136fd41b25364fffb0b9907c7f3e320c",
"version_major": 2,
"version_minor": 0
},
"text/plain": [
"VBox(children=(FloatText(value=0.005, description='Contour min:', style=DescriptionStyle(description_width='in…"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"app.contours"
]
},
{
"cell_type": "markdown",
"id": "official-advice",
"metadata": {},
"source": [
"## Max Interpolation Distance\n",
"\n",
"Maximum 3D radial distance between the vertices/cells of the input and output objects. No triangulation occurs beyont that distance."
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "impossible-still",
"metadata": {},
"outputs": [
{
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "633ac06d95ea47fc9b51ba379176e451",
"version_major": 2,
"version_minor": 0
},
"text/plain": [
"FloatText(value=500.0, description='Max Interpolation Distance (m):', style=DescriptionStyle(description_width…"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"app.max_distance"
]
},
{
"cell_type": "markdown",
"id": "center-arizona",
"metadata": {},
"source": [
"## Base Grid Resolution\n",
"\n",
"For all object types other than `BlockModel`, data values are first interpolated using an inverse distance averaging onto a uniform 3D grid, as required by the [Scikit-Image.measure.marching_cube](https://scikit-image.org/docs/dev/api/skimage.measure.html?highlight=marching_cubes#marching-cubes) routine. The `resolution` parameter controls the cell size of the underlying grid used for the interpolation. \n",
"\n",
"For more control on the iso-surface creation, users have the option to interpolate the data to a `BlockModel` as a pre-processing step using the [Data Transfer](data_interpolation.ipynb) application."
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "palestinian-choir",
"metadata": {},
"outputs": [
{
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "e6b6fec4c78d444e8834b51c420225d8",
"version_major": 2,
"version_minor": 0
},
"text/plain": [
"FloatText(value=50.0, description='Base grid resolution (m):', style=DescriptionStyle(description_width='initi…"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"app.resolution"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "distributed-fortune",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Starting the isosurface creation.\n",
"Isosurface completed. -> 5 surface(s) created.\n"
]
},
{
"data": {
"text/html": [
"