print('With Shrink Filter (NS)')
# En 3D, le shrink ne doit pas se faire sur la geometrie car elle ne represente que la peau
shrink = svtk.svtkShrinkFilter()
shrink.SetInputConnection(geometry.GetOutputPort())
shrink.SetShrinkFactor(.8)
else:
print('No Shrink Filter (NS)')
shrink = geometry
# LookupTable
lut = svtk.svtkLookupTable()
lut.SetHueRange(0.66, 0)
lut.Build()
# Mappers
mapper = svtk.svtkDataSetMapper()
mapper.SetInputConnection(shrink.GetOutputPort())
shrink.Update()
dataRange = shrink.GetOutput().GetCellData().GetArray('scalar').GetRange()
print('DataRange (after shrink): ', dataRange)
mapper.SetLookupTable(lut)
mapper.SetColorModeToMapScalars()
mapper.SetScalarModeToUseCellFieldData()
mapper.SelectColorArray('scalar')
mapper.SetScalarRange(dataRange[0], dataRange[1])
# Actors
actor = svtk.svtkActor()
actor.SetMapper(mapper)
t1 = svtk.svtkMatrixToLinearTransform() m1 = svtk.svtkMatrix4x4() t1.SetInput(m1) m1.SetElement(0, 0, 1.127631) m1.SetElement(0, 1, 0.205212) m1.SetElement(0, 2, -0.355438) m1.SetElement(1, 0, 0.000000) m1.SetElement(1, 1, 0.692820) m1.SetElement(1, 2, 0.400000) m1.SetElement(2, 0, 0.200000) m1.SetElement(2, 1, -0.469846) m1.SetElement(2, 2, 0.813798) f11 = svtk.svtkTransformPolyDataFilter() f11.SetInputConnection(ap.GetOutputPort()) f11.SetTransform(t1) m11 = svtk.svtkDataSetMapper() m11.SetInputConnection(f11.GetOutputPort()) a11 = svtk.svtkActor() a11.SetMapper(m11) a11.GetProperty().SetColor(1, 0, 0) a11.GetProperty().SetRepresentationToWireframe() ren11 = svtk.svtkRenderer() ren11.SetViewport(0.0, 0.5, 0.25, 1.0) ren11.ResetCamera(-0.5, 0.5, -0.5, 0.5, -1, 1) ren11.AddActor(a11) renWin.AddRenderer(ren11) # inverse identity transform f12 = svtk.svtkTransformPolyDataFilter() f12.SetInputConnection(ap.GetOutputPort()) f12.SetTransform(t1.GetInverse()) m12 = svtk.svtkDataSetMapper()
points.InsertPoint(9, 1, 0.5, 0)
points.InsertPoint(10, 0.5, 1, 0)
points.InsertPoint(11, 0, 0.5, 0)
points.InsertPoint(12, 0.5, 0.5, 0)
i = 0
while i < 13:
aConvex.GetPointIds().InsertId(i, i)
i += 1
aConvexGrid = svtk.svtkUnstructuredGrid()
aConvexGrid.Allocate(1, 1)
aConvexGrid.InsertNextCell(aConvex.GetCellType(), aConvex.GetPointIds())
aConvexGrid.SetPoints(points)
# Display the cell
dsm = svtk.svtkDataSetMapper()
dsm.SetInputData(aConvexGrid)
a = svtk.svtkActor()
a.SetMapper(dsm)
a.GetProperty().SetColor(0, 1, 0)
# Contour and clip the cell with elevation scalars
ele = svtk.svtkElevationFilter()
ele.SetInputData(aConvexGrid)
ele.SetLowPoint(-1, -1, -1)
ele.SetHighPoint(1, 1, 1)
ele.SetScalarRange(-1, 1)
# Clip
#
clip = svtk.svtkClipDataSet()
aPolyVertexGrid.GetPointData().SetScalars(sizes) aPolyVertexGrid.GetPointData().AddArray(colors) sphere = svtk.svtkSphereSource() sphere.SetRadius(1.0) sphere.Update() glyphs = svtk.svtkGlyph3D() glyphs.ScalingOn() glyphs.SetColorModeToColorByScalar() glyphs.SetScaleModeToScaleByScalar() glyphs.SetScaleFactor(1) glyphs.SetInputData(aPolyVertexGrid) glyphs.SetSourceConnection(sphere.GetOutputPort()) glyphs.SetInputArrayToProcess(0, 0, 0, 0, "Sizes") glyphs.SetInputArrayToProcess(3, 0, 0, 0, "Colors") #puts [ glyphs GetScaleModeAsString] glyphsMapper = svtk.svtkDataSetMapper() glyphsMapper.SetInputConnection(glyphs.GetOutputPort()) glyphsActor = svtk.svtkActor() glyphsActor.SetMapper(glyphsMapper) glyphsActor.GetProperty().BackfaceCullingOn() ren1.SetBackground(.1, .2, .4) ren1.AddActor(glyphsActor) glyphsActor.GetProperty().SetDiffuseColor(1, 1, 1) ren1.ResetCamera() ren1.GetActiveCamera().Azimuth(30) ren1.GetActiveCamera().Elevation(20) ren1.GetActiveCamera().Dolly(1.25) ren1.ResetCameraClippingRange() renWin.Render() # render the image #
edgePoints.InsertPoint(2, 0.5, 0.25, 0) edgeScalars = svtk.svtkFloatArray() edgeScalars.SetNumberOfTuples(3) edgeScalars.InsertValue(0, 0.0) edgeScalars.InsertValue(1, 0.0) edgeScalars.InsertValue(2, 0.9) aEdge = svtk.svtkQuadraticEdge() aEdge.GetPointIds().SetId(0, 0) aEdge.GetPointIds().SetId(1, 1) aEdge.GetPointIds().SetId(2, 2) aEdgeGrid = svtk.svtkUnstructuredGrid() aEdgeGrid.Allocate(1, 1) aEdgeGrid.InsertNextCell(aEdge.GetCellType(), aEdge.GetPointIds()) aEdgeGrid.SetPoints(edgePoints) aEdgeGrid.GetPointData().SetScalars(edgeScalars) aEdgeMapper = svtk.svtkDataSetMapper() aEdgeMapper.SetInputData(aEdgeGrid) aEdgeMapper.ScalarVisibilityOff() aEdgeActor = svtk.svtkActor() aEdgeActor.SetMapper(aEdgeMapper) aEdgeActor.GetProperty().SetRepresentationToWireframe() aEdgeActor.GetProperty().SetAmbient(1.0) # Quadratic triangle triPoints = svtk.svtkPoints() triPoints.SetNumberOfPoints(6) triPoints.InsertPoint(0, 0.0, 0.0, 0.0) triPoints.InsertPoint(1, 1.0, 0.0, 0.0) triPoints.InsertPoint(2, 0.5, 0.8, 0.0) triPoints.InsertPoint(3, 0.5, 0.0, 0.0) triPoints.InsertPoint(4, 0.75, 0.4, 0.0)
def testQuadricCut(self):
solidTexture = (255, 255)
clearTexture = (255, 0)
edgeTexture = (0, 255)
def makeBooleanTexture(caseNumber, resolution, thickness):
#global solidTexture, clearTexture, edgeTexture
booleanTexturecaseNumber = svtk.svtkBooleanTexture()
booleanTexturecaseNumber.SetXSize(resolution)
booleanTexturecaseNumber.SetYSize(resolution)
booleanTexturecaseNumber.SetThickness(thickness)
if caseNumber == 0:
booleanTexturecaseNumber.SetInIn(solidTexture)
booleanTexturecaseNumber.SetOutIn(solidTexture)
booleanTexturecaseNumber.SetInOut(solidTexture)
booleanTexturecaseNumber.SetOutOut(solidTexture)
booleanTexturecaseNumber.SetOnOn(solidTexture)
booleanTexturecaseNumber.SetOnIn(solidTexture)
booleanTexturecaseNumber.SetOnOut(solidTexture)
booleanTexturecaseNumber.SetInOn(solidTexture)
booleanTexturecaseNumber.SetOutOn(solidTexture)
elif caseNumber == 1:
booleanTexturecaseNumber.SetInIn(clearTexture)
booleanTexturecaseNumber.SetOutIn(solidTexture)
booleanTexturecaseNumber.SetInOut(solidTexture)
booleanTexturecaseNumber.SetOutOut(solidTexture)
booleanTexturecaseNumber.SetOnOn(edgeTexture)
booleanTexturecaseNumber.SetOnIn(edgeTexture)
booleanTexturecaseNumber.SetOnOut(solidTexture)
booleanTexturecaseNumber.SetInOn(edgeTexture)
elif caseNumber == 2:
booleanTexturecaseNumber.SetInIn(solidTexture)
booleanTexturecaseNumber.SetOutIn(clearTexture)
booleanTexturecaseNumber.SetInOut(solidTexture)
booleanTexturecaseNumber.SetOutOut(solidTexture)
booleanTexturecaseNumber.SetOnOn(edgeTexture)
booleanTexturecaseNumber.SetOnIn(edgeTexture)
booleanTexturecaseNumber.SetOnOut(solidTexture)
booleanTexturecaseNumber.SetInOn(solidTexture)
booleanTexturecaseNumber.SetOutOn(edgeTexture)
elif caseNumber == 3:
booleanTexturecaseNumber.SetInIn(clearTexture)
booleanTexturecaseNumber.SetOutIn(clearTexture)
booleanTexturecaseNumber.SetInOut(solidTexture)
booleanTexturecaseNumber.SetOutOut(solidTexture)
booleanTexturecaseNumber.SetOnOn(edgeTexture)
booleanTexturecaseNumber.SetOnIn(clearTexture)
booleanTexturecaseNumber.SetOnOut(solidTexture)
booleanTexturecaseNumber.SetInOn(edgeTexture)
booleanTexturecaseNumber.SetOutOn(edgeTexture)
elif caseNumber == 4:
booleanTexturecaseNumber.SetInIn(solidTexture)
booleanTexturecaseNumber.SetOutIn(solidTexture)
booleanTexturecaseNumber.SetInOut(clearTexture)
booleanTexturecaseNumber.SetOutOut(solidTexture)
booleanTexturecaseNumber.SetOnOn(edgeTexture)
booleanTexturecaseNumber.SetOnIn(solidTexture)
booleanTexturecaseNumber.SetOnOut(edgeTexture)
booleanTexturecaseNumber.SetInOn(edgeTexture)
booleanTexturecaseNumber.SetOutOn(solidTexture)
elif caseNumber == 5:
booleanTexturecaseNumber.SetInIn(clearTexture)
booleanTexturecaseNumber.SetOutIn(solidTexture)
booleanTexturecaseNumber.SetInOut(clearTexture)
booleanTexturecaseNumber.SetOutOut(solidTexture)
booleanTexturecaseNumber.SetOnOn(edgeTexture)
booleanTexturecaseNumber.SetOnIn(edgeTexture)
booleanTexturecaseNumber.SetOnOut(edgeTexture)
booleanTexturecaseNumber.SetInOn(clearTexture)
booleanTexturecaseNumber.SetOutOn(solidTexture)
elif caseNumber == 6:
booleanTexturecaseNumber.SetInIn(solidTexture)
booleanTexturecaseNumber.SetOutIn(clearTexture)
booleanTexturecaseNumber.SetInOut(clearTexture)
booleanTexturecaseNumber.SetOutOut(solidTexture)
booleanTexturecaseNumber.SetOnOn(edgeTexture)
booleanTexturecaseNumber.SetOnIn(edgeTexture)
booleanTexturecaseNumber.SetOnOut(edgeTexture)
booleanTexturecaseNumber.SetInOn(edgeTexture)
booleanTexturecaseNumber.SetOutOn(edgeTexture)
elif caseNumber == 7:
booleanTexturecaseNumber.SetInIn(clearTexture)
booleanTexturecaseNumber.SetOutIn(clearTexture)
booleanTexturecaseNumber.SetInOut(clearTexture)
booleanTexturecaseNumber.SetOutOut(solidTexture)
booleanTexturecaseNumber.SetOnOn(edgeTexture)
booleanTexturecaseNumber.SetOnIn(clearTexture)
booleanTexturecaseNumber.SetOnOut(edgeTexture)
booleanTexturecaseNumber.SetInOn(clearTexture)
booleanTexturecaseNumber.SetOutOn(edgeTexture)
elif caseNumber == 8:
booleanTexturecaseNumber.SetInIn(solidTexture)
booleanTexturecaseNumber.SetOutIn(solidTexture)
booleanTexturecaseNumber.SetInOut(solidTexture)
booleanTexturecaseNumber.SetOutOut(clearTexture)
booleanTexturecaseNumber.SetOnOn(edgeTexture)
booleanTexturecaseNumber.SetOnIn(solidTexture)
booleanTexturecaseNumber.SetOnOut(edgeTexture)
booleanTexturecaseNumber.SetInOn(solidTexture)
booleanTexturecaseNumber.SetOutOn(edgeTexture)
elif caseNumber == 9:
booleanTexturecaseNumber.SetInIn(clearTexture)
booleanTexturecaseNumber.SetInOut(solidTexture)
booleanTexturecaseNumber.SetOutIn(solidTexture)
booleanTexturecaseNumber.SetOutOut(clearTexture)
booleanTexturecaseNumber.SetOnOn(edgeTexture)
booleanTexturecaseNumber.SetOnIn(edgeTexture)
booleanTexturecaseNumber.SetOnOut(edgeTexture)
booleanTexturecaseNumber.SetInOn(edgeTexture)
booleanTexturecaseNumber.SetOutOn(edgeTexture)
elif caseNumber == 10:
booleanTexturecaseNumber.SetInIn(solidTexture)
booleanTexturecaseNumber.SetInOut(solidTexture)
booleanTexturecaseNumber.SetOutIn(clearTexture)
booleanTexturecaseNumber.SetOutOut(clearTexture)
booleanTexturecaseNumber.SetOnOn(edgeTexture)
booleanTexturecaseNumber.SetOnIn(edgeTexture)
booleanTexturecaseNumber.SetOnOut(edgeTexture)
booleanTexturecaseNumber.SetInOn(solidTexture)
booleanTexturecaseNumber.SetOutOn(clearTexture)
elif caseNumber == 11:
booleanTexturecaseNumber.SetInIn(clearTexture)
booleanTexturecaseNumber.SetInOut(solidTexture)
booleanTexturecaseNumber.SetOutIn(clearTexture)
booleanTexturecaseNumber.SetOutOut(clearTexture)
booleanTexturecaseNumber.SetOnOn(edgeTexture)
booleanTexturecaseNumber.SetOnIn(clearTexture)
booleanTexturecaseNumber.SetOnOut(edgeTexture)
booleanTexturecaseNumber.SetInOn(edgeTexture)
booleanTexturecaseNumber.SetOutOn(clearTexture)
elif caseNumber == 12:
booleanTexturecaseNumber.SetInIn(solidTexture)
booleanTexturecaseNumber.SetInOut(clearTexture)
booleanTexturecaseNumber.SetOutIn(solidTexture)
booleanTexturecaseNumber.SetOutOut(clearTexture)
booleanTexturecaseNumber.SetOnOn(edgeTexture)
booleanTexturecaseNumber.SetOnIn(solidTexture)
booleanTexturecaseNumber.SetOnOut(clearTexture)
booleanTexturecaseNumber.SetInOn(edgeTexture)
booleanTexturecaseNumber.SetOutOn(edgeTexture)
elif caseNumber == 13:
booleanTexturecaseNumber.SetInIn(clearTexture)
booleanTexturecaseNumber.SetInOut(clearTexture)
booleanTexturecaseNumber.SetOutIn(solidTexture)
booleanTexturecaseNumber.SetOutOut(clearTexture)
booleanTexturecaseNumber.SetOnOn(edgeTexture)
booleanTexturecaseNumber.SetOnIn(edgeTexture)
booleanTexturecaseNumber.SetOnOut(clearTexture)
booleanTexturecaseNumber.SetInOn(clearTexture)
booleanTexturecaseNumber.SetOutOn(edgeTexture)
elif caseNumber == 14:
booleanTexturecaseNumber.SetInIn(solidTexture)
booleanTexturecaseNumber.SetInOut(clearTexture)
booleanTexturecaseNumber.SetOutIn(clearTexture)
booleanTexturecaseNumber.SetOutOut(clearTexture)
booleanTexturecaseNumber.SetOnOn(edgeTexture)
booleanTexturecaseNumber.SetOnIn(edgeTexture)
booleanTexturecaseNumber.SetOnOut(clearTexture)
booleanTexturecaseNumber.SetInOn(edgeTexture)
booleanTexturecaseNumber.SetOutOn(clearTexture)
elif caseNumber == 15:
booleanTexturecaseNumber.SetInIn(clearTexture)
booleanTexturecaseNumber.SetInOut(clearTexture)
booleanTexturecaseNumber.SetOutIn(clearTexture)
booleanTexturecaseNumber.SetOutOut(clearTexture)
booleanTexturecaseNumber.SetOnOn(clearTexture)
booleanTexturecaseNumber.SetOnIn(clearTexture)
booleanTexturecaseNumber.SetOnOut(clearTexture)
booleanTexturecaseNumber.SetInOn(clearTexture)
booleanTexturecaseNumber.SetOutOn(clearTexture)
booleanTexturecaseNumber.Update()
return booleanTexturecaseNumber
# A list of positions
positions = []
positions.append((-4, 4, 0))
positions.append((-2, 4, 0))
positions.append((0, 4, 0))
positions.append((2, 4, 0))
positions.append((-4, 2, 0))
positions.append((-2, 2, 0))
positions.append((0, 2, 0))
positions.append((2, 2, 0))
positions.append((-4, 0, 0))
positions.append((-2, 0, 0))
positions.append((0, 0, 0))
positions.append((2, 0, 0))
positions.append((-4, -2, 0))
positions.append((-2, -2, 0))
positions.append((0, -2, 0))
positions.append((2, -2, 0))
ren = svtk.svtkRenderer()
renWin = svtk.svtkRenderWindow()
renWin.AddRenderer(ren)
# define two elliptical cylinders
quadric1 = svtk.svtkQuadric()
quadric1.SetCoefficients(1, 2, 0, 0, 0, 0, 0, 0, 0, -.07)
quadric2 = svtk.svtkQuadric()
quadric2.SetCoefficients(2, 1, 0, 0, 0, 0, 0, 0, 0, -.07)
# create a sphere for all to use
aSphere = svtk.svtkSphereSource()
aSphere.SetPhiResolution(50)
aSphere.SetThetaResolution(50)
# create texture coordinates for all
tcoords = svtk.svtkImplicitTextureCoords()
tcoords.SetInputConnection(aSphere.GetOutputPort())
tcoords.SetRFunction(quadric1)
tcoords.SetSFunction(quadric2)
aMapper = svtk.svtkDataSetMapper()
aMapper.SetInputConnection(tcoords.GetOutputPort())
# create a mapper, sphere and texture map for each case
aTexture = []
anActor = []
for i in range(0, 16):
aTexture.append(svtk.svtkTexture())
aTexture[i].SetInputData(makeBooleanTexture(i, 256, 1).GetOutput())
aTexture[i].InterpolateOff()
aTexture[i].RepeatOff()
anActor.append(svtk.svtkActor())
anActor[i].SetMapper(aMapper)
anActor[i].SetTexture(aTexture[i])
anActor[i].SetPosition(positions[i])
anActor[i].SetScale(2.0, 2.0, 2.0)
ren.AddActor(anActor[i])
ren.SetBackground(0.4392, 0.5020, 0.5647)
ren.ResetCamera()
ren.GetActiveCamera().Zoom(1.4)
renWin.SetSize(500, 500)
# render and interact with data
iRen = svtk.svtkRenderWindowInteractor()
iRen.SetRenderWindow(renWin)
renWin.Render()
img_file = "quadricCut.png"
svtk.test.Testing.compareImage(
iRen.GetRenderWindow(),
svtk.test.Testing.getAbsImagePath(img_file),
threshold=25)
svtk.test.Testing.interact()