lassoan

# This script can be used with this example scene:
# https://github.com/lassoan/PublicTestingData/releases/download/data/ColoredVolumeRenderingScene.mrb

volumeNode = getNode('Panoramix-cropped')
segmentationNode = getNode('Panoramix-cropped segmentation')
boostSegmentOpacity = False  # Increase voxel value inside segments to allow making them them more opaque

volumesLogic = slicer.modules.volumes.logic()

segmentIds = vtk.vtkStringArray()

lassoan

# Convert volume nodes of 3D Slicer to/from xarray and read/write into zarr format

def xarrayFromVolume(volumeNode) -> "xr.DataArray":
    """Convert an volume node to an xarray.DataArray.
    Origin and spacing metadata is preserved in the xarray's coords. The
    Direction is set in the `direction` attribute.
    Dims are labeled as `x`, `y`, `z`, `t`, and `c`.
    This interface is and behavior is experimental and is subject to possible
    future changes."""
    import xarray as xr

lassoan

# Create a markups fiducial node and place 6 points then copy paste this code into 3D Slicer's Python console.
# Two circles will be created (one for each 3 points) and the axis will be computed by connecting their centers.
# The circles and the axis is updated as the fiducial point positions are adjusted.

def sphereFrom3Points(markupsNode, startPointIndex):
    """Compute center and radius of 3-point sphere from 3 fiducial points
    source: https://stackoverflow.com/questions/20314306/find-arc-circle-equation-given-three-points-in-space-3d
    """
    import numpy as np
    A = np.zeros(3)

lassoan

The example was moved to Sandbox extension so that users can more easily install it:

https://github.com/PerkLab/SlicerSandbox/blob/master/LoadRemoteFile/LoadRemoteFile.py

More information: https://github.com/PerkLab/SlicerSandbox/blob/master/README.md#loadremotefile

lassoan

segmentationNode = slicer.mrmlScene.GetFirstNodeByClass('vtkMRMLSegmentationNode')
# set value to the size of the larges cracks in the segment surfaces
maximumHoleSizeMm = 2.0

############
masterVolumeNode = segmentationNode.GetNodeReference(segmentationNode.GetReferenceImageGeometryReferenceRole())

# Create segment editor to get access to effects
segmentEditorWidget = slicer.qMRMLSegmentEditorWidget()
# To show segment editor widget (useful for debugging):

lassoan

# This script randomly warps a 3D volume and adds random translations, rotations,
# and save each resulting 3D volume (and a screenshot for quick overview)
#
# The script can be executed by copy-pasting into 3D Slicer's Python console
# or in a Jupyter notebook running 3D Slicer kernel (provided by SlicerJupyter extension).
#
# Prerequisites:
# - Recent Slicer-4.11 version
# - SlicerIGT extension installed (for random deformations)

lassoan

# This gist contained a script that automatically removes table from a CT volume.

# The code was improved and turned into a Slicer module, available in the Sandbox extension:
# https://github.com/PerkLab/SlicerSandbox#remove-ct-table
#
# Source code:
# https://github.com/PerkLab/SlicerSandbox/blob/master/RemoveCtTable/RemoveCtTable.py#L251

lassoan

# Source data can be obtained from http://lifesciencedb.jp/bp3d/ (tested with partof_BP3D_4.0_obj_99)
# These input file must be loaded into the scene:
# - partof_inclusion_relation_list.txt and partof_element_parts.txt (loaded as Table nodes)
# - OBJ files

shNode = slicer.mrmlScene.GetSubjectHierarchyNode()
sceneItemID = shNode.GetSceneItemID()

def getItemParentsFmaIds(shNode, itemShItemId):
    existingParentShItemId = shNode.GetItemParent(itemShItemId)

lassoan

# Load dry bone CT of skull into the scene and run this script to automatically segment endocranium

masterVolumeNode = slicer.mrmlScene.GetFirstNodeByClass("vtkMRMLScalarVolumeNode")
smoothingKernelSizeMm = 3.0  # this is used for closing small holes in the se

# Compute bone threshold value automatically
import vtkITK
thresholdCalculator = vtkITK.vtkITKImageThresholdCalculator()
thresholdCalculator.SetInputData(masterVolumeNode.GetImageData())
thresholdCalculator.SetMethodToOtsu()

lassoan

# Load/generate input data
################################################

# Load master volume
import SampleData
sampleDataLogic = SampleData.SampleDataLogic()
masterVolumeNode = sampleDataLogic.downloadMRBrainTumor1()

# Define boundary points
import numpy as np