def pointIsNear(locations, distance, inputs):
array = vtkDoubleArray()
array.SetNumberOfComponents(3)
array.SetNumberOfTuples(len(locations))
for i in range(len(locations)):
array.SetTuple(i, locations[i])
node = vtkSelectionNode()
node.SetFieldType(vtkSelectionNode.POINT)
node.SetContentType(vtkSelectionNode.LOCATIONS)
node.GetProperties().Set(vtkSelectionNode.EPSILON(), distance)
node.SetSelectionList(array)
from paraview.vtk.vtkFiltersExtraction import vtkLocationSelector
selector = vtkLocationSelector()
selector.Initialize(node, "vtkInsidedness")
inputDO = inputs[0].VTKObject
outputDO = inputDO.NewInstance()
outputDO.CopyStructure(inputDO)
output = dsa.WrapDataObject(outputDO)
if outputDO.IsA('vtkCompositeDataSet'):
it = inputDO.NewIterator()
it.InitTraversal()
while not it.IsDoneWithTraversal():
outputDO.SetDataSet(it, inputDO.GetDataSet(it).NewInstance())
it.GoToNextItem()
selector.ComputeSelectedElements(inputDO, outputDO)
return output.PointData.GetArray('vtkInsidedness')
def pointIsNear(locations, distance, inputs):
array = vtkDoubleArray()
array.SetNumberOfComponents(3)
array.SetNumberOfTuples(len(locations))
for i in range(len(locations)):
array.SetTuple(i, locations[i])
node = vtkSelectionNode()
node.SetFieldType(vtkSelectionNode.POINT)
node.SetContentType(vtkSelectionNode.LOCATIONS)
node.GetProperties().Set(vtkSelectionNode.EPSILON(), distance)
node.SetSelectionList(array)
from paraview.vtk.vtkFiltersExtraction import vtkLocationSelector
selector = vtkLocationSelector()
selector.SetInsidednessArrayName("vtkInsidedness")
selector.Initialize(node)
inputDO = inputs[0].VTKObject
outputDO = inputDO.NewInstance()
outputDO.CopyStructure(inputDO)
output = dsa.WrapDataObject(outputDO)
if outputDO.IsA('vtkCompositeDataSet'):
it = inputDO.NewIterator()
it.InitTraversal()
while not it.IsDoneWithTraversal():
outputDO.SetDataSet(it, inputDO.GetDataSet(it).NewInstance())
it.GoToNextItem()
selector.Execute(inputDO, outputDO)
return output.PointData.GetArray('vtkInsidedness')
def cellContainsPoint(inputs, locations):
array = vtkDoubleArray()
array.SetNumberOfComponents(3)
array.SetNumberOfTuples(len(locations))
for i in range(len(locations)):
array.SetTuple(i, locations[i])
node = vtkSelectionNode()
node.SetFieldType(vtkSelectionNode.CELL)
node.SetContentType(vtkSelectionNode.LOCATIONS)
node.SetSelectionList(array)
from paraview.vtk.vtkFiltersExtraction import vtkLocationSelector
selector = vtkLocationSelector()
selector.Initialize(node, "vtkInsidedness")
inputDO = inputs[0].VTKObject
outputDO = inputDO.NewInstance()
outputDO.CopyStructure(inputDO)
output = dsa.WrapDataObject(outputDO)
if outputDO.IsA('vtkCompositeDataSet'):
it = inputDO.NewIterator()
it.InitTraversal()
while not it.IsDoneWithTraversal():
outputDO.SetDataSet(it, inputDO.GetDataSet(it).NewInstance())
it.GoToNextItem()
selector.ComputeSelectedElements(inputDO, outputDO)
return output.CellData.GetArray('vtkInsidedness')
