Extract values at voxels within a segment

[muratmaga]

If you mean the physical (RAS) coordinates, I think you can rework this example:

Convert between RAS and numpy array indices

First, download MRHead data from Sample Data module import numpy as np volumeArray = slicer.util.array("MRHead") # Get the image data in numpy array format volumeArray[60, 127, 150] = 0 # Make one voxel black, numpy array is in KJI, not IJK. volumeNode = slicer.mrmlScene.GetFirstNodeByName("MRHead") volumeNode.Modified() # I think this is needed to let Slicer know the image was modified. ijkToRasMatrix = vtk.vtkMatrix4x4() volumeNode.GetIJKToRASMatrix(ijkToRasMatrix) ijkPoint = np.array([150…