def InterpolatePatchCenterlines(patchCenterlines, parentCenterlines,
additionalPoint, lower, version):
if additionalPoint is not None:
additionalPointIds = []
for i in range(parentCenterlines.GetNumberOfCells()):
line = ExtractSingleLine(parentCenterlines, i)
additionalPointIds.append(line.FindPoint(additionalPoint))
else:
additionalPointIds = [
"" for i in range(parentCenterlines.GetNumberOfCells())
]
interpolatedLines = vtk.vtkPolyData()
interpolatedPoints = vtk.vtkPoints()
interpolatedCellArray = vtk.vtkCellArray()
pointsInserted = 0
interpolatedCellArray.Initialize()
for i in range(parentCenterlines.GetNumberOfCells()):
startingCell = vtk.vtkGenericCell()
endingCell = vtk.vtkGenericCell()
numberOfInterpolationPoints = parentCenterlines.GetCell(
i).GetNumberOfPoints()
patchCenterlines.GetCell(0, startingCell)
patchCenterlines.GetCell(i + 1, endingCell)
if version:
splinePoints = InterpolateSpline(startingCell, endingCell,
additionalPoint)
else:
splinePoints = InterpolateTwoCells(startingCell, endingCell, \
numberOfInterpolationPoints, \
additionalPointIds[i],
additionalPoint, lower)
interpolatedCellArray.InsertNextCell(splinePoints.GetNumberOfPoints())
for j in range(splinePoints.GetNumberOfPoints()):
interpolatedPoints.InsertNextPoint(splinePoints.GetPoint(j))
interpolatedCellArray.InsertCellPoint(pointsInserted + j)
pointsInserted += splinePoints.GetNumberOfPoints()
interpolatedLines.SetPoints(interpolatedPoints)
interpolatedLines.SetLines(interpolatedCellArray)
attributeFilter = vtkvmtk.vtkvmtkCenterlineAttributesFilter()
attributeFilter.SetInput(interpolatedLines)
attributeFilter.SetAbscissasArrayName(AbscissasArrayName)
attributeFilter.SetParallelTransportNormalsArrayName(
parallelTransportNormalsArrayName)
attributeFilter.Update()
attributeInterpolatedLines = attributeFilter.GetOutput()
return attributeInterpolatedLines
def InterpolatePatchCenterlines(patchCenterlines, parentCenterlines):
additionalPoint = [-1.0, -1.0, -1.0]
additionalPointIds = []
if (useAdditionalInterpolationPoint == 1):
additionalPoint = divergingPoints.GetPoint(0)
line1 = ExtractSingleLine(parentCenterlines, 0)
line2 = ExtractSingleLine(parentCenterlines, 1)
additionalPointIds.append(line1.FindPoint(additionalPoint))
additionalPointIds.append(line2.FindPoint(additionalPoint))
else:
for i in range(parentCenterlines.GetNumberOfCells()):
additionalPoint = clippingPoints.GetPoint(0)
line1 = ExtractSingleLine(parentCenterlines, 0)
additionalPointIds.append(line1.FindPoint(additionalPoint))
interpolatedLines = vtk.vtkPolyData()
interpolatedPoints = vtk.vtkPoints()
interpolatedCellArray = vtk.vtkCellArray()
pointsInserted = 0
interpolatedCellArray.Initialize()
for i in range(parentCenterlines.GetNumberOfCells()):
startingCell = vtk.vtkGenericCell()
endingCell = vtk.vtkGenericCell()
numberOfInterpolationPoints = parentCenterlines.GetCell(
i).GetNumberOfPoints()
patchCenterlines.GetCell(0, startingCell)
patchCenterlines.GetCell(i + 1, endingCell)
splinePoints = InterpolateTwoCells(startingCell, endingCell,
numberOfInterpolationPoints,
additionalPointIds[i],
additionalPoint)
interpolatedCellArray.InsertNextCell(splinePoints.GetNumberOfPoints())
for j in range(splinePoints.GetNumberOfPoints()):
interpolatedPoints.InsertNextPoint(splinePoints.GetPoint(j))
interpolatedCellArray.InsertCellPoint(pointsInserted + j)
pointsInserted += splinePoints.GetNumberOfPoints()
interpolatedLines.SetPoints(interpolatedPoints)
interpolatedLines.SetLines(interpolatedCellArray)
attributeFilter = vtkvmtk.vtkvmtkCenterlineAttributesFilter()
attributeFilter.SetInputData(interpolatedLines)
attributeFilter.SetAbscissasArrayName(abscissasArrayName)
attributeFilter.SetParallelTransportNormalsArrayName(
parallelTransportNormalsArrayName)
attributeFilter.Update()
attributeInterpolatedLines = attributeFilter.GetOutput()
return attributeInterpolatedLines
def InterpolatePatchCenterlines(patchCenterlines, parentCenterlines):
additionalPoint = [-1.0, -1.0, -1.0]
additionalPointIds = []
if useAdditionalInterpolationPoint == 1:
additionalPoint = divergingPoints.GetPoint(0)
line1 = ExtractSingleLine(parentCenterlines, 0)
line2 = ExtractSingleLine(parentCenterlines, 1)
additionalPointIds.append(line1.FindPoint(additionalPoint))
additionalPointIds.append(line2.FindPoint(additionalPoint))
else:
for i in range(parentCenterlines.GetNumberOfCells()):
additionalPoint = clippingPoints.GetPoint(0)
line1 = ExtractSingleLine(parentCenterlines, 0)
additionalPointIds.append(line1.FindPoint(additionalPoint))
interpolatedLines = vtk.vtkPolyData()
interpolatedPoints = vtk.vtkPoints()
interpolatedCellArray = vtk.vtkCellArray()
pointsInserted = 0
interpolatedCellArray.Initialize()
for i in range(parentCenterlines.GetNumberOfCells()):
startingCell = vtk.vtkGenericCell()
endingCell = vtk.vtkGenericCell()
numberOfInterpolationPoints = parentCenterlines.GetCell(i).GetNumberOfPoints()
patchCenterlines.GetCell(0, startingCell)
patchCenterlines.GetCell(i + 1, endingCell)
splinePoints = InterpolateTwoCells(
startingCell, endingCell, numberOfInterpolationPoints, additionalPointIds[i], additionalPoint
)
interpolatedCellArray.InsertNextCell(splinePoints.GetNumberOfPoints())
for j in range(splinePoints.GetNumberOfPoints()):
interpolatedPoints.InsertNextPoint(splinePoints.GetPoint(j))
interpolatedCellArray.InsertCellPoint(pointsInserted + j)
pointsInserted += splinePoints.GetNumberOfPoints()
interpolatedLines.SetPoints(interpolatedPoints)
interpolatedLines.SetLines(interpolatedCellArray)
attributeFilter = vtkvmtk.vtkvmtkCenterlineAttributesFilter()
attributeFilter.SetInput(interpolatedLines)
attributeFilter.SetAbscissasArrayName(abscissasArrayName)
attributeFilter.SetParallelTransportNormalsArrayName(parallelTransportNormalsArrayName)
attributeFilter.Update()
attributeInterpolatedLines = attributeFilter.GetOutput()
return attributeInterpolatedLines
def Execute(self):
if self.Centerlines == None:
self.PrintError('Error: No input centerlines.')
centerlineAttributes = vtkvmtk.vtkvmtkCenterlineAttributesFilter()
centerlineAttributes.SetInputData(self.Centerlines)
centerlineAttributes.SetAbscissasArrayName(self.AbscissasArrayName)
centerlineAttributes.SetParallelTransportNormalsArrayName(self.NormalsArrayName)
centerlineAttributes.Update()
self.Centerlines = centerlineAttributes.GetOutput()
def Execute(self):
if self.Centerlines == None:
self.PrintError('Error: No input centerlines.')
centerlineAttributes = vtkvmtk.vtkvmtkCenterlineAttributesFilter()
centerlineAttributes.SetInputData(self.Centerlines)
centerlineAttributes.SetAbscissasArrayName(self.AbscissasArrayName)
centerlineAttributes.SetParallelTransportNormalsArrayName(self.NormalsArrayName)
centerlineAttributes.Update()
self.Centerlines = centerlineAttributes.GetOutput()